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.     

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.     

Simulation modeling to understand how selective foraging by beaver can drive the structure and function of a willow community

Beaver–willow (Castor-Salix) communities are a unique and vital component of healthy wetlands throughout the Holarctic region. Beaver selectively forage willow to provide fresh food, stored winter food, and construction material. The effects of this complex foraging behavior on the structure and function of willow communities is poorly understood. Simulation modeling may help ecologists understand these complex interactions. In this study, a modified version of the SAVANNA ecosystem model was developed to better understand how beaver foraging affects the structure and function of a willow community in a simulated riparian ecosystem in Rocky Mountain National Park, Colorado (RMNP). The model represents willow in terms of plant and stem dynamics and beaver foraging in terms of the quantity and quality of stems cut to meet the energetic and life history requirements of beaver. Given a site where all stems were equally available, the model suggested a simulated beaver family of 2 adults, 2 yearlings, and 2 kits required a minimum of 4 ha of willow (containing about10 stems m⁻²) to persist in a steady-state condition. Beaver created a willow community where the annual net primary productivity (ANPP) was 2 times higher and plant architecture was more diverse than the willow community without beaver. Beaver foraging created a plant architecture dominated by medium size willow plants, which likely explains how beaver can increase ANPP. Long-term simulations suggested that woody biomass stabilized at similar values even though availability differed greatly at initial condition. Simulations also suggested that willow ANPP increased across a range of beaver densities until beaver became food limited. Thus, selective foraging by beaver increased productivity, decreased biomass, and increased structural heterogeneity in a simulated willow community.     

Multi-metric evaluation of the models WARM,CropSyst, and WOFOST for rice

WARM (Water Accounting Rice Model) simulates paddy rice (Oryza sativa L.), based on temperature-driven development and radiation-driven crop growth. It also simulates: biomass partitioning, floodwater effect on temperature, spikelet sterility, floodwater and chemicals management, and soil hydrology. Biomass estimates from WARM were evaluated and compared with the ones from two generic crop models (CropSyst, WOFOST). The test-area was the Po Valley (Italy). Data collected at six sites from 1989 to 2004 from rice crops grown under flooded and non-limiting conditions were split into a calibration (to estimate some model parameters) and a validation set. For model evaluation, a fuzzy-logic based multiple-metrics indicator (MQI) was used: 0 (best) ≤ MQI ≤ 1 (worst). WARM estimates compared well with the actual data (mean MQI = 0.037 against 0.167 and 0.173 with CropSyst and WOFOST, respectively). On an average, the three models performed similarly for individual validation metrics such as modelling efficiency (EF > 0.90) and correlation coefficient (R > 0.98). WARM performed best in a weighed measure of the Akaike Information Criterion: (worst) 0<wk<1$0<w_k<1$ (best), considering estimation accuracy and number of parameters required to achieve it (mean wk=0.983$w_k=0.983$ against 0.007 and ∼0.000 with CropSyst and WOFOST, respectively). WARM results were sensitive to 30% of the model parameters (ratio being lower with both CropSyst, <10%, and WOFOST, <20%), but appeared the easiest model to use because of the lowest number of crop parameters required (10 against 15 and 34 with CropSyst and WOFOST, respectively). This study provides a concrete example of the possibilities offered using a range of assessment metrics to evaluate model estimates, predictive capabilities, and complexity.     

Modelling llama population development under environmental and market constraints in the Bolivian highlands

Changes in the size of animal populations over time are mainly determined by demographic and environmental factors. Livestock population dynamics are additionally influenced by harvesting decisions taken by herders. In Bolivia, not much is known about current llama husbandry and the main influencing factors determining population sizes. We collected data on demography, environmental factors and market values affecting the current and future llama population in three different regions in Bolivia. We generated a population model and assessed the future development of the llama population including environmental factors (rangeland carrying capacity, disturbance phenomena), herd structure and dynamics, and economic market demands. We calibrated and validated the llama model on the basis of 20-year data sets of the regions of Oruro, Potosi and La Paz, respectively. Model calibration by means of the Gauss-Marquardt-Levenberg algorithm yielded a model efficiency of 0.94. For model validation, however, the simulation slightly overestimated the observed llama population yielding model efficiencies of 0.91 and 0.87 for Potosi and La Paz, respectively. Model outcomes were most sensitive to death and birth rates of juveniles and death rate of females compared to environmental or other demographic factors. Population trajectories approached an overall carrying capacity for Oruro, Potosi and La Paz of 8.8 × 10⁵, 9.1 × 10⁵, and 9.0 × 10⁵ llama individuals after 100 years of simulation. Hence, detailed monitoring of demographic, environmental, and economic factors can improve predictions of llama population development over time. Further management should focus on improving birth rates and lowering female mortality through providing supplemental food and shelters against the harsh environmental conditions of the Andean highlands.     

GIS-based modeling of the geographic distribution of Quercus emoryi Torr. (Fagaceae) in México and identification of significant environmental factors influencing the species’ distribution

The greatest concentration of oak species in the world is believed to be found in Mexico. These species are potentially useful for reforestation because of their capacity to adapt to diverse environments. Knowledge of their geographic distribution and of species–environment relations is essential for decision-making in the management and conservation of natural resources. The objectives of this study were to develop a model of the distribution of Quercus emoryi Torr. in Mexico, using geographic information systems and data layers of climatic and other variables, and to determine the variables that significantly influence the distribution of the species. The study consisted of the following steps: (A) selection of the target species from a botanical scientific collection, (B) characterization of the collecting sites using images with values or categories of the variables, (C) model building with the overlay of images that meet the habitat conditions determined from the characterization of sites, (D) model validation with independent data in order to determine the precision of the model, (E) model calibration through adjustment of the intervals of some variables, and (F) sensitivity analysis using precision and concordance non-parametric statistics applied to pairs of images. Results show that the intervals of the variables that best describe the species’ habitat are the following: altitude from 1650 to 2750 amsl, slope from 0 to 66°; average minimum temperature of January from −12 to −3 °C; mean temperature of June from 11 to 25 °C; mean annual precipitation from 218 to 1225 mm; soil units: lithosol, eutric cambisol, haplic phaeozem, chromic luvisol, rendzina, luvic xerosol, mollic planosol, pellic vertisol, eutric regosol; type of vegetation: oak forest, oak–pine forest, pine forest, pine–oak forest, juniperus forest, low open forest, natural grassland and chaparral. The resulting model of the geographic distribution of Quercus emoryi in Mexico had the following values for non-parametric statistics of precision and agreement: Kappa index of 0.613 and 0.788, overall accuracy of 0.806 and 0.894, sensitivity of 0.650 and 0.825, specificity of 0.963, positive predictive value of 0.945 and 0.957 and negative predictive value of 0.733 and 0.846. Results indicate that the variable average minimum temperature of January, with a maximum value of −3 °C, is an important factor in limiting the species’ distribution.     

Interpreting spatial heterogeneity of crop yield with a process model and remote sensing

A process-based crop growth model (Vegetation Interface Processes (VIP) model) is used to estimate crop yield with remote sensing over the North China Plain. Spatial pattern of the key parameter—maximum catalytic capacity of Rubisco (V_cmax) for assimilation is retrieved from Normalized Difference of Vegetation Index (NDVI) from Terra-MODIS and statistical yield records. The regional simulation shows that the agreements between the simulated winter wheat yields and census data at county-level are quite well with R² being 0.41-0.50 during 2001-2005. Spatial variability of photosynthetic capacity and yield in irrigated regions depend greatly on nitrogen input. Due to the heavy soil salinity, the photosynthetic capacity and yield in coastal region is less than 50 μmol C m⁻² s⁻¹ and 3000 kg ha⁻¹, respectively, which are much lower than that in non-salinized region, 84.5 μmol C m⁻² s⁻¹ and 5700 kg ha⁻¹. The predicted yield for irrigated wheat ranges from 4000 to 7800 kg ha⁻¹, which is significantly larger than that of rainfed, 1500-3000 kg ha⁻¹. According to the path coefficient analysis, nitrogen significantly affects yield, by which water exerts noticeably indirect influences on yield. The effect of water on yield is regulated, to a certain extent, by crop photosynthetic capacity and nitrogen application. It is believed that photosynthetic parameters retrieved from remote sensing are reliable for regional production prediction with a process-based model.     

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

Individual-based modeling of flooding and salinity effects on a coastal swamp forest

Coastal swamps are among the rapidly vanishing wetland habitats in Louisiana. Increased flooding, nutrient and sediment deprivation, and salt-water intrusion have been implicated as probable causes of the decline of coastal swamps. We developed a two-species individual-based forest succession model to compare the growth and composition of a cypress-tupelo swamp under various combinations of flooding intensity and salinity levels, using historical time-series of stage and salinity data as inputs. Our model simulates forest succession over 500 years by representing the growth, mortality, and reproduction of individual Taxodium distichum (baldcypress) and Nyssa aquatica (water tupelo) trees in a 1-km² spatial grid of 10 m × 10 m cells that vary in water levels and salinity through differences in elevation. We independently adjusted the elevations of each cell to obtain different grid-wide mean elevations and standard deviations of elevation; this affected the temporal and spatial pattern of flooding. We calibrated the model by adjusting selected parameters until averaged basal area, stem density and wood production rates under two different mean elevations (partially versus highly flooded) were qualitatively similar to comparable values reported for swamps in the literature. Corroboration involved comparing model predictions to four well-monitored contrasting habitat sites within the Maurepas Basin, Louisiana, USA. Model predictions of both species combined showed the same patterns among sites as the data, but the model overestimated wood production and the dominance of T. distichum. Exploratory simulations predicted that increased flooding leads to swamps with reduced basal areas and stem densities, while increased salinity resulted in lower basal areas at low salinity concentration (∼1-3 psu) and complete tree mortality at higher salinity concentrations (∼2-6 psu). Our model can provide insight into the succession dynamics of coastal swamps and information for the effective design of restoration actions.     

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.     

Application of a spatial meta-population model with stochastic parameters to the management of the invasive grass Nassella trichotoma in North Canterbury, New Zealand

Optimising the management of invasive plants requires the identification of the population size outcomes for alternative management strategies. Mathematical models can be useful tools for making such management strategy comparisons. In this paper we develop a generic landscape meta-population model and apply it to the weedy grass, Nassella trichotoma, an invasive species occupying approximately 800 land parcels, predominantly pastoral farms, in the Hurunui district, North Canterbury, New Zealand. Empirical evidence reveals that this meta-population is currently stable (at a median density of 6 plants ha⁻¹) under a community strategy requiring manual removal (termed ‘grubbing’) of plants annually from all land parcels. Reduction in population size requires an alternative management strategy. Field data, collected over a 12 year period, were used to provide stochastic parameter values for land parcel size, carrying capacity, rates of local population growth and grubbing.The model reveals that at steady state, the most significant contribution to population growth on a land parcel comes from within the land parcel itself; the expected annual per capita growth on an individual land parcel is 4 orders of magnitude greater than the expected annual contribution from plants arising from other land parcels. This result implies that many of the farms currently supporting N. trichotoma may pose little or no threat to, nor are threatened themselves by, other farms infested by the weed. However, the steady state distribution (of the weed's population density) was sensitive to the spread rate, revealing a need for data on this process. It was also sensitive to how any increase in the grubbing rate is distributed; increasing it via a uniform distribution U(0, 1) where all rates between 0 and 100% year⁻¹ are equally probable did not affect the steady state, whereas increasing the rates via the uniform distribution U(0.25, 0.75) resulted in fewer farms with high population densities. In general the model provides a basis for exploring the effects of a wide range of alternative grubbing strategies on population growth in N. trichotoma.     

Evaluating weather effects on interannual variation in net ecosystem productivity of a coastal temperate forest landscape: A model intercomparison

Forest productivity is strongly affected by seasonal weather patterns and by natural or anthropogenic disturbances. However weather effects on forest productivity are not currently represented in inventory-based models such as CBM-CFS3 used in national forest C accounting programs. To evaluate different approaches to modelling these effects, a model intercomparison was conducted among CBM-CFS3 and four process models (ecosys, CN-CLASS, Can-IBIS and 3PG) over a 2500 ha landscape in the Oyster River (OR) area of British Columbia, Canada. The process models used local weather data to simulate net primary productivity (NPP), net ecosystem productivity (NEP) and net biome productivity (NBP) from 1920 to 2005. Other inputs used by the process and inventory models were generated from soil, land cover and disturbance records. During a period of intense disturbance from 1928 to 1943, simulated NBP diverged considerably among the models. This divergence was attributed to differences among models in the sizes of detrital and humus C stocks in different soil layers to which a uniform set of soil C transformation coefficients was applied during disturbances. After the disturbance period, divergence in modelled NBP among models was much smaller, and attributed mainly to differences in simulated NPP caused by different approaches to modelling weather effects on productivity. In spite of these differences, age-detrended variation in annual NPP and NEP of closed canopy forest stands was negatively correlated with mean daily maximum air temperature during July-September (T_amax) in all process models (R² = 0.4-0.6), indicating that these correlations were robust. The negative correlation between T_amax and NEP was attributed to different processes in different models, which were tested by comparing CO₂ fluxes from these models with those measured by eddy covariance (EC) under contrasting air temperatures (T_a). The general agreement in sensitivity of annual NPP to T_amax among the process models led to the development of a generalized algorithm for weather effects on NPP of coastal temperate coniferous forests for use in inventory-based models such as CBM-CFS3: NPP′ = NPP − 57.1 (T_amax − 18.6), where NPP and NPP′ are the current and temperature-adjusted annual NPP estimates from the inventory-based model, 18.6 is the long-term mean daily maximum air temperature during July-September, and T_amax is the mean value for the current year. Our analysis indicated that the sensitivity of NPP to T_amax was nonlinear, so that this algorithm should not be extrapolated beyond the conditions of this study. However the process-based methodology to estimate weather effects on NPP and NEP developed in this study is widely applicable to other forest types and may be adopted for other inventory based forest carbon cycle models.     

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.     

Modelling coupled peach tree-aphid population dynamics and their control by winter pruning and nitrogen fertilization

Nitrogen fertilization and winter pruning are commonly used to control crop production in peach [Prunus persica (L.) Batsch] orchards. They are also known to affect the dynamics of Myzus persicae (Sulzer) (Homoptera: Aphididae) aphid populations via bottom-up regulation processes. Interactions between crops and pests can cause complex system behaviour in response to management practices. An integrated approach will therefore improve the understanding of the effects of these two cultural practices on aphid and peach performances.We developed a simulation model that describes the cultural control of interacting peach tree and aphid population dynamics. It uses the principles of common trophic models while gathering available knowledge and explicit assumptions on peach and aphid functioning and the effects of cultural practices.The model was able to qualitatively reproduce the system behaviour observed in the field. It accounted for actions and feedback such as stimulation of foliar growth by winter pruning, consecutive aphid population increase, subsequent damage to foliage, and partial compensatory growth of foliage. The model also reproduced low losses in fruit production due to aphid infestations. However, it called for further integration of ‘long-term’ effects. Analysis of the model showed the complexity of peach tree and aphid responses to leaf N × winter pruning interactions. Simulations indicated that fruit production losses remained low within a range of realistic values of leaf N and pruning intensity, whereas manipulating peach and aphid dynamics, their interactions and their relationships to practices could result in higher losses.The model is useful to evaluate the relevance of cultural practices for a bottom-up regulation of aphid dynamics in crop-pest management. After considering other control methods and fruit quality, it can be used to find a combination of practices that optimises trade-offs between fruit production and environmental conservation goals. A modelling approach that links crop growth and pest population dynamics and integrates management practice effects has strong potential for improving crop-pest management in an integrated crop production context.     

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.     

The effect of estimated PAR uncertainties on the physiological processes of biosphere models

Photosynthetically active radiation (PAR) energy reaching on the vegetated surface is a key determinant of plant physiological processes. Most of biosphere or crop models use the ratio of PAR to incoming solar radiation (R_s), PAR/R_s, to convert R_s into PAR in order to reduce weather data-input requirements. Several existing models simply specify a constant ratio, PAR/R_s = 0.5. However, some field experiments have reported that the ratio PAR/R_s may not be constant. Previous empirical equations of PAR/R_s were derived based on the data of monthly or daily timescales collected from only a few measurement sites, hence they may not be appropriate to be used in current global biosphere models usually with hourly simulation time steps. Here, we represent the exponential correlation between PAR/R_s and sky clearness index (0-1) using hourly data from 54 Ameriflux measurement sites. It is found that PAR/R_s increases up to 0.6 in cloudy conditions when the clearness index (CI) is below ∼0.2, whereas it is nearly constant at ∼0.42 when CI is above 0.2. When the identified empirical equation is used in the model simulation, it results in −4 to 2% difference in the stomatal conductance compared to that using the constant ratio PAR/R_s = 0.5.     

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.     

Identification of important factors for water vapor flux and CO2 exchange in a cropland

Water vapor flux and carbon dioxide (CO₂) exchange in croplands are crucial to water and carbon cycle research as well as to global warming evaluation. In this study, a standard three-layer feed-forward back propagation neural network technique associated with the Bayesian technique of automatic relevance determination (ARD) was employed to investigate water vapor and CO₂ exchange between the canopy of summer maize and atmosphere in responses to variations of environmental and physiological factors. These factors, namely the photosynthetically active radiation (PAR), air temperature (T), vapor pressure deficient (VPD), leaf-area index (LAI), soil water content in root zone (W), and friction velocity (U*), were used as inputs in neural network analysis. Results showed that PAR, VPD, T and LAI were the primary factors regulating both water vapor and CO₂ fluxes with VPD and W more critical to water vapor flux and PAR and T more crucial to CO₂ exchange. Furthermore, two time variables “day of the year (DOY)” and “time of the day (TOD)” could also improve the simulation results of neural network analysis. The important factors identified by the neural network technique used in this study were in the order of PAR > T > VPD > LAI > U* > TOD for water vapor flux and in the order of VPD > W > LAI > T > PAR > DOY for CO₂ exchange. This study suggests that neural network technique associated with ARD could be a useful tool for identifying important factors regulating water vapor and CO₂ fluxes in terrestrial ecosystem.     

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.