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.     

Biogeographical differences between the two Capra pyrenaica subspecies, C. p. victoriae and C. p. hispanica, inhabiting the Iberian Peninsula: Implications for conservation

The Iberian wild goat (Capra pyrenaica), an endemic species to the Iberian Peninsula, was represented in the past by four subspecies. Currently, the conservation status of the recognized subspecies - C. p. hispanica and C. p. victoriae - is quite different, being the former broadly distributed while the latter inhabits fewer and smaller areas. In this global context, our specific aims were to assess the biogeographical relationships, similarities and differences between the Iberian wild goat subspecies, and to evaluate their conservation status from a biogeographical perspective. We modelled the ecogeographical favourability for each subspecies, and also for the whole species, and we comparatively analyzed the spatial relationships among models in order to detect differential biogeographical patterns. Our results showed a differential pattern for each Iberian wild goat subspecies, and thus each one should be independently considered when designing management strategies. The distribution of C. p. victoriae was closely related to climatic conditions, thus it is highly sensitive to variations in climatic characteristics. Finally, from the perspective of metapopulation theory, and according to the geographical structure of their favourable areas, C. p. victoriae is at higher risk than C. p. hispanica. We discuss the implications of our results for species/subspecies conservation, and highlight the importance of conservation biogeography.     

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.     

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.     

Effects of temperature variation on TSD in turtle (C. picta) populations

We formulate a two-sex model of temperature-dependent sex determination (TSD) for a freshwater turtle (C. picta) population. The aim is to understand how environmental temperature variations and nest heat conduction properties affect the long term dynamics of the population. This is a key to understanding how global temperature changes may affect their survival. With stochastic inputs of ambient temperature and solar radiation, the model uses the heat equation to determine the temperature in the egg layer in the nest; in turn, this determines the sex ratio in the egg clutch using a variable degree-day model. Finally, a nonlinear Leslie type, stage-based, two-sex model, is used to determine the long term male and female populations. A two-sex model is required because of different development rates for males and females. The model is flexible enough to enable other researchers to examine the effects of temperature variation variations on other species with TSD, e.g., crocodilians, reptilians, as well as other turtle species. It can be adapted to study effects of nest location, soil type, rain events, different incubation periods, and density effects, for example, the dependence of the mating function on the ratio of males to females and each’s contribution to the sex of hatchlings. Modifications can be easily made to fit a specific life history traits. The model is a beginning step in understanding the long term, high fitness shown by many reptile species with TSD, and it may suggest to experimentalists what data may be relevant to these issues; it can also be useful to wildlife managers in developing strategies for intervention if needed. Among the principal findings are that temperature variability and detailed nest heat conduction properties may buffer projected negative effects on a population.     

Ecological niche similarities of Peponapis bees and non-domesticated Cucurbita species

Peponapis bees are considered specialized pollinators of Cucurbita flowers, a genus that presents several species of economic value (squashes and pumpkins). Both genera originated in the Americas, and their diversity dispersion center is in Mexico. Ten species of Peponapis and ten species of Cucurbita (only non-domesticated species) were analyzed considering the similarity of their ecological niche characteristics with respect to climatic conditions of their occurrence areas (abiotic variables) and interactions between species (biotic variables). The similarity of climatic conditions (temperature and precipitation) was estimated through cluster analyses. The areas of potential occurrence of the most similar species were obtained through ecological niche modeling and summed with geographic information system tools. Three main clusters were obtained: one with species that shared potential occurrence areas mainly in deserts (P. pruinosa, P. timberlakei, C. digitata, C. palmata, C. foetidissima), another in moist forests (P. limitaris, P. atrata, C. lundelliana, C. o. martinezii) and a third mainly in dry forests (C. a. sororia, C. radicans, C. pedatifolia, P. azteca, P. smithi, P. crassidentata, P. utahensis). Some species with similar ecological niche presented potential shared areas that are also similar to their geographical distribution, like those occurring predominantly on deserts. However, some clustered species presented larger geographical areas, such as P. pruinosa and C. foetidissima suggesting other drivers than climatic conditions to shape their distributions. The domestication of Cucurbita and also the natural history of both genera were considered also as important factors.     

A stage-structured, Aedes albopictus population model

Aedes albopictus has been the fastest spreading invasive animal species in the world from the mid-1980s until the mid-2000s. In areas it infests, it disrupts native mosquito ecology and can potentially vector up to 21 viruses. To better understand the population dynamics of this species, we created a temperature dependent population model. A stage-structured model was chosen to allow each life-stage to have different temperature dependent mortality and development rates, and each stage was modeled with an ordinary differential equation. Model parameters and distributions were based upon literature values. Initially, a basic model was constructed. This model then had parameters that were forced based upon daily average temperatures. Several criteria were used to evaluate the model, including a comparison to field data from Lubbock, TX. In a stochastic version of the model, a 95% confidence limit contained 70.7% of the field data points. Based upon these results, we feel reasonably confident that we have captured the role of temperature in driving the population dynamics of Ae. albopictus.     

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.     

Introduction to the special issue of Ecological Modelling: “Advances in Modeling Estuarine and Coastal Ecosystems: Approaches, Validation, and Applications”

The development and application of ecosystem models in estuarine and coastal systems has grown exponentially over the past four decades. Models have become ensconced as major tools for both heuristic study of ecosystem structure and function as well as for informing management decisions, particularly with respect to cultural eutrophication. In recent years an ever-expanding toolbox of modeling approaches is being offered to complement traditional methods. This expansion of modeling in estuarine and coastal science was exemplified by four sessions devoted to modeling at the 2007 biennial conference of the Estuarine Research Federation in Providence, RI. We felt the time was right to propose a special session of Ecological Modelling to synthesize talks from these sessions to present the state of the art in coastal and estuarine modeling. The collection of papers contained in this special issue presents a diversity of traditional and novel modeling approaches, methods for assessing model validity and predictability, and the utility of models in management applications. We believe that together these papers provide an excellent overview of current approaches to modeling estuarine hydrodynamics, water quality, and ecosystem/food web dynamics, applications of complex and relatively simple modeling approaches, applications in both deep and shallow coastal systems, goals relevant for both heuristic and management applications, and perspectives based on traditional mechanistic model development as well as more recent alternative approaches.     

A stochastic simulation model of brown shrimp (Penaeus aztecus ives) growth,movement, and survival in Galveston Bay,Texas

A stochastic simulation model of brown shrimp (Penaeus aztecus Ives) population dynamics in Galveston Bay, Texas, is described, validated, and used to evaluate the effects of management alternatives and changing environmental conditions on shrimp dynamics. The model is composed of submodels representing: (1) recruitment, (2) growth, (3) natural mortality, (4) fishing mortality, and (5) emigration of brown shrimp. The model predicts significant changes in total annual harvest from the food shrimp, bait, and recreational fisheries resulting from (1) closure of the bay system to all fishing except during the spring and fall open seasons, (2) two-week postponement of the opening and closing of the open seasons for the food shrimp fishery, (3) a 2.5°C increase and (4) a 2.5°C decrease in mean water temperature, (5) an 80% increase and (6) an 80% decrease in fishing effort. No significant change in the total annual harvest is predicted when the food shrimp fishing season is extended from May 15 through December 15. Sensitivity analysis suggests that field experimentation designed specifically to test the hypothesis of a 60-day time lag between brown shrimp recruitment into the bays and exposure to the fishery should receive high priority. Simulation results are discussed within a management framework.     

Identification of Illex coindetii, I. illecebrosus and I. argentinus (Cephalopoda: Ommastrephidae) throughout the Atlantic Ocean; by body and beak characters

Shortfinned squid species of the genus Illex support commercial fisheries throughout the Atlantic Ocean and Mediterranean Sea. Previous identification of interspecific and intraspecific populations by morphological and size-at-maturity studies have not provided conclusive results. We analysed morphometric body and beak variables (24 characters) in three species of the genus (I. coindetii, I. illecebrosus and I. argentinus), using a geographic and seasonal series of 33 populations for 1,500 specimens of I. coindetii, I. illecebrosus and I. argentinus. Residuals of the regression between each morphometric body and beak variable and mantle length were used as input in a stepwise discriminant analysis. Species discrimination by body and hectocotylus characters required at least eight variables and resulted in high correct-classification percentages for I. coindetii and I. argentinus (75% and 90%, respectively), whereas the best identification resulted from beak characters (83% correctly classified). Size of the suckerless basal arm, sucker-bearing length and beak lateral wall discriminated best among I. coindetii from northern Iberia, northwest Iberia (year-1996) and Ireland in the Atlantic and western Mediterranean versus middle and eastern Mediterranean samples. Canadian shelf and American samples were discriminated from Canadian slope I. illecebrosus. Winter/shelf and winter/slope samples of I. argentinus seemed to form a single biological group separated from Falkland Island, 46°S/autumn spawners and 46°S/1996 specimens along the Patagonian Shelf. No significant sexual or maturity polymorphism was obtained. Discriminant analysis optimised population diagnosis on a morphometric basis of interest in fisheries strategies. Electronic supplementary material to this paper can be obtained by using the Springer LINK server located at http://dx.doi.org/10.1007/s00227-002-0796-7.     

Predicting functional response and size selectivity of juvenile Notonecta maculata foraging on Daphnia magna

A mechanistic model was developed to assess the impact of predation of juvenile Notonecta maculata on size structured Daphnia magna populations and to provide a framework for quantifying the backswimmers uptake of food. Results of experiments and model predictions clearly demonstrate selective predation of backswimmers when fed with a choice of daphnid size classes, with patterns of selectivity differing across N. maculata instars. The model describes the foraging process empirically on the base of a general predation cycle including four conditional events instead of using classic functional response curves. For model parameterisation components of predation, namely probability of encounter, attack and success as well as time spent on handling prey was directly observed by means of video tracking experiments. Since attack rate, capture success and handling time appeared to be a function of prey size differing between Notonecta instars, the model takes into account ontogenic changes in both predator and prey characteristics. Independent data of functional response and size selectivity experiments were used for model validation and proved the model outcome to be consistent with observations.     

Agent-based simulation of effects of stress on forest tent caterpillar (Malacosoma disstria Hübner) population dynamics

The forest tent caterpillar (Malacosoma disstria Hübner) (FTC) has an outbreak cycle of approximately 10 years; however, smaller spatial scale analyses show some regions have longer or more frequent periods of high defoliation. This may be a result of local forest fragmentation, pollution or other sources of stress that may affect FTC directly or indirectly through stress on their hosts or parasitoids. Population dynamics of FTC were examined to investigate how stress may alter the severity and frequency of defoliation. We developed a spatially explicit agent-based model to simulate the host-parasitoid dynamics of FTC. Theoretical and empirically derived parameters were established using past literature and over 50 years of population data of FTC from Ontario, Canada. We find that increasing FTC fecundity, FTC dispersal or parasitoid mortality resulted in more severe outbreaks while a decrease in parasitoid fecundity or searching efficiency resulted in an overall elevation of defoliation. Parasitoid efficiency was the most effective parameter for altering the FTC defoliation. Since plant stress has been shown to alter several of these parameters in nature due to changes in food quality, habitat suitability, and chemical cue interference, our results suggest that forests affected by stressors such as climate change and pollution will have more severe and frequent defoliation from these insects than surrounding unaffected forests. As stressors such as drought and pollution emissions are predicted to increase in frequency or intensity over the next few decades, understanding how they may affect the outbreak cycle of a forest defoliator can aid in planning strategies to reduce the detrimental effects of this insect.     

Integration of distance, direction and habitat into a predictive migratory movement model for blue-winged teal (Anas discors)

Historically, the migration of birds has been poorly understood in comparison to other life stages during the annual cycle. The goal of our research is to present a novel approach to predict the migratory movement of birds. Using a blue-winged teal case study, our process incorporates not only constraints on habitat (temperature, precipitation, elevation, and depth to water table), but also approximates the likely bearing and distance traveled from a starting location. The method allows for movement predictions to be made from unsampled areas across large spatial scales. We used USGS’ Bird Banding Laboratory database as the source of banding and recovery locations. We used recovery locations from banding sites with multiple within-30-day recoveries were used to build core maximum entropy models. Because the core models encompass information regarding likely habitat, distance, and bearing, we used core models to project (or forecast) probability of movement from starting locations that lacked sufficient data for independent predictions. The final model for an unsampled area was based on an inverse-distance weighted averaged prediction from the three nearest core models. To illustrate this approach, three unsampled locations were selected to probabilistically predict where migratory blue-wing teals would stopover. These locations, despite having little or none data, are assumed to have populations. For the blue-winged teal case study, 104 suitable locations were identified to generate core models. These locations ranged from 20 to 228 within-30-day recoveries, and all core models had AUC scores greater than 0.80. We can infer based on model performance assessment, that our novel approach to predicting migratory movement is well-grounded and provides a reasonable approximation of migratory movement.     

Flow and transport characteristics at an Epinephelus guttatus (red hind grouper) spawning aggregation site in St. Thomas (US Virgin Islands)

This study characterizes the flow field at a spawning site located at the shelf break of a Caribbean island for the Epinephelus guttatus (red hind grouper) in relation to this species spawning events. In order to understand the oceanographic dynamics targeted by the fishes, current measurements were profiled throughout the water column for almost a year at the spawning site. The characteristics of the flow field and its evolution after spawning were investigated by using a numerical ocean model that resolved the observed tide and simulated the island scale flow where passive, neutrally buoyant virtual particles were released for 10 days to trace the flow pathways.Observed currents during the spawning period revealed that the flow was vertically sheared, to the south and weakest at the bottom, and to the west or east at the surface. The tidal analysis revealed that the flow at the time of spawning was directed across and on-shelf, although weaker close to the bottom. The model showed that the initial on-shelf transport was counteracted by the bottom flow directed to the shelf break, where virtual particles were entrained by the downwelling flow. A significant percent of particles resided less than two hundred meters deep, in the vicinity of the chlorophyll maximum and returned to the shelf break, close to the release location within 8-10 days. This journey was largely controlled by the timing between downwelling at the spawning site and upwelling further east at the shelf break, which was driven by the coupling between wind and tide induced vertical movements at the shelf break and deeper. The release location, vertical rotation of its flow field, and its transport properties were shown to be relatively resilient to the passage of transient sub-mesoscale eddies as well as to acute mesoscale flow reversals, suggesting that physical retention is maximized in the area surrounding the spawning site.     

Identifying key life-traits for the dynamics and gene flow in a weedy crop relative: Sensitivity analysis of the GeneSys simulation model for weed beet (Beta vulgaris ssp. vulgaris)

The benefits of genetically modified herbicide-tolerant (GMHT) sugar beet (Beta vulgaris) varieties stem from their presumed ability to improve weed control and reduce its cost, particularly targeting weed beet, a harmful annual weedy form of the genus Beta (i.e. B. vulgaris ssp. vulgaris) frequent in sugar beet fields. As weed beet is totally interfertile with sugar beet, it is thus likely to inherit the herbicide-tolerance transgene through pollen-mediated gene flow. Hence, the foreseeable advent of HT weed beet populations is a serious threat to the sustainability of GM sugar beet cropping systems. For studying and quantifying the long-term effects of cropping system components (crop succession and cultivation techniques) on weed beet population dynamics and gene flow, we developed a biophysical process-based model called GeneSys-Beet in a previous study. In the present paper, the model was employed to identify and rank the weed life-traits as function of their effect on weed beet densities and genotypes, using a global sensitivity analysis to model parameters. Monte Carlo simulations with simultaneous randomization of all life-trait parameters were carried out in three cropping systems contrasting for their risk for infestation by HT weed beets. Simulated weed plants and bolters (i.e. beet plants with flowering and seed-producing stems) were then analysed with regression models as a function of model parameters to rank processes and life-traits and quantify their effects. Key parameters were those determining the timing and success of growth, development, seed maturation and the physiological end of seed production. Timing parameters were usually more important than success parameters, showing for instance that optimal timing of weed management operations is more important than its exact efficacy. The ranking of life-traits though depended on the cropping system and, to a lesser extent, on the target variable (i.e. GM weeds vs. total weed population). For instance, post-emergence parameters were crucial in rotations with frequent sugar beet crops whereas pre-emergence parameters were most important when sugar beet was rare. In the rotations with frequent sugar beet and insufficient weed control, interactions between traits were small, indicating diverse populations with contrasted traits could prosper. Conversely, when sugar beet was rare and weed control optimal, traits had little impact individually, indicating that a small number of optimal combinations of traits would be successful. Based on the analysis of sugar beet parameters and genetic traits, advice for the future selection of sugar beet varieties was also given. In climatic conditions similar to those used here, the priority should be given to limiting the presence of hybrid seeds in seed lots rather than decreasing varietal sensitivity to vernalization.     

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.