Saltwater Forest is a Dacrydium cupressinum-dominated lowland forest covering 9000 ha in south Westland, South Island, New Zealand. Four thousand hectares is managed for sustainable production of indigenous timber. The aim of this study was to provide an integrated analysis of soils, soil-landform relationships, and soil-vegetation relationships at broad and detailed scales. The broad scale understandings provide a framework in which existing or future studies can be placed and the detailed studies elucidate sources of soil and forest variability. Glacial landforms dominate. They include late Pleistocene lateral, terminal and ablation moraines, and outwash aggradation and degradation terraces. Deposits and landforms from six glacial advances have been recognised ranging from latest Last (Otira) Glaciation to Penultimate (Waimea) Glaciation. The absolute ages of landforms were established by analysis of the thickness and soil stratigraphy of loess coverbeds, augmented with radiocarbon dating and phytolith and pollen analysis. In the prevailing high rainfall of Westland soil formation is rapid. The rate of loess accretion in Saltwater Forest (ca. 30 mm ka⁻¹) has been low enough that soil formation and loess accretion took place contemporaneously. Soils formed in this manner are known as upbuilding soils. The significant difference between upbuilding pedogenesis and pedogenesis in a topdown sense into an existing sediment body is that each subsoil increment of an upbuilding soil has experienced processes of all horizons above. In Saltwater Forest subsoils of upbuilding soils are strongly altered because they have experienced the extremely acid environment of the soil surface at some earlier time. Some soil chronosequence studies in Westland have included upbuilding soils formed in loess as the older members of the sequence. Rates and types of processes inferred from these soils should be reviewed because upbuilding is a different pedogenic pathway to topdown pedogenesis. Landform age and morphology were used as a primary stratification for a study of the soil pattern and nature of soil variability in the 4000 ha production area of Saltwater Forest. The age of landforms (> 14 ka) and rapid soil formation mean that soils are uniformly strongly weathered and leached. Soils include Humic Organic Soils, Perch-gley Podzols, Acid Gley Soils, Allophanic Brown Soils, and Orthic or Pan Podzols. The major influence on the nature of soils is site hydrology which is determined by macroscale features of landforms (slope, relief, drainage density), mesoscale effects related to position on landforms, and microscale influences determined by microtopography and individual tree effects. Much of the soil variability arises at microscales so that it is not possible to map areas of uniform soils at practical map scales. The distribution of soil variability across spatial scales, in relation to the intensity of forest management, dictates that it is most appropriate to map soil complexes with boundaries coinciding with landforms. Disturbance of canopy trees is an important agent in forest dynamics. The frequency of forest disturbance in the production area of Saltwater Forest varies in a systematic way among landforms in accord with changes in abundance of different soils. The frequency of forest turnover is highest on landforms with the greatest abundance of extremely poorly-drained Organic Soils. As the abundance of better-drained soils increases the frequency of forest turnover declines. Changes in turnover frequency are reflected in the mean size and density of canopy trees (Dacrydium cupressinum) among landforms. Terrace and ablation moraine landforms with the greatest abundance of extremely poorly-drained soils have on average the smallest trees growing most densely. The steep lateral moraines, characterised by well drained soils, have fewer, larger trees. The changes manifested at the landform scale are an integration of processes operating over much shorter range as a result of short-range soil variability. The systematic changes in forest structure and turnover frequency among landforms and soils have important implications for sustainable forest management.

Saltwater Forest is a Dacrydium cupressinum-dominated lowland forest covering 9000 ha in south Westland, South Island, New Zealand. Four thousand hectares is managed for sustainable production of indigenous timber. The aim of this study was to provide an integrated analysis of soils, soil-landform relationships, and soil-vegetation relationships at broad and detailed scales. The broad scale understandings provide a framework in which existing or future studies can be placed and the detailed studies elucidate sources of soil and forest variability. Glacial landforms dominate. They include late Pleistocene lateral, terminal and ablation moraines, and outwash aggradation and degradation terraces. Deposits and landforms from six glacial advances have been recognised ranging from latest Last (Otira) Glaciation to Penultimate (Waimea) Glaciation. The absolute ages of landforms were established by analysis of the thickness and soil stratigraphy of loess coverbeds, augmented with radiocarbon dating and phytolith and pollen analysis. In the prevailing high rainfall of Westland soil formation is rapid. The rate of loess accretion in Saltwater Forest (ca. 30 mm ka⁻¹) has been low enough that soil formation and loess accretion took place contemporaneously. Soils formed in this manner are known as upbuilding soils. The significant difference between upbuilding pedogenesis and pedogenesis in a topdown sense into an existing sediment body is that each subsoil increment of an upbuilding soil has experienced processes of all horizons above. In Saltwater Forest subsoils of upbuilding soils are strongly altered because they have experienced the extremely acid environment of the soil surface at some earlier time. Some soil chronosequence studies in Westland have included upbuilding soils formed in loess as the older members of the sequence. Rates and types of processes inferred from these soils should be reviewed because upbuilding is a different pedogenic pathway to topdown pedogenesis. Landform age and morphology were used as a primary stratification for a study of the soil pattern and nature of soil variability in the 4000 ha production area of Saltwater Forest. The age of landforms (> 14 ka) and rapid soil formation mean that soils are uniformly strongly weathered and leached. Soils include Humic Organic Soils, Perch-gley Podzols, Acid Gley Soils, Allophanic Brown Soils, and Orthic or Pan Podzols. The major influence on the nature of soils is site hydrology which is determined by macroscale features of landforms (slope, relief, drainage density), mesoscale effects related to position on landforms, and microscale influences determined by microtopography and individual tree effects. Much of the soil variability arises at microscales so that it is not possible to map areas of uniform soils at practical map scales. The distribution of soil variability across spatial scales, in relation to the intensity of forest management, dictates that it is most appropriate to map soil complexes with boundaries coinciding with landforms. Disturbance of canopy trees is an important agent in forest dynamics. The frequency of forest disturbance in the production area of Saltwater Forest varies in a systematic way among landforms in accord with changes in abundance of different soils. The frequency of forest turnover is highest on landforms with the greatest abundance of extremely poorly-drained Organic Soils. As the abundance of better-drained soils increases the frequency of forest turnover declines. Changes in turnover frequency are reflected in the mean size and density of canopy trees (Dacrydium cupressinum) among landforms. Terrace and ablation moraine landforms with the greatest abundance of extremely poorly-drained soils have on average the smallest trees growing most densely. The steep lateral moraines, characterised by well drained soils, have fewer, larger trees. The changes manifested at the landform scale are an integration of processes operating over much shorter range as a result of short-range soil variability. The systematic changes in forest structure and turnover frequency among landforms and soils have important implications for sustainable forest management.

Worldwide there is considerable concern over the effects of human activities on the quantity and quality of freshwater. Measurement of infiltration behaviour will be important for improving freshwater management. This study identifies that New Zealand has a sporadic history of measuring soil water movement attributes on a limited number of soil types, although the current practical demand should be large for management of irrigation, dairy farm effluent disposal, as well as municipal / domestic waste- and storm-water disposal. Previous research has demonstrated that infiltration behaviour is governed by the interplay between numerous mechanisms including hydrophobicity and preferential flow, the latter being an important mechanism of contaminant leaching for many NZ soils. Future characterisation will need to recognise the dynamic nature of these interactions, and be able to reliably characterise the key infiltration mechanisms. Since macropores are responsible for preferential flow, it is critical that infiltration studies use a representative sample of the macropore network. The aim of this project was to study the mechanisms governing the infiltration behaviour of a layered soil in large (50 x 70 cm) monolith lysimeters, where the connectivity of the macropore network remains undisturbed. Four lysimeters of the Gorge silt loam were collected, a structured soil with four distinct layers. On each lysimeter there were four separate infiltration experiments, with water applied under suctions of 0, 0.5, 1, and 1.5 kPa by a custom-built tension infiltrometer. Each lysimeter was instrumented with 30 tensiometers, located in arrays at the layer boundaries. There was also a field experiment using ponded dye infiltration to visually define preferential flowpaths. Analysis of dye patterns, temporal variability in soil matric potential (Ψm), and solute breakthrough curves all show that preferential flow is an important infiltration mechanism. Preferential flowpaths were activated when Ψm was above -1.5 kPa. During saturated infiltration, at least 97% of drainage was through the ‘mobile’ pore volume of the lysimeter (θm), estimated among the lysimeters at 5.4 – 8.7 % of the lysimeter volume. Early-time infiltration behaviour did not show the classical square-root of time behaviour, indicating sorptivity was not the governing mechanism. This was consistent across the four lysimeters, and during infiltration under different surface imposed suctions. The most likely mechanism restricting sorptivity is weak hydrophobicity, which appears to restrict infiltration for the first 5 – 10 mm of infiltration. Overall, the Gorge soil’s early-time infiltration behaviour is governed by the dynamic interaction between sorptivity, hydrophobicity, the network of air-filled pores, preferential flow and air encapsulation. Long-time infiltration behaviour was intimately linked to the temporal dynamics of Ψm, which was in turn controlled by preferential flow and soil layer interactions. Preferential flowpaths created strong inter-layer connectivity by allowing an irregular wetting front to reach lower layers within 2 – 15 mm of infiltration. Thereafter, layer interactions dominate infiltration for long-time periods, as Ψm in soil layers with different K(Ψm) relationships self-adjusts to try to maintain a constant Darcy velocity. An important finding was that Ψm rarely attained the value set by the tension infiltrometer during unsaturated infiltration. The results show that ‘true’ steady-state infiltration is unlikely to occur in layered soils. A quasi-steady state was identified once the whole column had fully wet and layer interactions had settled to where Ψm changes occurred in unison through each soil layer. Quasi-steady state was difficult to identify from just the cumulative infiltration curve, but more robustly identified as when infiltration matched drainage, and Ψm measurements showed each layer had a stable hydraulic gradient. I conclude that the in-situ hydraulic conductivity, K(Ψm), of individual soil layers can be accurately and meaningfully determined from lysimeter-scale infiltration experiments. My results show that K(Ψm) is different for each soil layer, and that differences are consistent among the four lysimeters. Under saturated flow the subsoil had the lowest conductivity, and was the restricting layer. Most interestingly this pattern reversed during unsaturated flow. As Ψm decreased below -0.5 to -1 kPa, the subsoil was markedly more conductive, and the topsoil layers became the restricting layers. All four soil layers demonstrate a sharp decline in K(Ψm) as Ψm decreases, with a break in slope at ~ -1 kPa indicating the dual-permeability nature of all layers.