In this paper, we construct a hierarchical model for spatial compositional
data, which is used to reconstruct past land-cover compositions (in terms of
coniferous forest, broadleaved forest, and unforested/open land) for five time
periods during the past $6\,000$ years over Europe. The model consists of a
Gaussian Markov Random Field (GMRF) with Dirichlet observations. A block
updated Markov chain Monte Carlo (MCMC), including an adaptive Metropolis
adjusted Langevin step, is used to estimate model parameters. The sparse
precision matrix in the GMRF provides computational advantages leading to a
fast MCMC algorithm. Reconstructions are obtained by combining pollen-based
estimates of vegetation cover at a limited number of locations with scenarios
of past deforestation and output from a dynamic vegetation model. To evaluate
uncertainties in the predictions a novel way of constructing joint confidence
regions for the entire composition at each prediction location is proposed. The
hierarchical model's ability to reconstruct past land cover is evaluated
through cross validation for all time periods, and by comparing reconstructions
for the recent past to a present day European forest map. The evaluation
results are promising and the model is able to capture known structures in past
land-cover compositions.
In this paper, we construct a hierarchical model for spatial compositional data, which is used to reconstruct past land-cover compositions (in terms of coniferous forest, broadleaved forest, and unforested/open land) for five time periods during the past $6\,000$ years over Europe. The model consists of a Gaussian Markov Random Field (GMRF) with Dirichlet observations. A block updated Markov chain Monte Carlo (MCMC), including an adaptive Metropolis adjusted Langevin step, is used to estimate model parameters. The sparse precision matrix in the GMRF provides computational advantages leading to a fast MCMC algorithm. Reconstructions are obtained by combining pollen-based estimates of vegetation cover at a limited number of locations with scenarios of past deforestation and output from a dynamic vegetation model. To evaluate uncertainties in the predictions a novel way of constructing joint confidence regions for the entire composition at each prediction location is proposed. The hierarchical model's ability to reconstruct past land cover is evaluated through cross validation for all time periods, and by comparing reconstructions for the recent past to a present day European forest map. The evaluation results are promising and the model is able to capture known structures in past land-cover compositions.
We look at the operational semantics of languages with interactive I/O through the glasses of constructive type theory. Following on from our earlier work on coinductive trace-based semantics for While, we define several big-step semantics for While with interactive I/O, based on resumptions and termination-sensitive weak bisimilarity. These require nesting inductive definitions in coinductive definitions, which is interesting both mathematically and from the point-of-view of implementation in a proof assistant. After first defining a basic semantics of statements in terms of resumptions with explicit internal actions (delays), we introduce a semantics in terms of delay-free resumptions that essentially removes finite sequences of delays on the fly from those resumptions that are responsive. Finally, we also look at a semantics in terms of delay-free resumptions supplemented with a silent divergence option. This semantics hinges on decisions between convergence and divergence and is only equivalent to the basic one classically. We have fully formalized our development in Coq.