Local-scale topoclimate effects on treeline elevations: a country-wide investigation of New Zealand’s southern beech treelines

Although treeline elevations are limited globally by growing season temperature, at regional scales treelines frequently deviate below their climatic limit. The cause of these deviations relate to a host of climatic, disturbance, and geomorphic factors that operate at multiple scales. The ability to disentangle the relative effects of these factors is currently hampered by the lack of reliable topoclimatic data, which describe how regional climatic characteristics are modified by topographic effects in mountain areas. In this study we present an analysis of the combined effects of local- and regional-scale factors on southern beech treeline elevation variability at 28 study areas across New Zealand. We apply a mesoscale atmospheric model to generate local-scale (200 m) meteorological data at these treelines and, from these data, we derive a set of topoclimatic indices that reflect possible detrimental and ameliorative influences on tree physiological functioning. Principal components analysis of meteorological data revealed geographic structure in how study areas were situated in multivariate space along gradients of topoclimate. Random forest and conditional inference tree modelling enabled us to tease apart the relative effects of 17 explanatory factors on local-scale treeline elevation variability. Overall, modelling explained about 50% of the variation in treeline elevation variability across the 28 study areas, with local landform and topoclimatic effects generally outweighing those from regional-scale factors across the 28 study areas. Further, the nature of the relationships between treeline elevation variability and the explanatory variables were complex, frequently non-linear, and consistent with the treeline literature. To our knowledge, this is the first study where model-generated meteorological data, and derived topoclimatic indices, have been developed and applied to explain treeline variation. Our results demonstrate the potential of such an approach for ecological research in mountainous environments.