Landscape disturbance history in Westland, New Zealand: The importance of …

Research Papers, Lincoln University

This thesis investigates landscape disturbance history in Westland since 1350 AD. Specifically, I test the hypothesis that large-magnitude regional episodes of natural disturbance have periodically devastated portions of the landscape and forest, and that these were caused by infrequent earthquakes along the Alpine Fault. Forest stand history reconstruction was used to determine the timing and extent of erosion and sedimentation events that initiated new forest cohorts in a 1412 ha study area in the Karangarua River catchment, south Westland. Over 85 % of the study area was disturbed sufficiently by erosion/sedimentation since 1350 AD to initiate new forest cohorts. During this time four episodes of catchment-wide disturbance impacted the study area, and these took place about 1825 AD ± 5 years (Ruera episode), 1715 AD ± 5 years (Sparkling episode), 1615 AD ± 5 years (McTaggart episode), and 1445 AD ± 15 years (Junction episode). The three most recent episodes disturbed 10 %, 35-40 % and 32-50 % respectively of the study area. The Junction episode disturbed at least 6 % of the study area, but elimination of evidence by more recent disturbances prevented an upper limit being defined. The three earliest episodes correspond to the date-ranges for three Alpine Fault earthquakes from geological data, and are the only episodes of disturbance within each date-range. An earthquake cause is also consistent with features of the disturbance record: large portions of the study area were disturbed, disturbance occurred on all types 'of landforms, and terrace surfaces were abandoned upstream of the Alpine Fault. On this basis erosion/sedimentation induced by Alpine Fault earthquakes has disturbed 14-20 % of the land surface in the study area per century. Storms and other non-seismic erosional processes have disturbed 3-4 % per century. To examine the importance of the Alpine Fault earthquakes to forest disturbance throughout Westland, I collated all available data on conifer stand age structures in the region and identified dates of disturbance events from 55 even-aged cohorts of trees. Three region-wide episodes of forest disturbance since 1350 AD were found in this sample, and these matched the three Alpine Fault earthquake-caused episodes found in the Karangarua. Forest disturbance at these times was widespread across Westland over at least 200 km from Paringa to Hokitika, and originated from both tree fall and erosion processes. This disturbance history can explain the long-observed regional conifer forest pattern in Westland, of a predominance of similar-sized stands of trees and a relative lack of small-sized (young) stands. The many similar-sized stands are a consequence of synchronous forest disturbance and re-establishment accompanying the infrequent Alpine Fault earthquakes, while the dominance of mature stands of trees and relative lack of young small-sized trees in stands is explained by the long lapsed time since the last Alpine Fault earthquake (c. 280 years). I applied the landscape disturbance history information to the existing geological data to reconstruct the paleoseismicity of the Alpine Fault since 1350 AD. Best estimates for the timing of the most recent three rupture events from these data are 1715 AD ± 5 years, 1615 AD ± 5 years and 1445 AD ± 15 years. Earthquake recurrence intervals were variable, ranging from about 100 years to at least 280 years (the lapsed time since the last event). All three events caused forest and geomorphic disturbance over at least a 200 km section of Fault between the Karangarua and Hokitika Rivers, and were probably single rupture events. Suppressions in cross dated tree-ring chronologies in the western South Island suggest that the last rupture occurred in 1717 AD, and extended as a single rupture from Haupiri to Fiordland, a distance along the Fault of 375 km.

Multi-scale influences of abiotic factors on alpine treelines

Research Papers, Lincoln University

Globally, the maximum elevations at which treelines are observed to occur coincide with a 6.4 °C soil isotherm. However, when observed at finer scales, treelines display a considerable degree of spatial complexity in their patterns across the landscape and are often found occurring at lower elevations than expected relative to the global-scale pattern. There is still a lack of understanding of how the abiotic environment imposes constraints on treeline patterns, the scales at which different effects are acting, and how these effects vary over large spatial extents. In this thesis, I examined abrupt Nothofagus treelines across seven degrees of latitude in New Zealand in order to investigate two broad questions: (1) What is the nature and extent of spatial variability in Nothofagus treelines across the country? (2) How is this variation associated with abiotic variation at different spatial scales? A range of GIS, statistical, and atmospheric modelling methods were applied to address these two questions. First, I characterised Nothofagus treeline patterns at a 15x15km scale across New Zealand using a set of seven, GIS-derived, quantitative metrics that describe different aspects of treeline position, shape, spatial configuration, and relationships with adjacent vegetation. Multivariate clustering of these metrics revealed distinct treeline types that showed strong spatial aggregation across the country. This suggests a strong spatial structuring of the abiotic environment which, in turn, drives treeline patterns. About half of the multivariate treeline metric variation was explained by patterns of climate, substrate, topographic and disturbance variability; on the whole, climatic and disturbance factors were most influential. Second, I developed a conceptual model that describes how treeline elevation may vary at different scales according to three categories of effects: thermal modifying effects, physiological stressors, and disturbance effects. I tested the relevance of this model for Nothofagus treelines by investigating treeline elevation variation at five nested scales (regional to local) using a hierarchical design based on nested river catchments. Hierarchical linear modelling revealed that the majority of the variation in treeline elevation resided at the broadest, regional scale, which was best explained by the thermal modifying effects of solar radiation, mountain mass, and differences in the potential for cold air ponding. Nonetheless, at finer scales, physiological and disturbance effects were important and acted to modify the regional trend at these scales. These results suggest that variation in abrupt treeline elevations are due to both broad-scale temperature-based growth limitation processes and finer-scale stress- and disturbance-related effects on seedling establishment. Third, I explored the applicability of a meso-scale atmospheric model, The Air Pollution Model (TAPM), for generating 200 m resolution, hourly topoclimatic data for temperature, incoming and outgoing radiation, relative humidity, and wind speeds. Initial assessments of TAPM outputs against data from two climate station locations over seven years showed that the model could generate predictions with a consistent level of accuracy for both sites, and which agreed with other evaluations in the literature. TAPM was then used to generate data at 28, 7x7 km Nothofagus treeline zones across New Zealand for January (summer) and July (winter) 2002. Using mixed-effects linear models, I determined that both site-level factors (mean growing season temperature, mountain mass, precipitation, earthquake intensity) and local-level landform (slope and convexity) and topoclimatic factors (solar radiation, photoinhibition index, frost index, desiccation index) were influential in explaining variation in treeline elevation within and among these sites. Treelines were generally closer to their site-level maxima in regions with higher mean growing season temperatures, larger mountains, and lower levels of precipitation. Within sites, higher treelines were associated with higher solar radiation, and lower photoinhibition and desiccation index values, in January, and lower desiccation index values in July. Higher treelines were also significantly associated with steeper, more convex landforms. Overall, this thesis shows that investigating treelines across extensive areas at multiple study scales enables the development of a more comprehensive understanding of treeline variability and underlying environmental constraints. These results can be used to formulate new hypotheses regarding the mechanisms driving treeline formation and to guide the optimal choice of field sites at which to test these hypotheses.