Impact of Land Surface Initialization Approach on Subseasonal Forecast Skill: a Regional Analysis in the Southern Hemisphere

The authors use a sophisticated coupled land-atmosphere modeling system for a Southern Hemisphere subdomain centered over southeastern Australia to evaluate differences in simulation skill from two different land surface initialization approaches. The first approach uses equilibrated land surface states obtained from offline simulations of the land surface model, and the second uses land surface states obtained from reanalyses. The authors find that land surface initialization using prior offline simulations contribute to relative gains in subseasonal forecast skill. In particular, relative gains in forecast skill for temperature of 10%-20% within the first 30 days of the forecast can be attributed to the land surface initialization method using offline states. For precipitation there is no distinct preference for the land surface initialization method, with limited gains in forecast skill irrespective of the lead time. The authors evaluated the asymmetry between maximum and minimum temperatures and found that maximum temperatures had the largest gains in relative forecast skill, exceeding 20% in some regions. These results were statistically significant at the 98% confidence level at up to 60 days into the forecast period. For minimum temperature, using reanalyses to initialize the land surface contributed to relative gains in forecast skill, reaching 40% in parts of the domain that were statistically significant at the 98% confidence level. The contrasting impact of the land surface initialization method between maximum and minimum temperature was associated with different soil moisture coupling mechanisms. Therefore, land surface initialization from prior offline simulations does improve predictability for temperature, particularly maximum temperature, but with less obvious improvements for precipitation and minimum temperature over southeastern Australia

A stand-alone tree demography and landscape structure module for Earth system models

We propose and demonstrate a new approach for the simulation of woody ecosystem stand dynamics, demography, and disturbance-mediated heterogeneity suitable for continental to global applications and designed for coupling to the terrestrial ecosystem component of any earth system model. The approach is encoded in a model called Populations-Order-Physiology (POP). We demonstrate the behavior and performance of POP coupled to the Community Atmosphere Biosphere Land Exchange model (CABLE) applied along the Northern Australian Tropical Transect, featuring gradients in rainfall and fire disturbance. The model is able to simultaneously reproduce observation-based estimates of key functional and structural variables along the transect, namely gross primary production, tree foliage projective cover, basal area, and maximum tree height. Prospects for the use of POP to address current vegetation dynamic deficiencies in earth system modeling are discussed

A vertically discretised canopy description for ORCHIDEE (SVN r2290) and the modifications to the energy, water and carbon fluxes

Since 70% of global forests are managed and forests impact the global carbon cycle and the energy exchange with the overlying atmosphere, forest management has the potential to mitigate climate change. Yet, none of the land surface models used in Earth system models, and therefore none of today’s predictions of future climate, account for the interactions between climate and forest management. We addressed this gap in modelling capability by developing and parametrizing a version of the land surface model ORCHIDEE to simulate the biogeochemical and biophysical effects of forest management. The most significant changes between the new branch called ORCHIDEE-CAN (SVN r2290) and the trunk version of ORCHIDEE (SVN r2243) are the allometric-based allocation of carbon to leaf, root, wood, fruit and reserve pools; the transmittance, absorbance and reflectance of radiation within the canopy; and the vertical discretisation of the energy budget calculations. In addition, conceptual changes were introduced towards a better process representation for the interaction of radiation with snow, the hydraulic architecture of plants, the representation of forest management and a numerical solution for the photosynthesis formalism of Farquhar, von Caemmerer and Berry. For consistency reasons, these changes were extensively linked throughout the code. Parametrization was revisited after introducing twelve new parameter sets that represent specific tree species or genera rather than a group of often distantly related or even unrelated species, as is the case in widely used plant functional types. Performance of the new model was compared against the trunk and validated against independent spatially explicit data for basal area, tree height, canopy strucure, GPP, albedo and evapotranspiration over Europe. For all tested variables ORCHIDEE-CAN outperformed the trunk regarding its ability to reproduce large-scale spatial patterns as well as their inter-annual variability over Europe. Depending on the data stream, ORCHIDEE-CAN had a 67% to 92% chance to reproduce the spatial and temporal variability of the validation data.JRC.H.5-Land Resources Managemen

Multi-decadal increase of forest burned area in Australia is linked to climate change

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

Sensitivity of isoprene emissions to drought over south-eastern Australia: Integrating models and satellite observations of soil moisture

South-east Australia, characterized by arid and semi-arid climate, has experienced large-scale rainfall reductions in recent decades. Larger temporal and spatial drought conditions are predicted in future. The temperate south east coastal zone is characterized by dense forests of Eucalyptus. Drought conditions have implications for the functioning of these indigenous ecosystems, and for emissions of reactive gases that are upwind of major metropolitan regions along the eastern coast. Here, we focus on the impact of drought on the emission of isoprene, a volatile organic compound emitted by a range of trees and shrubs. Previous model calculations grossly overestimate observations of the isoprene mixing ratio, potentially due to overestimated emission factors for native vegetation, but could also be due to drought-induced isoprene emission reductions. We develop the implementation of the Model of Emissions of Gases and Aerosols from Nature (MEGAN) within the CSIRO Chemical Transport Model to include a parameterization of drought using soil moisture. We test this parameterization using two approaches. First, we drive MEGAN using soil moisture fields from two Australian land surface models achieving reductions in isoprene emissions of 24–52% during summer. Second, we use a simple statistical approach to nudge model soil moisture towards satellite observations from the Soil Moisture and Ocean Salinity (SMOS) instrument. This work is the first application of SMOS towards isoprene emission modelling, providing a constraint on surface soil moisture which extends to global scales. Applying these soil moisture approaches, domain average isoprene emissions reduce by 38–58% in summer. Using these results, errors in basal emissions are likely in the region of 40%. Comparison of modelled soil moisture to observations using root density weighted averages across depths of 4 m showed minor differences of up to 0.04 m3 m−3. We find that the choice of land surface model used in the SMOS assimilations has a greater impact on isoprene emissions than adjusting either the nudging strength or the number of soil levels these satellite data may influence. However there are only small differences of 3% when using hourly or 24-hourly soil moisture input data to drive the emission calculations, suggesting that the lower temporal availability of satellite data does not reduce model quality. As the spatial resolution of satellite observations of atmospheric composition and land surface properties begin to approach the resolution of regulatory air quality models, we anticipate that these data will improve model predictive skills further

Evaluation of global terrestrial evapotranspiration using state-of-the-art approaches in remote sensing, machine learning and land surface modeling

Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm yr−1 (6.56×104 km3 yr−1) to 617.1 mm yr−1 (6.87×104 km3 yr−1). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm yr−2 with a significance level of p0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm yr−1. Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET

Understanding the uncertainty in global forest carbon turnover

The length of time that carbon remains in forest biomass is one of the largest uncertainties in the global carbon cycle, with both recent historical baselines and future responses to environmental change poorly constrained by available observations. In the absence of large-scale observations, models used for global assessments tend to fall back on simplified assumptions of the turnover rates of biomass and soil carbon pools. In this study, the biomass carbon turnover times calculated by an ensemble of contemporary terrestrial biosphere models (TBMs) are analysed to assess their current capability to accurately estimate biomass carbon turnover times in forests and how these times are anticipated to change in the future. Modelled baseline 1985-2014 global average forest biomass turnover times vary from 12.2 to 23.5 years between TBMs. TBM differences in phenological processes, which control allocation to, and turnover rate of, leaves and fine roots, are as important as tree mortality with regard to explaining the variation in total turnover among TBMs. The different governing mechanisms exhibited by each TBM result in a wide range of plausible turnover time projections for the end of the century. Based on these simulations, it is not possible to draw robust conclusions regarding likely future changes in turnover time, and thus biomass change, for different regions. Both spatial and temporal uncertainty in turnover time are strongly linked to model assumptions concerning plant functional type distributions and their controls. Thirteen model-based hypotheses of controls on turnover time are identified, along with recommendations for pragmatic steps to test them using existing and novel observations. Efforts to resolve uncertainty in turnover time, and thus its impacts on the future evolution of biomass carbon stocks across the world\u27s forests, will need to address both mortality and establishment components of forest demography, as well as allocation of carbon to woody versus non-woody biomass growth