A comparison of four process-based models and a statistical regression model to predict growth of Eucalyptus globulus plantations

In forest management and ecological research, consideration of the impacts and risks of climate change or management optimisation is complex. Computer models have long been applied as tools for these tasks. Process-based forest growth models claim to overcome the limitations of empirical statistical models, but the capacity of different process-based models and modelling approaches have rarely been compared directly. This study evaluates stepwise multiple regression models in comparison to four process-based modelling approaches (3-PG, 3-PG+, CABALA and Forest-DNDC) for greenfield predictions of Eucalyptus globulus plantation growth from 2 to 8 years after planting throughout southern Australia.     

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.     

A discrete model for estimating the development time from egg to infecting larva of Ostertagia ostertagi parametrized using a fuzzy rule-based system

Ostertagia ostertagi is a nematode, predominantly affecting cattle in the Pampean region of Argentina. A mathematical model parametrized using fuzzy rule-based systems of the Takagi-Sugeno-Kant type (FTSK) for estimating the development time from egg to infecting larval stage L3 of the gastrointestinal parasite O. ostertagi is here proposed. The estimation of development time of O. ostertagi is essential for the generation of appropriate control mechanisms, since this provides information about the time when parasites are ready to migrate to pastures. For the purpose of reflecting the natural environmental conditions, the mean daily temperature is taken as the main and only regulator of the development time. Humidity conditions are considered to be sufficient for the normal development of the larvae. Hence the individual's daily growth is a function of its length and the mean temperature recorded on the previous day. It is expressed in terms of a difference equation with fuzzy parameters, which are defined using laboratory data. Model outputs are tested against results of field experiments. Simulation results are very satisfactory, yielding a mean estimation error (MEE) of 0.64 weeks, with variance 0.34, and a determination coefficient R² = 0.74. The model clearly exhibits an inverse relationship between development time and temperature both in controlled and in field conditions. It also exhibits a very sensitive response both to the order in which the temperature sequence occurs, - reproducing the differences observed between spring and autumn - and to the amplitude of the temperature range.     

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.     

Estimating Plecoglossus altivelis altivelis migration using a mass balance model expressed by hydrological distribution parameters in a major limpid river basin in Japan

We examined the influence of several hydrological and meteorological parameters on the migratory movements of ayu Plecoglossus altivelis altivelis in central Japan. When comprehensively evaluating rivers and ayu behaviour on a catchment scale, the subjects of analysis typically include human activities and hydrological and meteorological phenomena. However, limiting analyses to such factors may be too restrictive when human activities are being conducted. Accordingly, we incorporated a biological viewpoint into the evaluation method, analysing hydrological data (river discharge, river water temperature, sea water temperature) to determine watershed characteristics and examining the relationship between these characteristics and the habitat conditions of ayu. Then we constructed a numerical model for ayu migratory runs that incorporated ayu ecology and watershed characteristics. Analyses of ayu movements from a lower estuarine dam demonstrated that downstream displacements were associated with high water flows of more than 200 m³ s⁻¹ at the beginning of summer. We conclude that it is important to consider the effects of environmental parameters on the movements of different fish species to understand the causes of spatial variation in fish distribution in lowland rivers.     

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.     

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.     

Computer simulation applied to the biological control of the insect Aphis gossypii for the parasitoid Lysiphlebus testaceipes

In integrated pest management (IPM), biological control is one of the possible options for the prevention or remediation of an unacceptable pest activity or damage. The success of forecast models in IPM depends, among other factors, on the knowledge of temperature effect over pests and its natural enemies. In this work, we simulated the effects of parasitism of Lysiphlebus testaceipes (Cresson, 1880) (Hymenoptera: Aphidiidae) on Aphis gossypii (Glover, 1877) (Hemiptera: Aphididae), a pest that is associated to crops of great economic importance in several parts of the world. We made use of experimental data relative to the host and its parasitoid at different temperatures. Age structure was incorporated into the dynamics through the Penna model. The results obtained showed that simulation, as a forecast model, can be a useful tool for biological control programs.     

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.     

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.     

Sea level rise and eelgrass (Zostera marina) production: A spatially explicit relative elevation model for Padilla Bay, WA

The dynamics that govern the elevation of a coastal wetland relative to sea level are complex, involving non-linear feedbacks among opposing processes. Changes in the balance between these processes can result in significant alterations to vegetation communities that are adapted to a specific range of water levels. Given that current sedimentation rates in Padilla Bay, Washington are likely less than historical levels and that eustatic sea level rise is accelerating, the extensive Zostera marina (eelgrass) meadows in the bay may be at risk of eventual submergence. We developed a spatially explicit relative elevation model and used it to project changes in the productivity and distribution of eelgrass in Padilla Bay over the next century. The model is mechanistic and incorporates many of the processes and feedbacks that govern coastal wetland elevation change. Accretion estimates made using ²¹⁰Pb dating of sediment cores, sediment characteristics measured within cores, and eelgrass productivity and decomposition data were used to initialize and calibrate the model. Validation was performed using an elevation change rate measured with a network of surface elevation tables. Both the field data and model simulations revealed a net accretion deficit for the bay. Simulations using current rates of sea level rise indicated an overall expansion of eelgrass within Padilla Bay over the next century as it migrates from the center of the bay shoreward.     

On the sudden disappearance of Egeria densa from a Ramsar wetland site of Southern Chile: A climatic event trigger model

Contemporary shallow lakes theory proposes that these ecosystems may experience abrupt regime shifts due to small changes in controlling variables or triggers. So far, these triggers have been related mostly to nutrients as the immediate driver. During May 2004 the río Cruces wetland, a Ramsar site located in Southern Chile, underwent a major regime shift, from a clear water state, vastly dominated by the invasive macrophyte Egeria densa, to a turbid water state. In this article we show, through the analysis of long-term meteorological data that late fall 2004 was anomalous due to the presence of a high-pressure cell that persisted most of the month of May over Southern Chile. This climatic event caused an almost complete absence of precipitations and lower temperatures during this period, including several freezing nights. Eco-physiological experiments showed that 6 h exposure to desiccation kill the macrophyte. We developed a simple-biology dynamic model, under Stella Research 9.1, to show that the climatic anomaly of May 2004, plus the increased sedimentation of the wetland's floodplains, and the associated response of E. densa, explains its sudden disappearance from río Cruces wetland.     

Allowing macroalgae growth forms to emerge: Use of an agent-based model to understand the growth and spread of macroalgae in Florida coral reefs,with emphasis on Halimeda tuna

The growth patterns of macroalgae in three-dimensional space can provide important information regarding the environments in which they live, and insights into changes that may occur when those environments change due to anthropogenic and/or natural causes. To decipher these patterns and their attendant mechanisms and influencing factors, a spatially explicit model has been developed. The model SPREAD (SPatially-explicit Reef Algae Dynamics), which incorporates the key morphogenetic characteristics of clonality and morphological plasticity, is used to investigate the influences of light, temperature, nutrients and disturbance on the growth and spatial occupancy of dominant macroalgae in the Florida Reef Tract. The model species, Halimeda and Dictyota spp., are modular organisms, with an “individual” being made up of repeating structures. These species can also propagate asexually through clonal fragmentation. These traits lead to potentially indefinite growth and plastic morphology that can respond to environmental conditions in various ways. The growth of an individual is modeled as the iteration of discrete macroalgal modules whose dynamics are affected by the light, temperature, and nutrient regimes. Fragmentation is included as a source of asexual reproduction and/or mortality. Model outputs are the same metrics that are obtained in the field, thus allowing for easy comparison. The performance of SPREAD was tested through sensitivity analysis and comparison with independent field data from four study sites in the Florida Reef Tract. Halimeda tuna was selected for initial model comparisons because the relatively untangled growth form permits detailed characterization in the field. Differences in the growth patterns of H. tuna were observed among these reefs. SPREAD was able to closely reproduce these variations, and indicate the potential importance of light and nutrient variations in producing these patterns.     

Modeling the role of macroalgae in a shallow sub-estuary of Narragansett Bay, RI (USA)

Proliferation of macroalgal mats is a frequent consequence of nutrient-driven eutrophication in shallow, photic coastal marine ecosystems. These macroalgae have the potential to significantly modify water quality, plankton productivity, nutrient cycling, and dissolved oxygen dynamics. We developed a model for Ulva lactuca and Gracilaria tikvahiae in Greenwich Bay, RI (USA), a shallow sub-estuary of Narragansett Bay, as part of a larger estuarine ecosystem model. The model predicts the biomass of both species in units of carbon, nitrogen, and phosphorus as a function of primary production, respiration, grazing, decay, and physical exchange, with particular attention to the effects of biomass layering on light attenuation and suppression of metabolic rates. The model successfully reproduced the magnitude and seasonal cycle of area-weighted and peak biomass in Greenwich Bay along with tissue C:N ratios, and highlighted the importance of grazing and inclusion of self-limitation primarily in the form of self-shading to overcome an order of magnitude difference in rates of production and respiration. Inclusion of luxury nutrient uptake demonstrated the importance of internal nutrient storage in fueling production when nutrients are limiting. Macroalgae were predicted to contribute a small fraction of total system primary production and their removal had little effect on predicted water quality. Despite a lack of data for calibration and a fair amount of sensitivity to individual parameter values, which highlights the need for further autecological studies to constrain formulations, the model successfully predicted macroalgal biomass dynamics and their role in ecosystem functioning. Our formulations should be exportable to other temperate systems where macroalgae occur in abundance.     

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.     

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.     

An individual-based model to study the reproduction of egg bearing copepods: Application to Eurytemora affinis (Copepoda Calanoida) from the Seine estuary,France

Limited empirical studies have elucidated the daily egg production and associated reproductive processes of egg bearing copepod. Herein, we present an individual-based model which constitutes a realistic representation of the reproduction in egg bearing copepods. The model has been parameterized using an extensive set of experimental data obtained from the literature and from the laboratory and field experiments on the estuarine copepod Eurytemora affinis. The proposed model takes into account the adult female longevity, the clutch size and interclutch duration, which is a function of egg maturation time and latency time required by the female after egg hatching to produce a new clutch. The embryonic development time and hatching success are also taken into account. The effect of temperature on the means and variances of above-mentioned reproductive parameters has been also incorporated. A multi agent system based generic platform “Mobidyc” has been used to generate and calibrate the model. The model demonstrates the reproductive parameters of females of E. affinis which is validated through individual based experiments. Temperature specific simulations provide a dynamical explanation of temperature effect on the cumulative egg production. The daily survival principally affects the number of clutches produced per female during its life span. The results obtained in the present study by combining temperature and survival effects reveal the relatively greater importance of the first factor on the daily egg production of egg-carrying copepods. The present model is generic and hence easily applicable to other animals with comparable reproductive strategy.     

A system dynamics model for the management of the Manila clam, Ruditapes philippinarum (Adams and Reeve, 1850) in the Bay of Arcachon (France)

The Manila clam Ruditapes philippinarum (Adams and Reeve, 1850) is one of the mollusc species that, driven mainly by the shellfish market industry, has extended throughout the world, far beyond the limits of its original habitat. The Manila clam was introduced into France for aquaculture purposes, between 1972 and 1975. In France, this venerid culture became increasingly widespread and, since 1988, this species has colonised most of the embayments along the French Atlantic coast. In 2004, this development resulted in a fishery of ca. 520 t in Arcachon Bay.     

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.