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 Community Atmosphere–Biosphere Land Exchange model (CABLE) is a land surface model (LSM) that can be applied stand-alone and provides the land surface–atmosphere exchange within the Australian Community Climate and Earth System Simulator (ACCESS). We describe new developments that extend the applicability of CABLE for regional and global carbon–climate simulations, accounting for vegetation responses to biophysical and anthropogenic forcings. A land use and land cover change module driven by gross land use transitions and wood harvest area was implemented, tailored to the needs of the Coupled Model Intercomparison Project 6 (CMIP6). Novel aspects include the treatment of secondary woody vegetation, which benefits from a tight coupling between the land use module and the Population Orders Physiology (POP) module for woody demography and disturbance-mediated landscape heterogeneity. Land use transitions and harvest associated with secondary forest tiles modify the annually resolved patch age distribution within secondary vegetated tiles, in turn affecting biomass accumulation and turnover rates and hence the magnitude of the secondary forest sink. Additionally, we implemented a novel approach to constrain modelled GPP consistent with the coordination hypothesis and predicted by evolutionary theory, which suggests that electron-transport- and Rubisco-limited rates adjust seasonally and across biomes to be co-limiting. We show that the default prior assumption – common to CABLE and other LSMs – of a fixed ratio of electron transport to carboxylation capacity at standard temperature (J_max,0∕V_cmax,0) is at odds with this hypothesis; we implement an alternative algorithm for dynamic optimisation of this ratio such that coordination is achieved as an outcome of fitness maximisation. The results have significant implications for the magnitude of the simulated CO₂ fertilisation effect on photosynthesis in comparison to alternative estimates and observational proxies. These new developments enhance CABLE's capability for use within an Earth system model and in stand-alone applications to attribute trends and variability in the terrestrial carbon cycle to regions, processes and drivers. Model evaluation shows that the new model version satisfies several key observational constraints: (i) trend and interannual variations in the global land carbon sink, including sensitivities of interannual variations to global precipitation and temperature anomalies; (ii) centennial trends in global GPP; (iii) coordination of Rubisco-limited and electron-transport-limited photosynthesis; (iv) spatial distributions of global ET, GPP, biomass and soil carbon; and (v) age-dependent rates of biomass accumulation in boreal, temperate and tropical secondary forests. CABLE simulations agree with recent independent assessments of the global land–atmosphere flux partition that use a combination of atmospheric inversions and bottom-up constraints. In particular, there is agreement that the strong CO₂-driven sink in the tropics is largely cancelled by net deforestation and forest degradation emissions, leaving the Northern Hemisphere (NH) extratropics as the dominant contributor to the net land sink.

Several lines of evidence point to an increase in the activity of the terrestrial biosphere over recent decades, impacting the global net land carbon sink (NLS) and its control on the growth of atmospheric carbon dioxide (c_a). Global terrestrial gross primary production (GPP)—the rate of carbon fixation by photosynthesis—is estimated to have risen by (31 ± 5)% since 1900, but the relative contributions of different putative drivers to this increase are not well known. Here we identify the rising atmospheric CO₂ concentration as the dominant driver. We reconcile leaf‐level and global atmospheric constraints on trends in modeled biospheric activity to reveal a global CO₂ fertilization effect on photosynthesis of 30% since 1900, or 47% for a doubling of c_a above the pre‐industrial level. Our historic value is nearly twice as high as current estimates (17 ± 4)% that do not use the full range of available constraints. Consequently, under a future low‐emission scenario, we project a land carbon sink (174 PgC, 2006–2099) that is 57 PgC larger than if a lower CO₂ fertilization effect comparable with current estimates is assumed. These findings suggest a larger beneficial role of the land carbon sink in modulating future excess anthropogenic CO₂ consistent with the target of the Paris Agreement to stay below 2°C warming, and underscore the importance of preserving terrestrial carbon sinks.

The relative influences of Indian and Pacific Ocean modes of variability on Australian rainfall and soil moisture are investigated for seasonal, interannual, and decadal time scales. For the period 1900–2006, observations, reanalysis products, and hindcasts of soil moisture during the cool season (June–October) are used to assess the impacts of El Niño–Southern Oscillation (ENSO) and the Indian Ocean dipole (IOD) on southeastern Australia and the Murray–Darling Basin, two regions that have recently suffered severe droughts. A distinct asymmetry is found in the impacts of the opposite phases of both ENSO and IOD on Australian rainfall and soil moisture. There are significant differences between the dominant drivers of drought at interannual and decadal time scales. On interannual time scales, both ENSO and the IOD modify southeastern Australian soil moisture, with the driest (wettest) conditions over the southeast and more broadly over large parts of Australia occurring during years when an El Niño and a positive IOD event (La Niña and a negative IOD event) co-occur. The atmospheric circulation associated with these responses is discussed. Lower-frequency variability over southeastern Australia, however, including multiyear drought periods, seems to be more robustly related to Indian Ocean temperatures than Pacific conditions. The frequencies of both positive and negative IOD events are significantly different during periods of prolonged drought compared to extended periods of “normal” rainfall. In contrast, the frequency of ENSO events remains largely unchanged during prolonged dry and wet periods. For the Murray–Darling Basin, there appears to be a significant influence by La Niña and both positive and negative IOD events. In particular, La Niña plays a much more prominent role than for more southern regions, especially on interannual time scales and during prolonged wet periods. For prolonged dry (wet) periods, positive IOD events also occur in unusually high (low) numbers.