Bridge to the future: Important lessons from 20 years of ecosystem observations made by the OzFlux network.

In 2020, the Australian and New Zealand flux research and monitoring network, OzFlux, celebrated its 20th anniversary by reflecting on the lessons learned through two decades of ecosystem studies on global change biology. OzFlux is a network not only for ecosystem researchers, but also for those 'next users' of the knowledge, information and data that such networks provide. Here, we focus on eight lessons across topics of climate change and variability, disturbance and resilience, drought and heat stress and synergies with remote sensing and modelling. In distilling the key lessons learned, we also identify where further research is needed to fill knowledge gaps and improve the utility and relevance of the outputs from OzFlux. Extreme climate variability across Australia and New Zealand (droughts and flooding rains) provides a natural laboratory for a global understanding of ecosystems in this time of accelerating climate change. As evidence of worsening global fire risk emerges, the natural ability of these ecosystems to recover from disturbances, such as fire and cyclones, provides lessons on adaptation and resilience to disturbance. Drought and heatwaves are common occurrences across large parts of the region and can tip an ecosystem's carbon budget from a net CO2 sink to a net CO2 source. Despite such responses to stress, ecosystems at OzFlux sites show their resilience to climate variability by rapidly pivoting back to a strong carbon sink upon the return of favourable conditions. Located in under-represented areas, OzFlux data have the potential for reducing uncertainties in global remote sensing products, and these data provide several opportunities to develop new theories and improve our ecosystem models. The accumulated impacts of these lessons over the last 20 years highlights the value of long-term flux observations for natural and managed systems. A future vision for OzFlux includes ongoing and newly developed synergies with ecophysiologists, ecologists, geologists, remote sensors and modellers.

Could patterns of animal behaviour cause the observed differences in soil carbon between adjacent irrigated and unirrigated pastures?

Previous soil sampling from grazed pastures in New Zealand compared the changes of soil organic carbon (SOC) in adjacent irrigated and unirrigated portions of the same paddocks. It showed that irrigated portions had lower SOC stocks than unirrigated portions, with an average difference of 7.0 tC ha-1 or 0.6 tC ha-1 yr-1. These findings have formed the basis of an assessment for the net effect of conversion of New Zealand's grazed pastures to irrigation. However, since cattle could move freely between irrigated and unirrigated portions of the studied paddocks, there could have been different grazing intensities and/or excreta transfer between the irrigated and unirrigated portions of the same paddocks. Both these factors could have affected SOC stocks. In this study, we used the process-based model, CenW, to simulate the consequences of this possible carbon transfer via animal excreta and different grazing intensities. We found that the observed increase of 0.6 tC ha-1 yr-1 in SOC stock in the unirrigated portions could result from a transfer of 20% excreta from the irrigated to unirrigated portions (with an area ratio of 6:1) of a paddock and with the unirrigated portions being grazed only lightly with 2.0 tDM ha-1 in foliage biomass residuals remaining after grazing. That means that the observed higher SOC stocks in the unirrigated portions could potentially be attributable to the behaviour of grazing animals. We suggest that a realistic extent of carbon transfer and/or differences in grazing intensities could be sufficient to account for the observed differences in SOC stocks even if irrigation per se caused no differences in carbon stocks. It is therefore inappropriate to ascribe the change of SOC to irrigation effects based on experimental findings where SOC changes can be affected by the behaviour of grazing animals.

Irrigation increases forage production of newly established lucerne but enhances net ecosystem carbon losses.

In New Zealand, dairy farming faces increasing scrutiny for its environmental impacts, including those on soil carbon (C) stocks; hence, alternative management practices are required. One such practice is usage of deep-rooting forage, such as lucerne (Medicago sativa L.). We measured the C and water exchange of two neighbouring lucerne fields on stony, well-drained soil for 3 years, following conversion from grassland. One field received irrigation and effluent; the other received neither. Net CO2 exchange and evaporation were measured by eddy covariance, drainage and leaching with lysimeters, and water inputs with rain gauges. Biomass removal from harvesting and grazing was recorded by direct sampling. In the conversion year, irrigated lucerne was C-neutral despite two harvests and losses from the conversion process. In the 2nd and 3rd years combined, the biomass-C removal exceeded net CO2 uptake, causing net losses of 450 g C m-2 and 210 g C m-2 for irrigated and non-irrigated lucerne, respectively. Leaching losses accounted for 1 to 9 % of annual net C uptake from the atmosphere. The ratio of ecosystem respiration to gross photosynthetic productivity (GPP) increased from <0.7 in spring to ≈ 1 in autumn. Consequently, the net C balance for both lucerne crops showed gains in the first two growth periods of each year and losses in the subsequent two to four growth periods. Irrigation made no difference to the photosynthetic water-use efficiency at field scale (GPP/evaporation), but enhanced production water-use efficiency (biomass/water input). Irrigation increased both the absolute amount of drainage and the fraction of water inputs lost by drainage. In one year, significant summer drainage occurred for the irrigated lucerne. To prevent that, soil-water content should be kept well below field capacity but above the crop's water-stress level. Such practice would likely also help retain soil carbon.

Nitrous Oxide Fluxes, Soil Oxygen, and Denitrification Potential of Urine- and Non-Urine-Treated Soil under Different Irrigation Frequencies.

Despite increased use of irrigation to improve forage quality and quantity for grazing cattle ( Linnaeus), there is a lack of data that assess how irrigation practices influence nitrous oxide (NO) emissions from urine-affected soils. Irrigation effects on soil oxygen (O) availability, a primary controller of NO fluxes, is poorly understood. It was hypothesized that increased irrigation frequency would result in lower NO emissions by increasing soil moisture and decreasing soil O concentrations. This would favor more NO reduction to dinitrogen (N). We examined effects of high (3-d) versus low (6-d) irrigation frequency with and without bovine urine addition to pasture. Nitrous oxide fluxes were measured daily for 35 d. Soil O, temperature, and water content were continuously measured at multiple depths. Inorganic nitrogen, organic carbon, and soil pH were measured at 6-d intervals. Measurements of denitrification enzyme activity with and without acetylene inhibition were used to infer the NO/(NO + N) ratio. The NO/(NO + N) ratio was lower under high- compared with low-frequency irrigation, suggesting greater potential for NO reduction to N with more frequent irrigation. Although NO fluxes were increased by urine addition, they were not affected by irrigation frequency. Soil O decreased temporarily after urine deposition, but O dynamics did not explain NO dynamics. Relative soil gas diffusivity (/) was a better predictor of NO fluxes than O concentration. On a free-draining soil, increasing irrigation frequency while providing the same total water volume did not enhance NO emissions under ruminant urine patches in a grazed pasture.

Footprint methods to separate N2O emission rates from adjacent paddock areas.

Using micrometeorological techniques to measure greenhouse gas emissions from differently treated adjacent plots is a promising avenue to verify the effect of mitigation strategies at the field scale. In pursuing such an approach, it is crucial to accurately characterize the source area of the fluxes measured at each sampling point. Hence, a comprehensive footprint analysis method is required so that emission rates can be obtained for a specific field within a biochemically heterogeneous area. In this study, a footprint analysis method is developed to estimate the emission for an experiment where the flux of N2O is measured from several control and treated plots. The emission rate of an individual plot is estimated using an inverse footprint fraction approach where the footprint fractions are obtained from an analytical footprint model. A numerical solution for obtaining the background flux for such a multiplot measurement system is also provided. Results of the footprint analysis method are assessed, first, by comparing footprint fractions obtained from both an analytical footprint model and a "forward" simulation of a backward Lagrangian stochastic (bLs) model; and second, by comparing the emission rates of a control plot obtained from the footprint analysis method and from the "backward" simulation of the bLs model. It is found that the analytical footprint fractions compare well with the values obtained from the bLs model (correlation coefficient of 0.58 and 0.66 within p value <0.001). An average of 4.3 % of the measured fluxes is found to be contributed by sources outside the measured area and, excluding this outside area contribution to the measured flux, footprint corrected emission rates within the defined domain are found to increase by 2.1 to 5.8 % of the measured flux. Also, the proposed method of emission rate estimation is found to work well under a wide range of atmospheric stability.

Rapid changes in δ¹³C of ecosystem-respired CO2 after sunset are consistent with transient ¹³C enrichment of leaf respired CO2.

The CO2 respired by darkened, light-adapted, leaves is enriched in ¹³C during the first minutes, and this effect may be related to rapid changes in leaf respiratory biochemistry upon darkening. We hypothesized that this effect would be evident at the ecosystem scale. High temporal resolution measurements of the carbon isotope composition of ecosystem respiration were made over 28 diel periods in an abandoned temperate pasture, and were compared with leaf-level measurements at differing levels of pre-illumination. At the leaf level, CO2 respired by darkened leaves that had been preadapted to high light was strongly enriched in ¹³C, but such a ¹³C-enrichment rapidly declined over 60-100 min. The ¹³C-enrichment was less pronounced when leaves were preadapted to low light. These leaf-level responses were mirrored at the ecosystem scale; after sunset following clear, sunny days respired CO2 was first ¹³C enriched, but the ¹³C-enrichment rapidly declined over 60-100 min. Further, this response was less pronounced following cloudy days. We conclude that the dynamics of leaf respiratory isotopic signal caused variations in ecosystem-scale ¹²CO2/¹³) CO2 exchange. Such rapid isotope kinetics should be considered when applying ¹³C-based techniques to elucidate ecosystem carbon cycling.