Australia's Tinderbox Drought: An extreme natural event likely worsened by human-caused climate change.

We examine the characteristics and causes of southeast Australia's Tinderbox Drought (2017 to 2019) that preceded the Black Summer fire disaster. The Tinderbox Drought was characterized by cool season rainfall deficits of around -50% in three consecutive years, which was exceptionally unlikely in the context of natural variability alone. The precipitation deficits were initiated and sustained by an anomalous atmospheric circulation that diverted oceanic moisture away from the region, despite traditional indicators of drought risk in southeast Australia generally being in neutral states. Moisture deficits were intensified by unusually high temperatures, high vapor pressure deficits, and sustained reductions in terrestrial water availability. Anthropogenic forcing intensified the rainfall deficits of the Tinderbox Drought by around 18% with an interquartile range of 34.9 to -13.3% highlighting the considerable uncertainty in attributing droughts of this kind to human activity. Skillful predictability of this drought was possible by incorporating multiple remote and local predictors through machine learning, providing prospects for improving forecasting of droughts.

Predicting resilience through the lens of competing adjustments to vegetation function.

There is a pressing need to better understand ecosystem resilience to droughts and heatwaves. Eco-evolutionary optimization approaches have been proposed as means to build this understanding in land surface models and improve their predictive capability, but competing approaches are yet to be tested together. Here, we coupled approaches that optimize canopy gas exchange and leaf nitrogen investment, respectively, extending both approaches to account for hydraulic impairment. We assessed model predictions using observations from a native Eucalyptus woodland that experienced repeated droughts and heatwaves between 2013 and 2020, whilst exposed to an elevated [CO2 ] treatment. Our combined approaches improved predictions of transpiration and enhanced the simulated magnitude of the CO2 fertilization effect on gross primary productivity. The competing approaches also worked consistently along axes of change in soil moisture, leaf area, and [CO2 ]. Despite predictions of a significant percentage loss of hydraulic conductivity due to embolism (PLC) in 2013, 2014, 2016, and 2017 (99th percentile PLC > 45%), simulated hydraulic legacy effects were small and short-lived (2 months). Our analysis suggests that leaf shedding and/or suppressed foliage growth formed a strategy to mitigate drought risk. Accounting for foliage responses to water availability has the potential to improve model predictions of ecosystem resilience.

Annual precipitation explains variability in dryland vegetation greenness globally but not locally.

Dryland vegetation productivity is strongly modulated by water availability. As precipitation patterns and variability are altered by climate change, there is a pressing need to better understand vegetation responses to precipitation variability in these ecologically fragile regions. Here we present a global analysis of dryland sensitivity to annual precipitation variations using long-term records of normalized difference vegetation index (NDVI). We show that while precipitation explains 66% of spatial gradients in NDVI across dryland regions, precipitation only accounts for <26% of temporal NDVI variability over most (>75%) dryland regions. We observed this weaker temporal relative to spatial relationship between NDVI and precipitation across all global drylands. We confirmed this result using three alternative water availability metrics that account for water loss to evaporation, and growing season and precipitation timing. This suggests that predicting vegetation responses to future rainfall using space-for-time substitution will strongly overestimate precipitation control on interannual variability in aboveground growth. We explore multiple mechanisms to explain the discrepancy between spatial and temporal responses and find contributions from multiple factors including local-scale vegetation characteristics, climate and soil properties. Earth system models (ESMs) from the latest Coupled Model Intercomparison Project overestimate the observed vegetation sensitivity to precipitation variability up to threefold, particularly during dry years. Given projections of increasing meteorological drought, ESMs are likely to overestimate the impacts of future drought on dryland vegetation with observations suggesting that dryland vegetation is more resistant to annual precipitation variations than ESMs project.

Global hotspots for the occurrence of compound events.

Compound events (CEs) are weather and climate events that result from multiple hazards or drivers with the potential to cause severe socio-economic impacts. Compared with isolated hazards, the multiple hazards/drivers associated with CEs can lead to higher economic losses and death tolls. Here, we provide the first analysis of multiple multivariate CEs potentially causing high-impact floods, droughts, and fires. Using observations and reanalysis data during 1980-2014, we analyse 27 hazard pairs and provide the first spatial estimates of their occurrences on the global scale. We identify hotspots of multivariate CEs including many socio-economically important regions such as North America, Russia and western Europe. We analyse the relative importance of different multivariate CEs in six continental regions to highlight CEs posing the highest risk. Our results provide initial guidance to assess the regional risk of CE events and an observationally-based dataset to aid evaluation of climate models for simulating multivariate CEs.

Plant profit maximization improves predictions of European forest responses to drought.

Knowledge of how water stress impacts the carbon and water cycles is a key uncertainty in terrestrial biosphere models. We tested a new profit maximization model, where photosynthetic uptake of CO2 is optimally traded against plant hydraulic function, as an alternative to the empirical functions commonly used in models to regulate gas exchange during periods of water stress. We conducted a multi-site evaluation of this model at the ecosystem scale, before and during major droughts in Europe. Additionally, we asked whether the maximum hydraulic conductance in the soil-plant continuum kmax (a key model parameter which is not commonly measured) could be predicted from long-term site climate. Compared with a control model with an empirical soil moisture function, the profit maximization model improved the simulation of evapotranspiration during the growing season, reducing the normalized mean square error by c. 63%, across mesic and xeric sites. We also showed that kmax could be estimated from long-term climate, with improvements in the simulation of evapotranspiration at eight out of the 10 forest sites during drought. Although the generalization of this approach is contingent upon determining kmax , it presents a mechanistic trait-based alternative to regulate canopy gas exchange in global models.

Impact of the representation of stomatal conductance on model projections of heatwave intensity.

Stomatal conductance links plant water use and carbon uptake, and is a critical process for the land surface component of climate models. However, stomatal conductance schemes commonly assume that all vegetation with the same photosynthetic pathway use identical plant water use strategies whereas observations indicate otherwise. Here, we implement a new stomatal scheme derived from optimal stomatal theory and constrained by a recent global synthesis of stomatal conductance measurements from 314 species, across 56 field sites. Using this new stomatal scheme, within a global climate model, subtantially increases the intensity of future heatwaves across Northern Eurasia. This indicates that our climate model has previously been under-predicting heatwave intensity. Our results have widespread implications for other climate models, many of which do not account for differences in stomatal water-use across different plant functional types, and hence, are also likely under projecting heatwave intensity in the future.

Allowable CO2 emissions based on regional and impact-related climate targets.

Global temperature targets, such as the widely accepted limit of an increase above pre-industrial temperatures of two degrees Celsius, may fail to communicate the urgency of reducing carbon dioxide (CO2) emissions. The translation of CO2 emissions into regional- and impact-related climate targets could be more powerful because such targets are more directly aligned with individual national interests. We illustrate this approach using regional changes in extreme temperatures and precipitation. These scale robustly with global temperature across scenarios, and thus with cumulative CO2 emissions. This is particularly relevant for changes in regional extreme temperatures on land, which are much greater than changes in the associated global mean.

The plumbing of land surface models: is poor performance a result of methodology or data quality?

The 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, surface air temperature and relative humidity. These results are explored here in greater detail and possible causes are investigated. We examine 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. We demonstrate that energy conservation in the observational data is not responsible for these results. We also show 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, we present evidence suggesting that the nature of this partitioning problem is likely shared among all contributing LSMs. While we do not find a single candidate explanation for why land surface models perform poorly relative to empirical benchmarks in PLUMBER, we do exclude multiple possible explanations and provide guidance on where future research should focus.

Factors that shape seed mass evolution.

We used correlated divergence analysis to determine which factors have been most closely associated with changes in seed mass during seed plant evolution. We found that divergences in seed mass have been more consistently associated with divergences in growth form than with divergences in any other variable. This finding is consistent with the strong relationship between seed mass and growth form across present-day species and with the available data from the paleobotanical literature. Divergences in seed mass have also been associated with divergences in latitude, net primary productivity, temperature, precipitation, and leaf area index. However, these environmental variables had much less explanatory power than did plant traits such as seed dispersal syndrome and plant growth form.