Foliage profiles of individual trees determine competition,self-thinning,biomass and NPP of a Cryptomeria japonica forest stand: A simulation study based on a stand-scale process-based forest model

A simulation study was carried out to investigate simultaneously the effects of eco-physiological parameters on competitive asymmetry, self-thinning, stand biomass and NPP in a temperate forest using an atmosphere–vegetation dynamics interactive model (MINoSGI). In this study, we selected three eco-physiological relevant parameters as foliage profiles (i.e. vertical distribution of leaf area density) of individual trees (distribution pattern is described by the parameter η), biomass allocation pattern in individual tree growth (χ) and the maximum carboxylation velocity (V_max). The position of the maximal leaf area density shifts upward in the canopy with increasing η. For scenarios with η < 4 (foliage concentrated in the lowest canopy layer) or η > 12 (foliage concentrated in the uppermost canopy layer), a low degree of competitive asymmetry was produced. These scenarios resulted in the survival of subordinate trees due to a brighter lower canopy environment when η < 4 or the generation of spatially separated foliage profiles between dominant and subordinate trees when η > 12. In contrast, competition between trees was most asymmetric when 4 ≤ η ≤ 12 (vertically widespread foliage profile in the canopy), especially when η = 8. In such cases, vertically widespread foliage of dominant trees lowered the opportunity of light acquisition for subordinate trees and reduced their carbon gain. The resulting reduction in carbon gain of subordinate trees yielded a higher degree of competitive asymmetry and ultimately higher mortality of subordinate trees. It was also shown that 4 ≤ η ≤ 12 generated higher self-thinning speed, smaller accumulated NPP, litter-fall and potential stand biomass as compared with the scenarios with η < 4 or η > 12. In contrast, our simulation revealed small effects of χ or V_max on the above-mentioned variables as compared with those of η. In particular, it is notable that greater V_max would not produce greater potential stand biomass and accumulated NPP although it has been thought that physiological parameters relevant to photosynthesis such as V_max influence dynamic changes in forest stand biomass and NPP (e.g. the greater the V_max, the greater the NPP). Overall, it is suggested that foliage profiles rather than biomass allocation or maximum carboxylation velocity greatly govern forest dynamics, stand biomass, NPP and litter-fall.     

Projecting population trends of endangered amphibian species in the face of uncertainty: A pattern-oriented approach

Amphibian populations have been declining worldwide for the last three decades. Determining the risk of extinction is one of the major goals of amphibian conservation, yet few quantitative models have been developed for amphibian populations. Like most rare or threatened populations, there is a paucity of life history data available for most amphibian populations. Data on the critical juvenile life stage are particularly lacking. Pattern oriented modeling (POM) has been used successfully to estimate life history parameters indirectly when critical data lacking, but has not been applied to amphibian populations. We describe a spatially explicit, individual-based, stochastic simulation model developed to project population dynamics of pond-breeding amphibian populations. We parameterized the model with life history and habitat data collected for the endangered Houston toad (Bufohoustonensis), a species for which there is a high degree of uncertainty for juvenile and adult male survival. During model evaluation, we focused on explicitly reducing this uncertainty, evaluating 16 different versions of the model that represented the range of parametric uncertainty for juvenile and adult male survival. Following POM protocol, we compared simulation results to four population-level patterns observed in the field: population size, adult sex ratio, proportion of toads returning to their natal pond, and mean maximum distance moved. Based on these comparisons, we rejected 11 of the 16 model versions. Results of the remaining versions confirmed that population persistence depends heavily on juvenile survival, and further suggested that probability of juvenile survival is likely between 0.0075 and 0.015 (previous estimates ranged from 0.003 to 0.02), and that annual male survival is near 0.15 (previous estimates ranged up to 0.43).     

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).     

Sanitary felling of Norway spruce due to spruce bark beetles in Slovenia: A model and projections for various climate change scenarios

A model is presented to predict sanitary felling of Norway spruce (Picea abies) due to spruce bark beetles (Ips typographus, Pityogenes chalcographus) in Slovenia according to different climate change scenarios. The model incorporates 21 variables that are directly or indirectly related to the dependent variable, and that can be arranged into five groups: climate, forest, landscape, topography, and soil. The soil properties are represented by 8 variables, 4 variables define the topography, 4 describe the climate, 4 define the landscape, and one additional variable provides the quantity of Norway spruce present in the model cell. The model was developed using the M5′ model tree. The basic spatial unit of the model is 1 km², and the time resolution is 1 year. The model evaluation was performed by three different measures: (1) the correlation coefficient (51.9%), (2) the Theil's inequality coefficient (0.49) and (3) the modelling efficiency (0.32). Validation of the model was carried out by 10-fold cross-validation. The model tree consists of 28 linear models, and model was calculated for three different climate change scenarios extending over a period until 2100, in 10-year intervals. The model is valid for the entire area of Slovenia; however, climate change projections were made only for the Maribor region (596 km²). The model assumes that relationships among the incorporated factors will remain unchanged under climate change, and the influence of humans was not taken into account. The structure of the model reveals the great importance of landscape variables, which proved to be positively correlated with the dependent variable. Variables that describe the water regime in the model cell were also highly correlated with the dependent variable, with evapotranspiration and parent material being of particular importance. The results of the model support the hypothesis that bark beetles do greater damage to Norway spruce artificially planted out of its native range in Slovenia, i.e., lowlands and soils rich in N, P, and K. The model calculation for climate change scenarios in the Maribor region shows an increase in sanitary felling of Norway spruce due to spruce bark beetles, for all scenarios. The model provides a path towards better understanding of the complex ecological interactions involved in bark beetle outbreaks. Potential application of the results in forest management and planning is discussed.     

A simulation model of rhizome networks for Fallopia japonica (Japanese knotweed) in the United Kingdom

Fallopia japonica (Japanese knotweed) is an aggressively invasive herbaceous perennial that causes substantial economic and environmental damage in the United Kingdom (UK). As such, it is of considerable concern to councils, environmental groups, private landowners and property developers. We construct a 3D correlated random walk model of the development of the subterranean rhizome network for a single stand of F. japonica. The formulation of this model uses detailed knowledge of the morphology and physiology of the plant, both of which differ in the UK to that of its native habitat due to factors including a lack of predation and competition, longer growth seasons and favourable environmental conditions in the UK. Field data obtained as a part of this study are discussed and used in the model for parameterisation and validation. The simulation captures the field data well and predicts, for example, quadratic growth in time for the stand area. Furthermore, the role of a selection of parameters on long-term stand development are discussed, highlighting some key factors affecting vegetative spread rates.     

A distributed-delay model to predict ageing and survival rates of adults of Cryptolestes ferrugineus (Stephens) (Coleoptera: Laemophloeidae) in granaries filled with wheat

A time-varying distributed-delay model simulating effects of multifactors was developed. Prediction of the ageing rate and survival distribution of adults of the rusty grain beetle, Cryptolestes ferrugineus (Stephens) (Coleoptera: Laemophloeidae) in various environments found in wheat-filled granaries was conducted as an example to illustrate the application of this developed model. Published adult mortalities, determined at different temperatures, relative humidities, and food sources, were directly used to find the average ageing rate and family of cumulative function of adult mortality. The developed model could predict the adult survival rate at constant or transient temperatures with different relative humidities. This model could also simulate the effect of adult acclimation to their environment when they experience temperature and moisture fluctuations inside granaries. To validate the developed model, the simulation results were compared with available experimental data from the literature. There was no difference between predicted and measured mortalities in two granaries in which the mortalities were determined in a 4-month experiment.     

A concept model to estimate the potential distribution of the Asiatic citrus psyllid (Diaphorina citri Kuwayama) in Australia under climate change—A means for assessing biosecurity risk

Increasing global temperatures as a result of climate change are widely considered inevitable for Australia. Despite this, the specific effects of climate change on Australian agriculture are little studied and the effects on agricultural pests and diseases are virtually unknown. In this paper we consider the impact of climate change on the Asiatic citrus psyllid (Diaphorina citri Kuwayama [Hemiptera: Psyllidae]); one of two known vectors of huanglongbing (citrus greening); a debilitating disease which is caused in Asia by a phloem-limited bacterium ‘Candidatus Liberibacter asiaticus’ (α-Proteobacteria). D. citri does not occur in Australia, but if introduced would pose a major threat to the viability of the Australian citrus industry and to native Citrus species. This paper presents an approach developed to understand how climate change may influence the behaviour, distribution and breeding potential of D. citri. Here we developed and describe an initial dynamic point model of D. citri biology in relation to its citrus host and applied it to a scenario of increasing temperatures, as indicators of climate change, on a continental scale. A comparison between model outputs for the three time frames considered (1990, 2030 and 2070) confirms that increasing temperatures projected under climate change will affect the timing and duration of new citrus growth (flush) necessary for psyllid development throughout Australia. Flushing will start progressively earlier as the temperature increases and be of shorter duration. There will also be a gradual southward expansion of shorter durations of the occurrence of flush. Increasing temperatures will impact on D. citri both directly through alteration of its temperature dependant development cycle and indirectly through the impact on the host flushing cycle. For the whole of Australia, a comparison between model outputs for the three scenarios considered indicates the seasonality of D. citri development will change to match changes in citrus flush initiation. Results indicate that the risk of establishment by D. citri is projected to decrease under increasing temperatures, mainly due to shortened intervals when it can feed on new leaf flushes of the host. However, the spatially heterogeneous results also suggest that regions located on the southern coastline of Australia could become more suitable for D. citri than projected under current temperatures. These results confirm the value of a linked host-pest approach as based on D. citri climatic requirements alone the model would have accounted only for shorter development periods and predicted an increased risk of potential distribution.     

A population model for the peach fruit moth, Carposina sasakii Matsumura (Lepidoptera: Carposinidae), in a Korean orchard system

A population model for the peach fruit moth, Carposina sasakii Matsumura, was constructed to understand the population dynamics of this pest species and to develop an effective management strategy for various orchard (apple, peach, apple + peach) systems. The model was structured by the five developmental stages of C. sasakii: egg, larva, pupa, larval-cocoon (overwintering larva), and adult. The model consisted of a series of component models: (1) a bimodal spring adult emergence model, (2) an adult oviposition model, (3) stage emergence models of eggs, larvae, and pupae, (4) a larval survival rate model in fruits, (5) a larval-cocoon formation model, and (6) an insecticide effect model. Simulations using the model described the typical patterns of C. sasakii adult abundance in various orchard systems well, and was specific to the composition of host plants: three adult abundance peaks (first peak, mid-season peak, and last peak) a year with decreased peaks after the first peak in monoculture orchards of late apple, two adult peaks a year with a much higher last peak in monoculture orchards of early peach, and three adult peaks a year with much higher later peaks in mixed orchards of late apple and early peach. The average deviation between model outputs and actual records for first and second adult peak dates was 2.8 and 3.9 d, respectively, in simulations without an insecticide effect. The deviation decreased when insecticide effects were incorporated into the model. We also performed a sensitivity analysis of our model, and suggest possible applications of the model.     

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.