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.     

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.     

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

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

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

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

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.     

Modelling habitat overlap among sympatric mesocarnivores in southern Illinois,USA

Few researchers have developed large-scale habitat models for sympatric carnivore species. We created habitat models for red foxes (Vulpes vulpes), coyotes (Canis latrans) and bobcats (Lynx rufus) in southern Illinois, USA, using the Penrose distance statistic, remotely sensed landscape data, and sighting location data within a GIS. Our objectives were to quantify and spatially model potential habitat differences among species. Habitat variables were quantified for 1-km² buffered areas around mesocarnivore sighting locations. Following variable reduction procedures, five habitat variables (percentage of grassland patches, interspersion–juxtaposition of forest patches, mean fractal dimension of wetland patches and the landscape, and road density) were used for analysis. Only one variable differed (P < 0.05) between red fox and coyote sighting areas (road density) and bobcat and coyote sighting areas (mean fractal dimension of the landscape). However, all five variables differed between red fox and bobcat sighting areas, indicating considerable differences in habitat affiliation between this pair-group. Compared to bobcats, red fox sightings were affiliated with more grassland cover and larger grassland patches, higher road densities, lower interspersion and juxtaposition of forest patches, and lower mean fractal dimension of wetland patches. These differences can be explained by different life history requirements relative to specific cover types. We then used the Penrose distance statistic to create habitat models for red foxes and bobcats, respectively, based on the five-variable dataset. An independent set of sighting locations were used to validate these models; model fit was good with 65% of mesocarnivore locations within the top 50% of Penrose distance values. In general, red foxes were affiliated with mixtures of agricultural and grassland cover, whereas bobcats were associated with a combination of grassland, wetland, and forest cover. The greatest habitat overlap between red foxes and bobcats was found at the interface between forested areas and more open cover types. Our study provides insight into habitat overlap among sympatric mesocarnivores, and the distance-based modelling approach we used has numerous applications for modelling wildlife–habitat relationships over large scales.     

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.     

Waterfowls habitat modeling: Simulation of nest site selection for the migratory Little Tern (Sterna albifrons) in the Nakdong estuary

This paper aims to find patterns in nest site selection by Little Terns Sterna albifrons, in the Nakdong estuary in South Korea. This estuary is important waterfowl stopover and breeding habitat, located in the middle of the East Asia-Australasian Flyway. The Little Tern is a common species easily observed near the seashore but their number is gradually declining around the world. We investigated their nests and eggs on a barrier islet in the Nakdong estuary during the breeding season (May to June, 2007), and a pattern for the nest site selection was identified using genetic programming (GP). The GP generated a predictive rule-set model for the number of Little Tern nests (training: R² = 0.48 and test: 0.46). The physical features of average elevation, variation of elevation, plant coverage, and average plant height were estimated to determine the influence on nest numbers for Little Tern. A series of sensitivity analyses stressed that mean elevation and vegetation played an important role in nest distribution for Little Tern. The influence of these two variables could be maximized when elevation changed moderately within the sampled quadrats. The study results are regarded as a good example of applying GP to vertebrate distribution patterning and prediction with several important advantages compared to conventional modeling techniques, and can help establish a management or restoration strategy for the species.     

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.     

Using classification trees to analyze the impact of exotic species on the ecological assessment of polder lakes in Flanders, Belgium

Polder lakes in Flanders are stagnant waters that were flooded by the sea in the past. Several of these systems are colonized by exotic species, but have hardly been studied until present. The aim of the present study was: (1) to assess the influence of exotic macrobenthic species on the outcome of the Multimetric Macroinvertebrate Index Flanders (MMIF) and (2) to use classification trees for evaluating to what extent physical-chemical characteristics affect the presence of exotic species.In total, 27 mollusc and 10 macro-crustacean species were present in the monitored lakes of which respectively five and four were exotic. The exclusion of the exotic species from the MMIF resulted in a significant decline of this ecological index (−0.03 ± 0.04; p = 0.00). This elimination often resulted into a lower ecological water quality class and more samples were classified into the bad and poor ecological water quality classes.Single-target classification trees for Gammarus tigrinus and Potamopyrgus antipodarum were constructed, relating environmental parameters and ecological status (MMIF) to the occurrence of both exotic invasive species. The major advantages of using single-target classification trees are the transparency of the rule sets and the possibility to use relatively small datasets. However, this classification technique only predicts a single-target attribute and the trees of the different species are often hard to integrate and use for water managers. As a solution, a multi-target approach was used in the present study. Exotic molluscs and crustaceans communities were modelled based on environmental parameters and the ecological status (MMIF) using multi-target classification trees. Multi-target classification trees can be used in management planning and investment decisions as they can lead to integrated decisions for the whole set of exotic species and avoid the construction of many models for each individual species. These trees provide general insights concerning the occurrence patterns of individual crustaceans and molluscs in an integrated way.     

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.     

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.     

Tropical amphibians in shifting thermal landscapes under land‐use and climate change

Land‐cover and climate change are both expected to alter species distributions and contribute to future biodiversity loss. However, the combined effects of land‐cover and climate change on assemblages, especially at the landscape scale, remain understudied. Lowland tropical amphibians may be particularly susceptible to changes in land cover and climate warming because many species have narrow thermal safety margins resulting from air and body temperatures that are close to their critical thermal maxima (CT_max). We examined how changing thermal landscapes may alter the area of thermally suitable habitat (TSH) for tropical amphibians. We measured microclimates in 6 land‐cover types and CT_max of 16 frog species in lowland northeastern Costa Rica. We used a biophysical model to estimate core body temperatures of frogs exposed to habitat‐specific microclimates while accounting for evaporative cooling and behavior. Thermally suitable habitat area was estimated as the portion of the landscape where species CT_max exceeded their habitat‐specific maximum body temperatures. We projected changes in TSH area 80 years into the future as a function of land‐cover change only, climate change only, and combinations of land‐cover and climate‐change scenarios representing low and moderate rates of change. Projected decreases in TSH area ranged from 16% under low emissions and reduced forest loss to 30% under moderate emissions and business‐as‐usual land‐cover change. Under a moderate emissions scenario (A1B), climate change alone contributed to 1.7‐ to 4.5‐fold greater losses in TSH area than land‐cover change only, suggesting that future decreases in TSH from climate change may outpace structural habitat loss. Forest‐restricted species had lower mean CT_max than species that occurred in altered habitats, indicating that thermal tolerances will likely shape assemblages in changing thermal landscapes. In the face of ongoing land‐cover and climate change, it will be critical to consider changing thermal landscapes in strategies to conserve ectotherm species.     

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.     

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

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

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.     

Predicting suitable habitat for the European lobster (Homarus gammarus), on the Basque continental shelf (Bay of Biscay), using Ecological-Niche Factor Analysis

Predicting species distribution and habitat suitability (HS) modelling, across broad spatial scales, is now a major challenge in marine ecology. The resulting knowledge is of considerable use in supporting the implementation of environmental legislation, integrated coastal zone management and ecosystem-based fisheries management. This contribution considers the identification of seafloor morphological characteristics, together with wave energy conditions, that determine the presence of European lobster (Homarus gammarus); and it predicts suitable habitats over the Basque continental shelf (Bay of Biscay), in summer. The results obtained, by applying Ecological-Niche Factor Analysis (ENFA), indicate that lobster habitat differs considerably from the mean environmental condition over the study area; likewise, that it is restrictive in terms of the range of conditions in which they dwell. The best of the environmental predictors found to be: distance to the rock substrate; Benthic Position Index; wave flux over the seafloor; and the underlying bathymetry. A habitat suitability map was produced, with a high model quality (Boyce index: 0.98 ± 0.06). The most suitable habitat for European lobster are locations at the boundary between sedimentary- and rocky-bottoms, coincident with seafloor depressions with a steep slope, with medium to high wave energy conditions, and located within a range of water depths of 35–40 m. This approach demonstrates the applicability of the method in case studies where only presence data are available, together with the inclusion of environmental variables obtained from different sources.     

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.     

Behaviour moderates climate warming vulnerability in high-rocky-shore snails: interactions of habitat use, energy consumption and environmental temperature

High-rocky-shore intertidal animals are predicted to be generally more vulnerable to climate warming than lower-shore species, because their thermal tolerances lie closer to maximum environmental temperatures (T _e). However, this prediction is based on taxonomically and ecologically limited information. The present study investigated the effect of habitat use on climate warming vulnerability of the tropical high-shore snail, Echinolittorina malaccana (from Brunei Darussalam, 5°N), which aestivates in sun-exposed or shaded habitats. The thermal regimes of these habitats differed vastly, but snails showed similar daily energy consumption in either habitat, due to temperature-insensitive metabolism (TIM) between 35 and 46 °C in the sun-resting snails. However, maximum T _e values in the shade and the sun were 35 and 46 °C, respectively, suggesting that sun-resting snails, which presently experience temperatures near the incipient lethal temperature range (46–56 °C), should be more threatened by further warming than shade-resting snails, which have an 11 °C ‘safety margin’. Thus, vulnerability of high-shore species to climate warming could be moderated by availability of shaded habitat, making predictions for these species more complex than previously realized.     

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.