Explore the vast range of titles available.

Fire activity in Australia is strongly affected by high inter-annual climate variability and extremes. Through changes in the climate, anthropogenic climate change has the potential to alter fire dynamics. Here we compile satellite (19 and 32 years) and ground-based (90 years) burned area datasets, climate and weather observations, and simulated fuel loads for Australian forests. Burned area in Australia’s forests shows a linear positive annual trend but an exponential increase during autumn and winter. The mean number of years since the last fire has decreased consecutively in each of the past four decades, while the frequency of forest megafire years (>1 Mha burned) has markedly increased since 2000. The increase in forest burned area is consistent with increasingly more dangerous fire weather conditions, increased risk factors associated with pyroconvection, including fire-generated thunderstorms, and increased ignitions from dry lightning, all associated to varying degrees with anthropogenic climate change. The degree to which wildfire activity in Australia is affected by climate change is not well quantified. Here, the authors show that the frequency of forest fires and the area burned have increased significantly over recent decades, mainly due to an increase in dangerous fire weather conditions through warmer temperature and circulation changes.

Although the existence of a large carbon sink in terrestrial ecosystems is well-established, the drivers of this sink remain uncertain. It has been suggested that perturbations to forest demography caused by past land-use change, management, and natural disturbances may be causing a large component of current carbon uptake. Here we use a global compilation of forest age observations, combined with a terrestrial biosphere model with explicit modeling of forest regrowth, to partition the global forest carbon sink between old-growth and regrowth stands over the period 1981–2010. For 2001–2010 we find a carbon sink of 0.85 (0.66–0.96) Pg year−1 located in intact old-growth forest, primarily in the moist tropics and boreal Siberia, and 1.30 (1.03–1.96) Pg year−1 located in stands regrowing after past disturbance. Approaching half of the sink in regrowth stands would have occurred from demographic changes alone, in the absence of other environmental changes. These age-constrained results show consistency with those simulated using an ensemble of demographically-enabled terrestrial biosphere models following an independent reconstruction of historical land use and management. We estimate that forests will accumulate an additional 69 (44–131) Pg C in live biomass from changes in demography alone if natural disturbances, wood harvest, and reforestation continue at rates comparable to those during 1981–2010. Our results confirm that it is not possible to understand the current global terrestrial carbon sink without accounting for the sizeable sink due to forest demography. They also imply that a large portion of the current terrestrial carbon sink is strictly transient in nature.

Summary Atmospheric carbon dioxide concentration ([CO2]) is increasing, which increases leaf‐scale photosynthesis and intrinsic water‐use efficiency. These direct responses have the potential to increase plant growth, vegetation biomass, and soil organic matter; transferring carbon from the atmosphere into terrestrial ecosystems (a carbon sink). A substantial global terrestrial carbon sink would slow the rate of [CO2] increase and thus climate change. However, ecosystem CO2 responses are complex or confounded by concurrent changes in multiple agents of global change and evidence for a [CO2]‐driven terrestrial carbon sink can appear contradictory. Here we synthesize theory and broad, multidisciplinary evidence for the effects of increasing [CO2] (iCO2) on the global terrestrial carbon sink. Evidence suggests a substantial increase in global photosynthesis since pre‐industrial times. Established theory, supported by experiments, indicates that iCO2 is likely responsible for about half of the increase. Global carbon budgeting, atmospheric data, and forest inventories indicate a historical carbon sink, and these apparent iCO2 responses are high in comparison to experiments and predictions from theory. Plant mortality and soil carbon iCO2 responses are highly uncertain. In conclusion, a range of evidence supports a positive terrestrial carbon sink in response to iCO2, albeit with uncertain magnitude and strong suggestion of a role for additional agents of global change. See also the Commentary on this article by Way et al., 229: 2383–2385.

The model Soil‐Litter‐Iso (SLI) calculates coupled heat and water transport in soil. It was recently implemented into the Australian land surface model CABLE, which is the land component of the Australian Community Climate and Earth System Simulator (ACCESS). Here we extended SLI to include accurate freeze‐thaw processes in the soil and snow. SLI provides thence an implicit solution of the energy and water balances of soil and snow as a standalone model and within CABLE. The enhanced SLI was tested extensively against theoretical formulations, laboratory experiments, field data, and satellite retrievals. The model performed well for all experiments at wide‐ranging temporal and spatial scales. SLI melts snow faster at the end of the cold season compared to observations though because there is no subgrid variability within SLI given by the implicit, coupled solution of energy and water. Combined CABLE‐SLI shows very realistic dynamics and extent of permafrost on the Northern hemisphere. It illustrated, however, also the limits of possible comparisons between large‐scale land surface models and local permafrost observations. CABLE‐SLI exhibits the same patterns of snow depth and snow water equivalent on the Northern hemisphere compared to satellite‐derived observations but quantitative comparisons depend largely on the given meteorological input fields. Further extension of CABLE‐SLI with depth‐dependence of soil carbon will allow realistic projections of the development of permafrost and frozen carbon stocks in a changing climate. Key Points Freeze‐thaw and snow implemented into soil model Soil‐Litter‐Iso: stand‐alone and part of the land surface scheme CABLE Excels against theoretical formulations, laboratory experiments, field data, and satellite observations Very realistic dynamics and extent of permafrost as well as of Northern hemispheric snow

Twenty-seven models participated in the Earth System Model–Snow Model Intercomparison Project (ESM-SnowMIP), the most data-rich MIP dedicated to snow modeling. Our findings do not support the hypothesis advanced by previous snow MIPs: evaluating models against more variables and providing evaluation datasets extended temporally and spatially does not facilitate identification of key new processes requiring improvement to model snow mass and energy budgets, even at point scales. In fact, the same modeling issues identified by previous snow MIPs arose: albedo is a major source of uncertainty, surface exchange parameterizations are problematic, and individual model performance is inconsistent. This lack of progress is attributed partly to the large number of human errors that led to anomalous model behavior and to numerous resubmissions. It is unclear how widespread such errors are in our field and others; dedicated time and resources will be needed to tackle this issue to prevent highly sophisticated models and their research outputs from being vulnerable because of avoidable human mistakes. The design of and the data available to successive snow MIPs were also questioned. Evaluation of models against bulk snow properties was found to be sufficient for some but inappropriate for more complex snow models whose skills at simulating internal snow properties remained untested. Discussions between the authors of this paper on the purpose of MIPs revealed varied, and sometimes contradictory, motivations behind their participation. These findings started a collaborative effort to adapt future snow MIPs to respond to the diverse needs of the community.

OzFlux is the regional Australian and New Zealand flux tower network that aims to provide a continental-scale national research facility to monitor and assess trends, and improve predictions, of Australia's terrestrial biosphere and climate. This paper describes the evolution, design, and current status of OzFlux as well as provides an overview of data processing. We analyse measurements from all sites within the Australian portion of the OzFlux network and two sites from New Zealand. The response of the Australian biomes to climate was largely consistent with global studies except that Australian systems had a lower ecosystem water-use efficiency. Australian semi-arid/arid ecosystems are important because of their huge extent (70 %) and they have evolved with common moisture limitations. We also found that Australian ecosystems had a similar radiation-use efficiency per unit leaf area compared to global values that indicates a convergence toward a similar biochemical efficiency. The two New Zealand sites represented extremes in productivity for a moist temperate climate zone, with the grazed dairy farm site having the highest GPP of any OzFlux site (2620 gC m−2 yr−1) and the natural raised peat bog site having a very low GPP (820 gC m−2 yr−1). The paper discusses the utility of the flux data and the synergies between flux, remote sensing, and modelling. Lastly, the paper looks ahead at the future direction of the network and concludes that there has been a substantial contribution by OzFlux, and considerable opportunities remain to further advance our understanding of ecosystem response to disturbances, including drought, fire, land-use and land-cover change, land management, and climate change, which are relevant both nationally and internationally. It is suggested that a synergistic approach is required to address all of the spatial, ecological, human, and cultural challenges of managing the delicately balanced ecosystems in Australasia.

New evidence is emerging that semi-arid ecosystems dominate interannual variability (IAV) of the global carbon cycle, largely via fluctuating water availability associated with El Niño Southern Oscillation. Recent evidence from global terrestrial biosphere modelling and satellite-based inversion of atmospheric CO2 point to a large role of Australian ecosystems in global carbon cycle variability, including a large contribution from Australia to the record land sink of 2011. However the specific mechanisms governing this variability, and their bioclimatic distribution within Australia, have not been identified. Here we provide a regional assessment, based on best available observational data, of IAV in the Australian terrestrial carbon cycle and the role of Australia in the record land sink anomaly of 2011. We find that IAV in Australian net carbon uptake is dominated by semi-arid ecosystems in the east of the continent, whereas the 2011 anomaly was more uniformly spread across most of the continent. Further, and in contrast to global modelling results suggesting that IAV in Australian net carbon uptake is amplified by lags between production and decomposition, we find that, at continental scale, annual variations in production are dampened by annual variations in decomposition, with both fluxes responding positively to precipitation anomalies.

Light use efficiency (LUE) defines the vegetation efficiency of converting radiative energy into biochemical energy through photosynthesis. Estimating the maximum LUE (ε max) is critical yet challenging for quantifying gross primary production (GPP) using LUE-based models. This study describes an analytical method for estimating ε max based on water use efficiency (WUE) as determined by plant water use and carbon gain. Unlike other complex parameterization schemes, this WUE-based method is simple and requires four variables relatively easy to acquire. The WUE-based ε max estimates compare favorably well with values based on traditional curve fitting method and that reported in the literature, and clearly distinguished ε max between C3 (1.48 ± 0.33 \\({\\text{g C M}}{{\\text{J}}^{ - 1}}{\\text{ PAR}}\\)) and C4 (2.63 ± 0.21 \\({\\text{g C M}}{{\\text{J}}^{ - 1}}{\\text{ PAR}}\\)) dominated ecosystems. The range in ε max estimates was narrow across different years and sites within a biome. The WUE-based ε max estimate is theoretically constrained by vegetation water use and can be directly incorporated into LUE models for GPP estimation across ecosystems.

The Protocol for the Analysis of Land Surface Models (PALS) Land Surface Model Benchmarking Evaluation Project (PLUMBER) illustrated the value of prescribing a priori performance targets in model intercomparisons. It showed that the performance of turbulent energy flux predictions from different land surface models, at a broad range of flux tower sites using common evaluation metrics, was on average worse than relatively simple empirical models. For sensible heat fluxes, all land surface models were outperformed by a linear regression against downward shortwave radiation. For latent heat flux, all land surface models were outperformed by a regression against downward shortwave radiation, surface air temperature, and relative humidity. These results are explored here in greater detail and possible causes are investigated. It is examined whether particular metrics or sites unduly influence the collated results, whether results change according to time-scale aggregation, and whether a lack of energy conservation in flux tower data gives the empirical models an unfair advantage in the intercomparison. It is demonstrated that energy conservation in the observational data is not responsible for these results. It is also shown that the partitioning between sensible and latent heat fluxes in LSMs, rather than the calculation of available energy, is the cause of the original findings. Finally, evidence is presented that suggests that the nature of this partitioning problem is likely shared among all contributing LSMs. While a single candidate explanation for why land surface models perform poorly relative to empirical benchmarks in PLUMBER could not be found, multiple possible explanations are excluded and guidance is provided on where future research should focus.

Year-to-year variability in CO 2 fluxes can yield insight into climate-carbon cycle relationships, a fundamental yet uncertain aspect of the terrestrial carbon cycle. In this study, we use global observations from NASA’s Orbiting Carbon Observatory-2 (OCO-2) satellite for years 2015–2019 and a geostatistical inverse model to evaluate 5 years of interannual variability (IAV) in CO 2 fluxes and its relationships with environmental drivers. OCO-2 launched in late 2014, and we specifically evaluate IAV during the time period when OCO-2 observations are available. We then compare inferences from OCO-2 with state-of-the-art process-based models (terrestrial biosphere model, TBMs). Results from OCO-2 suggest that the tropical grasslands biome (including grasslands, savanna, and agricultural lands within the tropics) makes contributions to global IAV during the 5 year study period that are comparable to tropical forests, a result that differs from a majority of TBMs. Furthermore, existing studies disagree on the environmental variables that drive IAV during this time period, and the analysis using OCO-2 suggests that both temperature and precipitation make comparable contributions. TBMs, by contrast, tend to estimate larger IAV during this time and usually estimate larger relative contributions from the extra-tropics. With that said, TBMs show little consensus on both the magnitude and the contributions of different regions to IAV. We further find that TBMs show a wide range of responses on the relationships of CO 2 fluxes with annual anomalies in temperature and precipitation, and these relationships across most of the TBMs have a larger magnitude than inferred from OCO-2. Overall, the findings of this study highlight large uncertainties in process-based estimates of IAV during recent years and provide an avenue for evaluating these processes against inferences from OCO-2.