Reducing over-reporting of deterministic co-occurrence patterns in biotic communities

Null models of species co-occurrence are widely used to infer the existence of various ecological processes. Here we investigate the susceptibility of the most commonly used of these models (the C-score in conjunction with the sequential swap algorithm) to type 1 and type 2 errors. To do this we use simulated datasets with a range of numbers of sites, species and coefficients of variation (CV) in species abundance. We find that this model is particularly susceptible to type 1 errors when applied to large matrices and those with low CV in species abundance. As expected, type 2 error rates decrease with increasing numbers of sites and species, although they increase with increasing CV in species abundance. Despite this, power remains acceptable over a wide range of parameter combinations. The susceptibility of this analytical method to type 1 errors indicates that many previous studies may have incorrectly reported the existence of deterministic patterns of species co-occurrence. We demonstrate that in order to overcome the problem of high type 1 error rates, the number of swaps used to generate null distributions for smaller matrices needs to be increased to over 50,000 swaps (well beyond the 5000 commonly used in published analyses and the 30,000 suggested by Lehsten and Harmand, 2006). We also show that this approach reduces type 1 error rates in real datasets. However, even using this solution, larger datasets still suffer from high type 1 error rates. Such datasets therefore require the use of very large numbers of swaps, which calls for improvements in the most commonly used software. In general, users of this powerful analytical method must be aware that they need surprisingly large numbers of swaps to obtain unbiased estimates of structuring in biotic communities.     

A Bayesian approach to evaluating site impairment

In the United States, each state is required to list water resources that are declared to be impaired under guidelines set by the Clean Water Act. Measurements are typically collected on a number of chemical constituents and compared with a standard. If there are too many measurements exceeding the standard, then the site is declared impaired. The approach is non-statistical but similar to a Binomial test. The Binomial approach would convert the measurements to binary data then test if the proportion exceeding the standard is excessive. Both methods convert measurements to binary values hence exclude potentially important information in the data. We present a statistical approach using a Bayesian model that uses the raw data instead of the binary transformed data. The population distribution of a family of location-scale parameter models is studied under the model. Posterior distributions from the Bayesian analysis are used in the decision-making process and error probabilities for the Bayesian and the Binomial approaches are compared for a normal population.     

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.     

Spatial error propagation when computing linear combinations of spectral bands: The case of vegetation indices

The collection of accurate and timely information on land use, crops, forest and vegetation are increasingly based on remote sensing spectral measurements produced by satellites. The most recent spacecrafts like the Earth Observing 1 (EO-1) produce a rich source of information being endowed with hyperspectral sensors that can provide up to 200 or more channels. In many instances such a multivariate signal has to be reduced to just one single value per pixel representing a particular characteristic of land. Linear combinations of bands are the general form of many indices. Since each individual image used to construct indices contains errors, when combined they produce a propagation of errors, a process that can distort a final output map. In this paper we measure the extent of error propagation when building linear vegetation indices. We consider three types of indices: the difference vegetation index (DVI), selected Kauth-Thomas indices (SBI, GVI and WET), and principal components, using benchmarking examples taken from the remote sensing literature. The main implication emerging from these examples is that the SBI and the first principal component are the indices more prone to error propagation. The formalization presented here allows a user to derive measures of error propagation in cases where technical characteristics of a sensor and physical characteristics of a landscape are known. These results can help a user to choose between alternative vegetation indices, and to associate a measure of reliability with such indices.     

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.     

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.     

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.     

Simulating the temporal pattern of waste production in farmed gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and Atlantic bluefin tuna (Thunnus thynnus)

Most fish farming waste output models provide gross waste rates as a function of stocked or produced biomass for a year or total culture cycle, but without contemplating the temporality of the discharges. This work aims to ascertain the temporal pattern of waste loads by coupling available growth and waste production models and developing simulation under real production rearing conditions, considering the overlapping of batches and management of stocks for three widely cultured species in the Mediterranean Sea: gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and Atlantic bluefin tuna (Thunnus thynnus). For a similar annual biomass production, the simulations showed that waste output and temporal dumping patterns differ between the three species as a result of the disparities in growth velocity, nutrient digestibility, maintenance metabolic budget and husbandry. The simulations allowed the temporal patterns including the periods of maximum discharge and the dissolved and particulate nitrogen and phosphorus content in the wastes released to be determined, both of which were seen to be species-specific.     

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

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.     

Measuring complexity in three-trophic level systems

This essay is a trial on measuring complexity in a three-trophic level system by using a convex function of the informational entropy. The complexity measure defined here is compatible with the fact that real complexity lies between ordered and disordered states. Applying this measure to the data collected for two three-trophic level systems some hints about their organization are obtained.     

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.     

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.     

Monitoring, modeling, and management impacts of bivalve filter feeders in the oligohaline and tidal fresh regions of the Chesapeake Bay system

Populations of bivalve filter feeders are distributed throughout the oligohaline waters of the Chesapeake Bay system and, to a lesser extent, in tidal fresh waters as well. Previous studies indicate these bivalves significantly diminish phytoplankton concentrations in one major tributary, the Potomac River, and observed chlorophyll concentrations suggest bivalve influence on phytoplankton in other oligohaline reaches. We incorporated a model of these bivalves into an existing eutrophication model of the system. The model indicated that bivalves may reduce phytoplankton concentrations in oligohaline and tidal fresh waters throughout the system but the most significant effects were noted in the Potomac and Patuxent tributaries. Bivalve impacts were related to hydraulic residence time. The greatest phytoplankton reductions occurred in the regions with the longest residence time. Model carbon and nutrient budgets indicated bivalves removed 14% to 40% of the carbon load, 11% to 23% of the nitrogen load, and 37% to 84% of the phosphorus load to the regions where their impact on computed chlorophyll was greatest.     

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

Density estimation using ranked-set sampling data

The ranked-set sampling (RSS) is applicable in practical problems where the variable of interest for an observed item is costly or time-consuming but the ranking of a set of items according to the variable can be easily done without actual measurement. In the context of RSS, the need for density estimation arises in certain statistical procedures. The density estimation also has its own interest. In this article, we develop a method for the density estimation using RSS data. We derive the properties of the resulted density estimate and compare it with its counterpart in simple random sampling (SRS). It is shown that the density estimate using RSS data provides a better estimate of the density than the usual density estimate using SRS data. The density estimate developed in this article can well serve various purposes in the context of RSS.     

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.     

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.     

学校图书信息管理系统设计

图书管理系统是典型的信息管理系统（MIS），本文主要介绍后台数据库的建立、维护及前端应用程序开发两个方面．该系统采用BIS架构，将图书入库、查询、借阅、维护、用户信息等功能集成在一起，具有用户介面友好、操作简单、数据安全、功能强大、实用性强的特点．     

大田作物对UV-B辐射的响应及其研究方法现状与展望

目前紫外线对作物影响的研究主要从两个方面进行,即人为模拟增强UV-B辐射和去除UV-B.这两个方面采用了不同的实验方法,但都会在很多环节上产生系统误差,本文概述了误差产生的原因及多种控制措施,并讨论了作物对两种处理的响应差异.相对来说,大田植物对去除UV-B实验响应更大,其生长及生物量的变化更明显;而增强的UV-B辐射对农田或自然生态系统的生产力影响较小,但其在调节农田生态系统过程上有更重要的作用,比如种间竞争、物候期的改变以及由于次生代谢产物的变化而导致营养循环和动物采食发生变化等.为了能更加深入了解UV-B的作用机制,研究者应采用增强和减弱并用的方法进行UV-B辐射研究.参43