Modelling spatial and temporal variability in a zero-inflated variable: The case of stone pine (Pinus pinea L.) cone production

Modelling masting habit, i.e. the spatial synchronized annual variability in fruit production, is a huge task due to two main circumstances: (1) the identification of main ecological factors controlling fruiting processes, and (2) the common departure of fruit data series from the main basic statistical assumptions of normality and independence. Stone pine (Pinus pinea L.) is one of the main species in the Mediterranean basin that is able to grow under hard limiting conditions (sandy soils and extreme continental climate), and typically defined as a masting species. Considering the high economical value associated with edible nut production, the masting habit of stone pine has been a main concern for the forest management of the species. In the present work we have used annual fruit data series from 740 stone pine trees measured during a 13 years period (1996-2008) in order: (a) to verify our main hypothesis pointing out to the existence of a weather control of the fruiting process in limiting environments, rather than resource depletion or endogenous inherent cycles; (b) to identify those site factors, stand attributes and climate events affecting specific traits involved in fruiting process; and (c) to construct a model for predicting spatial and temporal patterns of variability in stone pine cone production at different spatial extents as region, stand and tree. Given the nature of the data, the model has been formulated as zero-inflated log-normal, incorporating random components to carry out with the observed lack of independence. This model attains efficiencies close to 70-80% in predicting temporal and spatial variability at regional scale. Though efficiencies are reduced according to the spatial extent of the model, it leads to unbiased estimates and efficiencies over 35-50% when predicting annual yields at tree or stand scale, respectively. In this sense, the proposed model is a main tool for facilitating decision making in some management aspects such as the quantification of total amount of cones annually supplied to nut industry, design of cone harvest programs or the optimal application of seedling felling.