Modeling Diapause Termination of Rhipicephalus Appendiculatus using Statistical Tools to Detect Sudden Behavioral Changes and Time Dependencies

This paper presents statistical methodology to analyze longitudinal binary responses for which a sudden change in the response occurs in time. Probability plots, transition matrices, and change-point models and more advanced techniques such as generalized auto-regression models and hidden Markov chains are presented and applied on a study on the activity of Rhipicephalus appendiculatus, the major vector of Theileria parva, a fatal disease in cattle. This study presents individual measurements on female R. appendiculatus, which are terminating their diapause (resting status) and become active. Comprehending activity patterns is very important to better understand the ecology of R. appendiculatus. The model indicates that activity and non-activity act in an absorbing way meaning that once a tick becomes active it shows a tendency to remain active. The change-point model estimates that the sudden change in activity happens on December 10. The reaction of ticks on acceleration and changes in rainfall and temperature indicates that ticks can sense climatic changes. The study revealed the underlying not visually observable states during diapause development of the adult tick of R. appendiculatus. These states could be related to phases during the dynamic event of diapause development and post-diapause activity in R. appendiculatus.     

活性污泥过程数学模型进展

介绍了活性污泥过程模型ASM1，ASM2及ASM3的建模方法，并分析了各模型的微生物学基础。     

On the latent state estimation of nonlinear population dynamics using Bayesian and non-Bayesian state-space models

Nonlinear state-space models have been increasingly applied to study population dynamics and data assimilation in environmental sciences. State-space models can account for process error and measurement error simultaneously to correct for the bias in the estimates of system state and model parameters. However, few studies have compared the performance of different nonlinear state-space models for reconstructing the state of population dynamics from noisy time series. This study compared the performance of the extended Kalman filter (EKF), unscented Kalman filter (UKF) and Bayesian nonlinear state-space models (BNSSM) through simulations. Synthetic population time series were generated using the theta logistic model with known parameters, and normally distributed process and measurement errors were introduced using the Monte Carlo simulations. At higher levels of nonlinearity, the UKF and BNSSM had lower root mean square error (RMSE) than the EKF. The BNSSM performed reliably across all levels of nonlinearity, whereas increased levels of nonlinearity resulted in higher RMSE of the EKF. The Metropolis–Hastings algorithm within the Gibbs algorithm was used to fit the theta logistic model to synthetic time series to estimate model parameters. The estimated posterior distribution of the parameter θ indicated that the 95% credible intervals included the true values of θ (=0.5 and 1.5), but did not include 1.0 and 0.0. Future studies need to incorporate the adaptive Metropolis algorithm to estimate unknown model parameters for broad applications of Bayesian nonlinear state-space models in ecological studies.     

Examination of model predictions at different horizontal grid resolutions

While fluctuations in meteorological and air quality variables occur on a continuum of spatial scales, the horizontal grid spacing of coupled meteorological and photochemical models sets a lower limit on the spatial scales that they can resolve. However, both computational costs and data requirements increase significantly with increasing grid resolution. Therefore, it is important to examine, for any given application, whether the expected benefit of increased grid resolution justifies the extra costs. In this study, we examine temperature and ozone observations and model predictions for three high ozone episodes that occurred over the northeastern United States during the summer of 1995. In the first set of simulations, the meteorological model RAMS4a was run with three two-way nested grids of 108/36/12 km grid spacing covering the United States and the photochemical model UAM-V was run with two grids of 36/12 km grid spacing covering the eastern United States. In the second set of simulations, RAMS4a was run with four two-way nested grids of 108/36/12/4 km grid spacing and UAM-V was run with three grids of 36/12/4 km grid spacing with the finest resolution covering the northeastern United States. Our analysis focuses on the comparison of model predictions for the finest grid domain of the simulations, namely, the region overlapping the 12 km and 4 km domains. A comparison of 12 km versus 4 km fields shows that the increased grid resolution leads to finer texture in the model predictions; however, comparisons of model predictions with observations do not reveal the expected improvement in the predictions. While high-resolution modeling has scientific merit and potential uses, the currently available monitoring networks, in conjunction with the scarceness of highly resolved spatial input data and the limitations of model formulation, do not allow confirmation of the expected superiority of the high-resolution model predictions.The U.S. Governments right to retain a non-exclusive royalty-free licence in and to any copyright is acknowledged.     

On the Estimation of Dispersal Kernels from Individual Mark-Recapture Data

We present a new method for estimating a distribution of dispersal displacements (a dispersal kernel) from mark-recapture data. One conventional method of calculating the dispersal kernel assumes that the distribution of displacements are Gaussian (e.g. resulting from a diffusion process) and that individuals remain within sampled areas. The first assumption prohibits an analysis of dispersal data that do not exhibit the Gaussian distribution (a common situation); the second assumption leads to underestimation of dispersal distance because individuals that disperse outside of sampling areas are never recaptured. Our method eliminates these two assumptions. In addition, the method can also accommodate mortality during a sampling period. This new method uses integrodifference equations to express the probability of spatial mark-recapture data; associated dispersal, survival, and recapture parameters are then estimated using a maximum likelihood method. We examined the accuracy of the estimators by applying the method to simulated data sets. Our method suggests designs for future mark-recapture experiments. Received: January 2004 / Revised: July 2005     

Probability distributions and statistical inference for axial data

Observations on axes which lack information on the direction of propagation are referred to as axial data. Such data are often encountered in enviromental sciences, e.g. observations on propagations of cracks or on faults in mining walls. Even though such observations are recorded as angles, circular probability models are inappropriate for such data since the constraint that observations lie only in [0, π) needs to be enforced. Probability models for such axial data are argued here to have a general structure stemming from that of wrapping a circular distribution on a semi-circle. In particular, we consider the most popular circular model, the von Mises or circular normal distribution, and derive the corresponding axial normal distribution. Certain properties of this distribution are established. Maximum likelihood estimation of its parameters are shown to be surprisingly, in contrast to trigonometric moment estimation, numerically quite appealing. Finally we illustrate our results by several real life axial data sets. Received: September 2004/ Revised: December 2004     

Efficient and sustainable management of complex forest ecosystems

A large range of models has been developed for the analysis of optimal forest management strategies, with the well-known Faustmann models dating back to the mid-19th century. To date, however, there has been relatively little attention for the implications of complex ecosystem dynamics for optimal forest management. This paper examines the implications of irreversible ecosystem responses for efficient and sustainable forest management. The paper is built around two forest models that comprise two ecosystem components, forest cover and topsoil, the interactions between these components, and the supply of the ecosystem services ‘wood’ and ‘erosion control’. The first model represents a forest that responds in a reversible way to overharvesting. In the second model, an additional ecological process has been included and the ecosystem irreversibly collapses below certain thresholds in forest cover and topsoil depth. The paper presents a general model, and demonstrates the implications of pursuing efficient as well as sustainable forest management for the two forest ecosystems. Both fixed and variable harvesting cycles are examined. Efficient and sustainable harvesting cycles are compared, and it is shown that irreversible ecosystem behaviour reduces the possibilities to reconcile efficient and sustainable forest management through a variable harvesting cycle.     

Studying crop sequences with CarrotAge,a HMM-based data mining software

We have developed a knowledge discovery system based on high-order hidden Markov models for analyzing spatio-temporal data bases. This system, named CarrotAge , takes as input an array of discrete data – the rows represent the spatial sites and the columns the time slots – and builds a partition together with its a posteriori probability. CarrotAge has been developed for studying the cropping patterns of a territory. It uses therefore an agricultural drench database, named Ter-Uti , which records every year the land-use category of a set of sites regularly spaced. The results of CarrotAge are interpreted by agronomists and used in research works linking agricultural land use and water management. Moreover, CarrotAge can be used to find out and study crop sequences in large territories, that is a main question for agricultural and environmental research, as discussed in this paper.     

Aphid population response to agricultural landscape change: A spatially explicit,individual-based model

A spatially explicit individual-based simulation model has been developed to represent aphid population dynamics in agricultural landscapes. The application of the model to Rhopalosiphum padi (L.) population dynamics is detailed, including an outline of the construction of the model, its parameterisation and validation. Over time, the aphids interact with the landscape and with one another. The landscape is modified by varying a simple pesticide regime, and the multi-scale spatial and temporal implications for a population of aphids is analysed. The results show that a spatial modelling approach that considers the effects on the individual of landscape properties and factors such as wind speed and wind direction provides novel insight into aphid population dynamics both spatially and temporally. This forms the basis for the development of further simulation models that can be used to analyse how changes in landscape structure impact upon important species distributions and population dynamics.     

Model-based combination of spatial information for stream networks

Evolutionary improvements in Geographic Information Systems (GIS) now routinely allow the management and mapping of spatial-temporal information. In response, the development of statistical models to combine information of different types and spatial support is of vital importance to environmental science. In this paper we develop a hierarchical spatial statistical model for environmental indicators of stream and river systems in the United States Mid-Atlantic Region by combining information from separate monitoring surveys, available contextual information on hydrologic units and remote sensing information. These models are used to estimate the indicators throughout the riverine system based on information from multiple sources and aggregate scales. The analysis is based on information underlying the Landscape Atlas of the mid-Atlantic region produced by the US Environmental Monitoring and Assessment Program (EMAP). We also combine information from two overlapping separate monitoring surveys, the EMAP Stream and River Survey and the Maryland Biological Streams Survey. We present a general framework for comparative distributional analysis based on the concept of a relative spatial distribution. As an application, the spatial model is used to predict spatial distributions and relative spatial distributions for a watershed.