Ecological relevance of performance criteria for species distribution models

Species distribution models have often been developed based on ecological data. To develop reliable data-driven models, however, a sound model training and evaluation procedures are needed. A crucial step in these procedures is the assessment of the model performance, with as key component the applied performance criterion. Therefore, we reviewed seven performance criteria commonly applied in presence-absence modelling (the correctly classified instances, Kappa, sensitivity, specificity, the normalised mutual information statistic, the true skill statistic and the odds ratio) and analysed their application in both the model training and evaluation process. Although estimates of predictive performance have been used widely to assess final model quality, a systematic overview was missing because most analyses of performance criteria have been empirical and only focused on specific aspects of the performance criteria. This paper provides such an overview showing that different performance criteria evaluate a model differently and that this difference may be explained by the dependency of these criteria on the prevalence of the validation set. We showed theoretically that these prevalence effects only occur if the data are inseparable by an n-dimensional hyperplane, n being the number of input variables. Given this inseparability, different performance criteria focus on different aspects of model performance during model training, such as sensitivity, specificity or predictive accuracy. These findings have important consequences for ecological modelling because ecological data are mostly inseparable due to data noise and the complexity of the studied system. Consequently, it should be very clear which aspect of the model performance is evaluated, and models should be evaluated consistently, that is, independent of, or taking into account, species prevalence. The practical implications of these findings are clear. They provide further insight into the evaluation of ecological presence/absence models and attempt to assist modellers in their choice of suitable performance criteria.     

Comparison of sensitivity analysis techniques: A case study with the rice model WARM

The considerable complexity often included in biophysical models leads to the need of specifying a large number of parameters and inputs, which are available with various levels of uncertainty. Also, models may behave counter-intuitively, particularly when there are nonlinearities in multiple input-output relationships. Quantitative knowledge of the sensitivity of models to changes in their parameters is hence a prerequisite for operational use of models. This can be achieved using sensitivity analysis (SA) via methods which differ for specific characteristics, including computational resources required to perform the analysis. Running SA on biophysical models across several contexts requires flexible and computationally efficient SA approaches, which must be able to account also for possible interactions among parameters. A number of SA experiments were performed on a crop model for the simulation of rice growth (Water Accounting Rice Model, WARM) in Northern Italy. SAs were carried out using the Morris method, three regression-based methods (Latin hypercube sampling, random and Quasi-Random, LpTau), and two methods based on variance decomposition: Extended Fourier Amplitude Sensitivity Test (E-FAST) and Sobol’, with the latter adopted as benchmark. Aboveground biomass at physiological maturity was selected as reference output to facilitate the comparison of alternative SA methods. Rankings of crop parameters (from the most to the least relevant) were generated according to sensitivity experiments using different SA methods and alternate parameterizations for each method, and calculating the top-down coefficient of concordance (TDCC) as measure of agreement between rankings. With few exceptions, significant TDCC values were obtained both for different parameterizations within each method and for the comparison of each method to the Sobol’ one. The substantial stability observed in the rankings seem to indicate that, for a crop model of average complexity such as WARM, resource intensive SA methods could not be needed to identify most relevant parameters. In fact, the simplest among the SA methods used (i.e., Morris method) produced results comparable to those obtained by methods more computationally expensive.     

Measuring q-bits in three-trophic level systems

The use of quantum information has been proposed as an approach to deal with biological data (Piqueira, J.R.C., Serboncini, F.A., Monteiro, L.H.A., 2006. Biological models: measuring variability with classical and quantum information. J. Theor. Biol. 242 (2), 309–313). Using three-trophic level systems as examples, we show how to model population data by expressing the system states with q-bits. The system time evolution is given by the state transition matrices which relate the states to successive time intervals. It is a complementary way of looking at the problem which is usually modeled with deterministic differential equations. This is possible because the dynamics of interacting populations in three-trophic level systems is a problem with several coupled variables and, consequently, complex dynamical behaviors seem to result. The non deterministic dynamics generated by the state transition matrices is supposed to model the biological system as a whole, with real data expressing even the global effects of small disturbances in the ecological parameters.     

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.     

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.     

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.     

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.     

A parsimonious optimal foraging model explaining mortality patterns in Serengeti wildebeest

Based on data collected over 24 years in the Serengeti in Tanzania, Sinclair and Arcese (1995) indicated that the sensitivity of blue wildebeest Connochaetes taurinus to predation risk by lions Panthera leo may cause them to change habitats between open (low risk) and wooded (risky) habitats. They found that, in poor rainfall years, predators kill wildebeest that are in better condition than those that die of natural causes. In good rainfall years, predators kill wildebeest that are in worse condition than those that die of natural causes. Sinclair and Arcese (1995) proposed the “predation-sensitive food” hypothesis. This hypothesis suggests that, as food becomes limiting, animals take greater risks to obtain more food, and some of these animals are killed. I propose a more parsimonious hypothesis based on the marginal value theorem that is consistent with the observations made by Sinclair and Arcese (1995). Wildebeest follow a single decision rule in good and poor rainfall years, viz. move when foraging elsewhere increases your rate of intake of nutritious food. Similarly, predators follow a single decision rule in good and poor rainfall years, viz. take the prey item that maximizes the intake of energy per unit effort expended. This parsimonious model does not require differences in predator sensitivity as required by Sinclair and Arcese's (1995) model. I indicate ways in which my model can be falsified.     

Individual heterogeneity in recapture probability and survival estimates in cheetah

Accurate estimates of demographic parameters are key for understanding and predicting population dynamics and for providing insights for effective wildlife management. Up until recently, no suitable methodology has been available to estimate survival probabilities of species with asynchronous reproduction and a high level of individual variation in capture probabilities. The present work develops a capture-mark-recapture model for cheetahs in the Serengeti National Park, Tanzania, which (a) deals with continuous reproduction, (b) takes into account the high level of individual heterogeneity in capture probabilities and (c) is spatially explicit. Results show that (1) our approach, which is an extensive modification of the Cormack-Jolly-Seber model, provides a lower female adult survival estimate and a higher male adolescent survival estimate than previous approaches to estimate cheetah survival in the area, (2) using sighting location alone is not sufficient to capture the individual variation in resighting probabilities for both sexes, and (3) precision in estimated survival probabilities is generally increased. Species which are individually recognizable, wide-ranging and/or where individuals differ substantially in sightability are particularly appropriate to our modelling approach, and our methodology would thus be appropriate for a wide number of species to provide more accurate estimates of survival.     

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.     

Assessing self-organization of plant communities—A thermodynamic approach

Thermodynamics is a powerful tool for the study of system development and has the potential to be applied to studies of ecological complexity. Here, we develop a set of thermodynamic indicators including energy capture and energy dissipation to quantify plant community self-organization. The study ecosystems included a tropical seasonal rainforest, an artificial tropical rainforest, a rubber plantation, and two Chromolaena odorata (L.) R.M. King & H. Robinson communities aged 13 years and 1 year. The communities represent a complexity transect from primary vegetation, to transitional community, economic plantation, and fallows and are typical for Xishuangbanna, southwestern China. The indicators of ecosystem self-organization are sensitive to plant community type and seasonality, and demonstrate that the tropical seasonal rainforest is highly self-organized and plays an important role in local environmental stability via the land surface thermal regulation. The rubber plantation is at a very low level of self-organization as quantified by the thermodynamic indicators, especially during the dry season. The expansion of the area of rubber plantation and shrinkage of tropical seasonal rainforest would likely induce local surface warming and a larger daily temperature range.     

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.     

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.     

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.     

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.     

Simulating regeneration and vegetation dynamics in Mediterranean coniferous forests

This study aims to provide a quantitative framework to model the dynamics of Mediterranean coniferous forests by integrating existing ecological data within a generic mathematical simulator. We developed an individual-based vegetation dynamics model, constrained on long-term field regeneration data, analyses of tree-rings and seed germination experiments. The simulator implements an asymmetric competition algorithm which is based on the location and size of each individual. Growth is parameterized through the analysis of tree-rings from more than thirty individuals of each of the three species of interest. A super-individual approach is implemented to simulate regeneration dynamics, constrained with available regeneration data across time-since-disturbance and light-availability gradients. The study concerns an insular population of an endemic to Greece Mediterranean fir (Abies cephalonica Loudon) on the island of Cephalonia (Ionian Sea) and two interacting populations of a Mediterranean pine (Pinus brutia Ten.) and a more temperate-oriented pine (Pinus nigra Arn. ssp. pallasiana) on the island of Lesbos (NE Aegean Sea), Greece. The model was validated against plot-level observations in terms of species standing biomass and regeneration vigour and adequately captured regeneration patterns and overall vegetation dynamics in both study sites. The potential effects of changing climatic patterns on the regeneration dynamics of the three species of interest were subsequently explored. With the assumption that a warmer future would probably cause changes in the duration of cold days, we tested how this change would affect the overall dynamics of the study sites, by focusing on the process of cold stratification upon seed germination. Following scenarios of a warmer future and under the current model parameterization, changes in the overall regeneration vigour controlled by a reduction in the amount of cold days, did not alter the overall dynamics in all plant populations studied. No changes were identified in the relative dominance of the interacting pine populations on Lesbos, while the observed reduction in the amount of emerging seedlings of A. cephalonica on Cephalonia did not affect biomass yield at later stages of stand development.     

New hydroepidemiological models of indicator organisms and zoonotic pathogens in agricultural watersheds

Simple analytical models are derived to assess how a series of cattle animal farms affect the transport and fate of an indicator organism (Escherichia coli) and a zoonotic pathogen (Campylobacter) in a stream. Separate steady-state mass-balance models are developed and solved for the ultimate minimum and maximum concentrations for the two organisms. The E. coli model assumes that the organism is ubiquitous and abundant in the animals’ digestive tracts. In contrast, a simple dose-response model is employed to relate the Campylobacter prevalence to drinking water drawn from the stream. Because faecal indicators are commonly employed to assess the efficacy of best management practice (BMP) interventions, we also employ the models to assess how BMPs impact pathogen levels. The model provides predictions of (a) the relative removal efficacy for Campylobacter and (b) the prevalence of Campylobacter infection among farm animals after implementation of BMPs. Dimensionless numbers and simple graphs are developed to assess how prevalence is influenced by a number of factors including animal density and farm spacing. A significant outcome of this model development is that the numerous dimensional input and parameter variables are reduced to a group of just four dimensionless Campylobacter-related quantities, characterizing: animal density; in-stream attenuation; animal-to-animal transmission; and infection recovery. Calculations reveal that for some constellations of these four quantities there can be a greater-than-expected benefit in that the proportional reduction of stream Campylobacter concentrations post-BMP can substantially exceed the proportional reduction of concentrations of E. coli in that stream. In addition, a criterion for system sterility (i.e., the conditions required for the farm infection rate to decrease with downstream distance) is derived.     

Modelling N mineralization from green manure and farmyard manure from a laboratory incubation study

Predicting N mineralization from organic manures like farmyard manure (FYM) is more difficult than from fresh organic materials like crop residues, as the manures vary greatly in composition. A laboratory incubation experiment was carried out for 98 days at 30 °C under aerobic conditions to study the effects on N dynamics of Gliricidia (Gliricidia sepium, Jacquin) and FYM application to soil at 5 and 10 g kg⁻¹. Application of Gliricidia induced N mineralization from the start of incubation period, with the amount of N mineralized increasing with rate of application. In contrast, application of FYM resulted in immobilization of mineral N in soil, irrespective of the rate of application. The initial net immobilization from FYM was limited by availability of N in the soil for the higher rate of application.We used the APSIM SoilN module to simulate these contrasting patterns of mineralization of N from Gliricidia and from FYM. The prediction of N mineralized from Gliricidia was better than FYM. The default model parameters specify that the fresh organic matter pools (FPOOL1, FPOOL2 and FPOOL3) have the same C:N ratio and this assumption was ineffective in predicting N mineralized from FYM. The predictive ability of the model improved when this default assumption was modified based on the size of the individual pools (FPOOL1, FPOOL2 and FPOOL3), and the pool's C:N ratios. The modelling efficiency, a measure of goodness of fit between the simulated and observed data, improved markedly for the modified model. The discrepancy between the modelled and observed data was a tendency for the model to underestimate the rate of re-mineralization at the lower rate of application of FYM in the later part of incubation. Unfortunately the appropriate modification to the size and C:N ratios of the FPOOLs could not be determined on the basis of chemical analysis alone. Thus, a true predictive application of the model to a new FYM material is not yet possible.     

Determining the community structure of the coral Seriatopora hystrix from hydrodynamic and genetic networks

The exchange of genetic information between coral reefs through the transport of larvae can be described in terms of networks that capture the linkages between distant populations. A key question arising from these networks is the determination of the highly connected modules (communities). Communities can be defined using genetic similarity or distance statistics between multiple samples but due to limited specimen sampling capacity the boundaries of the communities for the known coral reefs in the seascape remain unresolved. In this study we use the microsatellite composition of individual corals to compare sample populations using a genetic dissimilarity measure (F_ST) which is then used to create a complex network. This network involved sampling 1025 colonies from 22 collection sites and examining 10 microsatellites loci. The links between each sampling site were given a strength that was created from the pair wise F_ST values. The result is an undirected weighted network describing the genetic dissimilarity between each sampled population. From this network we then determined the community structure using a leading eigenvector algorithm within graph theory. However, given the relatively limited sampling conducted, the representation of the regional genetic structure was incomplete. To assist with defining the boundaries of the genetically based communities we also integrated the communities derived from a hydrodynamic and distance based networks. The hydrodynamic network, though more comprehensive, was of smaller spatial extent than our genetic sampling. A Bayesian Belief network was developed to integrate the overlapping communities. The results indicate the genetic population structure of the Great Barrier Reef and provide guidance on where future genetic sampling should take place to complete the genetic diversity mapping.     

Identifying landscape scale patterns from individual scale processes

Extrapolating across scales is a critical problem in ecology. Explicit mechanistic models of ecological systems provide a bridge from measurements of processes at small and short scales to larger scales; spatial patterns at large scales can be used to test the outcomes of these models. However, it is necessary to identify patterns that are not dependent on initial conditions, because small scale initial conditions will not normally be measured at large scales. We examined one possible pattern that could meet these conditions, the relationship between mean and variance in abundance of a parasitic tick in an individual based model of a lizard tick interaction. We scaled discrepancies between the observed and simulated patterns with a transformation of the variance–covariance matrix of the observed pattern to objectively identify patterns that are “close”.