Quantile regression for modelling the impact of climate in cork growth quantiles in Portugal

The annual growth and the thickness of cork are known to be highly variable between trees located in the same geographical location. Researching how climate variables affect different trees within the same site is a step forward for the management of cork production since current knowledge focusses only on the average tree response. Quantile regression methodology was applied for the first time to a large data set containing measurements of cork growth, sampled in 35 stands across the cork oak distribution area in Portugal. This methodology proved to be useful for testing the hypothesis raised: does climate affect differently the annual cork growth, and ultimately cork thickness of individual trees located in the same stand? Estimating the amount of cork produced by one stand that has the required thickness for the production of natural cork stoppers is essential to support cork oak management. However, no model, before this work, had been developed to provide managers with this information. A downward parabolic relationship between annual cork growth and annual precipitation was determined for all quantiles, with optimum annual average precipitation value ranging from 1103 to 1007 mm. April to August monthly temperatures, spring average temperature or summer average temperature, showed a negative relationship with annual cork growth, in particular for lower quantiles. Maximum annual temperature was shown to negatively affect annual cork thickness, in particular for the trees under the 6th quantile. The ratio between annual precipitation and average temperature, that define the Lang index (LI), showed a downward parabolic relationship with annual cork growth. Best cork growth conditions are found for Lang index values around 60, corresponding for the transition between semi-arid climate and humid climate. The application of the final model developed for estimating cork thickness of an eight years' cork growth period allowed the prediction and mapping of the percentage of cork suitable for natural cork stopper production. It showed that higher values are expected in the Southern and Central coastal regions and along the Tagus River basin. The Northern coastal and mountain regions, characterised by Lang index values higher to 60 (humid climates), present lower estimated values for the percentage of cork suitable for natural cork stopper production. The estimated values are expected to be reduced under climate change scenarios in the Southern and Central coastal regions.