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.     

COSMOS,a spatially explicit model to simulate the epidemiology of Cosmopolites sordidus in banana fields

A stochastic individual-based model called COSMOS was developed to simulate the epidemiology of banana weevil Cosmopolites sordidus, a major pest of banana fields. The model is based on simple rules of local movement of adults, egg laying of females, development and mortality, and infestation of larvae inside the banana plants. The biological parameters were estimated from the literature, and the model was validated at the small-plot scale. Simulated and observed distributions of attacks were similar except for five plots out of 18, using a Kolmogorov–Smirnov test. These exceptions may be explained by variation in predation of eggs and measurement error. An exhaustive sensitivity analysis using the Morris method showed that predation rate of eggs, demographic parameters of adults and mortality rate of larvae were the most influential parameters. COSMOS was therefore used to test different spatial arrangements of banana plants on the epidemiology of C. sordidus. Planting bananas in groups increased the time required to colonise plots but also the percentage of banana plants with severe attacks. Spatial heterogeneity of banana stages had no effect on time required to colonise plots but increased the mean level of attacks. Our model helps explain key factors of population dynamics and the epidemiology of this tropical pest.     

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.     

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.     

Analysis of annual fluctuations of C. nodosa in the Venice lagoon: Modeling approach

A simple simulation model was developed to describe the growth trends of Cymodocea nodosa (Ucria) Ascherson based on data sets from the Venice lagoon. The model reproduces the seasonal fluctuations in the above and belowground biomass and in shoot density. The modeling results are in good agreement with data on net production, growth rates and chemical–physical parameters of water. It was assumed that light and temperature are the most important factors controlling C. nodosa development, and that the growth was not limited by nutrient availability. The aim was to simulate biomass production as a function of external forcing variables (light, water temperature) and internal control (plant density). A series of simulation experiments were performed with the basic model showing that among the most important phenomena affecting C. nodosa growth are: (1) inhibition of production and recruitment of new shoots by high temperature and (2) light attenuation due to seasonal fluctuation.     

A system-dynamic model on the competitive growth between Potamogeton malaianus Miq. and Spirogyra sp.

A system-dynamic model has been built to evaluate the competition between submerged macrophytes Potamogeton malaianus Miq. (PM) and filamentous green algae Spirogyra sp. (SP). The data background is based on a spring–summer and an autumn–winter experiment carried out in artificial field ponds. The experiments had the aim to acquire a knowledge base necessary to a successful restoration of submerged macrophyte vegetation in Lake Taihu, China by use of P. malaianus Miq. The model mainly focuses on variations in water volume; biomass dynamics of P. malaianus Miq., Spirogyra sp. and zoobenthos; nutrients cycling between water column, P. malaianus Miq., Spirogyra sp., zoobenthos, detritus and sediment. Sixteen state variables are included in the model: biomass of P. malaianus Miq., Spirogyra sp. and zoobenthos; nitrogen in sediments, detritus, in P. malaianus Miq., in Spirogyra sp. and zoobenthos; total dissolved nitrogen; phosphorus in sediments, detritus, in P. malaianus Miq., in Spirogyra sp. and in zoobenthos; total dissolved phosphorus, and water volume of the experiment pond. The calibration and validation of the model show a good accordance with the results of the spring–summer experiment and the autumn–winter experiment.     

Exploring the influence of lake water chemistry on chlorophyll a: A multivariate statistical model analysis

A multivariate statistical approach integrating the absolute principal components score (APCS) and multivariate linear regression (APCS-MLR), along with structural equation modeling (SEM), was used to model the influence of water chemistry variables on chlorophyll a (Chl a) in Lake Qilu, a severely polluted lake in southwestern China. Water quality was surveyed monthly from 2000 to 2005. APCS-MLR was used to identify key water chemistry variables, mine data for SEM, and predict Chl a. Seven principal components (PCs) were determined as eigenvalues >1, which explained 68.67% of the original variance. Four PCs were selected to predict Chl a using APCS-MLR. The results showed a good fit between the observed data and modeled values, with R² = 0.80. For SEM, Chl a and eight variables were used: NH₄-N (ammonia-nitrogen), total phosphorus (TP), Secchi disc depth (SD), cyanide (CN), arsenic (As), cadmium (Cd), fluoride (F), and temperature (T). A conceptual model was established to describe the relationships among the water chemistry variables and Chl a. Four latent variables were also introduced: physical factors, nutrients, toxic substances, and phytoplankton. In general, the SEM demonstrated good agreement between the sample covariance matrix of observed variables and the model-implied covariance matrix. Among the water chemistry factors, T and TP had the greatest positive influence on Chl a, whereas SD had the largest negative influence. These results will help researchers and decision-makers to better understand the influence of water chemistry on phytoplankton and to manage eutrophication adaptively in Lake Qilu.     

Incorporating natural enemies in an economic threshold for dynamically optimal pest management

The control of pests by their natural enemies represents an important regulating ecosystem service that helps maintain the stability of crop ecosystems. These services, however, are often ignored in pest management decision making. In addition, the use of broad-spectrum insecticides can damage the populations of natural enemies, reducing the cost-effectiveness of insecticide investment if unaccounted for in treatment decisions.The existing literature on modeling of biological control of insect pests has generally focused on simulations of the population dynamics of pest and natural enemy species and the processes underlying pest control. But agriculture is a managed ecosystem where predator–prey relationships are heavily influenced by human managers. In modeling managerial choices, this study develops an intra-seasonal dynamic bioeconomic optimization model for insecticide-based pest management that explicitly takes into account both the biological control effect of natural enemies on pest density and the nontarget mortality effect of insecticides on the level of natural pest control supplied. The model captures predator–prey interactions, linking them to crop growth and yield damage functions, which in turn are evaluated in a dynamic optimization framework. We introduce a new decision rule for judicious insecticide decisions using a natural enemy-adjusted economic threshold. This threshold represents the pest population density at which insecticide control becomes optimal in spite of the opportunity cost of injury to natural enemies of the target pest. Using field data from Michigan, the model is applied to the case of soybean aphid (Aphis glycines, Matsumura), a recent invasive pest of soybean (Glycine max), whose management is of both economic and environmental importance to the North Central region of the United States. As illustrated by the numerical examples, such natural enemy-adjusted threshold is likely to lead to fewer recommendations for insecticide use than naïve models that ignore natural enemies, resulting in less insecticide use, while maintaining profitability for farmers that rely on chemical pest control methods.The bioeconomic model developed in this study can be used to conduct a wide variety of analyses such as identifying dynamically optimal spray strategies and estimating the implied economic value of natural control services. Furthermore, with the incorporation of inter-year carry-over factors, such as overwintering of pests and natural enemies, the current model can contribute to building multi-year models for studying long-term pest management.     

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.     

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 amphibian energetics, habitat suitability, and movements of western toads, Anaxyrus (=Bufo) boreas, across present and future landscapes

Effective conservation of amphibian populations requires the prediction of how amphibians use and move through a landscape. Amphibians are closely coupled to their physical environment. Thus an approach that uses the physiological attributes of amphibians, together with knowledge of their natural history, should be helpful. We used Niche Mapper™ to model the known movements and habitat use patterns of a population of Western toads (Anaxyrus (=Bufo) boreas) occupying forested habitats in southeastern Idaho. Niche Mapper uses first principles of environmental biophysics to combine features of topography, climate, land cover, and animal features to model microclimates and animal physiology and behavior across landscapes. Niche Mapper reproduced core body temperatures (T_c) and evaporation rates of live toads with average errors of 1.6 ± 0.4 °C and 0.8 ± 0.2 g/h, respectively. For four different habitat types, it reproduced similar mid-summer daily temperature patterns as those measured in the field and calculated evaporation rates (g/h) with an average error rate of 7.2 ± 5.5%. Sensitivity analyses indicate these errors do not significantly affect estimates of food consumption or activity. Using Niche Mapper we predicted the daily habitats used by free-ranging toads; our accuracy for female toads was greater than for male toads (74.2 ± 6.8% and 53.6 ± 15.8%, respectively), reflecting the stronger patterns of habitat selection among females. Using these changing to construct a cost surface, we also reconstructed movement paths that were consistent with field observations. The effect of climate warming on toads depends on the interaction of temperature and atmospheric moisture. If climate change occurs as predicted, results from Niche Mapper suggests that climate warming will increase the physiological cost of landscapes thereby limiting the activity for toads in different habitats.     

Influence of host migration between woodland and pasture on the population dynamics of the tick Ixodes ricinus: A modelling approach

Ticks act as vectors of pathogens that can be harmful to animals and/or humans. Epidemiological models can be useful tools to investigate the potential effects of control strategies on diseases such as tick-borne diseases. The modelling of tick population dynamics is a prerequisite to simulating tick-borne diseases and the corresponding spread of the pathogen. We have developed a dynamic model to simulate changes in tick density at different stages (egg, larva, nymph and adult) under the influence of temperature. We have focused on the tick Ixodes ricinus, which is widespread in Europe. The main processes governing the biological cycles of ticks were taken into account: egg laying, hatching, development, host (small, mainly rodents, or large, like deer and cattle, mammals) questing, feeding and mortality. This model was first applied to a homogeneous habitat, where simulations showed the ability of the model to reproduce the general patterns of tick population dynamics. We considered thereafter a multi-habitat model, where three different habitats (woodland, ecotone and meadow) were connected through host migration. Based on this second application, it appears that migration from woodland, via the ecotone, is necessary to sustain the presence of ticks in the meadow. Woodland can therefore be considered as a source of ticks for the meadow, which in turn can be regarded as a sink. The influence of woodland on surrounding tick densities increases in line with the area of this habitat before reaching a plateau. A sensitivity analysis to parameter values was carried out and demonstrated that demographic parameters (sex ratio, development, mortality during feeding and questing, host finding) played a crucial role in the determination of questing nymph densities. This type of modelling approach provides insight into the influence of spatial heterogeneity on tick population dynamics.     

Assessing the impact of biological control of Plutella xylostella through the application of Lotka–Volterra model

The Lotka–Volterra model was applied to the population densities of diamondback moth (DBM), Plutella xylostella (L.) and its exotic larval parasitoid Diadegma semiclausum (Hellen) data that was collected earlier by icipe's DBM biological control team. The collections were done for 15 months before the release and 36 months after release of the parasitoid in two areas; in Werugha, Coast Province of Kenya and Tharuni, Central Province of Kenya, respectively. For each area in pre- and post-release periods, we estimated Lotka–Volterra model parameters from the minimization of the loss function between the theoretical and experimental time-series datasets following the Nelder-Mead multidimensional method. The model estimated a reduction in the value of the steady state of DBM population from 4.86 to 2.17 in Werugha and from 6.11 to 3.76 and 3.45 (with and without exclusion of the time before D. semiclausum recovery) in Tharuni when transiting from the pre- and post-release periods, respectively. This change was a consequence of the newly introduced parasitoid, in the areas. The study presented a successful and detailed technique for non-linear model parameters restoration which was demonstrated by the correct mimicking of empirical datasets from the classical biological control with D. semiclausum, in different areas of Kenya. The applied model has measured the parasitoids impact on the DBM biological control through a quantitative estimate of the effectiveness of the newly introduced species D. semiclausum. These equations may therefore be used as tool for decision making in the implementation for such pests’ management system strategy.     

A process-based population dynamics model to explore target and non-target impacts of a biological control agent

The risks and benefits associated with efforts to control invasive alien species using classical biological control are being subjected to increasing scrutiny. A process-based population dynamics model was developed to explore the interactions between a folivorous biological control agent, Cleopus japonicus, and its plant host Buddleja davidii. The model revealed that climate could have a significant impact upon the interactions between B. davidii and C. japonicus. At the coolest sites, the impact of C. japonicus on B. davidii was slowed, but it was still eventually capable of controlling populations of B. davidii. At the warmer sites where both B. davidii and C. japonicus grew faster, B. davidii succumbed rapidly to weevil damage. We hypothesise that barring an encounter with a natural enemy, C. japonicus will eventually be able to provide sustained control B. davidii throughout the North Island of New Zealand. The model scenarios illustrate the potential for the C. japonicus population to attain high densities rapidly, and to defoliate patches of B. davidii, creating the potential for spill-over feeding on non-target plants. The potential magnitude of this threat will depend partly on the climate suitability for C. japonicus, the pattern by which it migrates in response to a reduction in the available leaf resource, and the suitability of non-target plants as hosts. In all migration scenarios considered, the pattern of population growth and resource consumption by C. japonicus was exponential, with a strong tendency toward complete utilisation of resource patches more quickly at the warmer compared to colder sites. In addition to providing some useful hypotheses about the effects of climate on the biological control system, and the non-target risks, it also provides some insight into the mechanisms by which climate affects the system.     

Evaluation of discrete host–parasitoid models for diamondback moth and Diadegma semiclausum field time population density series

The performance of discrete mathematical models to describe the population dynamics of diamondback moth (DBM) (Plutella xylostella L.) and its parasitoid Diadegma semiclausum was investigated. The parameter values for several well-known models (Nicholson–Bailey, Hassell and Varley, Beddington, Free and Lawton, May, Holling type 2, 3 and Getz and Mills functional responses) were estimated. The models were tested on 20 consecutive sets of time series data collected at 14 days interval for pest and parasitoid populations obtained from a highland cabbage growing area in eastern Kenya. Model parameters were estimated from minimized squared difference between the numerical solution of the model equations and the empirical data using Powell's method. Maximum calculated DBM growth rates varied between 0.02 and 0.07. The carrying capacity determined at 16.5 DBM/plant by the Beddington et al. model was within the range of field data. However, all the estimated parameter values relating to the parasitoid, including the instantaneous searching rate (0.07–0.28), per capita searching efficiency (0.20–0.27), search time (5.20–5.33), handling time (0.77–0.90), and parasitism aggregation index (0.33), were well outside the range encountered empirically. All models evaluated for DBM under Durbin–Watson criteria, except the May model, were not autocorrelated with respect to residuals. In contrast, the criteria applied to the parasitoid residuals showed strong autocorrelations. Thus, these models failed to estimate parasitoid dynamics. We conclude that the interactions of the DBM with its parasitoid cannot be explained by any of the models tested. Two factors may be associated with this failure. First, the parasitoid in this integrated biological control system may not be playing a major role in regulating DBM population. Second, and perhaps more likely, poor correlations reflect gross inadequacies in the theoretical assumptions that underlie the existing models.     

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.     

Epidemiological risk assessment using linked network and grid based modelling: Phytophthora ramorum and Phytophthora kernoviae in the UK

We developed a stochastic simulation model incorporating most processes likely to be important in the spread of Phytophthora ramorum and similar diseases across the British landscape (covering Rhododendron ponticum in woodland and nurseries, and Vaccinium myrtillus in heathland). The simulation allows for movements of diseased plants within a realistically modelled trade network and long-distance natural dispersal. A series of simulation experiments were run with the model, representing an experiment varying the epidemic pressure and linkage between natural vegetation and horticultural trade, with or without disease spread in commercial trade, and with or without inspections-with-eradication, to give a 2 × 2 × 2 × 2 factorial started at 10 arbitrary locations spread across England. Fifty replicate simulations were made at each set of parameter values. Individual epidemics varied dramatically in size due to stochastic effects throughout the model. Across a range of epidemic pressures, the size of the epidemic was 5–13 times larger when commercial movement of plants was included. A key unknown factor in the system is the area of susceptible habitat outside the nursery system. Inspections, with a probability of detection and efficiency of infected-plant removal of 80% and made at 90-day intervals, reduced the size of epidemics by about 60% across the three sectors with a density of 1% susceptible plants in broadleaf woodland and heathland. Reducing this density to 0.1% largely isolated the trade network, so that inspections reduced the final epidemic size by over 90%, and most epidemics ended without escape into nature. Even in this case, however, major wild epidemics developed in a few percent of cases. Provided the number of new introductions remains low, the current inspection policy will control most epidemics. However, as the rate of introduction increases, it can overwhelm any reasonable inspection regime, largely due to spread prior to detection.     

Relationships between simulated water stress and mortality and growth rates in underplanted Toona ciliata Roem. in subtropical Argentinean plantations

Toona ciliata Roem. (Australian red cedar) requires a nurse-tree overstory to prevent damage from drought and irradiation in some regions of north-eastern Argentina. T. ciliata was planted in the understory of Pinus taeda L. (625 stems/ha), Pinus elliottii Engelm. × Pinus caribaea Morelet (625 stems/ha), and Grevillea robusta A. Cunn. (833 stems/ha) nurse trees, which were thinned to 0, 25, 50, 75 and 100% of the initial densities. We measured initial T. ciliata mortality and growth as well as Leaf Area Index (LAI) based on light transmission. T. ciliata soil water availability and its effect on early growth and mortality were examined by modelling drought stress using the two-dimensional forest hydrology model ForWaDy. Simulated patterns in T. ciliata water stress for the different overstory treatments were consistent with observed patterns of mortality. Early mortality was lowest with a G. robusta overstory, with corresponding lowest drought stress values and high modelled soil water contents in the top soil layer in intermediate and high overstory densities. Mortality was highest with a P. elliottii × P. caribaea overstory in treatments with the highest modelled drought stress values in the most open treatments. The model supported our field observations by indicating that water stress was an important limitation to T. ciliata survival and growth on our study sites. The linkage between T. ciliata establishment success, early growth and soil water availability as indicated by ForWaDy, leads us to conclude that the model is a suitable stand management tool for guiding establishment of T. ciliata plantations.     

Application of a Random Forest algorithm to predict spatial distribution of the potential yield of Ruditapes philippinarum in the Venice lagoon, Italy

We present a modelling framework that combines machine learning techniques and Geographic Information Systems to support the management of an important aquaculture species, Manila clam (Ruditapes philippinarum). We use the Venice lagoon (Italy), the first site in Europe for the production of R. philippinarum, to illustrate the potential of this modelling approach. To investigate the relationship between the yield of R. philippinarum and a set of environmental factors, we used a Random Forest (RF) algorithm. The RF model was tuned with a large data set (n = 1698) and validated by an independent data set (n = 841). Overall, the model provided good predictions of site-specific yields and the analysis of marginal effect of predictors showed substantial agreement among the modelled responses and available ecological knowledge for R. philippinarum. The most influent environmental factors for yield estimation were percentage of sand in the sediment, salinity, and water depth. Our results agree with findings from other North Adriatic lagoons. The application of the fitted RF model to continuous maps of all the environmental variables allowed estimates of the potential yield for the whole basin. Such a spatial representation enabled site-specific estimates of yield in different farming areas within the lagoon. We present a possible management application of our model by estimating the potential yield under the current farming distribution and comparing it to a proposed re-organization of the farming areas. Our analysis suggests a reduction of total yield is likely to result from the proposed re-organization.     

A model for measuring leaf growth rate of vegetative shoot in the seagrass, Zostera marina

A new model for determining leaf growth in vegetative shoots of the seagrass Zostera marina (eelgrass) is described. This model requires the weights of individual mature and immature whole leaves and leaf plastochrone interval (P_L) as parameters, differing from the conventional leaf marking technique (CLM) that requires cutting and separation between new and old tissue of leaves. The techniques required for the model are the same as for the plastochrone method, but the parameters differ between both methods in use of the weight of individual immature leaves. In a mesocosm study, eelgrass growth was examined, and parameters for the new model and plastochrone method (the weights of individual mature and immature leaves and P_L) were measured. Leaf growth rate was measured using the CLM and determined by the new method and the plastochrone method. The results were then compared between the CLM, the new model, and the plastochrone method. The results obtained with the new model were similar to those obtained with the CLM. However, the results of the plastochrone method differed from those of the CLM, while the weight of immature leaves varied seasonally. The new model was also used to determine leaf growth in a natural eelgrass bed in Mikawa Bay, Japan, and revealed the growth rates in all shoots and those of different ages. This method would be advantageous as an accurate means of direct measurement in fieldwork, and should therefore be a useful tool for monitoring seagrass growth.