A statistical test to detect vortices in the current fields of bodies of water

Vortices could play an important role in the occurrence of certain biological phenomena, such as the massive proliferation of harmful algae in bodies of water. Many measures exist to detect vortices in fluids, but little is known about the stochastic behavior of these quantities with data that contain statistical noise. Consequently they do not provide control over the probability of false positives and give little information about the risk of false negatives. Obtaining such control requires a statistical testing procedure. In this paper, we develop a test for vortices in random current fields using only the directions of the current observed at points on a regular grid. We construct a change-point test for spatially ordered angular data to detect the presence of a local vortex. A global vortex detection procedure based on this test is developed and applied to a data set from a lagoon located in the south of France. It is shown that this procedure can detect the presence of multiple vortices with good accuracy.     

Fourth-corner correlation is a score test statistic in a log-linear trait–environment model that is useful in permutation testing

Ecologists wish to understand the role of traits of species in determining where each species occurs in the environment. For this, they wish to detect associations between species traits and environmental variables from three data tables, species count data from sites with associated environmental data and species trait data from data bases. These three tables leave a missing part, the fourth-corner. The fourth-corner correlations between quantitative traits and environmental variables, heuristically proposed 20 years ago, fill this corner. Generalized linear (mixed) models have been proposed more recently as a model-based alternative. This paper shows that the squared fourth-corner correlation times the total count is precisely the score test statistic for testing the linear-by-linear interaction in a Poisson log-linear model that also contains species and sites as main effects. For multiple traits and environmental variables, the score test statistic is proportional to the total inertia of a doubly constrained correspondence analysis. When the count data are over-dispersed compared to the Poisson or when there are other deviations from the model such as unobserved traits or environmental variables that interact with the observed ones, the score test statistic does not have the usual chi-square distribution. For these types of deviations, row- and column-based permutation methods (and their sequential combination) are proposed to control the type I error without undue loss of power (unless no deviation is present), as illustrated in a small simulation study. The issues for valid statistical testing are illustrated using the well-known Dutch Dune Meadow data set.     

A set of associated statistical tests for spatial clustering

The purpose of this paper is to develop a set of associated statistical tests for spatial clustering. In particular, a set of three associated tests will be developed; these will correspond to the three types of tests set out by Besag and Newell (general tests, focused tests, and tests for the detection of clustering). The associated tests draw primarily, though not exclusively, upon existing tests and results. The principal contributions are based upon the score statistic for focused tests, which has been an important approach to testing for clustering around environmental hazards. The first contribution consists of the formulation of a global statistic for general tests that corresponds to focused score statistics, along with an assessment of the distribution of the statistic under the null hypothesis of no raised incidence. The local score statistics used for focused tests will have the property of summing to the global statistic used for the corresponding general test. Attention is also given to the maximum local score statistic for the “test for the detection of clustering”. The critical values of this statistic which are required for testing the null hypothesis are described. Application of the methods is made to leukemia data for central New York State.     

Developments in remote sensing, dynamic modelling and GIS applications for integrated coastal zone management

The International Institute for Aerospace Survey and Earth Sciences (ITC) has a research programme that should result in an integrated environmental coastal zone management system through three subprojects. The programme aims to develop methodologies and tools for assessing coastal zone changes, and for the evaluation of scenarios for coastal zone management, based on a spatio-temporal Geographical Information System (GIS) working platform which integrates remote sensing data, physical-morphodynamic and eco-hydrologic modelling, and a decision support system. The first subproject develops methodologies for the generation of optimum Remote Sensing (RS) data sets, leading to better interpretation and complementary use of conventional and new remote sensing imagery. It also integrates RS, GIS, and modelling through hypothesis generation, parameter estimation, evaluation and validation. The second subproject facilitates qualitative and quantitative analysis and prediction of the physical aspects of coastal landscape development under the influence of natural processes and human impacts. This subproject is based on the application of remote sensing and dynamic modelling. The third subproject leads to a spatio-temporal working platform which supports data integration of RS and in-situ measurements, and qualitative and quantitative analysis for the prediction of coastal landscape development. Both support decision making in Integrated Coastal Zone Management.     

遥感与GIS支持下的城市绿地信息提取方法研究--以深圳市为例

城市绿地作为城市结构中的自然生产力主体，在城市系统中起着重要作用。传统的城市绿地调查多采用人工普查并结合统计学方法进行城市绿地调查，需要投入大量的人力和资金，但得到的数据精度和现时性都较差，而采用高分辨率遥感卫星数据与地理信息系统相结合的方法，可以达到对城市绿地遥感数据的快速特征提取和准确的空间统计分析。文章采用基于遥感和GIS相结合的方法，并以深圳市为例，进行了城市绿地特征提取方法的研究和应用分析。     

Functional regression on remote sensing data in oceanography

The aim of this study is to propose the use of a functional data analysis approach as an alternative to the classical statistical methods most commonly used in oceanography and water quality management. In particular we consider the prediction of total suspended solids (TSS) based on remote sensing (RS) data. For this purpose several functional linear regression models and classical non-functional regression models are applied to 10 years of RS data obtained from medium resolution imaging spectrometer sensor to predict the TSS concentration in the coastal zone of the Guadalquivir estuary. The results of functional and classical approaches are compared in terms of their mean square prediction error values and the superiority of the functional models is established. A simulation study has been designed in order to support these findings and to determine the best prediction model for the TSS parameter in more general contexts.     

Evaluating the environmental flows of China's Wolonghu wetland and land use changes using a hydrological model, a water balance model, and remote sensing

Environmental flows are critical to sustaining a variety of plant and animal communities in wetlands. However, evaluation of environmental flows is hampered by the problem of hydrological and ecological data shortage, especially in many developing countries such as China. Based on a hydrological model, a water balance model and remote sensing data, we assessed the environmental flows of China's Wolonghu wetland with limited data. The hydrological model provides input data for the water balance model of the wetland, and the remote sensing data can be used to assess land use changes. Integration of these two models with the remote sensing data revealed both the environmental flows of the Wolonghu wetland and the relationships between these environmental flows and land use changes. The results demonstrate that environmental flows have direct and indirect influences on the wetland ecosystem and should be linked to sustainable wetland management.     

Statistical inference using the g or K point pattern spatial statistics

Spatial point pattern analysis provides a statistical method to compare an observed spatial pattern against a hypothesized spatial process model. The G statistic, which considers the distribution of nearest neighbor distances, and the K statistic, which evaluates the distribution of all neighbor distances, are commonly used in such analyses. One method of employing these statistics involves building a simulation envelope from the result of many simulated patterns of the hypothesized model. Specifically, a simulation envelope is created by calculating, at every distance, the minimum and maximum results computed across the simulated patterns. A statistical test is performed by evaluating where the results from an observed pattern fall with respect to the simulation envelope. However, this method, which differs from P. Diggle's suggested approach, is invalid for inference because it violates the assumptions of Monte Carlo methods and results in incorrect type I error rate performance. Similarly, using the simulation envelope to estimate the range of distances over which an observed pattern deviates from the hypothesized model is also suspect. The technical details of why the simulation envelope provides incorrect type I error rate performance are described. A valid test is then proposed, and details about how the number of simulated patterns impacts the statistical significance are explained. Finally, an example of using the proposed test within an exploratory data analysis framework is provided.     

A weight-adjusted Peto’s test when cause of death is not assigned

A new statistical testing approach using a weighted logrank statistic is developed for rodent tumorigenicity assays that have a single terminal sacrifice but not cause-of-death data. Instead of using cause-of-death assignment by pathologists, the number of fatal tumors is estimated by a constrained nonparametric maximum likelihood estimation method. For data lacking cause-of-death information, the Peto test is modified with estimated numbers of fatal tumors and a Fleming–Harrington-type weight, which is based on an estimated tumor survival function. A bootstrap resampling method is used to estimate the weight function. The proposed testing method with the weight adjustment appears to improve the performance in various situations of single-sacrifice animal experiments. A Monte Carlo simulation study for the proposed test is conducted to assess size and power of the test. This testing approach is illustrated using a real data set.     

浅水流动有限体积法/Osher格式的二维水流-水质模拟

应用无结构网格有限体积法 /Osher格式的二维水流 -水质耦合数值模型 ,采用黎曼近似解计算模拟区域水量、动量及污染物输送通量。通过对长江南京八卦洲江段浓度场的模拟计算 ,其结果与遥感卫星图片确定的污染带分布相吻合 ,证实了模型的合理性和模拟能力     

基于MODIS数据的河南省冬小麦产量遥感估算模型

小麦是世界上最重要的粮食作物,小麦生产对中国的粮食保障起着十分重要的作用,及时、准确、大范围对小麦产量进行监测预报,对于农学经济发展和粮食政策制定具有极为重要的现实意义。对作物产量进行遥感监测的原理是建立在其遥感特征基础之上的,通过建立作物长势指标与遥感信息的定量关系,可实现对作物产量的监测预报。文章基于2009年MODIS遥感数据和气象数据,利用Arcgis和ENVI提取纯小麦像元,并提取纯小麦像元对应的NDVI、NPP和LAI,获取分县NDVI、NPP和LAI均值,利用统计软件对产量数据和分县遥感参数均值进行数据整理和分析,建立了河南省冬小麦产量估算模型。以往研究多采用遥感图像上某像元和地面调查点进行研究,具有很大的不确定性,文章以县为单位,对冬小麦平均单产和县域内冬小麦种植像元遥感参数的均值进行相关研究,提高了模型模拟精度。同时文章选用多种遥感参数和多项气象因子建立估产模型,避免了针对一个参数进行估产的局限性。在最佳时相的选择上,根据冯美辰（2010）以往的研究结果,从4月以后,5月8日和4月20Et植被指数和产量相关性最大,4月份之前冬小麦处于返青到拔节期,对产量来说还有很多不确定闪素,因此文章选用5月8El和4月20日进行冬小麦估产研究。结果表明,5月8日的估产模型优于4月20日,加入气象冈子的遥感气象估产模型优于只采用遥感参数进行估产的遥感模型。利用2010年产量数据对模型精度进行检验,遥感气象模型预测精度在70．2％N99．7％之间,平均精度为90．7％;遥感模型预测精度在68．1％到95．5％之间,平均精度为83．9％。表明遥感气象模型模拟精度更高,其精度可以满足大面积估产要求,可以对产量预报提供科学参考。     

Upper level set scan statistic for detecting arbitrarily shaped hotspots

A declared need is around for geoinformatic surveillance statistical science and software infrastructure for spatial and spatiotemporal hotspot detection. Hotspot means something unusual, anomaly, aberration, outbreak, elevated cluster, critical resource area, etc. The declared need may be for monitoring, etiology, management, or early warning. The responsible factors may be natural, accidental, or intentional. This proof-of-concept paper suggests methods and tools for hotspot detection across geographic regions and across networks. The investigation proposes development of statistical methods and tools that have immediate potential for use in critical societal areas, such as public health and disease surveillance, ecosystem health, water resources and water services, transportation networks, persistent poverty typologies and trajectories, environmental justice, biosurveillance and biosecurity, among others. We introduce, for multidisciplinary use, an innovation of the health-area-popular circle-based spatial and spatiotemporal scan statistic. Our innovation employs the notion of an upper level set, and is accordingly called the upper level set scan statistic, pointing to a sophisticated analytical and computational system as the next generation of the present day popular SaTScan. Success of surveillance rests on potential elevated cluster detection capability. But the clusters can be of any shape, and cannot be captured only by circles. This is likely to give more of false alarms and more of false sense of security. What we need is capability to detect arbitrarily shaped clusters. The proposed upper level set scan statistic innovation is expected to fill this need     

Distribution-free ranked-set sample procedures for umbrella alternatives in the m-sample setting

In settings where measurements are costly and/or difficult to obtain but ranking of the potential sample data is relatively easy and reliable, the use of statistical methods based on a ranked-set sampling approach can lead to substantial improvement over analogous methods associated with simple random samples. Previous nonparametric work in this area has been concentrated almost exclusively on the one- and two-sample location problems. In this paper we develop ranked-set sample procedures for the m-sample location setting where the treatment effect parameters follow a restricted umbrella pattern. Distribution-free testing procedures are developed for both the case where the peak of the umbrella is known and for the case where it is unknown. Small sample and asymptotic null distribution properties are provided for the peak-known test statistic.     

中国农田作物植被碳储量研究进展

作物植被碳储量是全球陆地生态系统碳库的重要组成部分。中国农田作物植被碳储量的估算主要采用参数估算法、遥感资料反演法和环境参数模型法。通过对中国近几十年来全国和区域尺度作物植被碳储量的估算研究,获得了一些作物的经济系数、含碳率和作物收获部分水分系数等估算参数值,探讨了遥感反演和环境参数模型方法,并提出加强农田基本建设、改进农业生产技术与管理、调整作物结构和加强作物秸秆利用等固碳措施。目前对中国农田作物植被碳储量的估算仍存在较大的不确定性,获取的估算参数尚不充分,估算方法和模型有待完善,对作物植被碳储量变化的源/汇效应尚未取得统一认识。虽然在农田生态系统中土壤碳储量（密度）普遍大于作物植被碳储量（密度）,但作物植被碳储量仍然是一个数量可观、并有增加潜力和可能的碳库,其大小及秸秆利用情况直接影响着土壤碳库。因此,对农田作物植被碳储量应分时段和区域具体分析,才能认识其源/汇效应。今后应在以下几方面进一步加强作物植被碳储量的研究：进一步完善和改进估算方法;加强作物植被碳储量估算及固碳措施的区域个例研究,探索不同空间尺度作物植被碳储量的尺度转换;开展作物碳储量动态及固碳机理的综合研究。此外,还应就气候变化与作物植被碳储量的相互耦合关系进行探讨。     

A permutation-based combination of sign tests for assessing habitat selection

The analysis of habitat selection in radio-tagged animals is approached by comparing the portions of use against the portions of availability observed for each habitat type. Since data are linearly dependent with singular variance-covariance matrices, standard multivariate statistical tests cannot be applied. To bypass the problem, compositional data analysis is customarily performed via log-ratio transform of sample observations. The procedure is criticized in this paper, emphasizing the several drawbacks which may arise from the use of compositional analysis. An alternative nonparametric solution is proposed in the framework of multiple testing. The habitat use is assessed separately for each habitat type by means of the sign test performed on the original observations. The resulting p values are combined in an overall test statistic whose significance is determined permuting sample observations. The theoretical findings of the paper are checked by simulation studies. Applications to case studies previously considered in literature are discussed.     

Improved testing of species traits--environment relationships in the fourth-corner problem

The fourth-corner problem entails estimation and statistical testing of the relationship between species traits and environmental variables from the analysis of three data tables. In a 2008 paper, S. Dray and P. Legendre proposed and evaluated five permutation methods for statistical significance testing, including a new two-step testing procedure. However, none of these attained the correct type I error in all cases of interest. We solve this problem by showing that a small modification of their two-step procedure controls the type I error in all cases. The modification consists of adjusting the significance level from mean square root of alpha to alpha or, equivalently, of reporting the maximum of the individual P. values as the final one. The test is also applicable to the three-table ordination method RLQ.     

A GIS-based approach to spatio-temporal analysis of environmental pollution in urban areas: A case study of Prague's environment extended by LIDAR data

Environmental pollution of urban areas is one of key factors that state authorities and local agencies have to consider in the decision-making process. To find a compromise among many criteria, spatial analysis extended by geostatistical methods and dynamic models has to be carried out. In this case, spatial analysis includes processing of a wide range of air, water and soil pollution data and possibly noise assessment and waste management data. Other spatial inputs consist of data from remote sensing and GPS field measurements. Integration and spatial data management are carried out within the framework of a geographic information system (GIS). From a modeling point of view, GIS is used mainly for the preprocessing and postprocessing of data to be displayed in digital map layers and visualized in 3D scenes. Moreover, for preprocessing and postprocessing, deterministic and geostatistical methods (IDW, ordinary kriging) are used for spatial interpolation; geoprocessing and raster algebra are used in multi-criteria evaluation and risk assessment methods. GIS is also used as a platform for spatio-temporal analyses or for building relationships between the GIS database and stand-alone modeling tools. A case study is presented illustrating the application of spatial analysis to the urban areas of Prague. This involved incorporating environmental data from monitoring networks and field measurements into digital map layers. Extra data inputs were used to represent the 3D concentration fields of air pollutants (ozone, NO₂) measured by differential absorption LIDAR. ArcGIS was used to provide spatial data management and analysis, extended by modeling tools developed internally in the ArcObjects environment and external modules developed with MapObjects. Ordinary kriging methods were employed to predict ozone concentrations in selected 3D locations together with estimates of variability. Higher ozone concentrations were found above crossroads with their heavy traffic than above the surrounding areas. Ozone concentrations also varied with height above the digital elevation model. Processed data, spatial analysis and models are integrated within the framework of the GIS project, providing an approach that state and local authorities can use to address environmental protection issues.     

Experimental and environmental factors affect spurious detection of ecological thresholds

Threshold detection methods are increasingly popular for assessing nonlinear responses to environmental change, but their statistical performance remains poorly understood. We simulated linear change in stream benthic macroinvertebrate communities and evaluated the performance of commonly used threshold detection methods based on model fitting (piecewise quantile regression [PQR]), data partitioning (nonparametric change point analysis [NCPA]), and a hybrid approach (significant zero crossings [SiZer]). We demonstrated that false detection of ecological thresholds (type I errors) and inferences on threshold locations are influenced by sample size, rate of linear change, and frequency of observations across the environmental gradient (i.e., sample-environment distribution, SED). However, the relative importance of these factors varied among statistical methods and between inference types. False detection rates were influenced primarily by user-selected parameters for PQR (tau) and SiZer (bandwidth) and secondarily by sample size (for PQR) and SED (for SiZer). In contrast, the location of reported thresholds was influenced primarily by SED. Bootstrapped confidence intervals for NCPA threshold locations revealed strong correspondence to SED. We conclude that the choice of statistical methods for threshold detection should be matched to experimental and environmental constraints to minimize false detection rates and avoid spurious inferences regarding threshold location.     

Effects of experimental protocol on global vegetation model accuracy: A comparison of simulated and observed vegetation patterns for Asia

Prognostic vegetation models have been widely used to study the interactions between environmental change and biological systems. This study examines the sensitivity of vegetation model simulations to: (i) the selection of input climatologies representing different time periods and their associated atmospheric CO₂ concentrations, (ii) the choice of observed vegetation data for evaluating the model results, and (iii) the methods used to compare simulated and observed vegetation. We use vegetation simulated for Asia by the equilibrium vegetation model BIOME4 as a typical example of vegetation model output. BIOME4 was run using 19 different climatologies and their associated atmospheric CO₂ concentrations. The Kappa statistic, Fuzzy Kappa statistic and a newly developed map-comparison method, the Nomad index, were used to quantify the agreement between the biomes simulated under each scenario and the observed vegetation from three different global land- and tree-cover data sets: the global Potential Natural Vegetation data set (PNV), the Global Land Cover Characteristics data set (GLCC), and the Global Land Cover Facility data set (GLCF). The results indicate that the 30-year mean climatology (and its associated atmospheric CO₂ concentration) for the time period immediately preceding the collection date of the observed vegetation data produce the most accurate vegetation simulations when compared with all three observed vegetation data sets. The study also indicates that the BIOME4-simulated vegetation for Asia more closely matches the PNV data than the other two observed vegetation data sets. Given the same observed data, the accuracy assessments of the BIOME4 simulations made using the Kappa, Fuzzy Kappa and Nomad index map-comparison methods agree well when the compared vegetation types consist of a large number of spatially continuous grid cells. The results of this analysis can assist model users in designing experimental protocols for simulating vegetation.     

Use of ambiguous detections to improve estimates from species distribution models

As large carnivores recover throughout Europe, their distribution needs to be studied to determine their conservation status and assess the potential for human-carnivore conflicts. However, efficient monitoring of many large carnivore species is challenging due to their rarity, elusive behavior, and large home ranges. Their monitoring can include opportunistic sightings from citizens in addition to designed surveys. Two types of detection errors may occur in such monitoring schemes: false negatives and false positives. False-negative detections can be accounted for in species distribution models (SDMs) that deal with imperfect detection. False-positive detections, due to species misidentification, have rarely been accounted for in SDMs. Generally, researchers use ad hoc data-filtering methods to discard ambiguous observations prior to analysis. These practices may discard valuable ecological information on the distribution of a species. We investigated the costs and benefits of including data types that may include false positives rather than discarding them for SDMs of large carnivores. We used a dynamic occupancy model that simultaneously accounts for false negatives and positives to jointly analyze data that included both unambiguous detections and ambiguous detections. We used simulations to compare the performances of our model with a model fitted on unambiguous data only. We tested the 2 models in 4 scenarios in which parameters that control false-positive detections and true detections varied. We applied our model to data from the monitoring of the Eurasian lynx (Lynx lynx) in the European Alps. The addition of ambiguous detections increased the precision of parameter estimates. For the Eurasian lynx, incorporating ambiguous detections produced more precise estimates of the ecological parameters and revealed additional occupied sites in areas where the species is likely expanding. Overall, we found that ambiguous data should be considered when studying the distribution of large carnivores through the use of dynamic occupancy models that account for misidentification.