Modelling the spatial structure of forest stands by multivariate point processes with hierarchical interactions

A stochastic model is applied to describe the spatial structure of a forest stand. We aim at quantifying the strength of the competition process between the trees in terms of interaction within and between different size classes of trees using multivariate Gibbs point processes with hierarchical interactions introduced in [Högmander, H., Särkkä, A., 1999. Multitype spatial point patterns with hierarchical interactions. Biometrics 55, 1051–1058]. The new model overcomes the main limitation of the traditional use of the Gibbs models allowing to describe systems with non-symmetric interactions between different objects. When analyzing interactions between neighbouring trees it is natural to assume that the size of a tree determines its hierarchical level: the largest trees are not influenced by any other trees than the trees in the same size class, while trees in the other size classes are influenced by the other trees in the same class as well as by all larger trees. In this paper, we describe a wide range of Gibbs models with both hierarchical and non-hierarchical interactions as well as a simulation algorithm and a parameter estimation procedure for the hierarchical models. We apply the hierarchical interaction model to the analysis of forest data consisting of locations and diameters of tree stems.     

Efficient Markov chain Monte Carlo sampling for hierarchical hidden Markov models

Traditional Markov chain Monte Carlo (MCMC) sampling of hidden Markov models (HMMs) involves latent states underlying an imperfect observation process, and generates posterior samples for top-level parameters concurrently with nuisance latent variables. When potentially many HMMs are embedded within a hierarchical model, this can result in prohibitively long MCMC runtimes. We study combinations of existing methods, which are shown to vastly improve computational efficiency for these hierarchical models while maintaining the modeling flexibility provided by embedded HMMs. The methods include discrete filtering of the HMM likelihood to remove latent states, reduced data representations, and a novel procedure for dynamic block sampling of posterior dimensions. The first two methods have been used in isolation in existing application-specific software, but are not generally available for incorporation in arbitrary model structures. Using the NIMBLE package for R, we develop and test combined computational approaches using three examples from ecological capture–recapture, although our methods are generally applicable to any embedded discrete HMMs. These combinations provide several orders of magnitude improvement in MCMC sampling efficiency, defined as the rate of generating effectively independent posterior samples. In addition to being computationally significant for this class of hierarchical models, this result underscores the potential for vast improvements to MCMC sampling efficiency which can result from combinations of known algorithms.     

Natural control mechanisms in models of aquatic ecosystems

Based on cybernetic categories of natural control mechanisms, four generations of ecosystem models are distinguished: feed-forward, feedback, self-adaptation and self-organization models. The analysis of the natural control mechanisms in aquatic ecosystems suggests that different processes are controlled in different ways, and, although the four mechanisms were identified in historical sequence, they all operate simultaneously. The concept of self-organization of an ecosystem is introduced and specified for a model of an aquatic pelagic ecosystem. The concept of the ecosystem as a multilayer, multigoal and multiechelon hierarchical system with hierarchy of the levels of biological organization is also introduced.     

Natural control mechanisms in models of aquatic ecosystems

Based on cybernetic categories of natural control mechanisms, four generations of ecosystem models are distinguished: feed-forward, feedback, self-adaptation and self-organization models. The analysis of the natural control mechanisms in aquatic ecosystems suggests that different processes are controlled in different ways, and, although the four mechanisms were identified in historical sequence, they all operate simultaneously. The concept of self-organization of an ecosystem is introduced and specified for a model of an aquatic pelagic ecosystem. The concept of the ecosystem as a multilayer, multigoal and multiechelon hierarchical system with hierarchy of the levels of biological organization is also introduced.     

Using the maximum principle of impulse control for ecology-economical models

In this work we present mathematical models for population of single cohort and homogeneous animals. Investigating these mathematical models, we determine structure of optimal impulsive control which used maximum principle for optimal processes with impulse control.     

Poor-data and data-poor species stock assessment using a Bayesian hierarchical approach

Appropriate inference for stocks or species with low-quality data (poor data) or limited data (data poor) is extremely important. Hierarchical Bayesian methods are especially applicable to small-area, small-sample-size estimation problems because they allow poor-data species to borrow strength from species with good-quality data. We used a hammerhead shark complex as an example to investigate the advantages of using hierarchical Bayesian models in assessing the status of poor-data and data-poor exploited species. The hammerhead shark complex (Sphyrna spp.) along the Atlantic and Gulf of Mexico coasts of the United States is composed of three species: the scalloped hammerhead (S. lewini), the great hammerhead (S. mokarran), and the smooth hammerhead (S. zygaena) sharks. The scalloped hammerhead comprises 70-80% of the catch and has catch and relative abundance data of good quality, whereas great and smooth hammerheads have relative abundance indices that are both limited and of low quality presumably because of low stock density and limited sampling. Four hierarchical Bayesian state-space surplus production models were developed to simulate variability in population growth rates, carrying capacity, and catchability of the species. The results from the hierarchical Bayesian models were considerably more robust than those of the nonhierarchical models. The hierarchical Bayesian approach represents an intermediate strategy between traditional models that assume different population parameters for each species and those that assume all species share identical parameters. Use of the hierarchical Bayesian approach is suggested for future hammerhead shark stock assessments and for modeling fish complexes with species-specific data, because the poor-data species can borrow strength from the species with good data, making the estimation more stable and robust.     

Optimal model averaging by data perturbation for spatial noisy data

Geostatistical models play an important role in spatial data analysis, in which model selection is inevitable. Model selection methods, such as AIC and BIC, are popular for selecting appropriate models. In recent years, some model averaging methods, such as smoothed AIC and smoothed BIC, are also applied to spatial data models. However, the corresponding averaging estimators are outperformed by optimal model averaging estimators (Hansen in Econometrica 75:1175–1189, 2007) for the ordinary linear models. Therefore, this paper focuses on the optimal model averaging method for geostatistical models. We propose a weight choice criterion for the model averaging estimator on the basis of the generalized degrees of freedom and data perturbation technique. We further theoretically prove the resultant estimator is asymptotically optimal in terms of the mean squared error, and numerically demonstrate its satisfactory performance. Finally, the proposed method is applied to a mercury data set.     

Hierarchical modeling of abundance in closed population capture–recapture models under heterogeneity

Hierarchical modeling of abundance in space or time using closed-population mark-recapture under heterogeneity (model \(\hbox {M}_{\text {h}}\) ) presents two challenges: (i) finding a flexible likelihood in which abundance appears as an explicit parameter and (ii) fitting the hierarchical model for abundance. The first challenge arises because abundance not only indexes the population size, it also determines the dimension of the capture probabilities in heterogeneity models. A common approach is to use data augmentation to include these capture probabilities directly into the likelihood and fit the model using Bayesian inference via Markov chain Monte Carlo (MCMC). Two such examples of this approach are (i) explicit trans-dimensional MCMC, and (ii) superpopulation data augmentation. The superpopulation approach has the advantage of simple specification that is easily implemented in BUGS and related software. However, it reparameterizes the model so that abundance is no longer included, except as a derived quantity. This is a drawback when hierarchical models for abundance, or related parameters, are desired. Here, we analytically compare the two approaches and show that they are more closely related than might appear superficially. We exploit this relationship to specify the model in a way that allows us to include abundance as a parameter and that facilitates hierarchical modeling using readily available software such as BUGS. We use this approach to model trends in grizzly bear abundance in Yellowstone National Park from 1986 to 1998.     

Zero-inflated models with application to spatial count data

Count data arises in many contexts. Here our concern is with spatial count data which exhibit an excessive number of zeros. Using the class of zero-inflated count models provides a flexible way to address this problem. Available covariate information suggests formulation of such modeling within a regression framework. We employ zero-inflated Poisson regression models. Spatial association is introduced through suitable random effects yielding a hierarchical model. We propose fitting this model within a Bayesian framework considering issues of posterior propriety, informative prior specification and well-behaved simulation based model fitting. Finally, we illustrate the model fitting with a data set involving counts of isopod nest burrows for 1649 pixels over a portion of the Negev desert in Israel.     

Clustering species using a model of population dynamics and aggregation theory

The high species diversity of some ecosystems like tropical rainforests goes in pair with the scarcity of data for most species. This hinders the development of models that require enough data for fitting. The solution commonly adopted by modellers consists in grouping species to form more sizeable data sets. Classical methods for grouping species such as hierarchical cluster analysis do not take account of the variability of the species characteristics used for clustering. In this study a clustering method based on aggregation theory is presented. It takes account of the variability of species characteristics by searching for the grouping that minimizes the quadratic error (square bias plus variance) of some model’s prediction. This method allows one to check whether the gain in variance brought by data pooling compensate for the bias that it introduces. This method was applied to a data set on 94 tree species in a tropical rainforest in French Guiana, using a Usher matrix model to predict species dynamics. An optimal trade-off between bias and variance was found when grouping species. Grouping species appeared to decrease the quadratic error, except when the number of groups was very small. This clustering method yielded species groups similar to those of the hierarchical cluster analysis using Ward’s method when variance was small, that is when the number of groups was small.     

Tropical deforestation in Madagascar: analysis using hierarchical,spatially explicit,Bayesian regression models

Establishing cause–effect relationships for deforestation at various scales has proven difficult even when rates of deforestation appear well documented. There is a need for better explanatory models, which also provide insight into the process of deforestation. We propose a novel hierarchical modeling specification incorporating spatial association. The hierarchical aspect allows us to accommodate misalignment between the land-use (response) data layer and explanatory data layers. Spatial structure seems appropriate due to the inherently spatial nature of land use and data layers explaining land use. Typically, there will be missing values or holes in the response data. To accommodate this we propose an imputation strategy. We apply our modeling approach to develop a novel deforestation model for the eastern wet forested zone of Madagascar, a global rain forest “hot spot”. Using five data layers created for this region, we fit a suitable spatial hierarchical model. Though fitting such models is computationally much more demanding than fitting more standard models, we show that the resulting interpretation is much richer. Also, we employ a model choice criterion to argue that our fully Bayesian model performs better than simpler ones. To the best of our knowledge, this is the first work that applies hierarchical Bayesian modeling techniques to study deforestation processes. We conclude with a discussion of our findings and an indication of the broader ecological applicability of our modeling style.     

Combining functional data with hierarchical Gaussian process models

Gaussian process models have been used in applications ranging from machine learning to fisheries management. In the Bayesian framework, the Gaussian process is used as a prior for unknown functions, allowing the data to drive the relationship between inputs and outputs. In our research, we consider a scenario in which response and input data are available from several similar, but not necessarily identical, sources. When little information is known about one or more of the populations it may be advantageous to model all populations together. We present a hierarchical Gaussian process model with a structure that allows distinct features for each source as well as shared underlying characteristics. Key features and properties of the model are discussed and demonstrated in a number of simulation examples. The model is then applied to a data set consisting of three populations of Rotifer Brachionus calyciflorus Pallas. Specifically, we model the log growth rate of the populations using a combination of lagged population sizes. The various lag combinations are formally compared to obtain the best model inputs. We then formally compare the leading hierarchical Gaussian process model with the inferential results obtained under the independent Gaussian process model.     

An ecological model of the habitat mosaic in estuarine nursery areas: Part I—Interaction of dispersal theory and habitat variability in describing juvenile fish distributions

Identification of critical habitat in estuarine nursery areas is an important conservation and management objective. Habitat can be viewed as a mosaic of both temporally variable environmental features and spatially variable structural features that combine to define optimal habitat. Effective models of juvenile distributions should account for individual movement, as well as the full suite of habitat variability including both spatial and temporal components. We have extended a terrestrial model of small-scale movement patterns to describe habitat choices of an index juvenile fish in an estuarine nursery system. Movement of small juvenile fishes was found to be influenced by both spatial and temporal patterns in habitat quality, and it was a balanced mix of both that resulted in an optimal distribution. Fishes that perceive habitat on a scale much smaller than the scale of spatial heterogeneity may respond to temporal change as a movement cue allowing for more deterministic outcomes at larger scales despite perceptual limitations. These model outcomes suggest a hierarchical approach is best for describing habitat choice in juvenile fishes and this approach will be used in the future to explore individual and population responses to predictable habitat change.     

Uncertainty, learning, and the optimal management of wildlife

Wildlife management is limited by uncontrolled and often unrecognized environmental variation, by limited capabilities to observe and control animal populations, and by a lack of understanding about the biological processes driving population dynamics. In this paper I describe a comprehensive framework for management that includes multiple models and likelihood values to account for structural uncertainty, along with stochastic factors to account for environmental variation, random sampling, and partial controllability. Adaptive optimization is developed in terms of the optimal control of incompletely understood populations, with the expected value of perfect information measuring the potential for improving control through learning. The framework for optimal adaptive control is generalized by including partial observability and non-adaptive, sample-based updating of model likelihoods. Passive adaptive management is derived as a special case of constrained adaptive optimization, representing a potentially efficient suboptimal alternative that nonetheless accounts for structural uncertainty.     

A Bayesian hierarchical framework for calibrating aquatic biogeochemical models

Model practitioners increasingly place emphasis on rigorous quantitative error analysis in aquatic biogeochemical models and the existing initiatives range from the development of alternative metrics for goodness of fit, to data assimilation into operational models, to parameter estimation techniques. However, the treatment of error in many of these efforts is arguably selective and/or ad hoc. A Bayesian hierarchical framework enables the development of robust probabilistic analysis of error and uncertainty in model predictions by explicitly accommodating measurement error, parameter uncertainty, and model structure imperfection. This paper presents a Bayesian hierarchical formulation for simultaneously calibrating aquatic biogeochemical models at multiple systems (or sites of the same system) with differences in their trophic conditions, prior precisions of model parameters, available information, measurement error or inter-annual variability. Our statistical formulation also explicitly considers the uncertainty in model inputs (model parameters, initial conditions), the analytical/sampling error associated with the field data, and the discrepancy between model structure and the natural system dynamics (e.g., missing key ecological processes, erroneous formulations, misspecified forcing functions). The comparison between observations and posterior predictive monthly distributions indicates that the plankton models calibrated under the Bayesian hierarchical scheme provided accurate system representations for all the scenarios examined. Our results also suggest that the Bayesian hierarchical approach allows overcoming problems of insufficient local data by “borrowing strength” from well-studied sites and this feature will be highly relevant to conservation practices of regions with a high number of freshwater resources for which complete data could never be practically collected. Finally, we discuss the prospect of extending this framework to spatially explicit biogeochemical models (e.g., more effectively connect inshore with offshore areas) along with the benefits for environmental management, such as the optimization of the sampling design of monitoring programs and the alignment with the policy practice of adaptive management.     

Managing economic and social profit of cooperative models in three-echelon reverse supply chain for waste electrical and electronic equipment

In addition to maximizing economic benefits, reverse supply chains should further seek to maximize social benefits by increasing the quantity of waste electrical and electronic equipment (WEEE). The paper investigates cooperative models with different parties in a three-echelon reverse supply chain for WEEE consisting of a single collector, a single remanufacturer, and two retailers based on complete information. In addition, the optimal decisions of four cooperative models and the effect of the market demand of remanufactured WEEE products and the market share of two retailers on the optimal decisions are discussed. The results indicate that optimal total channel profit and recycle quantity in a reverse supply chain are maximized in a centralized model. The optimal total channel profit and recycle quantity increase with an increase in the market demand of remanufactured WEEE products. The three-echelon reverse supply chain consisting of duopolistic retailers maximizes total channel profit and recycle quantity in a reverse supply chain for WEEE.

Open image in new window

     

Sport and commercial fishing conflicts: A theoretical analysis

A recreational sector is added to a standard commercial fishing optimal control model to identify public decision variables which should be considered when determining optimal population levels and allocating harvestable fish between sport and commercial users. Both linear and nonlinear models are presented. A predator-prey relationship is also considered. Results derived from the models indicate that shortcomings exist with current economic inputs to policy making because of failure to consider relevant bioeconomic relationships. Future research topics to remedy this are discussed.     

An estimating function approach to the inference of catch-effort models

A class of catch-effort models, which allows for heterogeneous removal probabilities, is proposed for closed populations. The model includes three types of removal probabilities: multiplicative, Poisson and logistic. The usual removal and generalized removal models then become special cases. The equivalence of the proposed model and a special type of capture-recapture model is discussed. A unified estimating function approach is used to estimate the initial population size. For the homogeneous model, the resulting population size estimator based on optimal estimating functions is asymptotically equivalent to the maximum likelihood estimator. One advantage for our approach is that it can be extended to handle the heterogeneous populations in which the maximum likelihood estimators do not exist. The bootstrap method is applied to construct variance estimators and confidence intervals. We illustrate the method by two real data examples. Results of a simulation study investigating the performance of the proposed estimation procedure are presented.     

Using land as a control variable in density-dependent bioeconomic models

The bioeconomic analysis of endangered species without consumptive values can be problematic when analysed with density-dependent models that assume a fixed environment size. Most bioeconomic models use harvest as a control variable, yet when modelling non-harvestable species, frequently the only variable under control of conservationists is the quantity of habitat to be made available. The authors explore the implications of this in a model developed to analyse the potential population recovery of New Zealand’s yellow-eyed penguin. The penguin faces severe competition with man for the terrestrial resources required for breeding and has declined in population to perilously low levels. The model was developed to estimate the land use required for recovery and preservation of the species and to compare the results to current tourism-driven conservation efforts. It is demonstrated that land may serve as a useful control variable in bioeconomic models and that such a model may be useful for determining whether sufficient incentives exist to preserve a species. However, the model may generate less useful results for providing a specific estimate of the optimal allocation of land to such a species.     

Hierarchical Bayesian space-time models

Space-time data are ubiquitous in the environmental sciences. Often, as is the case with atmo- spheric and oceanographic processes, these data contain many different scales of spatial and temporal variability. Such data are often non-stationary in space and time and may involve many observation/prediction locations. These factors can limit the effectiveness of traditional space- time statistical models and methods. In this article, we propose the use of hierarchical space-time models to achieve more flexible models and methods for the analysis of environmental data distributed in space and time. The first stage of the hierarchical model specifies a measurement- error process for the observational data in terms of some 'state' process. The second stage allows for site-specific time series models for this state variable. This stage includes large-scale (e.g. seasonal) variability plus a space-time dynamic process for the anomalies'. Much of our interest is with this anomaly proc ess. In the third stage, the parameters of these time series models, which are distributed in space, are themselves given a joint distribution with spatial dependence (Markov random fields). The Bayesian formulation is completed in the last two stages by speci- fying priors on parameters. We implement the model in a Markov chain Monte Carlo framework and apply it to an atmospheric data set of monthly maximum temperature.