Multiple drivers of large-scale lichen decline in boreal forest canopies.

Thin, hair-like lichens (Alectoria, Bryoria, Usnea) form conspicuous epiphyte communities across the boreal biome. These poikilohydric organisms provide important ecosystem functions and are useful indicators of global change. We analyse how environmental drivers influence changes in occurrence and length of these lichens on Norway spruce (Picea abies) over 10 years in managed forests in Sweden using data from >6000 trees. Alectoria and Usnea showed strong declines in southern-central regions, whereas Bryoria declined in northern regions. Overall, relative loss rates across the country ranged from 1.7% per year in Alectoria to 0.5% in Bryoria. These losses contrasted with increased length of Bryoria and Usnea in some regions. Occurrence trajectories (extinction, colonization, presence, absence) on remeasured trees correlated best with temperature, rain, nitrogen deposition, and stand age in multinomial logistic regression models. Our analysis strongly suggests that industrial forestry, in combination with nitrogen, is the main driver of lichen declines. Logging of forests with long continuity of tree cover, short rotation cycles, substrate limitation and low light in dense forests are harmful for lichens. Nitrogen deposition has decreased but is apparently still sufficiently high to prevent recovery. Warming correlated with occurrence trajectories of Alectoria and Bryoria, likely by altering hydration regimes and increasing respiration during autumn/winter. The large-scale lichen decline on an important host has cascading effects on biodiversity and function of boreal forest canopies. Forest management must apply a broad spectrum of methods, including uneven-aged continuous cover forestry and retention of large patches, to secure the ecosystem functions of these important canopy components under future climates. Our findings highlight interactions among drivers of lichen decline (forestry, nitrogen, climate), functional traits (dispersal, lichen colour, sensitivity to nitrogen, water storage), and population processes (extinction/colonization).

A comparison of two procedures to estimate three basic monitoring landscape metrics for monitoring.

An interesting alternative to wall-to-wall mapping approaches for the estimation of landscape metrics is to use sampling. Sample-based approaches are cost-efficient, and measurement errors can be reduced considerably. The previous efforts of sample-based estimation of landscape metrics have mainly been focused on data collection methods, but in this study, we consider two estimation procedures. First, landscape metrics of interest are calculated separately for each sampled image and then the image values are averaged to obtain an estimate of the entire landscape (separated procedure, SP). Second, metric components are calculated in all sampled images and then the aggregated values are inserted into the landscape metric formulas (aggregated procedure, AP). The national land cover map (NLCM) of Sweden, reflecting the status of land cover in the year 2000, was used to provide population information to investigate the statistical performance of the estimation procedures. For this purpose, sampling simulation with a large number of replications was used. For all three landscape metrics, the second procedure (AP) produced a lower relative RMSE and bias than the first one (SP). A smaller sample unit size (50 ha) produced larger bias than a larger one (100 ha), whereas a smaller sample unit size produced a lower variance than a larger sample unit. The efficiency of a metric estimator is highly related to the degree of landscape fragmentation and the selected procedure. Incorporating information from all of the sampled images into a single one (aggregated procedure, AP) is one way to improve the statistical performance of estimators.

Spatially balanced sampling through the pivotal method.

A simple method to select a spatially balanced sample using equal or unequal inclusion probabilities is presented. For populations with spatial trends in the variables of interest, the estimation can be much improved by selecting samples that are well spread over the population. The method can be used for any number of dimensions and can hence also select spatially balanced samples in a space spanned by several auxiliary variables. Analysis and examples indicate that the suggested method achieves a high degree of spatial balance and is therefore efficient for populations with trends.