The conversion of marginal pastoral land in New Zealand to higher biomass shrubland consisting of manuka (Leptospermum scoparium) and kanuka (Kunzea ericoides var. ericoides) offers opportunity for carbon (C) sequestration, with potential co-benefits of soil erosion control. We therefore selected two areas with different soils in different climatic regions to investigate and compare soil respiration rates, methane and nitrous oxide emission profiles, and key carbon exchange processes controlling carbon sequestration. In addition, two shrubland stands of different ages were selected in each area, providing four sites in total. Regular (almost monthly) soil respiration measurements were made over a 2-year period, with less frequent methane and nitrous oxide flux measurements, and soil sampling once at the end of the study. The cooler, wetter volcanic soils had higher total organic C (6.39 ± 0.12% v. 5.51 ± 0.17%), soil C : nitrogen (N) ratios (20.55 ± 0.20 v. 18.45 ± 0.23), and slightly lower mineral N (3.30 ± 0.74 v. 4.89 ± 0.57 mg/kg) and microbial biomass C (1131 ± 108 v. 1502 ± 37 mg/kg) than the more drought-prone, stony, sedimentary soils. Mineral-N contents at all sites indicated N-limited ecosystems for allocation of below- and above-ground C.
The estimated mean annual cumulative respiration rate recorded in the volcanic soil was 10.26 ± 7.45 t CO2-C/ha.year compared with 9.85 ± 8.63 t CO2-C/ha.year in the stony sedimentary soil for the 2 years of our study. Older shrubland stands had higher respiration rates than younger stands in both study areas. Methane oxidation was estimated to be higher in the volcanic soil (4.10 ± 2.13 kg CH4-C/ha.year) than the sedimentary soil sites (2.51 ± 2.48 kg CH4-C/ha.year). The measured natural background levels of nitrous oxide emissions from these shrubland soils ranged between negligible and 0.30 ± 0.20 kg N2O-N/ha.year. A strong climatic control (temperature and moisture) on gas fluxes was observed at all sites. Our sampling strategy at each of the four sites was to estimate the mean soil respiration rates (n = 25) from an 8 by 8 m sampling grid positioned into a representative location. Soil respiration rates were also measured (by additional, less frequent sampling) in two adjacent grids (1-m offset and 100-m distant grid) to test the validity of these representative mean values. The 1-m offset grid (n = 25) provided a statistically different soil respiration rate from the main grid (n = 25) in 25% of the 12 sampling events. The 100-m grid (n = 25) provided a statistically different respiration rate to the main grid in 38% of the 26 sampling events. These differences are attributed to the spatially variable and sporadic nature of gaseous emissions from soils. The grid analysis tested the prediction uncertainty and it provides evidence for strong spatial and temporal control by edaphic processes in micro-sites. A partial least-squares regression model was used to relate the 2009 annual cumulative soil respiration to site-specific edaphic characteristics, i.e. biomass, nutrient availability, porosity and bulk density, measured at the end of that year. The model explained ≥80% of the variance at three of the four sites.