Next Article in Journal
Autonomous Mobile Robot with Attached Multispectral Camera to Monitor the Development of Crops and Detect Nutrient and Water Deficiencies in Vertical Farms
Previous Article in Journal
Metabolome and Transcriptome Analysis Reveal the Accumulation Mechanism of Carotenoids and the Causes of Color Differences in Persimmon (Diospyros kaki Thunb.) Fruits
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ant Diversity Is Enhanced by Ecological Infrastructures in Agroecosystems: A Case Study in Irrigated Mediterranean Farmland

by
Vera Zina
1,2,*,
André Fonseca
1,2,
Gonçalo Duarte
1,2,
Sofia Conde
1,2,
Maria Rosário Fernandes
1,2,
Maria Teresa Ferreira
1,2 and
José Carlos Franco
1,2
1
Forest Research Centre, School of Agriculture, University of Lisbon, Tapada da Ajuda, 1349-017 Lisbon, Portugal
2
Laboratory for Sustainable Land Use and Ecosystem Services, School of Agriculture, University of Lisbon, Tapada da Ajuda, 1349-017 Lisbon, Portugal
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(11), 2690; https://doi.org/10.3390/agronomy12112690
Submission received: 25 September 2022 / Revised: 24 October 2022 / Accepted: 26 October 2022 / Published: 29 October 2022
(This article belongs to the Topic Insects in Sustainable Agroecosystems)

Abstract

:
We aimed at assessing the role of ecological infrastructures (EI) in promoting ant biodiversity in floodplain Mediterranean agricultural crops. We examined and compared ant communities at the interface between EI (remnant vegetation patches) and adjoining agricultural matrix (maize, rice, others) in irrigated farmland. The study was conducted in 2019, in two agricultural landscapes in the valleys of the rivers Tagus and Sorraia, Central Portugal. We used the Akaike information criterion for model selection and to distinguish among a set of possible models describing the relationship between: the ant richness in the agricultural matrix and drivers associated with the surrounding landscape and crop type; the ant richness in EI and the habitat quality of EI patches, the characteristics of the surrounding landscape, and the presence of invasive ant species. We found that: EI patches supported a higher ant diversity and an overall specialized ant community, distinctive from the agricultural matrix; location but not vegetation physiognomy influenced ant diversity; ant richness within the agricultural matrix decreased with the distance to the EI, and that this relationship was influenced by the crop type; and that ant richness in the EI was associated with the absence of the invasive Argentine ant and the area of terrestrial EI in the surrounding landscape.

1. Introduction

Agricultural land use is pointed out as one of the major drivers of land use change, representing around 40% of the land surface [1]. In this context, agricultural intensification, characterized by crop monocultures, high soil disturbance, and use of pesticides, is generally associated with landscape simplification and biodiversity losses, reducing ecosystem services on which agriculture depends [2,3,4]. To overcome these negative effects, a sound management plan should be adopted to enhance habitat heterogeneity and support biodiversity and ecosystem services by conserving and promoting ecological infrastructures (EI), including a network of natural, semi-natural, and restored patches at different spatial scales [5,6].
In agroecosystems, EI are considered of such ecological value that different eco-schemes (i.e., payment aids) were introduced in the European Union, through the Common Agriculture Policy (CAP), to encourage farmers to progressively dedicate a proportion of their arable land to non-productive biodiversity-friendly features [7,8]. One of the present Green Deal targets under the CAP strategic plan is to bring back at least 10% of agricultural areas under high-diversity landscape features by 2030 [8]. These may include hedges, rows of trees, field copses, ponds, or fallow land [7]. However, these ecological approaches have been focused mainly on quantitative aspects, with higher compensation for those who have larger areas devoted to EI. Above all, the landscape features that compose the EI must fulfill the purpose of providing suitable habitats (of good quality) for promoting biodiversity and ecosystem services (ES) [9]. Two of the most important factors influencing this are the landscape configuration of EI (i.e., the spatial arrangement of land use patches) and their composition (i.e., the relative proportion of habitat types) [10,11,12]. In the last decades, quality-based assessment tools have been developed to evaluate the potential value of a certain area to harbor biodiversity, such as the ‘Indice de Biodiversité Potentielle’ (IBP) [13] and the Habitat Ecological Infrastructure Diversity Index (HEIDI) [14]. Nevertheless, further knowledge is needed to guarantee successful green measures implementation, such as guidelines for selecting the best elements of the landscape (e.g., in terms of composition and configuration) and redesigning sustainable and resilient crop production systems [2].
In the Mediterranean region, drought conditions and water scarcity pushed land use activities to the vicinity of watercourses, and EI are usually distributed among cultivated plots, along roads, paths, or water features [15,16]. As a result, patches are few in number and present simple configurations [14]. Nevertheless, these patches are still high-value habitats playing a crucial role in supporting food, refuge, and breeding areas for several biological communities [17,18,19,20]. In this regard, there is a growing awareness that agroecosystems should be a priority in the biological conservation agenda because some agroecosystems are repositories of high levels of biodiversity, including ants [21]. Some studies have been carried out on ant communities in Mediterranean drylands (e.g., [22,23]). There is a lack of information concerning irrigated lands.
Ants (Hymenoptera, Formicidae) are a group of eusocial and colonial insects with more than 16,500 described species worldwide [24]. Eusociality confers marked advantages in terms of ecological dominance, resources foraging, defense against enemies, and may allow ants to adapt or tolerate future environmental change [25,26,27]. Ants have evolved into distinctive life strategies, resulting in the interaction with many different animal and plant communities [28,29,30]. They also carry out important ecological functions and services, such as nutrient cycling, decomposition, soil movement, seed dispersal, and pest regulation [31,32,33]. However, these services are dependent on ant biodiversity, which has been threatened worldwide by agricultural intensification [34,35,36]. For instance, the intensification of coffee plantations in the tropics (involving the conversion of rustic systems, with shaded trees, to unshaded monocultures) significantly reduced ant species richness [21]. In temperate regions, ant richness and abundance have been also affected by agriculture, and its impacts may vary depending on agricultural practices [34] and landscape components [37,38,39]. Although the Mediterranean basin is considered a hotspot for ant richness [40], intensive agricultural landscapes in this region exhibited low potential for ant biodiversity, when compared to more extensive agricultural systems [4,14]. Despite considerable efforts to protect ant biodiversity in human-disturbed landscapes, many of the existing refuges are small, fragmented, isolated, or of poor quality [14]. In addition, invasion by exotic species is often reported, and changes in competitive interactions or colonization processes may also affect ant assemblages in disturbed areas [41].
To broaden our understanding of the role of EI in ant biodiversity, particularly in annual cropping systems, we examined ant communities at the interface between EI and the adjoining agricultural matrix in two intensively irrigated agricultural areas located in the Sorraia and Tagus river valleys in Portugal. We aimed at assessing the role of EI in promoting ant biodiversity in floodplain Mediterranean agricultural crops. We made four hypotheses, as follows: (1) the reduced habitat structure and lower food resources availability in the agricultural matrix would likely affect ant richness and composition, in comparison with EI; (2) the effect of EI on ant diversity would differ with their typologies, namely, location and vegetation physiognomy; (3) the ant richness in the agricultural matrix is dependent on the type of crops produced, as well as on other drivers associated with the surrounding landscape; and (4) the ant richness in EI is influenced by the habitat quality of EI patches, but also by the characteristics of the surrounding landscape, and can be affected by the presence of invasive ant species.
The results from this study will enable a better understanding of the role of EI in ant biodiversity. They will provide valuable clues for the improvement of management plans when aiming at the conservation and restoration of Mediterranean agricultural floodplains.

2. Materials and Methods

2.1. Study Area

The study was conducted in 2019, in two agricultural landscapes in the valleys of the rivers Tagus and Sorraia, Central Portugal (Figure 1). Both landscapes comprise irrigated cropland, composed of annual crops, with a predominance of rice paddies (Oryza sativa L.) in the alluvial plains of the Sorraia river and maize fields (Zea mays L.) in the Tagus river valley.
The region is characterized by a mild climate with rainy winters (average rainfall of circa 900 mm) and hot, dry summers (average annual air temperature of around 15 °C), with values varying regularly throughout the year, with a maximum in August and a minimum in January [42]. The rainfall pattern exhibits strong seasonal and inter-annual variability, with high floods usually occurring in autumn or early winter, a gradual decline in discharge, and subsequent drying out during late spring and summer [43].

2.2. Sampling Design

Sampling was carried out in 130 sites, based on a balanced subset of randomly distributed points, including 59 sites within the agricultural matrix (27 and 32 in the Tagus and Sorraia basin, respectively) and 71 sites in the EI patches (37 and 34 in the Tagus and Sorraia basin, respectively). We used field data and a geographic information system (GIS) image-based approach to characterize the sampling sites. The agricultural matrix and EI patches were manually and individually digitized through a visual analysis of the Esri World Imagery layer (1:1000 scale) [44]. We selected a minimum mapping unit of 200 m², with a minimum width of 5 m and a minimum gap distance among patches of 10 m [45,46]. These mapping thresholds were established to represent the minimum patch size and minimum distance between patches that are considered ecologically meaningful for Mediterranean ant communities [47,48]. The EI vegetation patches were then classified according to the vegetation physiognomy, i.e., those dominated by trees and shrubs (hereafter referred to as “woody EI”), and those where trees were absent or rare, including open areas with scarce vegetation or dominated by herbaceous plants and low bushes (hereafter referred to as “herbaceous EI”). We also classified EI vegetation patches according to their location, i.e., those located contiguous to a watercourse (hereafter referred to as “riparian EI”) and those that were non-contiguous to a watercourse (hereafter referred to as “terrestrial EI”).
A description of the variables used in the study is summarized in Table 1. Variables related to proximity and area/density were calculated using the Patch Analyst Vector format (ArcGis10.6) extension [49]. A 200 m buffer radius was considered around each sampling site to measure the total area occupied by the EI patch types and the agricultural land [37,38]. Habitat quality variables included shrub richness, the occurrence of the invasive Argentine ant, Linepithema humile (Mayr), and a measure of the potential habitat quality of the EI–HEIDI index, sensu Fonseca et al. [14]. Through the information collected in the field, we used different metrics to calculate the HEIDI index, including: vegetation structure (e.g., native-invasive plant species, vertical strata); the presence of specific habitats (e.g., microhabitat at trees, dead wood trunks on the ground, leaf litter cover); vegetation management (e.g., understory clearing); and floristic suitability (i.e., the ecological value of plant taxa for the provision of ES such as their potential for myrmecochory) (see Fonseca et al. [14] for a detailed description of the HEIDI index development).

2.3. Ant Sampling and Identification

Ant sampling was carried out in early summer, between June and July, i.e., the period with the highest activity for Mediterranean ant communities [50]. The sampling was performed using pitfall traps, a method that has been used in many studies of ant communities worldwide (e.g., [51,52,53]). It is a simple, cost-effective method for collecting epigaeic ants that provides good results in assessing species richness and composition patterns while allowing for continuous day and night sampling [54,55,56,57].
The sampling design consisted of five pitfall traps per site, distributed in a linear transect, with the middle trap considered as the geo-referenced point. A similar number of pitfall traps per site has been used by different authors (e.g., [58,59,60,61]). Each trap was five meters apart from the neighboring ones. Pitfall traps consisted of 100 mL plastic containers placed flush with the ground, and partially filled with a solution of 30% propylene glycol and a few drops of detergent, to retain and preserve the intercepted ants. They were left in the field for about 48 h. This sampling time was shown to be more cost effective and robust enough to estimate ant diversity, in comparison with 14-day sampling [62]. The ants collected per site were put together for sorting and identification.
All the entomological material was sorted and identified by the first author, who has more than 10 years of experience in the study and taxonomic identification of ants. Identification was carried out at the species level (except for Solenopsis spp.), under a stereomicroscope, using taxonomic keys for Portugal and Iberian Peninsula and knowledge expertise [63,64,65]. All identified specimens were preserved in 96% alcohol and kept at the laboratory of entomology of Instituto Superior de Agronomia.

2.4. Statistical Analysis

For testing Hypotheses 1 and 2, we used the ant richness (i.e., number of observed species in each sampling site) and ant species composition occurring in the EI and within the agricultural matrix. A one-way ANOVA (aov function in R package stats [66]) was used to determine statistically significant differences between and among independent groups (EI vs. agricultural matrix, herbaceous EI vs. woody EI vs. agricultural matrix, and riparian EI vs. terrestrial EI vs. agricultural matrix). If the p-value was statistically significant (p < 0.05) a post-hoc for multiple comparisons was performed using Tukey’s HSD Test (TukeyHSD function in R package stats [66]). Boxplots were created using graphics and ggplot2 R packages [66,67].
The differences among the communities were investigated using a permutational multivariate analysis of variance (PERMANOVA) (adonis function in R package vegan [68]) based on a Bray–Curtis distance matrix with 999 permutations. We used a site-per-species matrix containing incidence data for ant species at each sampling site. Ordination plots were created using metaMDS and ordiellipse functions in R package vegan [68].
To test Hypothesis 3, we performed generalized linear mixed models (GLMMs) (glmmTMB function in R package glmmTMB [69]). We used ant species richness sampled within the agricultural matrix as a dependent variable, and the models were fitted with a negative binomial distribution to handle overdispersion. The predictors were included as fixed factors, and all continuous ones were standardized (i.e., rescaled to the same unit), enabling comparisons of effect magnitude. We used the Akaike information criterion corrected for small sample size (AICc) to test the statistical relevance of including crop types (maize fields, rice paddies, and other mixed crops) [70]. The Argentine ant occurrence (present, absence) and sampling month (June, July) alone or together were added as random effects to account for the sampling design. Of these, we only include crop types as a predictor. Validation of the full multivariate models was carried out using DHArMa and performance R packages [71,72] with the help of diagnostic plots (Figure A1). We generated all possible models based on the full one and performed model selection with the MuMIn package [73]. Models were selected based on Akaike weights and AIC differences (ΔAIC) from the best-fitted model and were considered to be equally supported if AIC was less than two units [70]. Based on model-averaging, we estimate the predicted responses of ant species richness with associated 95% confidence intervals [74]. We performed also GLMMs in a similar way to test Hypothesis 4, using ant species richness sampled in the EI, and further assessed whether the influence of different predictors (the % area of agricultural matrix, riparian EI, and terrestrial EI, in a 200 m buffer; shrub richness; the HEIDI quality index sensu Fonseca et al. [14]; and the Argentine ant presence) could shape ant communities in the EI.
Data were stored in spreadsheets Microsoft® Excel for Mac version 16.16.27 (201012), and all analyses were conducted in the R environment [66].

3. Results

3.1. Agricultural Matrix Versus Ecological Infrastructures (Hypothesis 1)

3.1.1. Species Richness

We recorded 17 072 ants belonging to 47 species from 20 genera. We observed higher species richness in EI than in the agricultural matrix (Figure 2). Mean species richness was significantly different between EI and the agricultural matrix (p < 0.01, 95% C.I. = [1.72, 3.51]; Figure 2a), riparian EI and the agricultural matrix (p = 0.02, 95% C.I. = [−2.72, −0.22]), and between terrestrial EI and the agricultural matrix (p < 0.01, 95% C.I. = [0.46, 3.13]; Figure 2c).

3.1.2. Community Composition

Ant communities in EIs were significantly different from those in the agricultural matrix (F1,117 = 5.22, p = 0.007; Figure 3a and Figure 4). All 20 species observed in the agricultural matrix were also present in the EI, except for Cardiocondyla mauritanica Forel, which was only identified in the agricultural matrix; 27 ant species were only observed in association with the EI. Linepithema humile and Tapinoma nigerrimum Nylander were the most frequent species in EI and the agricultural matrix, respectively. In most of the cases, ant species were more frequent in EI, in comparison with the agricultural matrix. However, the opposite was observed in the case of T. nigerrimum, Tetramorium forte Forel, Hypoponera eduardi Forel, Messor bouvieri Bondroit, and Formica cunicularia Latreille.

3.2. Effect of the Ecological Infrastructure Typology (Hypothesis 2)

3.2.1. Species Richness

Ant species richness in terrestrial EI was significantly higher than in riparian EI (p < 0.01, 95% C.I. = [1.81, 4.72]; Figure 2c). No significant differences in ant species richness were observed regarding vegetation physiognomy, i.e., woody EI vs. herbaceous EI (F1, 126 = 0.03, p = 0.86; Figure 2b). Nevertheless, the total number of ant species identified in woody EI (40 species) was higher than in herbaceous EI (33 species) (Figure 4).

3.2.2. Community Composition

Ant communities in riparian EI were significantly different from those in terrestrial EI (F1,117 = 4.18, p = 0.009; Figure 3c). Of the 27 ant species observed only in EI, 10 and 5 species were specific to terrestrial and riparian EI, respectively. Most of the ant species were more frequent in the terrestrial EI than in the riparian EI, except L. humile and Lasius grandis Forel, which were favored by riparian habitats (Figure 4).
No significant differences were found between ant communities of herbaceous EI and woody EI (F1,117 = 1.61, p = 0.150; Figure 3b). Nevertheless, most ant species were more frequent in woody EI, occupying a total area of 564.95 ha, in comparison with herbaceous EI, which occupied 117.54 ha. Some ant species were identified only in one type of EI. For example, Tapinoma erraticum (Latreille), Goniomma hispanicum (André), Crematogaster auberti Emery, Tetramorium caespitum (Linnaeus), Temnothorax tuberum (Fabricius), and Camponotus micans (Nylander) were observed only in herbaceous EI, whereas T. lichtensteini (Bondroit), T. angustulus (Nylander), Myrmica ruginodis Nylander, C. fallax (Nylander), T. unifasciatus (Latreille), Oxyopomyrmex saulcyi Emery, Lasius brunneus (Latreille), F. cunicularia, and Colobopsis truncata (Spinola) were found only in woody EI.

3.3. Drivers of Ant Richness in the Agricultural Matrix (Hypothesis 3)

Ant species richness within the agricultural matrix was significantly associated with the distance to the EI and crop type, but not with the distance to the river or to urban areas (Table 2, Figure 5).
The best-fitted model included two out of four variables, i.e., the distance to EI and crop type. This model performed better than the others as it carried 52% of the cumulative model weight and has the lowest AIC score (Table A1). In addition, the next-two-best models and the null model were more than two AIC units higher (2.4, 2.4, and 11.1, respectively) and carried only 15.7%, 15.6%, and 0.2%, respectively, of the cumulative model weight (Table A1).
Predicted responses of the ant species richness within the agricultural matrix showed a significant negative effect of the distance to the nearest EI, indicating that ant species richness decreases with increasing distance to EI (Table 1, Figure 5). The number of ant species is predicted to drop to half if the EI patch is at a 200 m distance and almost zero if at a 600 m distance. This trend effect is irrespective of the crop type, yet is more pronounced in crops with higher ant species richness, such as rice paddies (Figure 5). Our model suggested that this type of crop is predicted to harbor a significantly higher number of species when compared with maize fields (Table 2, Figure 5). Nevertheless, these are both low levels of ant richness and are predicted to aggravate without the presence of EI patches at a short distance (Figure 5).

3.4. Effect of Ecological Infrastructure Habitat Quality and Characteristics of the Surrounding Landscape (Hypothesis 4)

Ant species richness in the EI was significantly associated with the absence of Argentine ant and the area of terrestrial EI in the surrounding landscape but not with the area of riparian EI or agricultural land, shrub richness, or HEIDI quality index (Table A2, Figure 6). Three models showed a ΔAIC < 2, yet the most parsimonious model of the three appear to be the one including fewer parameters (two out of six), i.e., the Argentine ant occurrence and the area of terrestrial EI in the surrounding landscape (Table 3, Table A2). The model with the lowest AIC includes another parameter, i.e., the area of riparian EI in the surrounding landscape. These two models explained 37.4% of the cumulative weight (Table A2).
Predicted responses of the ant species richness in the EI, showed a significant positive effect on the area occupied by EI in the surrounding landscape (Figure 6). The opposite effect is predicted to occur with the area occupied by the riparian EI in the surrounding landscape although it was not statistically significant (Table 3, Figure 6).

4. Discussion

We are facing the challenge of managing trade-offs between food production and conserving biodiversity under intensive cropland systems [1]. Our results provide evidence that EI may have an important role in conserving ant biodiversity in agricultural landscapes and that its effect is influenced by its typology and habitat quality.
Focusing on ants, a key group in providing valuable insights into disturbed habitats and in promoting ecosystem functions and services across agroecosystems, we aimed at assessing the importance of landscape components in promoting ant biodiversity in a floodplain Mediterranean agricultural mosaic. We first hypothesized that the reduced habitat structure and lower food resource availability in the agricultural matrix would likely affect ant richness and community composition, in comparison with EI (Hypothesis 1). As expected we found higher ant species richness in EI (46 species) than in the agricultural matrix (20 species). Few other studies also provided support for this hypothesis. For example, Dauber and Wolters [35] registered higher ant species richness in more complex habitats including fallow lands (13 species) and meadows (11 species), than in arable land (8 species), in Central Europe. This is possibly related to the level of disturbance inherent to the agricultural production system and its clear differences in habitat complexity [35,75], in accordance with the “habitat heterogeneity hypothesis” [76,77], which predicts that more structurally complex habitats, with further vegetation layers, can provide additional ecological niches and resources, thus supporting more species. The unique requirements of specialized ant fauna [35], associated with life-history species traits, such as feeding habits, colony size, dispersal capabilities [78], and their ability to cope with disturbance [79], may also contribute to the observed differences in ant species communities between EI and agricultural land. In fact, ant species communities observed foraging in EI were significantly different from those in the agricultural matrix. We found 27 species only occurring in the EI patches, whereas all species observed in the agricultural matrix were also present in the EI, except C. mauritanica, which is an invasive species favored by semi-arid areas and anthropogenic disturbance [80]. It seems that the agricultural matrix is dominated by opportunist species, such as T. nigerrimum, T. forte, and other broadly adapted omnivorous (e.g., F. cunicularia) and generalist (e.g., Pheidole pallidula Nylander) species, with wide habitat tolerance and higher dispersal capabilities, able of exploiting disturbed areas [81,82]. Landscape conversion to monocultures has been reported to favor generalist ants with large colonies [78]. These ant species, which might be less affected by local or landscape heterogeneity, can play an important role in pest management, such as T. nigerrimum [83,84,85] and F. cunicularia [29], or in decomposition, such as T. nigerrimum, L. humile, Tetramorium spp., P. pallidula, and Cataglyphis spp. They are scavengers and active foragers, collecting and consuming living prey and dead (carrion) animal matter [32,86,87]. Nevertheless, other ecosystem services might be also influenced. For instance, weed control is a potential service mediated by harvester ants, e.g., Messor spp. [88], that might be valuable in crops that resort to the use of herbicide applications, such as maize [89]. On the other hand, some ant species may be associated with disservices [29,32]. For example, L. humile became a major pest in many areas around the world [90], and can affect pollination [91], natural ant seed dispersal [92], and biological control of pests in agroecosystems [93].
Regarding the effect of EI typology (Hypothesis 2), we found significant differences between terrestrial and riparian EI but not between herbaceous and woody EI. Ant species richness in terrestrial EI (41 species) was higher than in riparian EI (33 species). Most of the ant species occurring in both EI locations were more frequent in terrestrial than in riparian EI, except L. grandis and L. humile, which were favored by riparian habitats. These are numerically dominant species known to resist anthropogenic disturbances, and with efficient foraging recruitment mechanisms [94,95]. The lower species richness observed in riparian EI may be explained, at least in part, by the higher frequency of L. humile in these habitats, in comparison with terrestrial EI. In fact, the negative impact of this invasive species on the diversity of native ants has been documented in different ecosystems [96,97,98,99].
The presence of alien invasive plant species in the riparian EI may also influence ant species richness. In Mediterranean irrigated cropland, riparian habitats suffer from the introduction and establishment of invasive plant species, such as the giant reed Arundo donax L. [100], that might affect the availability of crucial habitat resources for ant species, provided by native vegetation [36].
On the other hand, ant species richness and species composition were not influenced by vegetation physiognomy, i.e., woody EI vs. herbaceous EI. We would expect a different result, since woody EI are more complex than herbaceous EI and thus are likely to provide more diverse and suitable habitats for ants, reported recently for intensive agricultural landscapes [14]. Nevertheless, most ant species were more frequent in woody than in herbaceous EI.
We further hypothesized that ant richness in the agricultural matrix would be dependent on the crop type, as well as on drivers associated with the surrounding landscape, including the distance to EI, urban areas, and water lines (Hypothesis 3). We found that ant species richness within the agricultural matrix decreased with the distance to the EI and that this relationship was influenced by the crop type. Similar results were reported by Armbrecht and Perfecto [101] in Mexican intensive farmland, which observed a dramatic reduction of ant species richness in function of the distance to the forest fragment. These results and, as mentioned before, the fact that almost all ant species observed in the agricultural matrix were also present in the EI suggest that ant assemblages in agricultural lands depend on the recruitment of ant species from EI. The EI may act as refuges, buffering the negative consequences of insecticide application in adjacent fields [102], and reservoirs with less disturbed habitats increasing the ant species richness of adjacent higher disturbed habitats [35]. The survival of specialist ant species relies on the existence of EI where they can find food, nesting, and foraging resources [41,103].
However, no association was found between ant richness in the agricultural matrix with distance to the river or urban areas, as hypothesized. We would expect a positive effect with the river proximity, since several studies suggested that species richness and diversity of invertebrate communities are positively correlated with the structure and composition of riparian vegetation patches [104,105,106]. The observed lack of relationship between ant richness and river proximity might be related to the relatively low quality of riparian patches in the studied agricultural systems [14], since well-preserved riparian corridors are essential for the long-term maintenance of myrmecofauna [107]. The higher frequency of the Argentine ant presence in riparian habitats, and the expected negative impact on ant richness might also have contributed to this lack of relationship. Regarding urban areas, studies examining the impacts of urbanization on ant richness report diversity loss [108] and changes in ecological interactions [109]. We also would expect similar results with a negative effect on ant richness in the proximity of urban areas. However, urban environments may still provide habitat heterogeneity for small organisms such as ants [110], and in particular, for dry-adapted, heat-tolerant ant species [111].
In terms of crop type, our model suggested that rice paddies are predicted to harbor a significantly higher number of ant species when compared with maize fields. This might be due to differences in higher land crop type occupation and intensification [112]. Maize is the major cereal crop in Portugal, corresponding to 66% of cereal-producing farms, in 2016, while rice represented only 1% of the farms, usually associated with areas of lower cultivation value [113]. Due to sampling conflicts inherent to the crop being flooded (rice paddies), the sampling was carried out exclusively at the edges of the plots in the dry bunds (i.e., levees). Nevertheless, edges between different land use types did not increase ant species richness at the landscape scale, nor were they unique habitats for specialized ant fauna [35]. Furthermore, rice paddies are likely to have more resources, since they harbor both aquatic and terrestrial organisms, including a high-rich arthropod fauna inhabiting the vegetation, water, and surrounding bunds [114,115].
Finally, we hypothesized that the ant richness in EI would be influenced by the habitat quality of EI patches, but also by the characteristics of the surrounding landscape, and could be affected by the presence of invasive ant species (Hypothesis 4). We found that ant species richness in the EI was associated with the absence of the invasive Argentine ant and the area of terrestrial EI in the surrounding landscape. This is consistent with other studies that found an increase in ant species richness with an increasing percentage cover of forest or fallow land [37,116]. Our predicted responses of the ant species richness in the EI showed a positive effect on the area occupied by EI in the surrounding landscape. However, the selected habitat quality variables were not significant contributors to explaining ant species richness. This contradicts our initial expectation that both habitat quality and landscape features should contribute to the prediction of ant species richness, as Mediterranean ants are strongly linked to both patch and landscape variables [37]. Nevertheless, the regression model in our study showed low values of explained variance, so the habitat quality variables used, i.e., shrub richness and HEIDI index, might not be suitable for ants. Furthermore, these are highly disturbed areas with low variability. In addition, the relative importance of local and landscape factors may depend on the landscape context. Local allocation of habitats such as field boundaries and management practices seem to have a higher impact in structurally simple landscapes than in complex ones, demanding an even greater effort [117]. In our study region, agricultural land dominates most of the landscape.
The other relevant factor affecting ant species richness was the presence of the Argentine ant in EI patches. In fact, this invasive species is known to thrive in riparian, urban, and agricultural habitats with Mediterranean climates [118], and to disturb native ant communities [96,97,98,99]. The dispersion of the Argentine ant in the Iberian Peninsula was predicted to be possible mostly along the coast, but also into inland areas along river valleys [119], as in the case of the study area.

5. Conclusions

Agricultural intensification has been responsible for reducing biodiversity and ecosystem services. The conservation and increment of EI within the agricultural landscape have thus been proposed as an approach to mitigate the negative ecological impacts of intensification, by promoting functional biodiversity. Nevertheless, there is a need to assess the quantity and quality of semi-natural habitats in agricultural farmlands, and to better understand the role of different types of EI in enhancing biodiversity and ES in agricultural systems [9]. Here we investigated the effect of EI and its typologies (i.e., location and vegetation physiognomy) on ant species richness and community composition in floodplain Mediterranean agricultural systems. Our results showed that ant diversity in agroecosystems is enhanced by EI and that ant assemblages present in agricultural lands are likely recruited from ant communities of EI. Therefore, EI may have a critical role in determining ant diversity and species composition of ant assemblages in agricultural farmlands, and consequently in the corresponding services and disservices. However, the effect of the EI on ant diversity was influenced by the EI location, the type of agricultural crops, as well as other drivers and characteristics of the surrounding landscape. In fact, whereas EI location was a significant factor influencing ant species richness and composition, with higher diversity registered in terrestrial than in riparian EI, vegetation physiognomy of EI seems to have no influence. On the other hand, ant diversity may be negatively affected by the presence of invasive species, such as the Argentine ant.
Our findings provide new insights into the role of EI in ant diversity conservation in agroecosystems and may help define habitat conservation and restoration guidelines. Nevertheless, management guidelines should also integrate knowledge about the different animal guilds and related services. For example, two other studies carried out in the same agricultural system showed that the EI location and vegetation physiognomy may influence differently different animal guilds and services [20,120]. These differences should be considered in the global management system.

Author Contributions

Conceptualization, V.Z., A.F., G.D., M.T.F., M.R.F. and J.C.F.; methodology, V.Z., A.F., G.D., M.R.F. and J.C.F.; formal analysis, V.Z., G.D., M.R.F. and J.C.F.; investigation, V.Z., A.F., G.D. and S.C.; data curation, V.Z., A.F. and G.D.; writing—original draft preparation, V.Z. and J.C.F.; writing—review and editing, V.Z., A.F., G.D., M.R.F., M.T.F. and J.C.F.; supervision, M.R.F. and J.C.F.; project administration, G.D., M.R.F. and M.T.F.; funding acquisition, M.T.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Optimus Prime (FCT-PTDC/ASP-AGR/29771/2017) and CERES Interreg IV-B SUDOE-SOE2/P5/F0551 and received backing from Forest Research Centre (CEF), a research unit funded by Fundação para a Ciência e a Tecnologia (FCT), Portugal, Grant Number UIDB/00239/2020 and Laboratory for Sustainable Land Use and Ecosystem Services–TERRA (LA/P/0092/2020). G.D., A.F., and VZ were supported by national funds via FCT: PTDC/ASP-AGR/29771/2017 and UIDP/00239/2020; PD/BD/142884/2018; and PD/BD/142882/2018, respectively.

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful to the stakeholders and the farmers associations (Agrotejo and Associação de Regantes e Beneficiários do Vale do Sorraia–ARBVS) for all the support during fieldwork. We thank António Faro and Susana Rocha for their help during fieldwork and Alice Le Révérend for assistance sorting the pitfall traps. We thank the four anonymous reviewers for their valuable comments and suggestions that helped improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Diagnostic plots showing no significant problems detected in (a) GLMMs relating the effects of different predictors on the ant species richness in the agricultural matrix, and in (b) GLMMs relating the effects of different predictors on the ant species richness in the ecological infrastructures.
Figure A1. Diagnostic plots showing no significant problems detected in (a) GLMMs relating the effects of different predictors on the ant species richness in the agricultural matrix, and in (b) GLMMs relating the effects of different predictors on the ant species richness in the ecological infrastructures.
Agronomy 12 02690 g0a1
Table A1. Description of the GLMMs relating the effects of different predictors on the ant species richness in the agricultural matrix. Models are ranked in ascending order of AIC values and the number of parameters in the model (K), log-likelihood (logLik), AICc, ΔAICc, and AICc weight (ωi) are given for each model. The “+” signal means the inclusion of the nominal variables in the models.
Table A1. Description of the GLMMs relating the effects of different predictors on the ant species richness in the agricultural matrix. Models are ranked in ascending order of AIC values and the number of parameters in the model (K), log-likelihood (logLik), AICc, ΔAICc, and AICc weight (ωi) are given for each model. The “+” signal means the inclusion of the nominal variables in the models.
Modelssnd(Int)dsp(Int)CropDist_EIDist_riverDist_urbanKlogLikAICc∆AICcωi
Richness~Crop + scale(Dist_EI)1.17++−0.39 5−1052220.00.518
Richness~Crop + scale(Dist_EI) + scale(Dist_river)1.15++−0.37−0.041 6−1052242.40.157
Richness~Crop + scale(Dist_EI) + scale(Dist_urban)1.17++−0.39 −0.0356−1052242.40.156
Richness~Crop1.18++ 4−1092274.70.049
Richness~Crop + scale(Dist_EI) + scale(Dist_river) + scale(Dist_urban)1.15++−0.38−0.032−0.0257−1052274.90.044
Richness~Crop + scale(Dist_river)1.11++ −0.124 5−1092286.40.021
Richness~Crop + scale(Dist_urban)1.18++ 0.0365−1092297.00.015
Richness~scale(Dist_EI) + scale(Dist_river)0.72+ −0.28−0.205 4−1102307.90.010
Richness~scale(Dist_EI)0.74+ −0.31 3−1122308.20.008
Richness~Crop + scale(Dist_river) + scale(Dist_urban)1.12++ −0.1380.0636−1082308.70.007
Richness~scale(Dist_EI) + scale(Dist_urban)0.73+ −0.33 −0.1164−1112319.60.004
Richness~scale(Dist_river)0.75+ −0.257 3−1132329.70.004
Richness~scale(Dist_EI) + scale(Dist_river) + scale(Dist_urban)0.71+ −0.29−0.186−0.0485−11023210.20.003
Richness~1 (Null model)0.77+ 2−11423311.10.002
Richness~scale(Dist_river) + scale(Dist_urban)0.75+ −0.255−0.00364−11423412.00.001
Richness~scale(Dist_urban)0.77+ −0.08143−11423512.90.001
EI—ecological infrastructure; Dist_EI—Distance to the closest EI patch; Dist_river—Distance to river; Dist_urban—Distance to urban area; Crop—Crop type (maize field, rice paddy, others—mixed types). Please see Table 1 for a detailed description.
Table A2. Description of the GLMMs relating the effects of different predictors on the ant species richness in the ecological infrastructure (EI) patches. Models are ranked in ascending order of AIC values and the number of parameters in the model (K), log-likelihood (logLik), AICc, ΔAICc, and AICc weight (ωi) are given for each model. The “+” signal means the inclusion of the nominal variables in the models.
Table A2. Description of the GLMMs relating the effects of different predictors on the ant species richness in the ecological infrastructure (EI) patches. Models are ranked in ascending order of AIC values and the number of parameters in the model (K), log-likelihood (logLik), AICc, ΔAICc, and AICc weight (ωi) are given for each model. The “+” signal means the inclusion of the nominal variables in the models.
Modelssnd(Int)dsp(Int)Arg_antMatrix_areaRiparian_areaShrub_
richness
Terrestrial_areaHEIDI_indexKlogLikAICc∆AICcωi
Richness~Arg_ant + scale(Riparian_area) +
scale(Terrestrial_area)
1.3++ −0.14 0.16 5−1653400.000.192
Richness~Arg_ant + scale(Terrestrial_area)1.3++ 0.20 4−1663410.110.182
Richness~Arg_ant + scale(Matrix_area) +
scale(Terrestrial_area)
1.3++0.06 0.21 5−1663421.790.079
Richness~Arg_ant + scale(Terrestrial_area) +
scale(HEIDI_index)
1.3++ 0.21−0.0385−1663432.130.066
Richness~Arg_ant + scale(Matrix_area) +
scale(Riparian_area) + scale(Terrestrial_area)
1.4++0.03−0.13 0.17 6−1653432.190.064
Richness~Arg_ant + scale(Shrub_richness) +
scale(Terrestrial_area)
1.3++ −0.030.21 5−1663432.210.064
Richness~Arg_ant + scale(Riparian_area) +
scale(Shrub_richness) + scale(Terrestrial_area)
1.4++ −0.14−0.0250.17 6−1653432.260.062
Richness~Arg_ant + scale(Riparian_area) +
scale(Terrestrial_area) + scale(HEIDI_index)
1.3++ −0.14 0.16−0.236−1653432.270.062
Richness~scale(Riparian_area) +
scale(Terrestrial_area)
1.5+ −0.20 0.16 4−1683443.210.039
Richness~Arg_ant + scale(Riparian_area)1.4++ −0.23 4−1683443.330.036
Richness~Arg_ant + scale(Matrix_area) +
scale(Terrestrial_area) + scale(HEIDI_index)
1.3++0.06 0.22−0.046−1653443.800.029
Richness~Arg_ant + scale(Matrix_area) +
scale(Shrub_richness)) + scale(Terrestrial_area)
1.3++0.06 −0.040.23 6−1663443.900.027
Richness~Arg_ant + scale(Shrub_richness)) +
scale(Terrestrial_area) + scale(HEIDI_index)
1.3++ −0.020.21−0.036−1663454.430.021
Richness~Arg_ant + scale(Matrix_area) +
scale(Riparian_area) + scale(Shrub_richness) +
scale(Terrestrial_area)
1.4++0.04−0.13−0.030.18 7−1653454.490.020
Richness~Arg_ant + scale(Matrix_area) +
scale(Riparian_area)) + scale(Terrestrial_area) +
scale(HEIDI_index)
1.4++0.04−0.12 0.17−0.037−1653454.500.020
Richness~Arg_ant + scale(Riparian_area) +
scale(Shrub_richness)) + scale(Terrestrial_area) +
scale(HEIDI_index)
1.4++ −0.14−0.020.17−0.027−1653454.670.019
Richness~Arg_ant + scale(Matrix_area) +
scale(Riparian_area)) + scale(Terrestrial_area)
1.5++0.07−0.17 0.17 5−1673454.790.018
EI—ecological infrastructure; Arg_ant—ecological infrastructure; Matrix area—Distance to the closest EI patch; Riparian area—Area of riparian EI; Terrestrial area—Area of terrestrial EI; Shrub_richness—Shrub richness; HEIDI_index—HEIDI quality index. Please see Table 1 for a detailed description.

References

  1. Foley, J.A.; DeFries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K.; et al. Global Consequences of Land Use. Science 2005, 309, 570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Landis, D.A. Designing agricultural landscapes for biodiversity-based ecosystem services. Basic Appl. Ecol. 2016, 18, 1–2. [Google Scholar] [CrossRef] [Green Version]
  3. Rusch, A.; Chaplin-Kramerc, R.; Gardinere, M.M.; Hawrof, V.; Holland, J.; Landis, D.; Thiesi, C.; Tscharntke, T.; Weisserj, W.W.; Winqvistk, C.; et al. Agricultural landscape simplification reduces natural pest control: A quantitative synthesis. Agric. Ecosyst. Environ. 2016, 221, 198–204. [Google Scholar] [CrossRef] [Green Version]
  4. Holland, J.M.; Bianchi, F.J.; Entling, M.H.; Moonen, A.-C.; Smith, B.M.; Jeanneret, P. Structure, function and management of semi-natural habitats for conservation biological control: A review of European studies. Pest Manag. Sci. 2017, 72, 1638–1651. [Google Scholar] [CrossRef]
  5. Boller, E.F.; Hani, F.; Poehling, H.M. Ecological infrastructures. In Ideabook on Functional Biodiversity at the Farm Level; Boller, E.F., Hani, F., Poehling, H.M., Eds.; Swiss Centre for Agricultural Extension and Rural Development (LBL): Lindau, Switzerland, 2004; Volume XIV, p. 212. [Google Scholar]
  6. Silva, J.M.C.; Wheeler, E. Ecosystems as infrastructure. Perspect. Ecol. Conserv. 2017, 15, 32–35. [Google Scholar] [CrossRef]
  7. European Commission. The Post-2020 Common Agricultural Policy: Environmental Benefits and Simplification. Agriculture and Rural Developmen; European Commission: Brussels, Belgium, 2019. [Google Scholar]
  8. European Commission. List of Potential Practices That Eco-Schemes Could Support; European Commission: Brussels, Belgium, 2021. [Google Scholar]
  9. Rotcheés-Ribalta, R.; Ruas, S.; Ahmed, K.D.; Gormally, M.; Moran, J.; Stout, J.; White, B.; ÓhUallacháin, D. Assessment of semi-natural habitats and landscape features on Irish farmland: New insights to inform EU Common Agricultural Policy implementation. Ambio 2021, 50, 346–359. [Google Scholar] [CrossRef]
  10. Brosi, B.J.; Armsworth, P.R.; Daily, G.C. Optimal design of agricultural landscapes for pollination services. Conserv. Lett. 2008, 1, 27–36. [Google Scholar] [CrossRef]
  11. Fahrig, L.; Baudry, J.; Brotons, L.; Burel, F.G.; Crist, T.O.; Fuller, R.J.; Sirami, C.; Siriwardena, G.M.; Martin, J.L. Functional landscape heterogeneity and animal biodiversity in agricultural landscapes. Ecol. Lett. 2011, 14, 101–112. [Google Scholar] [CrossRef]
  12. Duflot, R.; Ernoult, A.; Aviron, S.; Fahrig, L.; Burel, F. Relative effects of landscape composition and configuration on multi-habitat gamma diversity in agricultural landscapes. Agric. Ecosyst. Environ. 2017, 241, 62–69. [Google Scholar] [CrossRef]
  13. Larrieu, L.; Gonin, P. Ĺlndice de biodiversité potentielle (IBP): Une méthode simple et rapide pour évaluer la biodiversité 798 potentielle des peuplements forestiers. Rev. For. Fr. 2008, 60, 727–748. [Google Scholar] [CrossRef]
  14. Fonseca, A.; Zina, V.; Duarte, G.; Aguiar, F.C.; Rodríguez-González, P.M.; Ferreira, M.T.; Fernandes, M.R. Riparian Ecological Infrastructures: Potential for biodiversity-related ecosystem services in Mediterranean human-dominated landscapes. Sustainability 2021, 13, 10508. [Google Scholar] [CrossRef]
  15. Corbacho, C.; Sánchez, J.M.; Costillo, E. Patterns of structural complexity and human disturbance of riparian vegetation in agricultural landscapes of a Mediterranean area. Agric. Ecosyst. Environ. 2003, 95, 495–507. [Google Scholar] [CrossRef]
  16. Fernandes, M.R.; Aguiar, F.C.; Ferreira, M.T. Assessing riparian vegetation structure and the influence of land use using landscape metrics and geostatistical tools. Landsc. Urban. Plan. 2011, 99, 166–177. [Google Scholar] [CrossRef] [Green Version]
  17. Naiman, R.J.; Decamps, H.; Pollock, M. The role of riparian corridors in maintaining regional biodiversity. Ecol. Appl. 1993, 3, 209–212. [Google Scholar] [CrossRef] [PubMed]
  18. Santos, A.; Fernandes, M.R.; Aguiar, F.C.; Branco, M.; Ferreira, M.T. Effects of riverine landscape changes on pollination services: A case study on the River Minho, Portugal. Ecol. Ind. 2018, 89, 656–666. [Google Scholar] [CrossRef]
  19. Riis, T.; Kelly-Quinn, M.; Aguiar, F.C.; Manolaki, P.; Bruno, D.; Bejarano, M.D.; Clerici, N.; Fernandes, M.R.; Franco, J.C.; Pettit, N.; et al. Global Overview of Ecosystem Services Provided by Riparian Vegetation. BioScience 2020, 70, 501–514. [Google Scholar] [CrossRef]
  20. Froidevaux, J.S.P.; Duarte, G.; Fonseca, A.; Zina, V.; Conde, S.; Ferreira, M.T.; Fernandes, M.R. The location and vegetation physiognomy of ecological infrastructures determine bat activity in Mediterranean floodplain landscapes. Agric. Ecosyst. Environ. 2022, 332, 107929. [Google Scholar] [CrossRef]
  21. Perfecto, I.; Armbrecht, I.; Philpott, S.M.; Soto-Pinto, L.; Dietsch, T.V. Shaded coffee and the stability of rainforest margins in northern Latin America. In Stability of Tropical Rainforest Margins. Environmental Science and Engineering; Tscharntke, T., Leuschner, C., Zeller, M., Guhardja, E., Bidin, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 227–263. [Google Scholar] [CrossRef]
  22. Wendt, C.F.; Nunes, A.; Verble, R.; Santini, G.; Boeiro, M.; Branquinho, C. Using a space-for-time approach to select the best biodiversity-based indicators to assess the effects of aridity on Mediterranean drylands. Ecol. Indic. 2020, 113, 106250. [Google Scholar] [CrossRef]
  23. Wendt, C.F.; Ceia-Hasse, A.; Nunes, A.; Verble, R.; Santini, G.; Boeiro, M.; Branquinho, C. Local environmental variables are key drivers of ant taxonomic and functional beta-diversity in a Mediterranean dryland. Sci. Rep. 2021, 11, 2292. [Google Scholar] [CrossRef] [PubMed]
  24. AntWeb. Version 8.78. California Academy of Science. Available online: https://www.antweb.org. (accessed on 11 August 2022).
  25. Hölldobler, B.; Wilson, E.O. The Ants; Harvard University Press: Cambridge, MA, USA, 1990; p. 746. [Google Scholar]
  26. Wilson, E.O.; Hölldobler, B. Eusociality: Origin and consequences. Proc. Natl. Acad. Sci. USA 2005, 102, 13367–13371. [Google Scholar] [CrossRef] [Green Version]
  27. Parr, C.L.; Bishop, T.R. The response of ants to climate change. Glob. Change Biol. 2022, 28, 3188–3205. [Google Scholar] [CrossRef] [PubMed]
  28. Sanders, D.; Van Veen, F.J.F. Ecosystem engineering and predation: The multi-trophic impact of two ant species. J. Anim. Ecol. 2011, 80, 569–576. [Google Scholar] [CrossRef] [PubMed]
  29. Wills, B.D.; Landis, D.A. The role of ants in north temperate grasslands: A review. Oecologia 2018, 186, 323–338. [Google Scholar] [CrossRef] [Green Version]
  30. Parker, J.; Kronauer, D.J.C. How ants shape biodiversity. Curr. Biol. 2021, 31, R1141–R1224. [Google Scholar] [CrossRef]
  31. Folgarait, P.J. Ant biodiversity and its relationship to ecosystem functioning: A review. Biodivers. Conserv. 1998, 7, 1221–1244. [Google Scholar] [CrossRef]
  32. Del Toro, I.; Ribbons, R.R.; Pelini, S.L. The little things than run the world revisited: A review of the ant-mediated ecosystem services and disservices (Hymenoptera: Formicidae). Myrmecol. News 2012, 17, 133–146. [Google Scholar]
  33. Elizalde, L.; Arbetman, M.; Arnan, X.; Eggleton, P.; Leal, I.R.; Lescano, M.N.; Saez, A.; Werenkraut, V.; Pirk, G.I. The ecosystem services provided by social insects: Traits, management tools and knowledge gaps. Biol. Rev. 2020, 95, 1418–1441. [Google Scholar] [CrossRef]
  34. Peck, S.L.; Mcquaid, B.; Campbell, C.L. Using ant species (Hymenoptera: Formicidae) as a biological indicator of agroecosystem condition. Environ. Entomol. 1998, 27, 1102–1110. [Google Scholar] [CrossRef] [Green Version]
  35. Dauber, J.; Wolters, V. Edge effects on ant community structure and species richness in an agricultural landscape. Biodiv. Conserv. 2004, 13, 901–915. [Google Scholar] [CrossRef]
  36. Ng, K.; Nowrouzi, S.; Staunton, K.M.; Barton, P.; Driscoll, D.A. Ant community responses to farmland use and revegetation in a fragmented agricultural landscape. Agric. Ecosyst. Environ. 2021, 311, 107316. [Google Scholar] [CrossRef]
  37. Dauber, J.; Hirsch, M.; Simmering, D.; Waldhardt, R.; Otte, A.; Wolters, V. Landscape structure as an indicator of biodiversity: Matrix effects on species richness. Agric. Ecosyst. Environ. 2003, 98, 321–329. [Google Scholar] [CrossRef]
  38. García-Martínez, M.A.; Valenzuela-González, J.E.; Escobar-Sarria, F.; Castaño-Meneses, G. The surrounding landscape influences the diversity of leaf-litter ants in riparian cloud forest remnants. PLoS ONE 2017, 12, e0172464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Frizzo, T.L.M.; Souza, L.M.; Sujii, E.R.; Togni, P.H.B. Ants provide biological control on tropical organic farms influenced by local and landscape factors. Biol. Control 2020, 151, 104378. [Google Scholar] [CrossRef]
  40. Kass, J.M.; Guénard, B.; Dudley, K.L.; Jenkins, C.N.; Azuma, F.; Fisher, B.L.; Parr, C.L.; Gibb, H.; Longino, J.T.; Ward, P.S.; et al. The global distribution of known and undiscovered ant biodiversity. Sci. Adv. 2022, 8, eabp9908. [Google Scholar] [CrossRef]
  41. Philpott, S.M.; Perfecto, I.; Armbrecht, I.; Parr, C.L. Ant diversity and function in disturbed and changing habitats. In Ant Ecology; Lach, L., Parr, C.L., Abbott, K., Eds.; Oxford University Press: Oxford, UK, 2010; pp. 137–157. [Google Scholar]
  42. European Commission; European Environmental Agency. The European Climate Adaptation Platform Climate—ADAPT. Available online: https://climate-adapt.eea.europa.eu/countries-regions/countries/portugal (accessed on 2 December 2021).
  43. Aguiar, F.C.; Ferreira, M.T.; Albuquerque, A.; Moreira, I. Alien and endemic flora at reference and non-reference sites in Mediterranean-type streams in Portugal. Aquat. Conserv. Mar. Freshw. Ecosyst. 2007, 17, 335–347. [Google Scholar] [CrossRef]
  44. Fonseca, A.; Duarte, G.; Zina, V.; Ferreira, M.T.; Fernandes, M.R. Ecological infrastructure Layers—Optimus Prime Project; Open Science Framework: Charlottesville, VA, USA, 2020. [Google Scholar] [CrossRef]
  45. Ferreira, M.T.; Aguiar, F.C.; Nogueira, C. Changes in riparian woods over space and time: Influence of environment and land 816 use. For. Ecol. Manag. 2005, 212, 145–159. [Google Scholar] [CrossRef]
  46. Wasser, L.; Chasmer, L.; Day, R.; Taylor, A. Quantifying land use effects on forested riparian buffer vegetation structure using 818 LiDAR data. Ecosphere 2015, 6, 1–17. [Google Scholar] [CrossRef]
  47. Gómez, C.; Espadaler, X. Seed dispersal curve of a Mediterranean myrmecochore: Influence of ant size and the distance to nests. Ecol. Res. 1998, 13, 347–354. [Google Scholar] [CrossRef]
  48. Gómez, C.; Espadaler, X. An update of the world survey of myrmecochorous dispersal distances. Ecography 1998, 36, 1193–1201. [Google Scholar] [CrossRef]
  49. Fonseca, A.; Duarte, G.; Zina, V.; Ferreira, M.T.; Fernandes, M.R. Explanatory Variables Database—Optimus Prime Project. Open Science Framework: Charlottesville, VA, USA, 2021. [Google Scholar] [CrossRef]
  50. Cros, S.; Cerdá, X.; Retana, J. Spatial and temporal variations in the activity patterns of Mediterranean ant communities. Ecoscience 1997, 4, 269–278. [Google Scholar] [CrossRef] [Green Version]
  51. Retana, J.; Cerdá, X. Patterns of diversity and composition of Mediterranean ground ant communities tracking spatial and temporal variability in the thermal environment. Oecologia 2000, 123, 436–444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Gómez, C.; Casellas, D.; Oliveras, J.; Bas, J.M. Structure of ground-foraging ant assemblages in relation to land-use change in the northwestern Mediterranean region. Biodivers. Conserv. 2003, 12, 2135–2146. [Google Scholar] [CrossRef]
  53. Angulo, E.; Boulay, R.; Ruano, F.; Tinaut, A.; Cerda, X. Anthropogenic impacts in protected areas: Assessing the efficiency of conservation efforts using Mediterranean ant communities. PeerJ 2016, 4, e2773. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Majer, J.D. The use of pitfall traps for sampling ants—A critique. Mem. Mus. Vic. 1997, 56, 323–329. [Google Scholar] [CrossRef]
  55. Parr, C.L.; Chown, S.L. Inventory and bioindicator sampling: Testing pitfall and Winkler methods with ants in a South African savanna. J. Insect Conserv. 2001, 5, 27–36. [Google Scholar] [CrossRef]
  56. Underwood, E.C.; Fisher, B.L. The role of ants in conservation monitoring: If, when, and how. Biol. Conserv. 2006, 132, 166–182. [Google Scholar] [CrossRef]
  57. Tista, M.; Fiedler, K. How to evaluate and reduce sampling effort for ants. J. Insect Conserv. 2011, 15, 547–559. [Google Scholar] [CrossRef]
  58. Samways, M.J. Community structure of ants (Hymenoptera: Formicidae) in a series of habitats associated with citrus. J. Appl. Ecol. 1983, 20, 833–847. [Google Scholar] [CrossRef]
  59. Vele, A.; Holusa, J.; Frouz, J. Sampling for ants in different-aged spruce forests: A comparison of methods. Eur. J. Soil Biol. 2009, 45, 301–305. [Google Scholar] [CrossRef]
  60. Schmidt, F.A.; Ribas, C.R.; Schoereder, J.H. How predictable is the response of ant assemblages to natural forest recovery? Implications for their use as bioindicators. Ecol. Indic. 2013, 24, 158–166. [Google Scholar] [CrossRef]
  61. Johnson, J.T.; Adkins, J.K. Canopy vegetation influences ant (Hymenoptera: Formicidae) communities in headwater stream riparian zones of central Appalachia. J. Insect Sci. 2014, 14, 237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Sheikh, A.H.; Ganaie, G.A.; Thomas, M.; Bhandari, R.; Rather, Y.A. Ant pitfall trap sampling: An overview. J. Soil Biol. 2018, 45, 301–305. [Google Scholar] [CrossRef]
  63. Collingwood, C.A.; Prince, A. A guide to ants of continental Portugal. Bol. Soc. Por. Entomol. 1998, 5, 1–49. [Google Scholar]
  64. Gómez, K.; Espadaler, X. Hormigas Ibéricas. Available online: http://www.hormigas.org (accessed on 3 April 2020).
  65. Lebas, C.; Galkowski, C.; Blatrix, R.; Wegnez, P. Guia de Campo de las Hormigas de Europa Ocidental; Omega: Barcelona, Spain, 2017; p. 416. [Google Scholar]
  66. R Core Team 2020. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 1 February 2022).
  67. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  68. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R Package Version 2.5–7. 2020. Available online: https://CRAN.R-project.org/package=vegan (accessed on 15 June 2022).
  69. Brooks, M.E.; Kristensen, K.; van Benthem, K.J.; Magnusson, A.; Berg, C.W.; Nielsen, A.; Skaug, H.J.; Maechler, M.; Bolker, B.M. glmmTMB Balances Speed and Flexibility Among Packages for Zero-inflated Generalized Linear Mixed Modeling. R J. 2017, 9, 378–400. [Google Scholar] [CrossRef] [Green Version]
  70. Burnham, K.P.; Anderson, D.R. Information and Likelihood Theory: A Basis for Model Selection and Inference. In Model Selection and Multimodel Inference, Burnham, K.P., Anderson, D.R., Eds.; Springer: New York, NY, USA, 2002; pp. 49–97. [Google Scholar]
  71. Hartig, F. DHARMa: Residual Diagnostics for Hierarchical (Multi-Level/Mixed) Regression Models. R Package Version 0.4.1. 2021. Available online: https://CRAN.R-project.org/package=DHARMa (accessed on 15 June 2022).
  72. Lüdecke, D.; Ben-Shachar, M.S.; Patil, I.; Waggoner, P.; Makowski, D. performance: An R Package for Assessment, Comparison and Testing of Statistical Models. J. Open Source Softw. 2021, 6, 3139. [Google Scholar] [CrossRef]
  73. Bartoń, K. MuMIn: Multi-Model Inference. R Package Version 1.43.17. 2020. Available online: https://CRAN.R-project.org/package=MuMIn (accessed on 15 June 2022).
  74. Burnham, K.P.; Anderson, D.R. Formal Inference from More Than One Model: Multimodel Inference (Mmi). In Model Selection and Multimodel Inference; Burnham, K.P., Anderson, D.R., Eds.; Springer: New York, NY, USA, 2002; pp. 149–205. [Google Scholar]
  75. Nooten, S.S.; Schultheiss, P.; Rowe, R.C.; Facey, S.L.; Cook, J.M. Habitat complexity affects functional traits and diversity of ant assemblages in urban green spaces (Hymenoptera: Formicidae). Myrmecol. News 2019, 29, 67–77. [Google Scholar]
  76. MacArthur, R.H.; MacArthur, J.W. On bird species diversity. Ecology 1961, 42, 594–598. [Google Scholar] [CrossRef]
  77. MacArthur, R.H.; Wilson, E.O. The Theory of Island Biogeography; Princeton University Press: Princeton, NJ, USA, 1967. [Google Scholar]
  78. Kreider, J.J.; Chen, T.-W.; Hartke, T.R.; Buchori, D.; Hidayat, P.; Nazarreta, R.; Scheu, S.; Drescher, J. Rainforest conversion to monocultures favors generalist ants with large colonies. Ecosphere 2012, 12, e03717. [Google Scholar] [CrossRef]
  79. Andersen, A.N. Responses of ant communities to disturbance: Five principles for understanding the disturbance dynamics of a globally dominant faunal group. J. Anim. Ecol. 2019, 88, 350–362. [Google Scholar] [CrossRef]
  80. Wetterer, J.K. Worldwide Spread of the Moorish sneaking Ant, Cardiocondyla mauritanica (Hymenoptera: Formicidae). Sociobiology 2012, 59, 3. [Google Scholar]
  81. Andersen, A.N. A Classification of Australian Ant Communities, Based on Functional Groups Which Parallel Plant Life-Forms in Relation to Stress and Disturbance. J. Biogeogr. 1995, 22, 15–29. [Google Scholar] [CrossRef]
  82. Hoffmann, B.D.; Andersen, A.N. Responses of ants to disturbance in Australia, with particular reference to functional groups. Austral. Ecol. 2003, 28, 444–464. [Google Scholar] [CrossRef]
  83. Morris, T.I.; Symondson, W.O.C.; Kidd, N.A.C.; Jervis, M.A.; Campos, M. Are ants significant predators of the olive moth Prays oleae? Crop Prot. 1998, 17, 365–366. [Google Scholar] [CrossRef]
  84. Campolo, O.; Palmeri, V.; Malacrinò, A.; Laudani, F.; Castracani, C.; Mori, A.; Grasso, D.A. Interaction between ants and the Mediterranean fruit fly: New insights for biological control. Biol. Control 2015, 90, 120–127. [Google Scholar] [CrossRef]
  85. Martínez-Núñez, C.; Rey, P.J.; Salido, T.; Manzaneda, A.J.; Camacho, F.M.; Isla, J. Ant community potential for pest control in olive groves: Management and landscape effects. Agric. Ecosyst. Environ. 2021, 305, 107185. [Google Scholar] [CrossRef]
  86. Cerdá, X.; Dejean, A. Predation by ants on arthropods and other Animals. In Predation in the Hymenoptera: An Evolutionary Perspective; Polidori, C., Ed.; Transworld Research Network: Kerala, India, 2011; pp. 39–78. [Google Scholar]
  87. Holway, D.A.; Cameron, E.K. The importance of scavenging in ant invasions. Curr. Opin. Insect Sci. 2021, 46, 39–42. [Google Scholar] [CrossRef]
  88. Zumeaga, H.; Azcárate, F.M.; Concepción, E.D.; Hevia, V.; Díaz, M. Landscape and agri-environmental scheme effects on ant communities in cereal croplands of central Spain. Agric. Ecosyst. Environ. 2021, 312, 107345. [Google Scholar] [CrossRef]
  89. Meissle, M.; Mouron, P.; Musa, T.; Bigler, F.; Pons, X.; Vasileiadis, V.P.; Otto, S.; Antichi, D.; Kiss, J.; Pálinkás, Z.; et al. Pests, pesticide use and alternative options in European maize production: Current status and future prospects. J. Appl. Entomol. 2010, 134, 357–375. [Google Scholar] [CrossRef]
  90. Wetterer, J.K.; Wild, A.L.; Suarez, A.V.; Roura-Pascual, N.; Espadaler, X. Worldwide spread of the Argentine ant, Linepithema humile (Hymenoptera: Formicidae). Myrmecol. News 2009, 12, 187–194. [Google Scholar]
  91. Lach, L. Argentine ants displace floral arthropods in a bio-diversity hotspot. Divers. Distrib. 2008, 14, 281–290. [Google Scholar] [CrossRef]
  92. Gómez, C.; Oliveras, J. Can the Argentine ant (Linepithema humile Mayr) replace native ants in myrmecochory? Acta Oecologica 2003, 24, 47–53. [Google Scholar] [CrossRef]
  93. Mgocheki, N.; Addison, P. Interference of ants (Hymenoptera: Formicidae) with biological control of the vine mealybug Planococcus ficus (Signoret) (Hemiptera: Pseudococcidae). Biol. Control 2009, 49, 180–185. [Google Scholar] [CrossRef]
  94. Oliveras, J.; Bas, J.M.; Casellas, D.; Gómez, C. Numerical dominance of the Argentine ant vs native ants and consequences on soil resource searching in Mediterranean cork-oak forests (Hymenoptera: Formicidae). Sociobiology 2005, 45, 1–16. [Google Scholar]
  95. Kolay, S.; Boulay, R.; D’Ettorre, P. Regulation of ants foraging: A review of the role of information use and personality. Front. Psychol. 2020, 11, 734. [Google Scholar] [CrossRef]
  96. Walters, A.C. Invasion of the Argentine ants (Hymenoptera, Formicidae) in South Australia: Impacts on community composition and abundance of invertebrates in urban parklands. Austral Ecol. 2006, 31, 367–376. [Google Scholar] [CrossRef]
  97. Roura-Pascual, N.; Bas, J.M.; Hui, C. The spread of the Argentine ant: Environmental determinants and impacts on native ant communities. Biol. Invasions 2010, 12, 2399–2412. [Google Scholar] [CrossRef]
  98. Rowles, A.D.; O’Dowd, D.J. Impacts of the invasive Argentine ant on native ants and other invertebrates in coastal scrub in south-eastern Australia. Austral Ecol. 2009, 34, 239–248. [Google Scholar] [CrossRef]
  99. Zina, V.; Branco, M.; Franco, J.C. Impact of the invasive Argentine ant in citrus agroecosystems: Effects on the diversity and frequency of native ant species foraging on tree canopy. Insects 2020, 11, 785. [Google Scholar] [CrossRef]
  100. Aguiar, F.C.; Ferreira, M.T. Plant invasions in the rivers of the Iberian Peninsula, south-western Europe: A review. Int. J. Plant Biol. 2013, 147, 4, 1107–1119. [Google Scholar] [CrossRef]
  101. Armbrecht, I.; Perfecto, I. Litter-twig dwelling ant species richness and predation potential within a forest fragment and neighboring coffee plantations of contrasting habitat quality in Mexico. Agric. Ecosyst. Environ. 2003, 97, 107–115. [Google Scholar] [CrossRef]
  102. Lee, J.C.; Menalled, F.D.; Landis, D.A. Refuge habitats modify impact of insecticide disturbance on carabid beetle communities. J. Appl. Ecol. 2001, 38, 472–483. [Google Scholar] [CrossRef]
  103. Crist, T.O. Biodiversity, species interactions, and functional roles of ants (Hymenoptera: Formicidae) in fragmented landscapes: A review. Myrmecol. News 2009, 12, 3–13. [Google Scholar]
  104. Santos, M.J.; Rosalino, L.M.; Matos, H.M.; Santos-Reis, M. Riparian ecosystem configuration influences mesocarnivores presence in Mediterranean landscapes. Eur. J. Wildl. Res. 2016, 62, 251–261. [Google Scholar] [CrossRef]
  105. Forio, M.A.E.; De Troyer, N.; Lock, K.; Witing, F.; Baert, L.; Saeyer, N.D.; Rîșnoveanu, G.; Popescu, C.; Burdon, F.J.; Kupilas, B.; et al. Small patches of riparian woody vegetation enhance biodiversity of invertebrates. Water 2020, 12, 3070. [Google Scholar] [CrossRef]
  106. Popescu, C.; Oprina-Pavelescu, M.; Dinu, V.; Cazacu, C.; Burdon, F.J.; Forio, M.A.E.; Kupilas, B.; Friberg, N.; Goethals, P.; McKie, B.G.; et al. Riparian Vegetation Structure Influences Terrestrial Invertebrate Communities in an Agricultural Landscape. Water 2021, 13, 188. [Google Scholar] [CrossRef]
  107. García-Martínez, M.A.; Escobar-Sarria, F.; López-Barrera, F.; Castaño-Meneses, G.; Valenzuela-González, J.E. Value of riparian vegetation remnants for leaf-litter ants (Hymenoptera: Formicidae) in a human-dominated landscape in Central Veracruz, Mexico. Environ. Entomol. 2015, 44, 1488–1497. [Google Scholar] [CrossRef]
  108. Buczkowski, G.; Richmond, D.S. The effect of urbanization on ant abundance and diversity: A temporal examination of factors affecting biodiversity. PLoS ONE 2012, 7, e41729. [Google Scholar] [CrossRef] [Green Version]
  109. Rocha, E.A.; Fellowes, M.D.E. Urbanization alters ecological interactions: Ant mutualists increase and specialist insect predators decrease on an urban gradient. Sci. Rep. 2020, 10, 6406. [Google Scholar] [CrossRef] [Green Version]
  110. Guénard, B.; Cardinal-De Casas, A.; Dunn, R.R. High diversity in an urban habitat: Are some animal assemblages resilient to long-term anthropogenic change? Urban Ecosyst. 2014, 18, 449–463. [Google Scholar] [CrossRef]
  111. Menke, S.B.; Guénard, B.; Sexton, J.O.; Weiser, M.D.; Dunn, R.R.; Silverman, J. Urban areas may serve as habitat and corridors for dry-adapted, heat tolerant species; an example from ants. Urban Ecosyst. 2011, 14, 135–163. [Google Scholar] [CrossRef]
  112. Crenna, E.; Sinkko, T.; Sala, S. Biodiversity impacts due to food consumption in Europe. J. Clean. Prod. 2019, 227, 378–391. [Google Scholar] [CrossRef] [PubMed]
  113. GPP—Gabinete de Planeamento, Políticas e Administração Geral. Análise Sectorial CEREAIS; GPP—Gabinete de Planeamento: Lisbon, Portugal, 2020. [Google Scholar]
  114. Bambaradeniya, C.N.B.; Amerasinghe, F.P. Biodiversity Associated with Rice Field Agroecosystem in Asian Countries: A Brief Review; Working Paper 63; International Water Management Institute: Colombo, Sri Lanka, 2003. [Google Scholar]
  115. Edirisinghe, J.P.; Bambaradeniya, C.N.B. Rice fields: An ecosystem rich in biodiversity. J. Natl. Sci. Found. Sri Lanka 2016, 34, 57–59. [Google Scholar] [CrossRef]
  116. De la Mora, A.; Murnen, C.J.; Philpott, S.M. Local and landscape drivers of biodiversity of four groups of ants in coffee landscapes. Biodivers Conserv. 2013, 22, 871–888. [Google Scholar] [CrossRef]
  117. Tscharntke, T.; Klein, A.M.; Kruess, A.; Steffan-Dewenter, I.; Thies, C. Landscape perspectives on agricultural intensification and biodiversity—Ecosystem service management. Ecol. Lett. 2005, 8, 857–874. [Google Scholar] [CrossRef]
  118. Ward, P.S. Distribution of the introduced Argentine ant (Iridomyrmex humilis) in natural habitats of the lower Sacramento Valley and its effects on the indigenous ant fauna. Hilgardia 1987, 55, 1–16. [Google Scholar] [CrossRef] [Green Version]
  119. Roura-Pascual, N.; Brotons, L.; Peterson, A.T.; Thuiller, W. Consensual predictions of potential distributional areas for invasive species: A case study of Argentine ants in the Iberian Peninsula. Biol. Invas. 2009, 11, 1017–1031. [Google Scholar] [CrossRef]
  120. Franco, J.C.; Branco, M.; Conde, S.; Garcia, A.; Fernandes, M.R.; Lima Santos, J.; Messina, T.; Duarte, G.; Fonseca, A.; Zina, V.; et al. Ecological Infrastructures May Enhance Lepidopterans Predation in Irrigated Mediterranean Farmland, Depending on Their Typology and the Predator Guild. Sustainability 2022, 14, 3874. [Google Scholar] [CrossRef]
Figure 1. Geographic location of study areas and sampling sites. (a) Tagus and (b) Sorraia valleys.
Figure 1. Geographic location of study areas and sampling sites. (a) Tagus and (b) Sorraia valleys.
Agronomy 12 02690 g001
Figure 2. Boxplots represent the ant species richness among groups. (a) Ecological infrastructures (EI) vs. agricultural matrices; (b) woody EI vs. herbaceous EI vs. agricultural matrices; (c) terrestrial EI vs. riparian EI vs. agricultural matrices. Different letters show significant differences between and among groups (Tukey’s HSD test).
Figure 2. Boxplots represent the ant species richness among groups. (a) Ecological infrastructures (EI) vs. agricultural matrices; (b) woody EI vs. herbaceous EI vs. agricultural matrices; (c) terrestrial EI vs. riparian EI vs. agricultural matrices. Different letters show significant differences between and among groups (Tukey’s HSD test).
Agronomy 12 02690 g002
Figure 3. Ordination plots of species composition of ant communities in Ecological infrastructures (EI) and the agricultural matrix assemblages based on Bray–Curtis dissimilarities: (a) EI vs. agricultural matrices; (b) woody EI vs. herbaceous EI vs. agricultural matrices; (c) terrestrial EI vs. riparian EI vs. agricultural matrices Legend: EI—ecological infrastructures (green); MATRIX—agricultural matrix (black); W—woody EI (dark green); H—herbaceous EI (light green); T—terrestrial EI (green); R—riparian EI (blue).
Figure 3. Ordination plots of species composition of ant communities in Ecological infrastructures (EI) and the agricultural matrix assemblages based on Bray–Curtis dissimilarities: (a) EI vs. agricultural matrices; (b) woody EI vs. herbaceous EI vs. agricultural matrices; (c) terrestrial EI vs. riparian EI vs. agricultural matrices Legend: EI—ecological infrastructures (green); MATRIX—agricultural matrix (black); W—woody EI (dark green); H—herbaceous EI (light green); T—terrestrial EI (green); R—riparian EI (blue).
Agronomy 12 02690 g003
Figure 4. Frequency distribution of ant species occurrence (in percentage) as a function of the sites sampled in the ecological infrastructures (EI) (left blue-green bars) and the agricultural matrix (right gray bars).
Figure 4. Frequency distribution of ant species occurrence (in percentage) as a function of the sites sampled in the ecological infrastructures (EI) (left blue-green bars) and the agricultural matrix (right gray bars).
Agronomy 12 02690 g004
Figure 5. Predicted responses of the ant species richness within the agricultural matrix in relation to the distance to the nearest ecological infrastructure, (a) overall and (b) considering the crop type. Model predictions from GLMMs and associated 95% confidence intervals are represented by the solid lines and shaded areas, respectively.
Figure 5. Predicted responses of the ant species richness within the agricultural matrix in relation to the distance to the nearest ecological infrastructure, (a) overall and (b) considering the crop type. Model predictions from GLMMs and associated 95% confidence intervals are represented by the solid lines and shaded areas, respectively.
Agronomy 12 02690 g005
Figure 6. Predicted responses of the ant species richness in relation to the riparian and terrestrial ecological infrastructure (EI) proportional area in the 200 m buffer, overall (a,c) and considering the Argentine ant presence (b,d). Model predictions from GLMMs and associated 95% confidence intervals are represented by the solid lines and shaded areas, respectively.
Figure 6. Predicted responses of the ant species richness in relation to the riparian and terrestrial ecological infrastructure (EI) proportional area in the 200 m buffer, overall (a,c) and considering the Argentine ant presence (b,d). Model predictions from GLMMs and associated 95% confidence intervals are represented by the solid lines and shaded areas, respectively.
Agronomy 12 02690 g006
Table 1. Description of the variables used in the study.
Table 1. Description of the variables used in the study.
Applied To:Variable NameAbbreviationCategoryType, Units, and RangeDescription
Agricultural matrix dataDistance to the closest EI patchDist_EIProximityContinuous:
meters; [0, ∞]
Distance to the closest remnant habitat patch–EI edge
Distance to urban areaDist_urbanProximityContinuous:
meters; [0, ∞]
Distance to the closest urban area using Level 4 of COS 2018 layer (i.e., the Portuguese Land use and Occupancy Charter of 2018, www.dgterritorio.pt (accessed on 26 March 2021))
Distance to riverDist_riverProximityContinuous:
meters; [0, ∞]
Distance to the closest watercourse defined by the adapted HIDCOD layer (i.e., the Portuguese waterline layer based on the Digital Elevation Model of 25 m)
Crop typeCrop_typeHabitat qualityNominal:
maize field, rice paddy, others (mixed types)
Crop type characterization within the agricultural matrix
Ecological infrastructure (EI) dataAgricultural landAgricultural_matrixArea/densityContinuous:
ha; [0, ∞]
Area of the agricultural matrix in a 200 m buffer contained within the study area
Area of riparian EIRiparian_EI_areaArea/densityContinuous:
ha; [0, ∞]
Sum of the areas of riparian EI; basic statistics of the spatial configuration
Area of terrestrial EITerrestrial_EI_areaArea/densityContinuous:
ha; [0, ∞]
Sum of the areas of terrestrial EI; basic statistics of the spatial configuration.
Shrub richnessShrub_richnessHabitat qualityContinuous:
none; [0, ∞]
Number of shrub plant species in the EI patches
HEIDI quality indexHEIDI_indexHabitat qualityContinuous:
none; [0, ∞]
HEIDI value for short distance dispersers in the EI patches sensu Fonseca et al. [14]
Argentine ant occurrenceArgentine_antHabitat qualityNominal Bolean:
presence, absence
Argentine ant occurrence in the EI patches
Table 2. Summary of the GLMM best-fitted model compared with the Null and Full models, to test Hypothesis 3. Significant p-values (p < 0.05) are displayed in bold.
Table 2. Summary of the GLMM best-fitted model compared with the Null and Full models, to test Hypothesis 3. Significant p-values (p < 0.05) are displayed in bold.
Null ModelBest-Fitted ModelFull Model
PredictorsIncidence Rate Ratio
(95% Confidence Interval)
Statisticp-ValueIncidence Rate Ratio
(95% Confidence Interval)
Statisticp-ValueIncidence Rate Ratio
(95% Confidence Interval)
Statisticp-Value
(Intercept)2.17
(1.70–2.78)
6.16<0.0013.23
(2.35–4.44)
7.23<0.0013.16
(2.24–4.47)
6.52<0.001
Crop [Maize_field] 0.39
(0.24–0.65)
–3.68<0.0010.41
(0.23–0.71)
–3.170.002
Crop [Other_mixed_types] 0.69
(0.40–1.18)
–1.370.1720.69
(0.40–1.20)
–1.310.191
Dist_EI 0.68
(0.50–0.92)
–2.470.0130.68
(0.49–0.94)
–2.330.020
Dist_urban 0.98
(0.76–1.26)
–0.190.847
Dist_river 0.97
(0.73–1.29)
–0.220.828
Table 3. Summary of the GLMMs model with the lowest AIC, compared with the Null and Full models, to test Hypothesis 4. Significant p-values (p < 0.05) are displayed in bold.
Table 3. Summary of the GLMMs model with the lowest AIC, compared with the Null and Full models, to test Hypothesis 4. Significant p-values (p < 0.05) are displayed in bold.
Null modelBest-fitted modelFull model
PredictorsIncidence Rate Ratio
(95% Confidence Interval)
Statisticp-ValueIncidence Rate Ratio
(95% Confidence Interval)
Statisticp-ValueIncidence Rate Ratio
(95% Confidence Interval)
Statisticp-Value
(Intercept)4.79
(4.09–5.61)
19.41<0.0013.85
(3.13–4.74)
12.72<0.0013.88
(3.15–4.78)
12.71<0.001
Argentine ant [absence] 1.41
(1.06–1.87)
2.370.0181.38
(1.03–1.85)
2.180.030
Riparian area 0.87
(0.73–1.04)
–1.530.1270.88
(0.73–1.07)
–1.290.196
Terrestrial area 1.17
(1.04–1.33)
2.540.0111.20
(1.04–1.38)
2.520.012
Matrix area 1.04
(0.90–1.21)
0.510.613
Shrub richness 0.98
(0.84–1.13)
–0.300.763
HEIDI quality index 0.98
(0.85–1.13)
–0.280.780
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zina, V.; Fonseca, A.; Duarte, G.; Conde, S.; Fernandes, M.R.; Ferreira, M.T.; Franco, J.C. Ant Diversity Is Enhanced by Ecological Infrastructures in Agroecosystems: A Case Study in Irrigated Mediterranean Farmland. Agronomy 2022, 12, 2690. https://doi.org/10.3390/agronomy12112690

AMA Style

Zina V, Fonseca A, Duarte G, Conde S, Fernandes MR, Ferreira MT, Franco JC. Ant Diversity Is Enhanced by Ecological Infrastructures in Agroecosystems: A Case Study in Irrigated Mediterranean Farmland. Agronomy. 2022; 12(11):2690. https://doi.org/10.3390/agronomy12112690

Chicago/Turabian Style

Zina, Vera, André Fonseca, Gonçalo Duarte, Sofia Conde, Maria Rosário Fernandes, Maria Teresa Ferreira, and José Carlos Franco. 2022. "Ant Diversity Is Enhanced by Ecological Infrastructures in Agroecosystems: A Case Study in Irrigated Mediterranean Farmland" Agronomy 12, no. 11: 2690. https://doi.org/10.3390/agronomy12112690

APA Style

Zina, V., Fonseca, A., Duarte, G., Conde, S., Fernandes, M. R., Ferreira, M. T., & Franco, J. C. (2022). Ant Diversity Is Enhanced by Ecological Infrastructures in Agroecosystems: A Case Study in Irrigated Mediterranean Farmland. Agronomy, 12(11), 2690. https://doi.org/10.3390/agronomy12112690

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop