Modeling spatiotemporal forest health monitoring data

Forest health monitoring schemes were set Lip across Europe in the 1980s in response to concerns about air pollution-related forest dieback (Waldsterben) and have continued since then. Recent threats to forest health are climatic extremes likely due to global climate change and increased ground ozone levels and nitrogen deposition. We model yearly data on tree crown defoliation, an indicator of tree health, from a monitoring survey carried Out in Baden-Wurttemberg, Germany since 1983. On a changing irregular grid, defoliation and other sue-specific variables are recorded, In Baden-Wurttemberg, the temporal trend of defoliation differs among areas because of site characteristics and pollution levels, making it necessary to allow for space-time interaction in the model. For this purpose, we propose using generalized additive mixed models (GAMMs) incorporating scale-invariant tensor product smooths of the space-time dimensions. The space-time smoother allows separate smoothing parameters and penalties for the space and time dimensions and thus avoids the need to make arbitrary or ad hoc choices about the relative scaling of space and time. The approach of using a space-time smoother has intuitive appeal, making it easy to explain and interpret when communicating the results to nonstaticians, such as environmental policy makers. The model incorporates a nonlinear effect for mean tree age. the most important predictor. allowing the separation of trends in time. which may be pollution-related, from trends that relate purely to the aging of the survey population. In addition to a temporal trend due to site characteristics and other conditions modeled with the space-time Smooth, We account for random temporal correlation at site level by an autoregressive moving average (ARMA) process. Model selection is carried Out using the Bayes information criterion (BIC). and the adequacy of the assumed spatial and temporal error Structure is investigated with the empirical semivariogram and the empirical autocorrelation function.