Coming up with Bayesian models for spatial data is easy, but performing inference with them can be challenging. Writing fast inference code for a complex spatial model with realistically‐sized datasets from scratch is time‐consuming, and if changes are made to the model, there is little guarantee that the code performs well. The key advantages of R‐INLA are the ease with which complex models can be created and modified, without the need to write complex code, and the speed at which inference can be done even for spatial problems with hundreds of thousands of observations. R‐INLA handles latent Gaussian models, where fixed effects, structured and unstructured Gaussian random effects are combined linearly in a linear predictor, and the elements of the linear predictor are observed through one or more likelihoods. The structured random effects can be both standard areal model such as the Besag and the BYM models, and geostatistical models from a subset of the Matérn Gaussian random fields. In this review, we discuss the large success of spatial modeling with R‐INLA and the types of spatial models that can be fitted, we give an overview of recent developments for areal models, and we give an overview of the stochastic partial differential equation (SPDE) approach and some of the ways it can be extended beyond the assumptions of isotropy and separability. In particular, we describe how slight changes to the SPDE approach leads to straight‐forward approaches for nonstationary spatial models and nonseparable space–time models. This article is categorized under: Statistical and Graphical Methods of Data Analysis > Bayesian Methods and Theory Statistical Models > Bayesian Models Data: Types and Structure > Massive Data

No abstract

The classical tools in spatial statistics are stationary models, like the Matérn field. However, in some applications there are boundaries, holes, or physical barriers in the study area, e.g. a coastline, and stationary models will inappropriately smooth over these features, requiring the use of a non-stationary model.We propose a new model, the Barrier model, which is different from the established methods as it is not based on the shortest distance around the physical barrier, nor on boundary conditions. The Barrier model is based on viewing the Matérn correlation, not as a correlation function on the shortest distance between two points, but as a collection of paths through a Simultaneous Autoregressive (SAR) model. We then manipulate these local dependencies to cut off paths that are crossing the physical barriers. To make the new SAR well behaved, we formulate it as a stochastic partial differential equation (SPDE) that can be discretised to represent the Gaussian field, with a sparse precision matrix that is automatically positive definite.The main advantage with the Barrier model is that the computational cost is the same as for the stationary model. The model is easy to use, and can deal with both sparse data and very complex barriers, as shown in an application in the Finnish Archipelago Sea. Additionally, the Barrier model is better at reconstructing the modified Horseshoe test function than the standard models used in R-INLA.

This study evaluates the dynamics of soil organic carbon (SOC) under perennial crops across the globe. It quantifies the effect of change from annual to perennial crops and the subsequent temporal changes in SOC stocks during the perennial crop cycle. It also presents an empirical model to estimate changes in the SOC content under crops as a function of time, land use, and site characteristics. We used a harmonized global dataset containing paired‐comparison empirical values of SOC and different types of perennial crops (perennial grasses, palms, and woody plants) with different end uses: bioenergy, food, other bio‐products, and short rotation coppice. Salient outcomes include: a 20‐year period encompassing a change from annual to perennial crops led to an average 20% increase in SOC at 0–30 cm (6.0 ± 4.6 Mg/ha gain) and a total 10% increase over the 0–100 cm soil profile (5.7 ± 10.9 Mg/ha). A change from natural pasture to perennial crop decreased SOC stocks by 1% over 0–30 cm (−2.5 ± 4.2 Mg/ha) and 10% over 0–100 cm (−13.6 ± 8.9 Mg/ha). The effect of a land use change from forest to perennial crops did not show significant impacts, probably due to the limited number of plots; but the data indicated that while a 2% increase in SOC was observed at 0–30 cm (16.81 ± 55.1 Mg/ha), a decrease in 24% was observed at 30–100 cm (−40.1 ± 16.8 Mg/ha). Perennial crops generally accumulate SOC through time, especially woody crops; and temperature was the main driver explaining differences in SOC dynamics, followed by crop age, soil bulk density, clay content, and depth. We present empirical evidence showing that the FAO perennialization strategy is reasonable, underscoring the role of perennial crops as a useful component of climate change mitigation strategies.

No abstract