Search

Mapping in situ eddy covariance measurements of terrestrial land–atmosphere fluxes to the globe is a key method for diagnosing the Earth system from a data-driven perspective. We describe the first global products (called X-BASE) from a newly implemented upscaling framework, FLUXCOM-X, representing an advancement from the previous generation of FLUXCOM products in terms of flexibility and technical capabilities. The X-BASE products are comprised of estimates of CO2 net ecosystem exchange (NEE), gross primary productivity (GPP), evapotranspiration (ET), and for the first time a novel, fully data-driven global transpiration product (ETT), at high spatial (0.05°) and temporal (hourly) resolution. X-BASE estimates the global NEE at −5.75 ± 0.33 Pg C yr−1 for the period 2001–2020, showing a much higher consistency with independent atmospheric carbon cycle constraints compared to the previous versions of FLUXCOM. The improvement of global NEE was likely only possible thanks to the international effort to increase the precision and consistency of eddy covariance collection and processing pipelines, as well as to the extension of the measurements to more site years resulting in a wider coverage of bioclimatic conditions. However, X-BASE global net ecosystem exchange shows a very low interannual variability, which is common to state-of-the-art data-driven flux products and remains a scientific challenge. With 125 ± 2.1 Pg C yr−1 for the same period, X-BASE GPP is slightly higher than previous FLUXCOM estimates, mostly in temperate and boreal areas. X-BASE evapotranspiration amounts to 74.7×103 ± 0.9×103 km3 globally for the years 2001–2020 but exceeds precipitation in many dry areas, likely indicating overestimation in these regions. On average 57 % of evapotranspiration is estimated to be transpiration, in good agreement with isotope-based approaches, but higher than estimates from many land surface models. Despite considerable improvements to the previous upscaling products, many further opportunities for development exist. Pathways of exploration include methodological choices in the selection and processing of eddy covariance and satellite observations, their ingestion into the framework, and the configuration of machine learning methods. For this, the new FLUXCOM-X framework was specifically designed to have the necessary flexibility to experiment, diagnose, and converge to more accurate global flux estimates.

Reliable and robust monitoring tools are crucial to assess the effectiveness of land mitigation techniques (LMTs) in enhancing carbon uptake, enabling informed decision making by policymakers. This study, addressing one of the scientific goals of the EU Horizon 2020 Land Use Based Mitigation for Resilient Climate Pathways (LANDMARC) project, examines the feasibility of using satellite solar-induced fluorescence (SIF) in combination with other satellite data as a monitoring proxy to evaluate the effects of LMTs on carbon uptake. Two distinct cases are explored: (1) instantaneous vegetation destruction caused by a 2019 eucalyptus wildfire in southeast Australia and (2) gradual forest gain resulting from reforestation efforts in northern China in 2007–2012. The cases are monitored using SIF from the TROPOspheric Monitoring Instrument (TROPOMI) and Global Ozone Monitoring Experiment-2A (GOME-2A), respectively. Comparing the temporal variability in SIF across the affected areas and nearby reference areas reveals that vegetation dynamics changed as a consequence of the land-use changes in both cases. Specifically, in the Australia case, TROPOMI demonstrated an immediate reduction in the SIF signal of 0.6 mW m−2 sr−1 nm−1 (−72 %) over the eucalypt forest right after the fire. Exploiting the strong correspondence between TROPOMI SIF and gross primary productivity (GPP) at the nearby representative eddy covariance Tumbarumba site and through the FluxSat product, we estimate that the Australian fire led to a loss in GPP of 130–200 GgC in the first 8 months after the fire. Over the northern Chinese provinces of Gansu, Shaanxi, Sichuan, Chongqing, and Shanxi, we report an increase in GOME-2A summertime SIF of 0.1–0.2 mW m−2 sr−1 nm−1, coinciding with reforestation efforts between 2007 and 2012. This increase in the SIF signal is likely driven by a combination of increasingly favourable natural conditions and by the reforestation effort itself. A multivariate model that takes into account growth factors such as water availability and maximum temperature as well as satellite-derived forest-cover data explains the observed variability in GOME-2A SIF in the Chinese case reasonably well (R2=0.72). The model suggests that increases in both forest cover and soil moisture have led, in step, to the observed increase in vegetation activity over northern China. In that region, for every 100 km2 of additional forest cover, SIF increases by 0.1 mW m−2 sr−1 nm−1 between 2007 and 2012. Our study highlights that the use of satellite-based SIF, together with supporting in situ, modelled, and satellite data, allows us to monitor the impact of LMT implementation on regional carbon uptake as long as the scale of the LMT is of sufficient spatial extent.

One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Aerosol- and cloud-induced changes in diffuse light have important impacts on the global land carbon cycle, as they alter light distribution and photosynthesis in vegetation canopies. However, this effect remains poorly represented or evaluated in current land surface models. Here, we add a light partitioning module and a new canopy light transmission module to the ORCHIDEE (Organising Carbon and Hydrology In Dynamic Ecosystems) land surface model (trunk version, v5453) and use the revised model, ORCHIDEE_DF, to estimate the fraction of diffuse light and its effect on gross primary production (GPP) in a multilayer canopy. We evaluate the new parameterizations using flux observations from 159 eddy covariance sites over the globe. Our results show that, compared with the original model, ORCHIDEE_DF improves the GPP simulation under sunny conditions and captures the observed higher photosynthesis under cloudier conditions in most plant functional types (PFTs). Our results also indicate that the larger GPP under cloudy conditions compared with sunny conditions is mainly driven by increased diffuse light in the morning and in the afternoon as well as by a decreased vapor pressure deficit (VPD) and decreased air temperature at midday. The observations show that the strongest positive effects of diffuse light on photosynthesis are found in the range from 5 to 20 ∘C and at a VPD 2 uptake. ORCHIDEE_DF underestimates the diffuse light effect at low temperature in all PFTs and overestimates this effect at high temperature and at a high VPD in grasslands and croplands. The new model has the potential to better investigate the impact of large-scale aerosol changes and long-term changes in cloudiness on the terrestrial carbon budget, both in the historical period and in the context of future air quality policies and/or climate engineering.

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 (Jmax, 0∕Vcmax, 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 CO2 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 CO2-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.

This paper presents a methodology for the attribution and characterisation of Sclerophyll forested landscapes over large areas. First we define a set of woody vegetation data primitives (e.g. canopy cover, leaf area index (LAI), bole density, canopy height), which are then scaled-up using multiple remote sensing data sources to characterise and extract landscape woody vegetation features. The advantage of this approach is that vegetation landscape features can be described from composites of these data primitives. The proposed data primitives act as building blocks for the re-creation of past woody characterisation schemes as well as allowing for re-compilation to support present and future policy and management and decision making needs. Three main research sites were attributed; representative of different sclerophyll woody vegetated systems (Box Iron-bark forest; Mountain Ash forest; Mixed Species foothills forest). High resolution hyperspectral and full waveform LiDAR data was acquired over the three research sites. At the same time, land management agencies (Victorian Department of Environment, Land Water and Planning) and researchers (RMIT, CRC for Spatial Information and CSIRO) conducted fieldwork to collect structural and functional measurements of vegetation, using traditional forest mensuration transects and plots, terrestrial lidar scanning and high temporal resolution in-situ autonomous laser (VegNet) scanners. Results are presented of: 1) inter-comparisons of LAI estimations made using ground based hemispherical photography, LAI 2200 PCA, CI-110 and terrestrial and airborne laser scanners; 2) canopy height and vertical canopy complexity derived from airborne LiDAR validated using ground observations; and, 3) time-series characterisation of land cover features. 1. Accuracy targets for remotely sensed LAI products to match within ground based estimates are ± 0.5 LAI or a 20% maximum (CEOS/GCOS) with new aspirational targets of 5%). In this research we conducted a total of 67 ground-based method-to-method pairwise comparisons across 11 plots in five sites, incorporating the previously mentioned LAI methods. Out of the 67 comparisons, 29 had an RMSE ≥ 0.5 LAIe. This has important implications for the validation of remotely sensed products since ground based techniques themselves exhibit LAI variations greater than internationally recommended guidelines for satellite product accuracies. 2. Two methods of canopy height derivation are proposed and tested over a large area (4 Million Ha). 99th percentile maximum height achieved a RMSE of 6.6%, whilst 95th percentile dominant height a RMSE = 10.3%. Vertical canopy complexity (i.e. the number of forest layers of strata) was calculated as the local maxima of vegetation density within the LiDAR canopy profile and determined using a cubic spline smoothing of Pgap. This was then validated against in-situ and LiDAR observations of canopy strata with an RMSE 0.39 canopy layers. 3. Preliminary results are presented of landcover characterisation using LandTrendr analysis of Landsat LEDAPS data. kNN is then used to link these features to a dense network of 800 field plots sites.

No abstract available.

This study introduces a prototype end-to-end Simulator software tool for simulating two-dimensional satellite multispectral imagery for a variety of satellite instrument models in aquatic environments. Using case studies, the impact of variable sensor configurations on the performance of value-added products for challenging applications, such as coral reefs and cyanobacterial algal blooms, is assessed. This demonstrates how decisions regarding satellite sensor design, driven by cost constraints, directly influence the quality of value-added remote sensing products. Furthermore, the Simulator is used to identify situations where retrieval algorithms require further parameterization before application to unsimulated satellite data, where error sources cannot always be identified or isolated. The application of the Simulator can verify whether a given instrument design meets the performance requirements of end-users before build and launch, critically allowing for the justification of the cost and specifications for planned and future sensors. It is hoped that the Simulator will enable engineers and scientists to understand important design trade-offs in phase 0/A studies easily, quickly, reliably, and accurately in future Earth observation satellites and systems.

In 2020, the Australian and New Zealand flux research and monitoring network, OzFlux, celebrated its 20 th 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 CO 2 sink to a net CO 2 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., (© 2022 The Authors. Global Change Biology published by John Wiley & Sons Ltd.)

Gross primary productivity (GPP) of wooded ecosystems (forests and savannas) is central to the global carbon cycle, comprising 67%-75% of total global terrestrial GPP. Climate change may alter this flux by increasing the frequency of temperatures beyond the thermal optimum of GPP (T opt ). We examined the relationship between GPP and air temperature (Ta) in 17 wooded ecosystems dominated by a single plant functional type (broadleaf evergreen trees) occurring over a broad climatic gradient encompassing five ecoregions across Australia ranging from tropical in the north to Mediterranean and temperate in the south. We applied a novel boundary-line analysis to eddy covariance flux observations to (a) derive ecosystem GPP-Ta relationships and T opt (including seasonal analyses for five tropical savannas); (b) quantitatively and qualitatively assess GPP-Ta relationships within and among ecoregions; (c) examine the relationship between T opt and mean daytime air temperature (MDTa) across all ecosystems; and (d) examine how down-welling short-wave radiation (Fsd) and vapour pressure deficit (VPD) influence the GPP-Ta relationship. GPP-Ta relationships were convex parabolas with narrow curves in tropical forests, tropical savannas (wet season), and temperate forests, and wider curves in temperate woodlands, Mediterranean woodlands, and tropical savannas (dry season). Ecosystem T opt ranged from 15℃ (temperate forest) to 32℃ (tropical savanna-wet and dry seasons). The shape of GPP-Ta curves was largely determined by daytime Ta range, MDTa, and maximum GPP with the upslope influenced by Fsd and the downslope influenced by VPD. Across all ecosystems, there was a strong positive linear relationship between T opt and MDTa (Adjusted R 2 : 0.81; Slope: 1.08) with T opt exceeding MDTa by >1℃ at all but two sites. We conclude that ecosystem GPP has adjusted to local MDTa within Australian broadleaf evergreen forests and that GPP is buffered against small Ta increases in the majority of these ecosystems., (© 2021 John Wiley & Sons Ltd.)

Ecosystem respiration is a major component of the global terrestrial carbon cycle and is strongly influenced by temperature. The global extent of the temperature-ecosystem respiration relationship, however, has not been fully explored. Here, we test linear and threshold models of ecosystem respiration across 210 globally distributed eddy covariance sites over an extensive temperature range. We find thresholds to the global temperature-ecosystem respiration relationship at high and low air temperatures and mid soil temperatures, which represent transitions in the temperature dependence and sensitivity of ecosystem respiration. Annual ecosystem respiration rates show a markedly reduced temperature dependence and sensitivity compared to half-hourly rates, and a single mid-temperature threshold for both air and soil temperature. Our study indicates a distinction in the influence of environmental factors, including temperature, on ecosystem respiration between latitudinal and climate gradients at short (half-hourly) and long (annual) timescales. Such climatological differences in the temperature sensitivity of ecosystem respiration have important consequences for the terrestrial net carbon sink under ongoing climate change.

Background: Earth Observation 'EO' remote sensing technology development enables original insights into vegetation function and health at ever finer temporal, spectral and spatial resolution. Research sites equipped with monitoring infrastructure such as flux towers operate at a key bridging scale between satellite platform measurements and on-the-ground leaf-level processes., Results: This paper presents the technical details of the design and operation of a proximal observation system 'THEMS' that generates unattended long-term high quality thermal and hyperspectral images of a forest canopy on a short (sub-daily) timescale. The primary purpose of the system is to measure canopy temperature, spectral reflectance and radiance coincident with a highly instrumented flux tower site for benchmarking purposes. Basic system capability is demonstrated through low level data product descriptions of the high-resolution multi-angular imagery and ancillary data streams. The system has been successfully operational for more than 2 years with little to no intervention., Conclusions: These data can then be used to derive remotely sensed proxies of canopy and ecosystem function to study temporal forest dynamics over a wide range of wavelengths, spatial scales (individual trees to canopy), and temporal scales (minutes to multiple years). The multi-purpose system is intended to provide unprecedented spatio-temporal ecophysiological insight and to underpin upscaling of remotely sensed dynamic ecosystem water, CO 2 , and energy exchange processes., Competing Interests: Competing interestsThe authors declare that they have no competing interests., (© The Author(s) 2020.)

The FLUXNET2015 dataset provides ecosystem-scale data on CO 2 , water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

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 2 concentration as the dominant driver. We reconcile leaf-level and global atmospheric constraints on trends in modeled biospheric activity to reveal a global CO 2 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 2 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 2 consistent with the target of the Paris Agreement to stay below 2°C warming, and underscore the importance of preserving terrestrial carbon sinks., (© 2020 The Authors. Global Change Biology published by John Wiley & Sons Ltd.)