Grassy field margins as potential corridors for butterflies in agricultural landscapes: A simulation study

Over the last decades, agricultural intensification has caused a dramatic reduction of grassy habitats. This habitat loss has had a strong negative effect on many meadow-living insect populations, including butterflies. As a part of the cross-compliance measures of the Common Agricultural Policy of the European Union, subsidies for creation and maintenance of grassy field margins (GFM) have been launched. Among other environmental issues, they may serve as corridors for movement of various meadow-living species between individual meadows. Their role as corridors has, however, not yet been demonstrated at the landscape scale and their characteristics that most significantly increase landscape connectivity are unknown. Empirical data for such studies are missing, as the GFM subsidies were launched only 3 years ago. One possibility to get some predictions of their outcomes is provided by simulation models. Here we present our simulation results, using an extension of the model developed by Kindlmann et al. (2004) for the Meadow Brown butterfly, Maniola jurtina. The extension includes the probability to cross a boundary (Conradt and Roper, 2006) that negatively influences dispersal rates but increases sensitivity to the corridor effect. Our simulations show that GFMs increase the dispersal rates between habitat patches and we predict the optimal combinations of width and number of GFMs in the landscape. This way we provide a decision-making tool for increasing landscape connectivity for M. jurtina and similar species. Although our simulations are based on a particular species, they may be generalized because this species shows dispersal rates that are typical of butterfly metapopulations (Conradt et al., 2000), and a potentially widespread dispersal kernel (i.e. “foray search”) that has been reported in a wide variety of species (see Conradt et al., 2003 for a review).     

Simulating the temporal pattern of waste production in farmed gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and Atlantic bluefin tuna (Thunnus thynnus)

Most fish farming waste output models provide gross waste rates as a function of stocked or produced biomass for a year or total culture cycle, but without contemplating the temporality of the discharges. This work aims to ascertain the temporal pattern of waste loads by coupling available growth and waste production models and developing simulation under real production rearing conditions, considering the overlapping of batches and management of stocks for three widely cultured species in the Mediterranean Sea: gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and Atlantic bluefin tuna (Thunnus thynnus). For a similar annual biomass production, the simulations showed that waste output and temporal dumping patterns differ between the three species as a result of the disparities in growth velocity, nutrient digestibility, maintenance metabolic budget and husbandry. The simulations allowed the temporal patterns including the periods of maximum discharge and the dissolved and particulate nitrogen and phosphorus content in the wastes released to be determined, both of which were seen to be species-specific.     

Application of a temperature-dependent von Bertalanffy growth model to bullhead (Cottus gobio)

The most studied and commonly applied model of fish growth is the von Bertalanffy model. However, this model does not take water temperature into account, which is one of the most important environmental factors affecting the life cycle of fish, as many physiological processes that determine growth, e.g. metabolic rate and oxygen supply, are directly influenced by temperature. In the present study we propose a version of the von Bertalanffy growth model that includes mean annual water temperatures by correlating the growth coefficient, k, explicitly and the asymptotic length, L_∞, implicitly to water temperature. All relationships include parameters with an obvious biological relevance that makes them easier to identify. The model is used to fit growth data of bullhead (Cottus gobio) at different locations in the Bez River network (Drme, France). We show that temperature explains much of the growth variability at the different sampling sites of the network.     

Modeling the population dynamics of cotton bollworm Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) over a wide area in northern China

The cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) is one of the most serious crop pests in northern China, calling for accurate prediction of pest outbreaks and strategies for pest control. A computer model is developed to simulate the population dynamics of H. armigera over a wide area in northern China. The area considered covers 12 provinces where serious outbreaks of H. armigera have been observed. In this model, pest development is driven by local ambient temperature, and adults migrate long distances between regions and select preferred hosts for oviposition within a region. Six types of host including cotton, wheat, corn, peanut, soybean and a single category composed of all other minor hosts are considered in this model. Survival rates of eggs and larvae are based on life-table data, and simulated as a function of host type, host phenology and temperature. The incidence of diapause depends on temperature and photoperiod experienced during the larval stage. Survival rate of non-diapause pupae is a nonlinear function of rainfall, and overwinter survival rate is a nonlinear function of temperature. Insecticide is applied when population density exceeds the economic threshold on a host crop within a region. Comparisons of model output with light-trap data indicate that our model reflects the pest population dynamics over a wide area, and could potentially be used for testing novel pest control strategies in northern China.     

Integrating farming techniques in an ecological matrix model: Implementation on the primrose (Primula vulgaris)

Several studies have proven the importance of field margins in sustaining biodiversity and other work has been done on the effect of field management on field margin flora. However few models have been built to predict the effects of field management on the flora. Our project addresses this need for a model capable of predicting the effect of cropping techniques and their timing on the flora of field margins. Primula vulgaris is a biodiversity indicator, characteristic of undisturbed flora and found in field margins and woodlands: its population has been declining for several years. We created a temporal matrix model of P. vulgaris populations on field margins, taking into account the effects of field, field margin and roadside management based on literature and expert knowledge. We then analysed its sensitivity to demographic parameters by comparing lambda (growth rate) sensitivity and elasticity. We compared the management parameter effect using the relative growth rate of the population after 6 years of simulation. Sensitivity analysis to biological parameters showed the importance of adult survival and seed production and germination. Results show that P. vulgaris is particularly sensitive to broad-spectrum herbicides and that other management techniques like early mowing, scything and scrub-killer (diluted broad-spectrum herbicide or specific herbicide) are less aggressive. Our simulations show that management of cash crops in Brittany is too aggressive for P. vulgaris populations and that 4-5 years of grassland in the adjacent field are necessary to maintain populations.     

Agent-based simulation of effects of stress on forest tent caterpillar (Malacosoma disstria Hübner) population dynamics

The forest tent caterpillar (Malacosoma disstria Hübner) (FTC) has an outbreak cycle of approximately 10 years; however, smaller spatial scale analyses show some regions have longer or more frequent periods of high defoliation. This may be a result of local forest fragmentation, pollution or other sources of stress that may affect FTC directly or indirectly through stress on their hosts or parasitoids. Population dynamics of FTC were examined to investigate how stress may alter the severity and frequency of defoliation. We developed a spatially explicit agent-based model to simulate the host-parasitoid dynamics of FTC. Theoretical and empirically derived parameters were established using past literature and over 50 years of population data of FTC from Ontario, Canada. We find that increasing FTC fecundity, FTC dispersal or parasitoid mortality resulted in more severe outbreaks while a decrease in parasitoid fecundity or searching efficiency resulted in an overall elevation of defoliation. Parasitoid efficiency was the most effective parameter for altering the FTC defoliation. Since plant stress has been shown to alter several of these parameters in nature due to changes in food quality, habitat suitability, and chemical cue interference, our results suggest that forests affected by stressors such as climate change and pollution will have more severe and frequent defoliation from these insects than surrounding unaffected forests. As stressors such as drought and pollution emissions are predicted to increase in frequency or intensity over the next few decades, understanding how they may affect the outbreak cycle of a forest defoliator can aid in planning strategies to reduce the detrimental effects of this insect.     

Assessing environmental factors in red spruce (Picea rubens Sarg.) growth in the Great Smoky Mountains National Park, USA: From conceptual model, envirogram, to simulation model

This study provides a method for assessing a multiplicity of environmental factors in red spruce growth in the Great Smoky Mountains National Park (GSMNP) of Southeastern USA. Direct and indirect factors in the annual growth increment are first organized into a schematic input-output envirogram (ARIRS), and this information is then used to construct a simulation model (ARIM). The envirogram represents a structured conceptualization of most environmental factors involved in growth, as developed from relevant literature. This interdisciplinary synthesis distinguishes direct vs. indirect factors in growth and takes account of the systems ecology concept that indirect factors may be as important as or more important than direct ones in regulating growth. The ARIRS envirogram summarizes hierarchically organized, within- and cross-scale, local-to-global interactions, and its construction makes it obvious that growth is influenced by many cross-scale spatiotemporal interactions. More research on genecology is still needed to clarify the role of phenotypic plasticity and adaptive capacity in nutrient cycling, global change, and human disturbance.     

Identifying key life-traits for the dynamics and gene flow in a weedy crop relative: Sensitivity analysis of the GeneSys simulation model for weed beet (Beta vulgaris ssp. vulgaris)

The benefits of genetically modified herbicide-tolerant (GMHT) sugar beet (Beta vulgaris) varieties stem from their presumed ability to improve weed control and reduce its cost, particularly targeting weed beet, a harmful annual weedy form of the genus Beta (i.e. B. vulgaris ssp. vulgaris) frequent in sugar beet fields. As weed beet is totally interfertile with sugar beet, it is thus likely to inherit the herbicide-tolerance transgene through pollen-mediated gene flow. Hence, the foreseeable advent of HT weed beet populations is a serious threat to the sustainability of GM sugar beet cropping systems. For studying and quantifying the long-term effects of cropping system components (crop succession and cultivation techniques) on weed beet population dynamics and gene flow, we developed a biophysical process-based model called GeneSys-Beet in a previous study. In the present paper, the model was employed to identify and rank the weed life-traits as function of their effect on weed beet densities and genotypes, using a global sensitivity analysis to model parameters. Monte Carlo simulations with simultaneous randomization of all life-trait parameters were carried out in three cropping systems contrasting for their risk for infestation by HT weed beets. Simulated weed plants and bolters (i.e. beet plants with flowering and seed-producing stems) were then analysed with regression models as a function of model parameters to rank processes and life-traits and quantify their effects. Key parameters were those determining the timing and success of growth, development, seed maturation and the physiological end of seed production. Timing parameters were usually more important than success parameters, showing for instance that optimal timing of weed management operations is more important than its exact efficacy. The ranking of life-traits though depended on the cropping system and, to a lesser extent, on the target variable (i.e. GM weeds vs. total weed population). For instance, post-emergence parameters were crucial in rotations with frequent sugar beet crops whereas pre-emergence parameters were most important when sugar beet was rare. In the rotations with frequent sugar beet and insufficient weed control, interactions between traits were small, indicating diverse populations with contrasted traits could prosper. Conversely, when sugar beet was rare and weed control optimal, traits had little impact individually, indicating that a small number of optimal combinations of traits would be successful. Based on the analysis of sugar beet parameters and genetic traits, advice for the future selection of sugar beet varieties was also given. In climatic conditions similar to those used here, the priority should be given to limiting the presence of hybrid seeds in seed lots rather than decreasing varietal sensitivity to vernalization.