Content area

Abstract

Simple Summary

The agricultural intensification over the last 80 years has led to the creation of large-scale crop fields and the loss of ecological elements, affecting natural communities. The implementation of measures such as floral margins allows for the creation of refuges for insect communities. The present long-term study (2013–2022) demonstrates that the sustained implementation of floral margins can effectively protect insect communities in intensive agricultural areas, highlighting their importance as a tool for fostering insect biodiversity.

Full text

Turn on search term navigation

1. Introduction

For decades, the intensification of agriculture has been imposing increasing pressure on biodiversity in agroecosystems. The range and abundance of thousands of plant and animal species have been in serious decline [1,2,3]. Consequently, ecosystems services have been decreasing over time [3,4,5], ultimately resulting in unsustainable agricultural production and even soil degradation [6,7,8,9].

Today, agriculture is practiced on roughly 50% of the usable land of planet Earth, making it one of humanity’s large impacts on the environment [10]. It is evident that effectively safeguarding our environment and preserving biodiversity is nearly impossible without incorporating agricultural landscapes into conservation efforts. Supporting and promoting biodiversity within agricultural landscapes is crucial both for the conservation of ecosystems and related services such as pollination, predation, or carbon sequestration. A well-established method proposed to protect biodiversity in agricultural systems is sowing crop field margins with wildflower mixes [11,12,13,14,15,16,17].

The use of margins, whether natural or implemented, based on seed mixtures of autochthonous species, appears to function as ecological corridors, linking isolated habitat patches and reducing landscape fragmentation [18,19]. Moreover, increasing the abundance of wildflowers, insects, and birds has been highlighted as an important way of promoting ecosystem services and supporting biodiversity conservation [20,21,22,23,24,25,26].

However, most studies on agroecosystems are based on short-term observations (1–3 years) [27,28,29] and are mainly focused on bees or bumblebees [18,30,31,32,33,34,35,36,37,38], with only a few considering other groups of insects such as beetles, butterflies, or hoverflies [28,35,39,40,41].

This work assumes that biodiversity can be assessed by measuring the abundance of insects and their presence in different environments in the long term, evaluating whether changes induced by floral margins are permanent over time. This leads to the testing of two hypotheses. First, there is a benefit of integrating floral margins to protect biodiversity. Second, the use of floral margins improves biodiversity over time. These hypotheses were tested on an intensive vegetable farm in Spain.

2. Materials and Methods

2.1. Areas of Study

The study was carried out on one highly productive Spanish vegetable farm located in Águilas (Murcia; 37°25′01.6″ N 1°36′13.7″ W) (Figure 1). The location area has a semi-arid Mediterranean climate [42] with hot summers (27.7 °C) and mild winters (13.7 °C) and with an average annual rainfall of 303 mm. Appendix A summarises the annual rainfall and temperature data recorded in the study area.

During the sampling period, the crops were rotated successively, leaving a fallow period between July and October. Crop rotation included: beet (Amaranthaceae; Beta vulgaris L.), celery (Apiacea; Apium graveolens L.), lettuce “romana” (Asteraceae; Lactuca sativa L. var. longifolia), lettuce “iceberg” (Asteraceae; Lactuca sativa L. var. capitata), lettuce “mini romana” (Asteraceae; Lactuca sativa L.), lettuce “baby gem” (Asteraceae; Lactuca sativa L.), onion (Amaryllidaceae; Allium cepa L.), and triticale (Poaceae; Triticosecale Wittm. ex A. Camus). All crops were planted in a design where the planting distance was 26 cm between the rows and 20 cm between the plants of a same row, except for onion, where the distance was 15 cm between rows and 15 cm between plants, and triticale, where the distance was 15 cm between rows and 3–6 cm between plants. The field size was 7.5 ha.

During the study, the growers stuck to their preferred agricultural practices, such as tillage, sowing, and fertilisation. Moreover, they continued with their same phytosanitary treatments as before, applying the appropriate products according to pest and disease thresholds. Any management measures were confined to the crop to avoid interference with the multifunctional margin.

2.2. Floral Margins and Plant Mixture Selection

The selection of plant species was based on several fundamental criteria such as the strict use of native species, ensuring a smooth climatic adaptation; not become a potential weed for the crop; featuring easy maintenance and capacity for self-sowing, as well as staggered flowering phenologies; and finally, being attractive for pollinators and natural enemies.

We established a floral margin using a herbaceous mixture consisting of Borago officinalis L. (7%) (borage; Fam. Boraginaceae), Calendula officinalis L. (17.5%) (pot marigold; Fam. Asteraceae), Coriandrum sativum L. (10%) (coriander; Fam. Apiaceae), Diplotaxis catholica (L.) DC. (5%) (wall-rocket; Fam. Brassicaceae), Echium vulgare L. (5%) (viper’s bugloss; Fam. Boraginaceae), Lobularia maritima (L.) Desv. (5%) (sweet alyssum; Fam. Brassicaceae), Melilotus officinalis (L.) Pall. (12.5%) (sweet yellow clover; Fam. Fabaceae), Nigella damascena (L.) (5%) (love-in-a-mist; Fam. Ranunculaceae), Phlomis purpurea L. (3%) (Jerusalem sage; Fam. Lamiaceae), Salvia verbenaca L. (10%) (wild clary; Fam. Lamiaceae), Silene vulgaris (Moench) Garcke (10%) (bladder campion; Fam. Caryophyllaceae), and Vicia sativa L. (10%) (common vetch; Fam. Fabaceae). The floral field margin of 3 m width × 300 m length was sown next to the crop area at 3 m from the field to facilitate daily work in the crop. Sowing took place using a seed electric drill with air distribution (APV 100 pneumatic, APV Technische Produkte GmbH, Hötzelsdorf, Austria) after the soil had been prepared with a flail mower. The seeds were covered using a rake. The seed sowing rate applied was 15 kg/ha. The field margin was mowed in autumn and then left to regrow the following season. During the first two years, supplementary seeds were added annually in March and April to ensure a consistent plant emergence by following this planting schedule; irrigation or watering was not necessary as the rainfall provided favourable growth conditions, as these species are adapted to the climatic conditions of the area.

2.3. Experimental Design and Sampling

To investigate the dynamics of effects of floral margin on insect biodiversity, the experiment was conducted over a period of 10 years (2013–2022). Insect abundance was assessed visually (flower observation) and by using sweeping nets (observed and captured specimens were merged to perform the corresponding analyses). The samplings were carried out one day per month between March and July by moving in a zigzag along 4 fixed transects of 50 × 2 m during 15 min per line and 4 times per day to avoid the light and temperature gradient and obtain a more representative sample. Replication became unfeasible due to the impossibility of locating uniform fields as constant crop rotations were dictated by market demands.

The collected specimens were introduced in a bottle with small amount of cyanide to keep them intact and to avoid discoloration. All specimens were identified to the genus or species level using appropriate entomological literature (see [43,44,45,46,47,48,49,50,51,52,53,54,55,56,57]). Specimens are deposited in the entomological collection of the National Museum of Natural Sciences (Madrid, Spain; MNCN).

2.4. Statistical Data Analysis

We implemented an approach based on count data regression modelling to study the temporal dynamics of the number of species and number of insects as a measure of their diversity. To accomplish that, we initially perform an exploratory data analysis of the diversity indexes, Shannon’s H, species richness, and Pielou’s evenness, comparing them among orders and across years to study their behaviour and detect patterns. We secondly characterise the temporal autocorrelation structure of the number of species and insects by building the correlation matrix of the yearly numbers to identify temporal autoregressive effects. We then propose a generalised linear mixed model (GLMM) for count data to describe the number of species and the number of insects. Thus, the model assumes that the number of species or insects follows either a Poisson distribution when its conditional mean and conditional variance are equal or a negative binomial when its conditional variance is greater than its conditional mean (overdispersion). The model is specified for the number of species and for the number of insects in the Equations (1) and (2), respectively.

(1)μi,t=expα+βt+γi+γit+λyik,t1

(2)μij,t=expα+βt+γi+γit+δj+λyik,t1

In Equation (1), μi,t represents the conditional mean of the number of species for the ith Order for a specific year t. In Equation (2), μij,t represents the conditional mean of the number of insects for the ith Order, the jth specie, for year t. To account for temporal variation of the counts, the log-linear predictor incorporates an overall linear trend βt (fixed effect), a γi term to represent the Order (fixed effect), a δj term for the Specie (random effect), and a interaction γit Order–year (fixed effect). On the other hand, the component λ captures the impact of the past of the process yik,t1 in its future. The parameters of the models in Equations (1) and (2) are estimated via maximum likelihood and assuming two possible distributions for the response variable. The fitted models are compared to choose the best model to explain the variability of the counts by using likelihood measures and information criterions (AIC, BIC). All statistical data analyses are conducted in R statistical software by using the lm4 package (version 1.1-35.5).

The choice of generalised linear mixed models (GLMMs) for our analysis was driven by the complex nature of our ecological data and the specific research questions we aimed to address. GLMMs are particularly well suited for analysing count data in ecological studies as they can accommodate non-normal error distributions and account for both fixed and random effects. In our case, the use of GLMMs allowed us to model the discrete, non-negative nature of species and insect counts while accounting for the hierarchical structure of our data (species nested within orders).

The inclusion of temporal components and autoregressive terms in our models was crucial for capturing the dynamic nature of insect populations over time. Ecological systems often exhibit temporal dependencies, where the state of the system at one time point influences future states [58]. By incorporating a linear time trend (βt), we could model overall temporal changes in species and insect abundance. The interaction term between order and time (γit) allowed us to capture order-specific temporal trends, addressing potential differences in how various insect orders respond to environmental changes over time.

The autoregressive component λyik,t1 was included to account for temporal autocorrelation in our data. This term captures the influence of population sizes in the previous year on current year populations, a common phenomenon in population dynamics. While more complex time series models like ARIMA could have been considered, the relatively short duration of our study (10 years) limited their applicability. Our approach of incorporating autoregressive terms within the GLMM framework provides a robust alternative that accounts for temporal dependencies while allowing for the inclusion of other important predictors and random effects.

To address our research questions, we analysed the models described in Equations (1) and (2), with particular focus on the parameter β. The statistical significance and magnitude of β provide insights into the temporal effects on biodiversity. A statistically significant β indicates a meaningful change over time in either the number of species or the number of insects. Specifically, a positive β (β>0) suggests an improvement in biodiversity metrics over time in the presence of floral margins. However, to fully answer our research questions, we also considered the interaction between time and the presence of floral margins (γit). This approach allows us to assess both the overall temporal trends and the specific impact of floral margins on biodiversity enhancement over time. This modelling approach strikes a balance between model complexity and biological realism, allowing us to address our research questions while accounting for the inherent structure and temporal nature of our ecological data.

In order to complement the diversity analyses and inquire into community structure, log-series, log-normal, and broken-stick models were also applied [59]. The log-series model represents a community composed of a few abundant species and a high number of rare species. The broken-stick model refers to the maximum occupation of an environment with equitable sharing of resources between species. Finally, the log-normal model reflects an intermediate situation between the two [59]. Each of these models was applied to data obtained from the farm to calculate the expected number of species and the log2 grouping of species according to abundance [59,60,61]. To test the significance of the model outputs, the expected species values were compared with those of the observed species through a chi-square analysis [62].

3. Results

3.1. Exploratory Data Analysis

Table 1 summarises the identified number of species and insects and their rate of change over the years. A total of 172 species belonging to seven orders, Coleoptera (24), Diptera (31), Hemiptera (3), Hymenoptera (71), Lepidoptera (41), Neuroptera (1), and Odonata (1), were captured during the ten-year research program. Appendix B compiles the list and abundance of each species captured during the study. Figure 2 shows the scatterplot of the counts of species and insects across the years and by orders. These results show that while the number of species increases during the first three years and then stabilises, the number of individuals approximately follows a linear growth trend over the analysed period. The population dynamics also exhibits that the most frequent order is Hymenoptera. The number of species and insects of this order doubles over almost the whole period in comparison with other orders.

Figure 3 presents the temporal evolution of the indices of species diversity, -diversity, Shannon’s Index H, species richness, and Pielou’s evenness J, by orders. The estimated values for the -diversity indexes show that H-index and richness exhibit the same behaviour: For all orders examined, the values approximately increase during the first three years and then stabilise. However, the increase seems to be higher for Hymenoptera and Lepidoptera than for Coleoptera and Diptera. On the other hand, Pielou’s evenness J shows differences between the orders. While in Hymenoptera and Coleoptera, the J-index follows the same pattern, with values increasing during the first three years and then stabilising; it shows a roughly linear decrease trend over the years for Diptera and almost constant values for Lepidoptera during the entire study period.

Figure 4 shows lower triangular matrices with correlation coefficients for the number of species and insects between years. Additionally, Figure 4 only lists statistically significant correlations at the significance level of 5%. Where no value is displayed, the associated correlation coefficient is not significant. The temporal pattern of the number of species is more stable. There is less dynamics across the years because the correlation coefficients between consecutive years are low and closer to zero. This may reflect the fact that the number of identified species is approximately constant over the study period. On the other hand, the temporal pattern of the number of insects is more dynamic across the years because the correlation coefficients between years are high and all are statistically significant. Moreover, in most of the cases, the correlation coefficients for the number of insects are higher between successive years than between non-consecutive years. Finally, the insect population dynamics exhibits some degree of autoregressive effects, i.e., the past counts of insects explain how future counts will look like.

3.2. Statistical Modelling

For each response variable, i.e., the number of species or the number of insects, four models were fitted based on the Equations (1) and (2), respectively. The latter means that for each response variable, the fitted models were: (a) a full model if the associated counts follow a Poisson distribution, (b) a reduced model, dropping the interaction Order:time, assuming that the associated counts follow a Poisson distribution, (c) a full model, assuming that the associated counts follow a negative binomial distribution, and, (d) a reduced model, dropping the interaction Order:time, assuming that the associated counts follow a negative binomial distribution.

Table 2 presents the analysis of deviance table and the statistics of goodness of fit of the fitted models. The results show that in the case of the number of species, the best model is the reduced model, assuming a Poisson distribution for the response variable. It means that the count of the number of species does not show evidence that their mean is different than their variance, i.e., the population dynamics of the number of species across the years is stable and is explained by the order, which is the unique statistically significant parameter. On the other hand, in the case of the number of insects, the best model is the reduced model, assuming a negative binomial distribution for the response variable. It means that the abundance of insects is a process with high variability as the best model assumes that the variance of the number of insects depends on its mean. Additionally, for this model, it was identified that the count of insects in the previous year, the linear temporal trend, and the order are statistically significant parameters. Thus, the number of insects increases linearly across the study period with different starting points for the observed orders and depends on the previous state of the population.

Table 3 summarises the exponentiated estimated coefficients, relative risk (RR), and their confidence intervals for the two selected reduced models with regards to the number of species (the selected model for the number of species only includes the order as an independent variable) and the number of insects. For both models, the reference category was the order Coleoptera. In the case of the model for the number of species, all parameters are statistically significant, and the estimated relative risks are higher than one. This means that it is 31%, 194%, and 80% more likely to identify a specimen of the orders Diptera, Hymenoptera, and Lepidoptera than one of the order Coleoptera. In the same way, with regards to the model applied to the number of insects, it was found that having fixed all other independent variables, an additional year is associated with 15% more identified insects. On the other hand, fixing all other independent variables, an increase of one identified insect in a specific year implies an increase of 2% more identified insects in the following year.

3.3. Community Structure Models

Figure 5 illustrates the evolution of species abundance classes over the 10-year study period. Using a log2-based grouping method, we established seven distinct abundance classes based on the number of individuals per species. The results reveal that Class 1, comprising species with fewer than 2.5 individuals, initially increases but shows a declining trend over time. In contrast, the remaining classes exhibit progressive growth as the abundance of individuals per species increases. This growth is particularly pronounced in classes representing more than nine individuals per species.

Moreover, the data analysis shows a notable variation in the rate of change across different abundance classes, as illustrated by the dashed arrows in the Figure 5. This variation demonstrates three clear patterns: rapid flux in low abundance classes (Classes 2–3), representing species with relatively low abundance, which exhibits a markedly high rate of change; moderate transitions in middle classes (Classes 4–6), indicating a level of stability for species with intermediate abundance levels; and stability in the high abundance class (Class 7), comprising the most abundant species, which shows a very slow rate of change, meaning a high degree of persistence and stability for dominant species within the ecosystem.

The analysis of community structure models (Table 4) reveals that the observed community patterns align significantly with both log-series and log-normal distributions (p-value > 0.05). However, the community structure deviates significantly from the broken-stick model (p-value < 0.05). This pattern is indicative of an unstable community characterised by a small number of abundant species coexisting with a large number of rare species. Interestingly, these results suggest that habitat factors were not the primary determinants of community structure. This conclusion is supported by the observation that the sampling area exhibited highly specific floral and faunal compositions.

4. Discussion

The selection of the seed mixture for implementing floral margins plays a decisive role in the successful attraction of insects and speed of biodiversity protection [4,28,38,39,40,41,63,64,65,66,67,68]. According to our first working hypothesis, this significant association between plants and insects can increase the number of species and individuals over time. Several studies have highlighted the influence of floral margins on insect abundance and their role as conservation practices, mainly associated with pollinators such as bees, butterflies, or beetles [32,63,65,66,67]. Most of the research carried out so far focuses on studies covering just one year (growing season) [41,64,65,66,67,68,69,70] or at most three [38,41]. In contrast, the present work covers a period of 10 consecutive years. Comparing our results with those of three-year studies, we found similar growth patterns triggered by floral margins, with rising numbers for both the number of insect species and individuals. Looking at the data that our research generated in the first three years (2013–2015), there is both an increase in the number of species (116.67%) and an increase in the number of individuals (78.75%).

According to our second hypothesis, the use of floral margins improves biodiversity over time. However, the dynamics of change in the number of species tends to stabilise after the third year (the increase of species is 0.63% between 2016 and 2022, with a total increase of 138.80% between 2013 and 2022), while the number of individuals shows a linear growth trend over the 10-year period (the increase of individuals is 173.92% bet-ween 2016 and 2022, with a total increase of 403.33% between 2013 and 2022). These results fit very well with the dynamics of populations: The number of species remains flat after three or four years because no new or modified management practices were implemented at the level of both crop and floral composition due to the fact that the planted crops (celery, lettuce, onion) have practically very similar agricultural management and, moreover, none of them are insect-dependent [71,72,73]. In contrast, the analysis of population dynamics reveals a self-induced process: insect abundance over successive years is related to insect species presence and numbers in the previous year, having a consistent increase in the abundance of certain species. This gradient in change across abundance provides valuable insights into community dynamics. It suggests that rare species are more susceptible to fluctuations, possibly due to environmental changes or competitive pressures; intermediate abundance species show a balance between stability and responsiveness to ecological factors, while the most abundant species demonstrate resilience, maintaining their dominance over time. The only articles we have found are based on studies conducted outside the agricultural ecosystem, where several authors have observed similar compliance with log-series and log-normal models in cerambycids (Coleoptera) [74,75] and braconids (Hymenoptera) [76,77,78,79]. Additionally, Lima et al. [80] observed similar trends in their population dynamics and demographics of the northern short-tailed shrew year after year. These observations have important implications for understanding ecosystem stability, succession processes, and potential responses to environmental changes or management interventions.

Moreover, our study is the first long-term work to assess the effects of floral margins on insect diversity. However, extrapolating the information from Noordijk et al. [81], a study of the impact of flower margins’ age on different groups of ground-dwelling species, we can observe that there is a correlation with species abundance growth over time. Recently, Claire et al. [29] analysed the effects of flower margins’ age on pollinator abundance in Hungary and found that abundance was higher when margins were younger and lower in older margins. However, a comparison of our research with the former is not practicable as in our study the margins were partially re-sown to maintain a good floral diversity.

5. Conclusions

Our findings provide compelling evidence that the establishment of floral margins in agricultural landscapes significantly enhances biodiversity over time. Field margins sown with diverse plant mixtures serve a crucial dual purpose: they not only contribute to biodiversity conservation but also boost the abundance of both species and individuals. This positive impact is observed in the short term and, importantly, persists over longer periods.

Moreover, our research underscores the critical role of flower margins as an essential and enduring strategy for biological conservation and ecosystem enhancement. This approach is particularly valuable in intensively farmed areas, where biodiversity is often under significant pressure. The implementation of flower margins offers a practical and effective method to counterbalance the ecological impacts of intensive agriculture.

Author Contributions

Conceptualisation, F.J.P.-F. and F.S.; methodology, R.A.-S. and O.A.-M.; software, F.S.; validation, L.M.-B.; formal analysis, F.S.; investigation, F.J.P.-F., R.A.-S. and O.A.-M.; resources, L.M.-B.; data curation, O.A.-M., F.J.P.-F., F.S. and A.L.G.-Q.; writing—original draft preparation, F.J.P.-F., F.S. and A.L.G.-Q.; writing—review and editing, F.J.P.-F. and F.S.; visualisation, F.J.P.-F.; supervision, L.M.-B. and F.G.-V.; project administration, L.M.-B.; funding acquisition, F.G.-V. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Appendix B.

Acknowledgments

We are very grateful to the staff of G’s (Murcia) and especially to Fernando Yutera for their kindness and support during our research. Moreover, we want to thank German Canomanuel for their involvement and work in the launch and development of the study. This research was funded by the project Operation Pollinator (Syngenta). Finally, we are very thankful to Rudolf Gugger for providing a new view of the MS.

Conflicts of Interest

Author Francisco Javier Peris-Felipo, Fernando Santa, Ana Lia Gayan-Quijano, Luis Miranda-Barroso and Francisco Garcia-Verde were employed by the company Syngenta. The authors declare that this study received funding from Syngenta. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables
View Image - Figure 1. Floral margin on a farm in Águilas (Murcia) and its location in Spain.

Figure 1. Floral margin on a farm in Águilas (Murcia) and its location in Spain.

View Image - Figure 2. Scatterplot of the number of species and insects across the years. (a) Number of species. (b) Number of insects.

Figure 2. Scatterplot of the number of species and insects across the years. (a) Number of species. (b) Number of insects.

View Image - Figure 3. Scatterplot of the α-diversity indices: Shannon’s H (a), species richness (b), and Pielou’s evenness J across years (c).

Figure 3. Scatterplot of the α-diversity indices: Shannon’s H (a), species richness (b), and Pielou’s evenness J across years (c).

View Image - Figure 4. Correlation matrices of the number of species and insects between years. (a) Number of species. (b) Number of insects.

Figure 4. Correlation matrices of the number of species and insects between years. (a) Number of species. (b) Number of insects.

View Image - Figure 5. Changes of the species by abundance classes throughout the study (the arrows indicate the direction of change).

Figure 5. Changes of the species by abundance classes throughout the study (the arrows indicate the direction of change).

Number of species and insects and their rate of change across the years.

Year Number ofSpecies Rate ofChange (%) Number ofIndividuals Rate ofChange (%)
2013 67 241
2014 108 61.2 349 44.8
2015 150 38.8 443 26.9
2016 159 6.0 471 6.3
2017 140 −11.9 581 23.3
2018 139 −0.7 535 −7.9
2019 168 20.8 882 64.8
2020 167 −0.6 1038 17.7
2021 147 −11.9 907 −12.6
2022 160 8.8 1273 40.3

Analysis of deviance table and statistics of goodness of fit of the fitted generalised linear models and generalised linear mixed models for the number of species and insects, respectively (*** [0, 0.001]; ** [0.001, 0.01]; * [0.01, 0.05]).

Number of Species
Analysis of Deviance Table
Effect Poisson Poisson Negative Binomial Negative Binomial
Full Reduced Full Reduced
LR Chisq LR Chisq LR Chisq LR Chisq
Lag num. Species 1.13 1.30 1.13 1.2976
Order 30.76 *** 30.76 *** 30.76 *** 30.76 ***
Time 1.25 1.25 1.25 1.25
Order:Time 0.85 0.85
Statistics of Goodness of Fit
AIC 229.3 224.1 231.3 226.1
BIC 243.5 233.6 247.1 237.2
Log.Lik. −105.630 −106.057 −105.630 −106.057
RMSE 4.22 4.3 4.22 4.3
Number of Insects
Analysis of Deviance Table
Effect Poisson Poisson Negative Binomial Negative Binomial
Full Reduced Full Reduced
LR Chisq LR Chisq LR Chisq LR Chisq
Lag num. Insects 0.41 0.13 22.53 *** 23.30 ***
Order 10.91 * 10.90 * 11.71 ** 11.75 **
Time 804.31 *** 812.50 *** 238.09 *** 237.30 ***
Order:Time 31.56 *** 3.06
Statistics of Goodness of Fit
AIC 6946.6 6971.9 6107.1 6104.1
BIC 6999.5 7008.9 6165.2 6146.4
ICC 0.7 0.7 0.5 0.5
RMSE 4.85 4.95 11.48 12.58

Exponentiated estimated regression coefficients and 95% confidence of the fitted generalised linear models and generalised linear mixed models for the number of species and insects, respectively (n.a.—not included because it was not statistically significant).

Parameter Number of Species Number of Insects
RR (95% CI) RR (95% CI)
(Intercept) 18.9 (16.33–21.72) 1.56 (1.09–2.24)
Order Diptera 1.31 (1.08–1.58) 0.67 (0.43–1.06)
Order Hymenoptera 2.94 (2.5–3.47) 0.56 (0.38–0.84)
Order Lepidoptera 1.8 (1.51–2.15) 0.88 (0.57–1.36)
Lag count _ 1.02 (1.01–1.02)
Time _ 1.15 (1.13–1.17)

Analysis of the community structure according to abundance models (log-normal, log-series, and broken-stick) for the insect community (*** [0, 0.001]; ** [0.001, 0.01]; * [0.01, 0.05]).

Community Structure Models
Model Year 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022
Log-normal Chisq 6.81 3.42 8.20 10.12 * 6.39 6.41 8.46 7.49 5.88 6.77
p-value 0.234 0.635 0.145 0.034 * 0.380 0.260 0.132 0.277 0.436 0.342
Log-series Chisq 9.12 4.260 3.238 3.208 4.013 3.420 4.288 3.263 2.565 2.611
p-value 0.166 0.512 0.663 0.523 0.674 0.6355 0.508 0.775 0.861 0.855
Broken-stick Chisq 21.212 18.464 17.809 17.832 20.313 15.750 20.595 27.982 23.284 32.623
p-value 0.0007 *** 0.002 ** 0.003 ** 0.001 ** 0.002 ** 0.007 ** 0.0009 *** 0.0001 *** 0.0007 *** 0.0001 ***

Appendix A. Annual Rainfall (mm) and Temperature (°C) Data Recorded at Águilas Farm During the Study Period (2013–2022)

Year Total Rainfall (mm) Avg. Temperature (°C) Avg. Max. Temp. ( ° C) Avg. Min. Temp. ( ° C)
2013 200.3 19.8 24.5 15.2
2014 187.7 20.6 25.4 16.0
2015 285.6 20.4 25.2 16.0
2016 281.4 20.5 25.0 16.2
2017 234.3 20.3 25.0 15.7
2018 357.6 20.1 24.5 15.8
2019 419.1 20.1 24.8 15.6
2020 294.7 20.1 24.8 15.7
2021 333.7 20.2 24.7 16.1
2022 441.6 20.9 25.6 16.4

Appendix B. List of Species and Their Abundance During the Study Period (2013–2022)

Order Species 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022
Lepidoptera Acontia lucida (Hufnagel, 1776) 1 3 1 3 4 2
Coleoptera Adalia bipunctata (L., 1758) 1 2 2 6 21 1 33 39 21 55
Coleoptera Agapanthia asphodeli (Latreille, 1804) 1 2 1 3 1 2 2 1
Coleoptera Agapanthia cardui (L., 1767) 2 2 1 1 1 1 2 4 3 4
Lepidoptera Aglais urticae (L., 1758) 1 2 1 2 2
Hymenoptera Amegilla quadrifasciata (Villers, 1789) 2 3 1 2 2 4 3 2
Hymenoptera Amegilla sp. 1 1 1 1 1 1 1 1 1
Hymenoptera Ammophila heydeni Dalhbom, 1845 1 2 3 3 3 1 5 13 8 1
Hymenoptera Ammophila sabulosa (L., 1758) 2 1 4 1 2 2 1 1 1
Hymenoptera Ancistrocerus biphaleratus (Saussure, 1852) 2 2 1 1 1 1 1 1
Hymenoptera Andrena albopunctata (Rossi, 1792) ssp. melona Warncke, 1967 1 1 4 2 2 2 2 3 2 2
Hymenoptera Andrena angustior (Kirby, 1802) ssp. impressa Warncke,1967 5 1 2 1
Hymenoptera Andrena asperrima Pérez, 1895 1 1 1 1 1 2 1 2
Hymenoptera Andrena flavipes Panzer, 1799 4 5 5 6 5 9 5 4
Hymenoptera Andrena fuscosa (Erichson, 1835) 1 1 1 2 1 2 3 2
Hymenoptera Andrena hispania Warncke, 1967 1 1 3 2 3 2 1 3 1 2
Hymenoptera Andrena nilotica Warncke, 1967 1 1 3 2
Hymenoptera Andrena vetula Lepeletier, 1841 1 1 2 2 5 10 16 14 13 3
Coleoptera Anisoplia baetica (Erichson, 1847) 1 1 3 2 1 5 1
Hymenoptera Anthidiellum strigatum (Panzer, 1805) 1 1 1 2 1 3
Hymenoptera Anthidium florentinum (F., 1775) 1 1 1 2 1 1 2 1 2
Hymenoptera Anthidium manicatum (L., 1758) 2 2
Hymenoptera Anthophora atroalba Lepeletier, 1841 1 2 1 2 2 3 1
Hymenoptera Anthophora ferruginea (Lepeletier, 1841) 3 1 3 1 1 4
Diptera Anthrax sp. 1 1 1 1 1 1
Hymenoptera Apis mellifera L., 1758 70 90 56 20 97 106 42 42 149 112
Lepidoptera Aporia crataegi (L. 1758) 1 2 9 2 3 3 3 1 5
Lepidoptera Aricia cramera (Escholtz, 1821) 2 3 7 5 6 6 9 7 7
Hymenoptera Athalia nevadensis Lacourt, 1987 1 1 1 1 1 2 1 1 1 1
Hymenoptera Bembyx merceti J. Parker, 1904 1 1 1 1 2 1 1 5 2
Hymenoptera Bombus terrestris (L., 1758) 7 3 3 12 9 14 14 7 24
Diptera Bombylella atra (Scopoli, 1763) 1 1 2 3 5 3 7
Diptera Bombylius discolor Mikan, 1796 1 3 2 1 5 2 5 4
Diptera Bombylius sp. 1 1 1 1 2 1 2
Lepidoptera Brintesia circe (F. 1775) ssp. hispanica (Spuler, 1908) 1 1 1 2 1 1 1 4 3
Lepidoptera Cacyreus marshalli (Butler, 1898) 6 9 11 8 9 2 4 9
Diptera Calliphora vomitoria (L., 1758) 1 1 1 1 1 2 2 2 2
Lepidoptera Callophrys rubi (L., 1758) 2 2 4 2 5 5 2 2 7
Hymenoptera Camponotus cruentatus (Latreille, 1802) 3 12 22 12 22 22 11 1
Lepidoptera Carcharodes alceae (Esper, 1780) 1 1 2 1 1 2 2 1 1 1
Hymenoptera Ceratina chalcites Germar, 1839 1 3 2 1 1 1 12 12 11 7
Hymenoptera Cerceris arenaria (L., 1758) 1 1 4 1 1 1 1 1 3
Coleoptera Cerocoma schreberi F., 1781 7 9 11 1 1 3
Coleoptera Chlorophorus trifasciatus (F., 1781) 2 3 5 2 3 6 5 7 12
Hymenoptera Chrysis sp. 1 1 2 2 1 1 3 2
Neuroptera Chrysoperla carnea (Stephens, 1836) 5 4 7 11 24 49 27 47
Diptera Chrysotoxum festivum (L., 1758) 3 5 9 16 12 16 14
Coleoptera Clanoptilus sp. 1 2 1 3 1 2 3
Coleoptera Coccinella septempunctata (L., 1758) 34 27 19 25 12 3 36 12 11 34
Lepidoptera Coenonympha pamphilus (L., 1758) 1 2 1 4 5 7 7 12 12 18
Lepidoptera Colias alfacariensis (Ribbe, 1905) 2 5 4 1 4 4 9 4 7
Lepidoptera Colias crocea (Geoffroyi in Fou- reroy, 1785) 6 5 8 11 3 16 16 22 12 48
Hymenoptera Colletes similis Schenck, 1853 2 1 1 2 1 1 3 9
Hymenoptera Colpa sexmaculata (F., 1781) 2 2 1 3 2 2 2 1 5
Hemiptera Coranus subapterus (De Geer, 1773) 1 1 1 3 2 5
Diptera Cylindromyia bicolor (Olivier, 1812) 2 2 2 1 2 2 1 5
Diptera Cylyndromyia pilipes (Loew, 1844) 1 2 1 1 2 1 1
Diptera Dischistus biroi (Becker, 1906) 1 3 1 2 1 1
Hymenoptera Epeolus cruciger (Panzer, 1799) 1 1 1 1
Hymenoptera Epeolus fallax Morawitz, 1872 1 1
Diptera Episyrphus balteatus (De Geer, 1776) 3 4 11 21 34 9 41 82 41 11
Diptera Eristalis arbustorum (L., 1758) 2 1 7 1 9 3 9 7
Diptera Eristalix tenax (L., 1758) 2 2 4 1 4 11 4 3
Lepidoptera Erynnis tages (L., 1758) 1 1 1 1 1 1 1 2
Hymenoptera Eucera cineraria Eversmann, 1852 1 1 2 1 1 2 3 5 2 2
Hymenoptera Eucera elongatula Vachal, 1907 1 2 2 2 6 4 1 2 2 1
Hymenoptera Eucera notata Lepeletier, 1841 2 4 7 6 7 11 9 9 2 2
Lepidoptera Euchloe belemia (Esper, 1800) 6 9 10 7 18 28 16 2 20
Hymenoptera Eumenes mediterraneus Kriechbaumer, 1879 1 2 2 2 1 1 3
Coleoptera Exosoma lusitanicum (L., 1758) 6 7 9 11 2 1 4
Lepidoptera Glaucopsyche alexis (Poda, 1761) 1 2 3 1 3 2 5 4 1
Diptera Graphomya maculata (Scopoli, 1763) 2 12 1 6 5 3 7 3 7 2
Diptera Gymnosoma rotundatum (L., 1758) 3 6 9 11 9 3 12 5 12 5
Hymenoptera Halictus scabiosae (Rossi,1790) 1 1 2 1 2 1 3 5 3 5
Coleoptera Heliotaurus ruficollis (F., 1781) 7 9 11 12 5 4 15 21 13 58
Diptera Hemipenthes morio (L., 1758) 1 2 2 1 2 1 1 3 1 6
Diptera Heteralonia algira (F., 1794) 1 1 2 2 4 2 2
Lepidoptera Hipparchia semele (L. 1758) 1 2 1 4 4 1 2
Hymenoptera Hoplitis adunca (Panzer, 1798) 2 3 6 4 8 1 9 12 19 2
Hymenoptera Hoplitis cristata (Fonscolombe, 1846) 1 1 2 2 2 1 5 3 2 1
Coleoptera Hycleus duodecimpunctatus (Oli-vier, 1811) 6 5 4 7 11 2 2 1
Hymenoptera Hylaeus (Prosopis) pictus (Smith, 1853) 2
Hymenoptera Hylaeus (Prosopis) variegatus (F., 1798) 1 1 1 2 3 5 2 12 17 11
Hymenoptera Icneumon sp. 2 2 3 3 2 2 3
Lepidoptera Iphiclides feisthamelii (Duponchel 1832) 1 2 2 2 5 1
Odonata Ischnura graellsii (Rambur, 1842) 2 1 1 4 2 1 1
Lepidoptera Issoria lathonia (L., 1758) 4 3 5 3 1 7 15 12 3 18
Coleoptera Lachnaia pubescens (Dufour, 1820) 2 5 4 7 5 3 3
Coleoptera Lachnaia tristigma (Lacordairei, 1868) 3 1 1 3 1 1 1
Lepidoptera Lampides boeticus (L., 1767) 11 4 7 6 3 9 9 29 21 36
Hymenoptera Lasioglossum aegyptiellum (Strand, 1909) 1 1 1 1 1
Hymenoptera Lasioglossum malachurum (Kirby, 1802) 1 1 1 4 6 11 9 14 9 19
Hymenoptera Lasioglossum sp. 1 1 2 2 6 3 2 3 11
Lepidoptera Leptotes pirithous (L., 1767) 2 2 4 2 2 6 6 3 3 24
Hymenoptera Lestica clypeata (Schreber, 1759) 1 2 1 2 3 5 5 2
Diptera Lomatia lateralis (Meigen, 1820) 3 5 7 9 11 5 10 5 10 8
Lepidoptera Lycaena phlaeas (L., 1761) ssp. lusitanica (Bryk, 1940) 5 4 9 5 7 10 10 32 19 34
Lepidoptera Lysandra bellargus (Rottemburg, 1775) 2 1 5 5 1 7
Lepidoptera Macroglossum stellatarum (L., 1758) 9 3 4 2 4 5 8 12 5 14
Lepidoptera Maniola jurtina (Linneo 1758) ssp. hispulla (Esper 1805) 2 2 1 2 2 1 1 5
Hymenoptera Megascolia maculata (Drury, 1773) ssp.flavifrons (F., 1773) 1 1 1 1 1 1 1 1 2
Lepidoptera Melanargia lachesis (Hübner, 1790) 1 1 1 1 1 1 1 5
Lepidoptera Melitaea phoebe (Denis y Schiffermüller, 1775) 1 1 2 2 2 2 1 1 7
Hymenoptera Melitta murciana Warncke, 1973 5 1 1 1
Hymenoptera Meria tripunctata (Rossi, 1790) 1 3 1 1 1 2 2 1 9
Hymenoptera Mutilla europaea L., 1758 1 1 1 1 2 1 1 1 1
Diptera Myathropa florea (L., 1758) 1 1 1 1 1 2
Coleoptera Mylabris quadripunctata (L., 1767) 1 2 2 1 1 3 6 2 4
Coleoptera Mylabris variabilis (Pallas, 1781) 2 1 1 1 1 1 4 2 1 5
Hymenoptera Nomada agrestis F., 1787 1 1 2 3 1 2 1 3
Coleoptera Oedemera nobilis (Scopoli, 1763) 3 2 1 2 5 12 8 11
Hymenoptera Ophion luteus (L., 1758) 1 1 2 3 5 4 7 4 6 12
Coleoptera Opsilia coerulescens (Scopoli, 1763) 1 3 1 3 5
Hymenoptera Osmia aurulenta Panzer, 1799 1 2 1 1 1 1 2
Hymenoptera Osmia bicornis (L., 1758) 2 1 1 1 1 1 1 1 2
Hymenoptera Osmia caerulescens L., 1758 1 1 2 1 1 2 2 1
Hymenoptera Oxybelus quattuordecimnotatus Jurine, 1807 1 1 2 1 1 1 1
Coleoptera Oxythyrea funesta (Poda, 1761) 2 3 3 1 4 1 7 13 11 19
Lepidoptera Pandoriana pandora (Denis y Schiffermüller, 1775) 2 2 3 3 1 1 2 1
Hymenoptera Panurgus calcaratus (Scopoli, 1763) 2 3 5 1 4 1 2 5 3 2
Hymenoptera Panurgus cephalotes Latreille, 1811 1 1 2 1 2 2 2 1 3 1
Lepidoptera Papilio machaon (L., 1758) 1 1 2 2 2 1 1 4
Lepidoptera Pararge aegeria (L., 1767) 3 1 1 1 4 4 5 7 7
Hemiptera Peirates stridulus (F., 1787) 1 2 2 1 2 1 5 5 2
Diptera Peleteria meridionalis (Robineau-Desvoidy, 1830) 1 1 3 2 3 2
Hymenoptera Philanthus coronatus ibericus (Thunberg, 1794) 1 1 1 1 1 2 2 1
Hymenoptera Philanthus triangulum (F., 1775) 14 2 1 1 2 29 1
Diptera Phthiria pulicaria (Mikan, 1796) 1 1 1 1 1 2
Diptera Physocephala rufipes (Olivier, 1795) 1 2 1 1 3 1 2 2 2 3
Lepidoptera Pieris brassicae (L., 1758) 7 9 11 6 27 5 6 5 3 4
Lepidoptera Pieris rapae (L., 1758) 2 5 9 4 11 13 15 3 3 2
Hymenoptera Polistes biglumis (L., 1758) 2 3 2 2 2 1 5 17 17 7
Lepidoptera Polyonmatus icarus (Rottemburg, 1775) 2 1 3 3 1 3 3 12 11 19
Lepidoptera Pontia daplidice (L., 1758) 1 1 2 1 3 3 5 6 3
Hymenoptera Pseudopipona lativentris (Saussure, 1853) 1 2 1 2 3 2 1 4
Lepidoptera Pyronia bathseba (F., 1793) ssp. amyclas (Fruhstofer 1910) 1 1 1 1 1 1 1 2
Lepidoptera Pyronia tithonus (L., 1771) 1 1 1 1 1 1 2
Coleoptera Rhagonycha fulva (Scopoli, 1763) 5 7 1 2 1 1 3 19 29 49
Hemiptera Rhynocoris iracundus (Poda, 1761) 1 1 1 2 1 1 2
Hymenoptera Rodanthidium sticticum (F., 1787) 1 1 2 3 5 4 7 18 21 13
Diptera Sarcophaga sp. 1 1 1 1 2 1 1
Diptera Saropogon leucocephalus (Meigen, 1820) 1 1 1 2 1 3 9 3 1
Diptera Saropogon sp. 1 1 1 1 1 1
Diptera Scaeva pyrastri (L., 1758) 1 1 1 34 26 34 21
Hymenoptera Sceliphron destillatorium (Illiger, 1807) 1 1 1 1 1 1 1 1 2
Hymenoptera Scolia carbonaria (L., 1767) 4 1 1 1 1 2 2 4 1
Hymenoptera Scolia erythrocephala (F., 1798) 1 1 2 1 1 1 1 2 1 1
Diptera Sphaerophoria scripta (L., 1758) 2 3 3 4 8 11 22 11 67
Hymenoptera Sphex funerarius Gussakovskij, 1934 1 1 1 2 3 2 2 2 2 3
Lepidoptera Spialia sertorius (Hoffmannsegg 1804) 1 1 1 1 1 3
Coleoptera Stenopterus ater L., 1767 1 2 1 2 1 1 1 2
Coleoptera Stenurella melanura (L., 1758) 1 1 2 2 2 4 9
Hymenoptera Stilbum cyanurum (Forster, 1771) 1 1 1 1 1 1 1 1
Lepidoptera Syrichtus proto (Esper 1808) 1 1 1 1 3 2
Diptera Syritta pipiens (L., 1758) 4 5 11 2 9 3 9 7
Diptera Tachina fera (L., 1761) 1 1 2 1 1 1 1 1 1
Hymenoptera Tenthredo baetica Spinola, 1843 2 1 2 2 3 2 5 5 2 3
Hymenoptera Tetraloniella nana Morawitz, 1873 5
Lepidoptera Thymelicus lineola (Ochsenheimer, 1808) 2 1 1 1 1 2 4 5
Diptera Thyridanthrax fenestratus (Fallen, 1814) 1 1 2 1 1 1 1 2 1 2
Coleoptera Trichodes leucopsideus (Olivier, 1795) 2 1 1 1 2 1 2 8 2 5
Coleoptera Trichodes octopunctatus (F., 1787) 1 1 1 1 3 1 2 2 2 2
Diptera Trichopoda pennipes (F., 1781) 2 1 1 1 1 1 2 1 1 1
Coleoptera Tropinota squalida (Scopoli, 1783) 1 1 1 3 1 1 3 2 4 22
Lepidoptera Vanessa atalanta (L., 1758) 1 4 5 2 6 9 3 13 1
Lepidoptera Vanessa cardui (L., 1758) 2 2 4 11 3 6 6 7 9 22
Hymenoptera Vespula germanica (F., 1793) 1 1 2 2 1 2 1 2 1
Diptera Villa hottentota (L., 1758) 1 1 2 1 2 7
Hymenoptera Xylocopa iris (Christ, 1791) ssp. uclesiensis Pérez, 1901 1 1 1 3
Hymenoptera Xylocopa violacea (L., 1758) 3 1 1 1 1 2 7 12 9 5
Lepidoptera Zygaena filipendulae (L., 1758) 2 2 2 1 2 3
Lepidoptera Zygaena sarpedon (Hübner, 1790) 1 2 1 1 4 2 1 8

References

1. Renard, D.; Tilman, D. Cultivate biodiversity to harvest food security and sustainability. Curr. Biol.; 2021; 31, pp. R1154-R1158. [DOI: https://dx.doi.org/10.1016/j.cub.2021.06.082] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34637721]

2. Malo, M. Sustainable Agriculture: Need of the Hour. Agric. Food E Newsl.; 2020; 2, pp. 298-300.

3. Loreau, M.; Naeem, S.; Inchausti, P.; Bengtsson, J.; Grime, J.P.; Hector, A.; Hooper, D.U.; Huston, M.A.; Raffaelli, D.; Schmid, B. et al. Biodiversity and ecosystem functioning: Current knowledge and future challenges. Science; 2001; 294, pp. 804-809. [DOI: https://dx.doi.org/10.1126/science.1064088] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11679658]

4. Rodríguez-Gasol, N.; Avilla, J.; Aparicio, Y.; Arnó, J.; Gabarra, R.; Riudavets, J.; Alegre, S.; Lordan, J.; Alins, G. The contribution of surrounding margins in the promotion of natural enemies in Mediterranean apple orchards. Insects; 2019; 10, 148. [DOI: https://dx.doi.org/10.3390/insects10050148]

5. Campbell, A.J.; Wilby, A.; Sutton, P.; Wäckers, F. Getting more power from your flowers: Multi-functional flower strips enhance pollinators and pest control agents in apple orchards. Insects; 2017; 8, 101. [DOI: https://dx.doi.org/10.3390/insects8030101]

6. Meena, R.S.; Kumar, S.; Datta, R.; Lal, R.; Vijayakumar, V.; Brtnicky, M.; Sharma, M.P.; Yadav, G.S.; Jhariya, M.K.; Jangir, C.K. et al. Impact of Agrochemicals on Soil Microbiota and Management: A Review. Land; 2020; 9, 34. [DOI: https://dx.doi.org/10.3390/land9020034]

7. Wratten, S.D.; Gillespie, M.; Decourtye, A.; Mader, E.; Desneux, N. Pollinator habitat enhancement: Benefits to other ecosystem services. Agric. Ecosyst. Environ.; 2012; 159, pp. 112-122. [DOI: https://dx.doi.org/10.1016/j.agee.2012.06.020]

8. Mei, Z.; de Groot, G.A.; Kleijn, D.; Dimmers, W.; van Gils, S.; Lammertsma, D.; van Kats, R.; Scheper, J. Flower availability drives effects of wildflower strips on ground-dwelling natural enemies and crop yield. Agric. Ecosyst. Environ.; 2021; 319, 107570. [DOI: https://dx.doi.org/10.1016/j.agee.2021.107570]

9. Pollier, A.; Tricault, Y.; Plantegenest, M.; Bischoff, A. Sowing of margin strips rich in floral resources improves herbivore control in adjacent crop fields. Agric. Forest. Entomol.; 2019; 21, pp. 119-129. [DOI: https://dx.doi.org/10.1111/afe.12318]

10. Our World in Data. Available online: https://ourworldindata.org/global-land-for-agriculture (accessed on 5 September 2022).

11. Haaland, C.; Naisbit, R.E.; Bersier, L.F. Sown wildflower strips for insect conservation: A review. Insect Conserv. Divers.; 2011; 4, pp. 60-80. [DOI: https://dx.doi.org/10.1111/j.1752-4598.2010.00098.x]

12. Holzschuh, A.; Steffan-Dewenter, I.; Tscharntke, T. Grass strip corridors in agricultural landscapes enhance nest-site colonization by solitary wasps. Ecol. Appl.; 2009; 19, pp. 123-132. [DOI: https://dx.doi.org/10.1890/08-0384.1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19323177]

13. Thomas, C.F.G.; Cooke, H.; Bauly, J.; Marshall, E.J.P. Invertebrate colonisation of overwintering sites in different field boundary habitats. Asp. Appl. Biol.; 1994; 40, pp. 229-232.

14. Karamaouna, F.; Jaques, J.A.; Kati, V. Practices to Conserve Pollinators and Natural Enemies in Agro-Ecosystems. Insects; 2021; 12, 31. [DOI: https://dx.doi.org/10.3390/insects12010031]

15. Celeste, A.; Medina, P.; Adán, A.; Sánchez-Ramos, I.; del Estal, P.; Fereres, A.; Viñuela, E. The Role of Annual Flowering Plant Strips on a Melon Crop in Central Spain. Influence on Pollinators and Crop. Insects; 2020; 11, 66. [DOI: https://dx.doi.org/10.3390/insects11010066]

16. Herz, A.; Cahenzli, F.; Penvern, S.; Pfiner, L.; Tasin, M.; Sigsgaard, L. Managing Floral Resources in Apple Orchards for Pest Control: Ideas, Experiences and Future Directions. Insects; 2019; 10, 247. [DOI: https://dx.doi.org/10.3390/insects10080247]

17. Antkowiak, M.; Kowalska, J.; Trzcinski, P. Flower Strips as an Ecological Tool to Strengthen the Environmental Balance of Fields: Case Study of a National Park in Western Poland. Sustainability; 2024; 16, 1251. [DOI: https://dx.doi.org/10.3390/su16031251]

18. Purvis, E.; Meehan, M.; Lindo, Z. Agricultural field margins provide food and nesting resources to bumble bees (Bombus spp., Hymenoptera: Apidae) in Southwestern Ontario, Canada. Insect Conserv. Diver.; 2020; 13, pp. 219-228. [DOI: https://dx.doi.org/10.1111/icad.12381]

19. Mota, L.; Hevia, V.; Rad, C.; Alves, J.; Silva, A.; González, J.A.; Ortega-Marcos, J.; Aguado, O.; Alcorlo, P.; Azcárate, F.M. et al. Flower strips and remnant semi-natural vegetation have different impacts on pollination and productivity of sunflower crops. J. Appl. Ecol.; 2022; 59, pp. 2386-2397. [DOI: https://dx.doi.org/10.1111/1365-2664.14241]

20. Smith, J.; Potts, S.G.; Woodcock, B.A.; Eggleton, P. Can arable field margins be managed to enhance their biodiversity, conservation and functional value for soil macrofauna?. J. Appl. Ecol.; 2008; 45, pp. 269-278. [DOI: https://dx.doi.org/10.1111/j.1365-2664.2007.01433.x]

21. Nowakowski, M.; Pywell, R. Habitat Creation and Management for Pollinators; Centre for Ecology & Hydrology: Wallingford, UK, 2016.

22. Schumacher, W. Gefährdete Ackerwildkräuter können auf ungespritzten Feldrändern erhalten werden. Mitt. LÖLF; 1984; 9, pp. 14-20.

23. Jörg, E. Field Margin-strip programmes. Proceedings of the Technical Seminar of the Landesanstalt für Pflanzenbau und Pflanzschutz, Pretty Print; Mainz, Germany, 25–27 May 1994.

24. De Snoo, G.R. Unsprayed field margins: Effects on environment, biodiversity and agricultural practice. Landsc. Urb. Plan.; 1999; 46, pp. 151-160. [DOI: https://dx.doi.org/10.1016/S0169-2046(99)00039-0]

25. Marshall, E.J.P.; Moonen, A.C. Field margins in northern Europe: Their functions and interactions with agriculture. Agric. Ecosyst. Environ.; 2002; 89, pp. 5-21. [DOI: https://dx.doi.org/10.1016/S0167-8809(01)00315-2]

26. Jacot, K.; Eggenschwiler, L.; Junge, X.; Luka, H.; Bosshard, A. Improved field margins for a higher biodiversity in agricultural landscapes. Asp. Appl. Biol.; 2007; 87, pp. 277-283.

27. Albrecht, M.; Kleijn, D.; Williams, N.M.; Tschumi, M.; Blaauw, B.R.; Bommarco, R.; Campbell, A.J.; Dainese, M.; Drummond, F.A.; Entling, M.H. et al. The effectiveness of flower strips and hedgerows on pest control, pollination services and crop yield: A quantitative synthesis. Ecol. Lett.; 2020; 23, pp. 1488-1498. [DOI: https://dx.doi.org/10.1111/ele.13576]

28. Miranda-Barroso, L.; Aguado, O.; Falcó-Garí, J.V.; Lopez, D.; Schade, M.; Vasileiadis, V.; Peris-Felipo, F.J. Multifunctional areas as a tool to enhance biodiversity and promote conservation in alfalfa fields. J. Insect Biodivers. Syst.; 2021; 7, pp. 251-261. [DOI: https://dx.doi.org/10.52547/jibs.7.3.251]

29. Brittain, C.; Benke, S.; Pecze, R.; Potts, S.G.; Peris-Felipo, F.J.; Vasileiadis, V. Flower margins: Attractiveness over time for different pollinator groups. Land; 2022; 11, 1933. [DOI: https://dx.doi.org/10.3390/land11111933]

30. Sánchez, J.A.; Carrasco, A.; La Spina, M.; Pérez-Marcos, M.; Ortíz-Sanchez, F.J. How bees respond differently to field margins of shrubby and herbaceous plants in intensive agricultural crops of the Mediterranean area. Insects; 2020; 11, 26. [DOI: https://dx.doi.org/10.3390/insects11010026]

31. Hines, H.M.; Hendrix, S.D. Bumble bee (Hymenoptera: Apidae) diversity and abundance in tallgrass prairie patches: Effects of local and landscape floral resources. Environ. Entomol.; 2005; 34, pp. 1477-1484. [DOI: https://dx.doi.org/10.1603/0046-225X-34.6.1477]

32. Colla, S.R.; Taylor-Pindar, A. Recovery Strategy for the Rusty-Patched Bumble Bee (Bombus Affinis) in Ontario; Ontario Recovery Strategy Series; Queens Printer for Ontario: Peterborough, ON, Canada, 2011.

33. Morandin, L.A.; Kremen, C. Hedgerow restoration promotes pollinator populations and exports native bees to adjacent fields. Ecol. Appl.; 2013; 23, pp. 829-839. [DOI: https://dx.doi.org/10.1890/12-1051.1]

34. Gilpin, A.-M.; O’Brien, C.; Kobel, C.; Brettell, L.; Cook, J.; Power, S. Co-flowering plants support diverse pollinator populations and facilitate pollinator visitation to sweet cherry crops. Basic Appl. Ecol.; 2022; 63, pp. 36-48. [DOI: https://dx.doi.org/10.1016/j.baae.2022.05.005]

35. Zamorano, J.; Bartomeus, I.; Grez, A.A.; Garibaldi, L.A. Field margin floral enhancements increase pollinator diversity at the field edge but show no consistent spillover into the crop field: A meta-analysis. Insect Conserv. Divers.; 2020; 13, pp. 519-531. [DOI: https://dx.doi.org/10.1111/icad.12454]

36. Carvell, C.; Meek, W.R.; Pywell, R.F.; Goulson, D.; Nowakowski, M. Comparing the efficacy of agri-environment schemes to enhance bumble bee abundance and diversity on arable field margins. J. Appl. Ecol.; 2007; 44, pp. 29-40. [DOI: https://dx.doi.org/10.1111/j.1365-2664.2006.01249.x]

37. Geppert, C.; Hass, A.; Földesi, R.; Donkó, B.; Akter, A.; Tscharntke, T.; Batáry, P. Agri-environment schemes enhance pollinator richness and abundance but bumblebee reproduction depends on field size. J. Appl. Ecol.; 2020; 57, pp. 1818-1828. [DOI: https://dx.doi.org/10.1111/1365-2664.13682]

38. Santa, F.; Aguado, L.O.; Falcó-Garí, J.V.; Jiménez-Peydró, R.; Schade, M.; Vasileiadis, V.; Miranda-Barroso, L.; Peris-Felipo, F.J. Effectiveness of Multifunctional Margins in Insect Biodiversity Enhancement and RTE Species Conservation in Intensive Agricultural Landscapes. Agronomy; 2021; 11, 2093. [DOI: https://dx.doi.org/10.3390/agronomy11112093]

39. Castle, D.; Grass, I.; Westphal, C. Fruit quantity of strawberries benefit from enhanced pollinator abundance at hedgerows in agricultural landscapes. Agric. Ecosyst. Environ.; 2019; 275, pp. 14-22. [DOI: https://dx.doi.org/10.1016/j.agee.2019.01.003]

40. Meek, B.; Loxton, D.; Sparks, T.; Pywell, R.; Pickett, H.; Nowakowski, M. The effect of arable field margin composition on invertebrate biodiversity. Biol. Conserv.; 2002; 106, pp. 259-271. [DOI: https://dx.doi.org/10.1016/S0006-3207(01)00252-X]

41. Peris-Felipo, F.J.; Santa, F.; Aguado, O.; Falcó-Garí, J.V.; Iborra, A.; Schade, M.; Brittain, C.; Vasileiadis, V.; Miranda-Barroso, L. Enhancement of the diversity of pollinators and beneficial insects in intensively managed vineyards. Insects; 2021; 12, 740. [DOI: https://dx.doi.org/10.3390/insects12080740]

42. AEMET. Agencia Estatal de Meteorología. Available online: http://www.aemet.es/es/portada (accessed on 18 September 2022).

43. Aguado, L.O.; Viñuelas, E.; Ferreres, A. Guía de Polinizadores de la Península Ibérica y de Los Archipiélagos Balear y Canario; Ediciones Mundiprensa & Syngenta: Madrid, Spain, 2016.

44. Amiet, F.; Herrmann, M.; Müller, A.; Neumeyer, R. Fauna Helvetica APIDAE 3: Halictus, Lasioglossum; Centre Suisse de Cartographie de la Faune (CSCF Info Fauna) Schweizerische Entomologische Gesellschaft (SEG/SES): Bern, Switzerland, 2001.

45. Amiet, F.; Hermann, M.; Müller, A.; Neumeyer, R. Fauna Helvetica APIDAE 5: Ammobates, Ammobatoides, Anthophora, Biastes, Ceratina, Dasypoda, Epeoloides, Epeolus, Eucera, Macropis, Melecta, Melitta, Nomada, Pasites, Tetralonia, Thyreus, Xylocopa; Centre Suisse de Cartographie de la Faune (CSCF Info Fauna) Schweizerische Entomologische Gesellschaft (SEG/SES): Bern, Switzerland, 2007.

46. Amiet, F.; Hermann, M.; Müller, A.; Neumeyer, R. Fauna Helvetica APIDAE 6: Andrena, Meliturga, Panurginus, Panurgus; Centre Suisse de Cartographie de la Faune (CSCF Info Fauna) Schweizerische Entomologische Gesellschaft (SEG/SES): Bern, Switzerland, 2010.

47. Cobos, A. Fauna Ibérica de Coleópteros Buprestidae; Consejo Superior de Investigaciones Científicas: Madrid, Spain, 1986.

48. Eizaguirre, S. Coleóptera, Coccinellidae. Fauna Ibérica; Ramos, M.A. Museo Nacional de Ciencias Naturales (CSIC): Madrid, Spain, 2004; Volume 40.

49. Fernández-Rubio, F. Guía de las Mariposas Diurnas de la Península Ibérica, Baleares, Canarias, Azores y Madeira. 1; Libytheidae, Nymphalidae, Riodinidae y Lycaenidae Ediciones Pirámide, S.A.: Madrid, Spain, 1991.

50. Fernández-Rubio, F. Guía de las mariposas diurnas de la Península Ibérica, Baleares, Canarias, Azores y Madeira. 2; Papilionidae, Pieridae, Danaidae, Satyridae y Hesperiidae Ediciones Pirámide, S.A.: Madrid, Spain, 1991.

51. Fernández-Rubio, F. Clave para la determinación de las especies españolas del género Zygaena F., 1775 (Insecta: Lepidoptera). Graellsia; 2006; 62, pp. 3-12. [DOI: https://dx.doi.org/10.3989/graellsia.2006.v62.i1.24]

52. Martín-Piera, F.; López-Colón, J.I. Coleóptera Scarabaeoidea I. Fauna Ibérica; Ramos, M.A. Museo Nacional de Ciencias Naturales (CSIC): Madrid, Spain, 2000; Volume 14.

53. Oosterbroek, P. The European Families of the Diptera. Identification, Diagnosis, Biology; KNNV Publishing: Utrecht, The Netherlands, 2006.

54. Scheuchl, E. Illustrierte Bestimmungstabellen der Wildbienen. Band I: Anthophoridae; Preisinger KG: Landhut, Germany, 2000.

55. Scheuchl, E.; Willner, W. Taschenlexikon der Wildbienen Mitteleuropas: Alle Arten im Portrait; Quelle & Meyer GmbH: Wiebelsheim, Germany, 2016.

56. Séméria, Y.; Berland, L. Atlas des Néuroptères de France et d´Europe. Mégaloptères, Raphidioptères, Néuroptères Planipennes, Mécoptères; Société Nouvelle des Éditions Boubée: Paris, France, 1988.

57. Vives, E. Coleóptera, Cerambycidae. Fauna Ibérica; Ramos, M.A. Museo Nacional de Ciencias Naturales: Madrid, Spain, 2000; Volume 12.

58. Ives, A.R.; Zhu, J. Statistics for correlated data: Phylogenies, space and time. Ecol. Appl.; 2006; 16, pp. 20-32. [DOI: https://dx.doi.org/10.1890/04-0702]

59. Magurran, A.E. Ecological Diversity and Its Measurement; Chapman and Hall: London, UK, 1991.

60. Soares, S.A.; Antonialli-Junior, W.F.; Limajunior, S.E. Diversidade de formigas epigéicas (Hymenoptera, Formicidae) em dois ambientes no Centro-Oeste do Brasil. Rev. Bras. Entomol.; 2010; 54, pp. 76-81. [DOI: https://dx.doi.org/10.1590/S0085-56262010000100009]

61. Tokeshi, M. Species abundance patterns and community structure. Adv. Ecol. Research; 1993; 24, 111186.

62. Krebs, C.J. Ecological Methodology; Benjamin Cummings: Menlo Park, CA, USA, 1999.

63. Zar, J.H. Biostatistical Analysis; Prentice Hall: Upper Saddle River, NJ, USA, 1999.

64. Danne, A.; Thomson, L.J.; Sharley, D.J.; Penfold, C.M.; Hoffmann, A.A. Effects on native grass cover crops on beneficial and pest invertebrates in Australian vineyards. Pest Manag.; 2010; 39, pp. 970-978. [DOI: https://dx.doi.org/10.1603/EN09144] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20550812]

65. Liu, Y.; Duan, M.; Zhang, X.; Zhang, X.; Yu, Z.; Axmacher, J.C.; Stewart, A.; Ewers, R. Effects of plant diversity, habitat and agricultural landscape structure on the functional diversity of carabid assemblages in the North China Plain. Insect Conserv. Diver.; 2015; 8, pp. 163-176. [DOI: https://dx.doi.org/10.1111/icad.12096]

66. Eckert, M.; Mathulwe, L.L.; Gaigher, R.; Joubert-van der Merwe, L.; Pryke, J.S. Native cover crops enhance arthropod diversity in vineyards of the Cape Floristic Region. J. Ins. Conser.; 2020; 24, pp. 133-149. [DOI: https://dx.doi.org/10.1007/s10841-019-00196-0]

67. Sáenz-Romo, M.G.; Veas-Bernal, A.; Martínez-García, H.; Campos-Herrera, R.; Ibáñez-Pascual, S.; Martínez-Villar, E.; Pérez-Moreno, I.; Marco-Mancebón, V.S. Ground cover management in a Mediterranean vineyard: Impact on insect abundance and diversity. Agric. Ecosyst. Environ.; 2019; 283, 106571. [DOI: https://dx.doi.org/10.1016/j.agee.2019.106571]

68. Sáenz-Romo, M.G.; Veas-Bernal, A.; Martínez-García, H.; Campos-Herrera, R.; Ibáñez-Pascual, S.; Martínez-Villar, E.; Marco-Mancebón, V.S.; Pérez-Moreno, I. Effects of ground cover management on insect predators and pests in a Mediterranean vineyard. Insects; 2019; 10, 421. [DOI: https://dx.doi.org/10.3390/insects10120421]

69. Dassou, A.G.; Dépigny, S.; Canard, E.; Vinatier, F.; Carval, D.; Tixier, P. Contrasting effects of plant diversity across arthropod trophic groups in plantain-based agroecosystems. Basic Appl. Ecol.; 2016; 17, pp. 11-20. [DOI: https://dx.doi.org/10.1016/j.baae.2015.09.003]

70. Von Königslow, V.; Fornoff, F.; Klein, A.-M. Wild bees communities benefit from temporal complementary of hedges and flower strips in apple orchards. J. Appl. Ecol.; 2022; 59, pp. 2814-2824. [DOI: https://dx.doi.org/10.1111/1365-2664.14277]

71. Japon-Quintero, J. La lechuga; Publicaciones de Extension Agraria, Ministerio de Agricultura: Madrid, Spain, 1977.

72. Jil-Martinez, P.; Muñoz-Soto, M. Manejo Productivo de la Cebolla (Allium cepa L.); Instituto de Investigaciones Agropecuarias: Madrid, Spain, 2018; 191.

73. Tarin-Zahonero, J.; Segarra-Dalmau, J. El Cultivo del Apio; Publicaciones de Extension Agraria, Ministerio de Agricultura: Madrid, Spain, 1975.

74. Peris-Felipo, F.J.; Jiménez-Peydró, R. Cerambycidae (Coleoptera) richness in Mediterranean landscapes of Spain: Diversity and community structure analysis. Biodivesity J.; 2012; 3, pp. 59-68.

75. Peris-Felipo, F.J.; Jiménez-Peydró, R. Cerambycidae (Coleoptera) diversity and community structure in the Mediterranean forest of the Natural Park of Sierra Calderona (Spain). Frustula Entom.; 2012; 33, pp. 180-191.

76. Jiménez-Peydró, R.; Peris-Felipo, F.J. Diversity and community structure analysis of Opiinae (Hymenoptera: Braconidae) in the Forest Estate of Artikutza (Spain). Fla. Entomol.; 2011; 94, pp. 472-479. [DOI: https://dx.doi.org/10.1653/024.094.0311]

77. Pérez-Rodríguez, J.; Oltra-Moscardó, M.T.; Peris-Felipo, F.J.; Jiménez-Peydró, R. Microgastrinae (Hymenoptera: Braconidae) in the Forest State of Artikutza (Navarra: Spain): Diversity and community structure. Insects; 2013; 4, pp. 493-505. [DOI: https://dx.doi.org/10.3390/insects4030493] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26462432]

78. Falcó-Garí, J.V.; Peris-Felipo, F.J.; Jiménez-Peydró, R. Diversity and phenology of the braconid community (Hymenoptera: Braconidae) in the Mediterranean protected landscape of Sierra Calderona. Open J. Ecol.; 2014; 4, pp. 174-181. [DOI: https://dx.doi.org/10.4236/oje.2014.44018]

79. Jiménez-Peydró, R.; Peris-Felipo, F.J. Diversity and community structure of Opiinae (Hymenoptera: Braconidae) in Mediterranean landscapes of Spain. J. Entomol. Res. Soc.; 2014; 16, pp. 75-85.

80. Lima, M.; Merritt, J.F.; Bozinovic, F. Numerical fluctuations in the northern short-tailed shrew: Evidence of non-linear feedback signatures on population dynamics and demography. J. Anim. Ecol.; 2002; 71, pp. 159-172. [DOI: https://dx.doi.org/10.1046/j.1365-2656.2002.00597.x]

81. Noordijk, J.; Musters, C.J.M.; van Dijk, J.; de Snoo, G.R. Invertebrates in field margins: Taxonomic group diversity and functional group abundance in relation to age. Biodivers. Conserv.; 2010; 19, pp. 3255-3268. [DOI: https://dx.doi.org/10.1007/s10531-010-9890-1]

© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.