Group defense in prey and hunting cooperation in predators are two important ecological phenomena and can occur concurrently. In this article, we consider cooperative hunting in generalist predators and group defense in prey under a mathematical framework to comprehend the enormous diversity the model could capture. To do so, we consider a modified Holling-Tanner model where we implement Holling type IV functional response to characterize grazing pattern of predators where prey species exhibit group defense. Additionally, we allow a modification in the attack rate of predators to quantify the hunting cooperation among them. The model admits three boundary equilibria and up to three coexistence equilibrium points. The geometry of the nontrivial prey and predator nullclines and thus the number of coexistence equilibria primarily depends on a specific threshold of the availability of alternative food for predators. We use linear stability analysis to determine the types of hyperbolic equilibrium points and characterize the non-hyperbolic equilibrium points through normal form and center manifold theory. Change in the model parameters leading to the occurrences of a series of local bifurcations from non-hyperbolic equilibrium points, namely, transcritical, saddle-node, Hopf, cusp and Bogdanov-Takens bifurcation; there are also occurrences of global bifurcations such as homoclinic bifurcation and saddle-node bifurcation of limit cycles. We observe two interesting closed \'bubble\' form induced by global bifurcations due to change in the strength of hunting cooperation and the availability of alternative food for predators. A three dimensional bifurcation diagram, concerning the original system parameters, captures how the alternation in model formulation induces gradual changes in the bifurcation scenarios. Our model highlights the stabilizing effects of group or gregarious behaviour in both prey and predator, hence supporting the predator-herbivore regulation hypothesis. Additionally, our model highlights the occurrence of \"saltatory equilibria\" in ecological systems and capture the dynamics observed for lion-herbivore interactions.

Grasshoppers represent a significant biological challenge in Inner Mongolia\'s grasslands, severely affecting the region\'s animal husbandry. Thus, dynamic monitoring of grasshopper infestation risk is crucial for sustainable livestock farming. This study employed the Maxent model, along with remote sensing data, to forecast Oedaleus decorus asiaticus occurrence during the growing season, using grasshopper suitability habitats as a base. The Maxent model\'s predictive accuracy was high, with an AUC of 0.966. The most influential environmental variables for grasshopper distribution were suitable habitat data (34.27%), the temperature-vegetation dryness index during the spawning period (18.81%), and various other meteorological and vegetation factors. The risk index model was applied to calculate the grasshopper distribution across different risk levels for the years 2019-2022. The data indicated that the level 1 risk area primarily spans central, eastern, and southwestern Inner Mongolia. By examining the variable weights, the primary drivers of risk level fluctuation from 2019 to 2022 were identified as accumulated precipitation and land surface temperature anomalies during the overwintering period. This study offers valuable insights for future O. decorus asiaticus monitoring in Inner Mongolia.

Interspecific interactions can influence species\' activity and movement patterns. In particular, species may avoid or attract each other through reactive responses in space and/or time. However, data and methods to study such reactive interactions have remained scarce and were generally limited to two interacting species. At this time, the deployment of camera traps opens new opportunities but adapted statistical techniques are still required to analyze interaction patterns with such data. We present the multivariate Hawkes process (MHP) and show how it can be used to analyze interactions between several species using camera trap data. Hawkes processes use flexible pairwise interaction functions, allowing us to consider asymmetries and variations over time when depicting reactive temporal interactions. After describing the theoretical foundations of the MHP, we outline how its framework can be used to study interspecific interactions with camera trap data. We design a simulation study to evaluate the performance of the MHP and of another existing method to infer interactions from camera trap-like data. We also use the MHP to infer reactive interactions from real camera trap data for five species from South African savannas (impala Aepyceros melampus, greater kudu Tragelaphus strepsiceros, lion Panthera leo, blue wildebeest Connochaetes taurinus and Burchell\'s zebra Equus quagga burchelli). The simulation study shows that the MHP can be used as a tool to benchmark other methods of interspecific interaction inference and that this model can reliably infer interactions when enough data are considered. The analysis of real data highlights evidence of predator avoidance by prey and herbivore-herbivore attraction. Lastly, we present the advantages and limits of the MHP and discuss how it can be improved to infer attraction/avoidance patterns more reliably. As camera traps are increasingly used, the multivariate Hawkes process provides a promising framework to decipher the complexity of interactions structuring ecological communities.

Plant life-history traits, such as size and flowering, contribute to shaping variation in herbivore abundance. Although plant genes involved in physical and chemical traits have been well studied, less is known about the loci linking plant life-history traits and herbivore abundance. Here, we conducted a genome-wide association study (GWAS) of aphid abundance in a field population of Arabidopsis thaliana. This GWAS of aphid abundance detected a relatively rare but significant variant on the third chromosome of A. thaliana, which was also suggestively but non-significantly associated with the presence or absence of inflorescence. Out of candidate genes near this significant variant, a mutant of a ribosomal gene (AT3G13882) exhibited slower growth and later flowering than a wild type under laboratory conditions. A no-choice assay with the turnip aphid, Lipaphis erysimi, found that aphids were unable to successfully establish on the mutant. Our GWAS of aphid abundance unexpectedly found a locus affecting plant growth and flowering.

We used European geometrid moths (>630 species) as a model group to investigate how life history traits linked to larval host plant use (i.e., diet breadth and host-plant growth form) and seasonal life cycle (i.e., voltinism, overwintering stage and caterpillar phenology) are related to adult body size in holometabolous insect herbivores. To do so, we applied phylogenetic comparative methods to account for shared evolutionary history among herbivore species. We further categorized larval diet breadth based on the phylogenetic structure of utilized host plant genera. Our results indicate that species associated with woody plants are, on average, larger than herb feeders and increase in size with increasing diet breadth. Obligatorily univoltine species are larger than multivoltine species, and attain larger sizes when their larvae occur exclusively in the early season. Furthermore, the adult body size is significantly smaller in species that overwinter in the pupal stage compared to those that overwinter as eggs or caterpillars. In summary, our results indicate that the ecological niche of holometabolous insect herbivores is strongly interrelated with body size at maturity.

Agricultural productivity relies on plant resistance to insect pests, with silicon (Si) being increasingly recognized as an important anti-herbivore defense. However, the processes by which Si works to counteract the effects of insect injury are not completely understood. The role of Si in mitigating the adverse effects of herbivory has been mostly studied at the species level in various crops, ignoring the sensitivity and variability at the genotypic level. Understanding such variation across genotypes is important because Si-derived benefits are associated with the amount of Si accumulated in the plant. Therefore, the present investigation was pursued to study the effect of different Si concentrations (0, 125, and 250 mg L⁻1) on Si accumulation and plant growth using two wheat genotypes (WW-101 and SW-2) under grasshopper herbivory for 48 h. The higher Si absorption increased the concentration of leaf chlorophyll, carotenoids, soluble sugars, and proteins. Silicon application at higher concentrations increased the dry weight, antioxidant enzyme activity, total phenolics, flavonoids and shoot Si concentration, whereas it decreased the electrolyte leakage, hydrogen peroxide (H2O2) and malonaldehyde (MDA) levels, thereby preventing leaf damage. We infer that the higher Si concentration alleviates the adverse effects of herbivory in wheat by improving the accumulation of secondary metabolites and enhancing the antioxidant defense system. The effects were pronounced in the genotype \'WW-101\' compared to \'SW-2\' for most of the studied traits, indicating overall stress response to be genotype-dependent. Thus, Si acquisition efficiency of genotypes should be considered while developing efficient crop management strategies.

Insect herbivores are amongst the most destructive plant pests, damaging both naturally occurring and domesticated plants. As sessile organisms, plants make use of structural and chemical barriers to counteract herbivores. However, over 75% of herbivorous insect species are well adapted to their host\'s defenses and these specialists are generally difficult to ward off. By actively antagonizing the number of insect eggs deposited on plants, future damage by the herbivore\'s offspring can be limited. Therefore, it is important to understand which plant traits influence attractiveness for oviposition, especially for specialist insects that are well adapted to their host plants. In this study, we investigated the oviposition preference of Pieris butterflies (Lepidoptera: Pieridae) by offering them the choice between 350 different naturally occurring Arabidopsis accessions. Using a genome-wide association study of the oviposition data and subsequent fine mapping with full genome sequences of 164 accessions, we identified WRKY42 and AOC1 as candidate genes that are associated with the oviposition preference observed for Pieris butterflies. Host plant choice assays with Arabidopsis genotypes impaired in WRKY42 or AOC1 function confirmed a clear role for WRKY42 in oviposition preference of female Pieris butterflies, while for AOC1 the effect was mild. In contrast, WRKY42-impaired plants, which were preferred for oviposition by butterflies, negatively impacted offspring performance. These findings exemplify that plant genotype can have opposite effects on oviposition preference and caterpillar performance. This knowledge can be used for breeding trap crops or crops that are unattractive for oviposition by pest insects.

The coexistence of plant-herbivore populations in an ecological system is a fundamental topic of research in mathematical ecology. Plant-herbivore interactions are often described by using discrete-time models in the case of non-overlapping generations: such generations have some specific time interval of life and their old generations are replaced by new generations after some regular interval of time. Keeping in mind the dynamical reliability of continuous-time models we presented two discrete-time plant-herbivore models. Mainly, by applying Euler\'s forward method a discrete-time plant-herbivore model is obtained from a continuous-time plant-herbivore model. In addition, a dynamically consistent discrete-time plant-herbivore model is obtained by applying a nonstandard difference scheme. Moreover, local stability is discussed and the existence of bifurcation about positive equilibrium is shown under some mathematical conditions. To control bifurcation and chaos, a modified hybrid technique is developed. Finally, to support our theocratical results and to show the dynamical reliability of the nonstandard difference scheme some numerical examples are provided.

Under the Ecosystem Exploitation Hypothesis ecosystem productivity predicts trophic complexity, but it is unclear if spatial and temporal drivers of productivity have similar impacts. Long-term studies are necessary to capture temporal impacts on trophic structure in variable ecosystems such as deserts. We sampled ants and measured plant resources in the Simpson Desert, central Australia over a 22-year period, during which rainfall varied 10-fold. We sampled dune swales (higher nutrient) and crests (lower nutrient) to account for spatial variation in productivity. We asked how temporal and spatial variation in productivity affects the abundance of ant trophic guilds. Precipitation increased vegetation cover, with the difference more pronounced on dune crests; seeding and flowering also increased with precipitation. Generalist activity increased over time, irrespective of productivity. Predators were more active in more productive (swale) habitat, i.e. spatial impacts of productivity were greatest at the highest trophic level. By contrast, herbivores (seed harvesters and sugar feeders) increased with long-term rainfall; seed harvesters also increased as seeding increased. Temporal impacts of productivity were therefore greatest for low trophic levels. Whether productivity variation leads to top-down or bottom-up structured ecosystems thus depends on the scale and dimension (spatial or temporal) of productivity.

Rewilding is increasingly considered as an option for environmental regeneration, with potential for enhancing both biodiversity and ecosystem services. So far, however, there is little practical information on how to gauge the benefits and limitations of rewilding schemes on ecosystem composition, structure and functioning. To address this knowledge gap, we explored how satellite remote sensing can contribute to informing the monitoring and evaluation of rewilding projects, using the Knepp estate as a case study. To our knowledge, this study is the first to assess the impacts of rewilding as an ecological regeneration strategy on landscape structure and functioning over several decades. Results show significant changes in land cover distribution over the past 20 years inside rewilded areas in the Knepp estate, with a 41.4% decrease in areas with brown agriculture and grass, a roughly sixfold increase in areas covered with shrubs, and a 40.9% increase in areas with trees; vegetation in the rewilded areas also showed a widespread increase in annual primary productivity. Changes in land cover and primary productivity are particularly pronounced in the part of the estate that began its rewilding journey with a period of large herbivore absence. Altogether, our approach clearly demonstrates how freely available satellite data can (1) provide vital insights about long-term changes in ecosystem composition, structure and functioning, even for small, heterogeneous and relatively intensively used landscapes; and (2) help deepen our understanding of the impacts of rewilding on vegetation distribution and dynamics, in ways that complement existing ground-based studies on the impacts of this approach on ecological communities.