Flux Tower Data Research Articles

Are the ecosystem-level evaporative stress indices representative of evaporative stress of vegetation?

Evaporative Stress Index (ESI), also sometimes referred as Evaporative Stress Ratio (ESR), has been widely used as an indicator of vegetation evaporative stress, and is often used to track forest and agriculture droughts. Lower the stress, higher is the value of ESI or ESR. The goal of this study is to assess the suitability of these indices for tracking vegetation evaporative stress. As the dynamics of water loss from vegetation through transpiration (T) can be different than that of evapotranspiration (ET) from the ecosystem, it is hypothesized that ESI or ESR may not be sufficiently representative of the vegetation evaporative stress. Using eddy covariance flux tower data of 518 site years, distributed across 49-sites and 9 land covers globally, our findings reveal underestimation of vegetation evaporative stress by ESI during periods of high vapor pressure deficit (VPD) and overestimation during dry, low-VPD periods. The results highlight the need to improve representativeness of ESI for monitoring vegetation evaporative stress. Notably, this may entail accurate estimation of ecosystem T in systems lacking in-situ data, a challenge that warrants further attention.

Better representation of vegetation phenology improves estimations of annual gross primary productivity

Carbon uptake by vegetation plays a vital role in the global carbon cycle. Annual gross primary productivity (AGPP) represents the total amount of carbon compounds produced by vegetation photosynthesis over a year and is a crucial metric for quantifying vegetation carbon uptake. Although some theories have been developed to explain the spatiotemporal variation in AGPP, they often overlook the seasonal differences in photosynthesis across vegetation phenological periods. This gap highlights the need for a more reasonable representation of vegetation phenology in AGPP modelling. Therefore, we developed a novel theoretical AGPP model that decomposes AGPP into detailed vegetation phenological periods (green-up, maturation, and senescence) and investigated the effects of climatic factors on the seasonal components of AGPP. Compared with existing models, our model considers the length of the multiple vegetation phenological periods, rather than just the length of the carbon uptake period. When evaluated against flux tower data, our model outperformed the comparative model, with a higher determination coefficient and lower root mean square error. The analysis also demonstrated that the developed model can reproduce spatial and temporal variations in satellite-based global AGPP. In addition, we identified distinct differences in the responses of AGPP's seasonal components to climatic factors: shortwave radiation predominantly affected the seasonal component of AGPP during senescence, air temperature during green-up, and vapor pressure deficit during maturation. Our study proposes a potential mechanism for AGPP estimation and highlights the importance of accurately representing vegetation phenology in AGPP modelling.

Predicting time series of vegetation leaf area index across North America based on climate variables for land surface modeling using attention-enhanced LSTM

ABSTRACT Ecosystem process modeling is vital for a broad range of applications. It requires a key metric to characterize vegetation canopies: the leaf area index (LAI). Many land surface models utilize satellite-derived LAI for modeling ecosystem processes. The use of satellite-derived LAI constrains the models for predicting ecosystem processes when observational data are unavailable. Enriching prognostic models for predicting the LAI is favorable to land surface studies for forecasting vegetation-related processes. An attention-enhanced long and short memory (AELSTM) model was developed to predict the vegetation LAI time series based on climatic data. The AELSTM outperformed comparative machine learning models when assessed using satellite and flux tower data. The LAI predicted by AELSTM are temporally and spatially consistent with satellite-derived LAI. The R2 values attained 0.86–0.93 across biomes. AELSTM can be coupled with the land surface models to predict gross primary productivity. Our proposed model was demonstrated to be effective in predicting the LAI time series under different Shared Socioeconomic Pathways in CMIP6. Our study demonstrated that deep learning approaches are capable of modeling and characterizing the spatial and temporal patterns of key vegetation variables such as LAI. Moreover, the study has provided references for vegetation processes in land surface studies.

Causal Drivers of Land‐Atmosphere Carbon Fluxes From Machine Learning Models and Data

AbstractInteractions among atmospheric, root‐soil, and vegetation processes drive carbon dioxide fluxes (Fc) from land to atmosphere. Eddy covariance measurements are commonly used to measure Fc at sub‐daily timescales and validate process‐based and data‐driven models. However, these validations do not reveal process interactions, thresholds, and key differences in how models replicate them. We use information theory‐based measures to explore multivariate information flow pathways from forcing data to observed and modeled hourly Fc, using flux tower data sets in the Midwestern U.S. in intensively managed corn‐soybean landscapes. We compare multiple linear regressions, long‐short term memory, and random forests (RF), and evaluate how different model structures use information from combinations of sources to predict Fc. We extend a framework for model predictive and functional performance, which examines a suite of dependencies from all forcing variables to the observed or modeled target. Of the three model types, RF exhibited the highest functional and predictive performance, correctly capturing strong dependencies between radiation and temperature variables with Fc. Regionally trained models demonstrate lower predictive but higher functional performance compared to site‐specific models, suggesting superior reproduction of observed relationships at the expense of predictive accuracy. This study shows that some metrics of predictive performance encapsulate functional behaviors better than others, highlighting the need for multiple metrics of both types. This study improves our understanding of carbon fluxes in an intensively managed landscape, and more generally provides insight into how model structures and forcing variables translate to interactions that are well versus poorly captured in models.

Technical note: Flagging inconsistencies in flux tower data

Abstract. Global collections of synthesized flux tower data such as FLUXNET have accelerated scientific progress beyond the eddy covariance community. However, remaining data issues in FLUXNET data pose challenges for users, particularly for multi-site synthesis and modelling activities. Here, we present complementary consistency flags (C2Fs) for flux tower data, which rely on multiple indications of inconsistency among variables, along with a methodology to detect discontinuities in time series. The C2F relates to carbon and energy fluxes, as well as to core meteorological variables, and consists of the following: (1) flags for daily data values, (2) flags for entire-site variables, and (3) flags at time stamps that mark large discontinuities in the time series. The flagging is primarily based on combining outlier scores from a set of predefined relationships among variables. The methodology to detect break points in the time series is based on a non-parametric test for the difference in distributions of model residuals. Applying C2F to the FLUXNET 2015 dataset reveals the following: (1) among the considered variables, gross primary productivity and ecosystem respiration data were flagged most frequently, in particular during rain pulses under dry and hot conditions. This information is useful for modelling and analysing ecohydrological responses. (2) There are elevated flagging frequencies for radiation variables (shortwave, photosynthetically active, and net). This information can improve the interpretation and modelling of ecosystem fluxes with respect to issues in the driver. (3) The majority of long-term sites show temporal discontinuities in the time series of latent energy, net ecosystem exchange, and radiation variables. This should be useful for carefully assessing the results in terms of interannual variations in and trends of ecosystem fluxes. The C2F methodology is flexible for customizing and allows for varying the desired strictness of consistency. We discuss the limitations of the approach that can present starting points for future developments.

Exploring the Potential of Long Short‐Term Memory Networks for Predicting Net CO2 Exchange Across Various Ecosystems With Multi‐Source Data

AbstractUpscaling flux tower measurements based on machine learning (ML) algorithms is an essential approach for large‐scale net ecosystem CO2 exchange (NEE) estimation, but existing ML upscaling methods face some challenges, particularly in capturing NEE interannual variations (IAVs) that may relate to lagged effects. With the capacity to characterize temporal memory effects, the Long Short‐Term Memory (LSTM) networks are expected to help solve this problem. Here we explored the potential of LSTM for predicting NEE across various ecosystems using flux tower data over 82 sites in North America. The LSTM model with differentiated plant function types (PFTs) demonstrates the capability to explain 79.19% (R2 = 0.79) of the monthly variations in NEE within the testing set, with RMSE and Mean Absolute Error values of 0.89 and 0.57 g C m−2 d−1 respectively (r = 0.89, p < 0.001). Moreover, the LSTM model performed robustly in predicting cross‐site variability, with 67.19% of the sites that can be predicted by both LSTM models with and without distinguished PFTs showing improved predictive ability. Most importantly, the IAV of predicted NEE highly correlated with that in flux observations (r = 0.81, p < 0.001), clearly outperforming that by the random forest model (r = −0.21, p = 0.011). Among all nine PFTs, solar‐induced chlorophyll fluorescence, downward shortwave radiation, and leaf area index are the most important variables for explaining NEE variations, collectively accounting for approximately 54.01% in total. This study highlights the great potential of LSTM for improving carbon flux upscaling with multi‐source remote sensing data.

AgriCarbon-EO v1.0.1: large-scale and high-resolution simulation of carbon fluxes by assimilation of Sentinel-2 and Landsat-8 reflectances using a Bayesian approach

Abstract. Soil organic carbon storage is a well-identified climate change mitigation solution. Quantification of the soil carbon storage in cropland for agricultural policy and offset carbon markets using in situ sampling would be excessively costly, especially at the intrafield scale. For this reason, comprehensive monitoring, reporting, and verification (MRV) of soil carbon and its explanatory variables at a large scale need to rely on hybrid approaches that combine remote sensing and modelling tools to provide the carbon budget components with their associated uncertainties at intrafield scale. Here, we present AgriCarbon-EO v1.0.1: an end-to-end processing chain that enables the estimation of carbon budget components for major and cover crops at intrafield resolution (10 m) and regional extents (e.g. 10 000 km2) by assimilating remote sensing data (e.g. Sentinel-2 and Landsat8) in a physically based radiative transfer (PROSAIL) and agronomic models (SAFYE-CO2). The data assimilation in AgriCarbon-EO is based on a novel Bayesian approach that combines normalized importance sampling and look-up table generation. This approach propagates the uncertainties across the processing chain from the reflectances to the output variables. After a presentation of the chain, we demonstrate the accuracy of the estimates of AgriCarbon-EO through an application over winter wheat in the southwest of France during the cropping seasons from 2017 to 2019. We validate the outputs with flux tower data for net ecosystem exchange, biomass destructive samples, and combined harvester yield maps. Our results show that the scalability and uncertainty estimates proposed by the approach do not hinder the accuracy of the estimates (net ecosystem exchange, NEE: RMSE =1.68–2.38 gC m−2, R2=0.87–0.77; biomass: RMSE =11.34 g m−2, R2=0.94). We also show the added value of intrafield simulations for the carbon components through scenario testing of pixel and field simulations (biomass: bias =-47 g m−2, −39 % variability). Our overall analysis shows satisfying accuracy, but it also points out the need to represent more soil processes and include synthetic aperture radar data that would enable a larger coverage of AgriCarbon-EO. The paper's findings confirm the suitability of the choices made in building AgriCarbon-EO as a hybrid solution for an MRV scheme to diagnose agro-ecosystem carbon fluxes.

Simulation of gross primary productivity and impact of drought in Liulin watershed of Taihang mountains over 2000–2020

As the largest flux in the global terrestrial carbon cycle, Gross Primary Productivity (GPP) is a key factor for the terrestrial carbon budget. Based on the eddy flux tower data, remote sensing data and the meteorological precipitation data, we calibrated three parameters of maximum light use efficiency, optimum temperature, and maximum temperature in the Vegetation Photosynthesis Model (VPM) to reproduce the GPP spatiotemporal characteristic in Liulin watershed of Taihang Mountains. We further analyzed the time lag effect and cumulative effect of the standard precipitation index (SPI) on GPP by the Pearson correlation coefficients. Results show that VPM captures the GPP characteristics of eddy flux tower observation in Liulin watershed. During the calibration period (R2 = 0.94, RMSE = 0.47 g C m-2·d-1) and validation period (R2 = 0.91, RMSE = 0.60 g C m-2 d-1), the adjusted VPM model successfully reproduced the temporal dynamics of site-observed GPP. From 2000 to 2020, VPM model can effectively reproduce the long-term variations in GPP (R2 = 0.69 vs. NIRv_GPP), with the GPP ranging from 446.4 to 828.1 g C m-2 y-1 and an average value of 545.6 g C m-2 y-1 in Liulin watershed. The GPP exhibited an increasing trend from 2000 to 2020. A mutation in the GPP occurred in 2010, with an average annual increase of 13.9 g C m-2 y-1 from 2000 to 2010 (p < 0.05), followed by an average annual increase of 10.22 g C m-2 y-1 from 2011 to 2020 (p > 0.05). The SPI effectively monitored drought events in the Liulin watershed. With a cumulative effect of 9 months and a lag effect of 6–7 months, the time lag and cumulative effects of drought had a long-term effect on GPP. It is imperative to have long-term eddy flux tower data and meteorological data to better evaluate the spatiotemporal dynamics of GPP and the impact of drought.

Hydrometeorological Factors Affecting the Carbon Exchange of the Himalayan Pine-dominated Ecosystem

The net carbon uptake of the natural and dense Himalayan ecosystems is not well explored at a sub-daily scale due to limited data availability. Here, we present the first-ever analysis of the hydrometeorological controls of the sub-daily scale Net Ecosystem Exchange (NEE) in a Western Himalayan pine-dominated ecosystem (Pinus roxburghii). We used observed half-hourly flux tower data for 104 weeks from December 2014 to November 2016. We found that the sub-daily (6 h) scale variability of NEE contributes 20–40% to its within-week variability. We generated transfer entropy (TE)-based weekly process networks with a maximum memory of 6 h to understand sub-daily scale processes. Using weekly networks, we present consistent causal links on a seasonal scale. The hydrometeorological variables affecting the NEE and Ecosystem Respiration (RE) at a sub-daily scale are Net Solar Radiation (NSR), Relative Humidity (RH), Surface Air Temperature (TA), and Sensible Heat Flux (SH). These variables are consistently active in the networks with seasonally varying magnitudes of TE. We also found that for both monsoon and post-monsoon, NEE receives more incoming links from meteorological variables during dry weeks than wet weeks. Precipitation (P) did not send a direct causal link to NEE or RE within a sub-daily scale but strongly influenced the causal associations between hydrometeorological and carbon exchange variables. The network does not show any immediate (within 6 h) impact of P on NEE, which probably implies a delayed response of vegetation to P through soil moisture. At higher memory, P may act as an influential background variable.

Improving the ability of solar-induced chlorophyll fluorescence to track gross primary production through differentiating sunlit and shaded leaves

Recently, solar-induced chlorophyll fluorescence (SIF) is a promising tool to estimate gross primary production (GPP). Photosynthesis gradually saturates with the increasing light, but fluorescence tends to keep increasing, leading to a nonlinear SIF-GPP relationship. This nonlinearity occurs for sunlit leaves but not for shaded leaves for which photosynthesis is light-limited. However, the separation of sunlit and shaded SIF has not been systematically investigated when estimating GPP from SIF. Therefore, it is promising to develop a model for GPP estimation considering such differences. This study proposed an approach to separate the total canopy SIF emission (SIFtotal) from TROPOspheric Monitoring Instrument (TROPOMI) SIF into their sunlit and shaded components (SIFsun and SIFshade). The nonlinearity and linearity in SIF-GPP relationships for sunlit and shaded leaves were incorporated into a two-leaf hybrid model, which was fitted using flux tower data and then evaluated using leave-one-site-out crossing validation. We also elucidated the distinct SIF-GPP relationships between sunlit and shaded leaves using the Soil-Canopy-Observation of Photosynthesis and the Energy balance (SCOPE) model simulation. Compared to previously used linear (R2 = 0.68, RMSE = 2.13 gC·m−2·d−1) or hyperbolic (R2 = 0.72, RMSE = 2.01 gC·m−2·d−1) model based on the big-leaf assumption, our proposed two-leaf hybrid model has the best performance on GPP estimation (R2 = 0.77, RMSE = 1.79 gC·m−2·d−1). We also applied this two-leaf hybrid model to estimate the global GPP during the main growing season in Northern Hemisphere, which were highly correlated with several existing GPP products, with R2 ranging from 0.79 to 0.88. These results will improve our understanding of the relationship between SIF and GPP for sunlit and shaded leaves and will advance application of satellite SIF data to GPP estimation.

A soil-air temperature model to determine the start of season phenology of deciduous forests

Forest ecosystems play a major role in sequestering atmospheric carbon dioxide, which can help offset the detrimental effects of anthropogenic carbon emissions. However, climate change has and will continue to affect the phenology of forest ecosystems’ carbon uptake, changing both the “carbon uptake transition date” - when forests shift from being a net carbon source to sink - and the “green-up date” reflecting the onset of bud burst. Previous studies have shown that a forest's carbon uptake transition date correlates to the date when soil temperature warms enough to surpass mean annual air temperature (soil-air temperature model). However, we still don't know if this simple relationship holds across different sites or over longer time periods. In this study, we explore the relationship between climate and both types of phenological transition dates using over 200 site years of data between 1997 and 2022. Using flux tower data from 18 sites across North America and Europe, we derive three potential carbon uptake transition dates corresponding to the dates when 10%, 25%, and 50% of seasonal net ecosystem exchange (NEE) amplitude is reached. Using PhenoCam data, we then derive three potential green-up dates corresponding to when 10%, 25%, and 50% of total seasonal green chromatic coordinate (GCC) is reached (the greenness model). We evaluate our model estimates using concordance coefficients, a metric of agreement between two measures, to determine which process, carbon uptake or budburst, is best predicted by the soil-air temperature model and to what extent. We find that variation in phenological relationships can be attributed to different regional and bioclimatic groups, highlighting potential biome-specific strengths and limitations of the soil-air temperature model. This model offers a simple approach to better understand phenological transitions and identify potential and limitations for a simple universal SOS prediction approach in deciduous forests.

Satellite-based assessment of water use and leaf area efficiencies of dryland conifer forests along an aridity gradient

Dryland forests worldwide are increasingly threatened by drought stress due to climate change. Understanding the relationships between forest structure and function is essential for managing dryland forests to adapt to these changes. We investigated the structure-function relationships in four dryland conifer forests distributed along a semiarid to subhumid climatic aridity gradient. Forest structure was represented by leaf area index (LAI) and function by gross primary productivity (GPP), evapotranspiration (ET), and the derived efficiencies of water use (WUE = GPP/ET) and leaf area (LAE = GPP/LAI). Estimates of GPP and ET were based on the observed relationships between high-resolution vegetation indices from VENμS and Sentinel-2A satellites and flux data from three eddy covariance towers in the study regions between November 2015 to October 2018. The red-edge-based MERIS Terrestrial Chlorophyll Index (MTCI) from VENμS and Sentinel-2A showed strong correlations to flux tower GPP and ET measurements for the three sites (R2cal > 0.91, R2val > 0.84). Using our approach, we showed that as LAI decreased with decreasing aridity index (AI) (i.e., dryer conditions), estimated GPP and ET decreased (R2 > 0.8 to LAI), while WUE (R2 = 0.68 to LAI) and LAE increased. The observed global-scale patterns are associated with a variety of forest vegetation characteristics, at the local scale, such as tree species composition and density. However, our results point towards a canopy-level mechanism, where the ecosystem-LAI and resultant proportion of sun-exposed vs. shaded leaves are primary determinants of WUE and LAE along the studied climatic aridity gradient. This work demonstrates the importance of high-resolution (spatially and spectrally) remote sensing data conjugated with flux tower data for monitoring dryland forests and understanding the intricate structure-function interactions in their response to drying conditions.

The effects of flash drought on the terrestrial ecosystem in Korea

Typified by rapid intensification in a short time, flash droughts (FDs) have attracted increasing attention because of their effects on water resources and ecosystems. Water use efficiency (WUE) is a key indicator for understanding the interaction between the carbon dynamic and water cycle. However, the effects of FDs on terrestrial ecosystems in Korea remain mostly unknown. Therefore, we identified FD events in Korea using the soil moisture (SM)-based rate of intensification (RI) and the evapotranspiration (ET)-based Standardized Evaporative Stress Ratio (SESR) with flux tower data. Continuous satellite data from the Moderate Resolution Imaging Spectroradiometer (MODIS) were used to examine spatiotemporal patterns and trends in FD-related factors, Potential ET (PET), Actual ET (AET), Gross Primary Productivity (GPP), Evaporative Stress Ratio (ESR), and WUE from 2001 to 2020. The frequency, duration, and intensity of FDs in Korea were investigated and responses in WUE to FDs in different landcover types were examined. The SESR and RI had a similar performance in capturing FD events in flux tower site data. At the regional scale, a lower ESR was mainly observed in the western and southern parts of the study area, indicating progressively drier climate conditions. The mean frequency percentage and number of Korean FD events in the past 20 years were 61.3% and 16, respectively. In addition, more FD events occurred in the northern and eastern regions of the study area, where the main landcover type is forest, and northeastern Korea suffered from long FD events (over 30 days). FDs caused a positive response of WUE in most of the study area. However, the change in WUE after FD events varied by FD severity and landcover type. The changes in WUE after moderate FD events were largely driven by the change in GPP, and the highest sensitivity of WUE to FD intensity was found in croplands. These findings highlight the characteristics and effects of FDs on WUE across landcover types, which will be helpful in managing potential environmental risks in Korea.

A Hybrid Framework for Simulating Actual Evapotranspiration in Data-Deficient Areas: A Case Study of the Inner Mongolia Section of the Yellow River Basin

Evapotranspiration (ET) plays an important role in transferring water and converting energy in the land–atmosphere system. Accurately estimating ET is crucial for understanding global climate change, ecological environmental problems, the water cycle, and hydrological processes. Machine learning (ML) algorithms have been considered as a promising method for estimating ET in recent years. However, due to the limitations associated with the spatial–temporal resolution of the flux tower data commonly used as the target set in ML algorithms, the ability of ML to discover the inherent laws within the data is reduced. In this study, a hybrid framework was established to simulate ET in data-deficient areas. ET simulation results of a coupled model comprising the Budyko function and complementary principle (BC2021) were used as the target set of the random forest model, instead of using the flux station observation data. By combining meteorological and hydrological data, the monthly ET of the Inner Mongolia section of the Yellow River Basin (IMSYRB) was simulated from 1982 to 2020, and good results were obtained (R2 = 0.94, MAE = 3.82 mm/mon, RMSE = 5.07 mm/mon). Furthermore, the temporal and spatial variations in ET and the influencing factors were analysed. In the past 40 years, annual ET in the IMSYRB ranged between 241.38 mm and 326.37 mm, showing a fluctuating growth trend (slope = 0.80 mm/yr), and the summer ET accounted for the highest proportion in the year. Spatially, ET in the IMSYRB showed a regular distribution of high ET in the eastern region and low ET in the western area. The high ET value areas gradually expanded from east to west over time, and the area increased continuously, with the largest increase observed in the 1980s. Temperature, precipitation, and normalized difference vegetation index (NDVI) were found to be the most important factors affecting ET in the region and play a positive role in promoting ET changes. These results provide an excellent example of long-term and large-scale accurate ET simulations in an area with sparse flux stations.

Meteorological Influences on Short-Term Carbon-Water Relationships in Two Forests in Subtropical China

Carbon and water fluxes in ecosystems are tightly coupled by gas diffusion through stomata. However, carbon–water (C–W) relationships vary largely across time scales, vegetation types, and regions. Subtropical forests in China play an important role in the global carbon and water cycles, yet studies of C–W relationships in this region remain limited. Here, we investigated summer-time C–W relationships in this region at two subtropical sites: the evergreen broadleaved forest at Dinghushan (23.17° N, 112.53° E, 300 m) and the evergreen coniferous forest at Qianyanzhou (26.74° N, 115.06° W, 106 m), using the flux tower data from the FLUXNET2015. The C–W relationship was examined using two measures. The first was daily water use efficiency (WUE), which is the ratio of daily gross primary productivity (GPP) to evapotranspiration (ET). The second was the correlation coefficient (r) of hourly GPP and ET. Our analysis showed that the daily WUE in the two forests ranged over 4–14 mg CO2 per g H2O, higher in the coniferous forest than in the broadleaved forest. The mean values of r for hourly C–W coupling were similar at the two forests, being 0.5–0.6, which suggests asynchronous diurnal variations in GPP and ET. Both daily WUE and r were modulated by meteorological conditions. In general, high radiation, air temperature, and humidity can reduce WUE at both sites. For the broadleaved forest, the most influential factor on WUE was VPD, followed by radiation, while in the coniferous forest, VPD, air temperature, and radiation were almost equally important. For hourly C–W coupling, VPD plays a significant role. The drier the air is, the weaker the coupling in the two forests. The daily WUE and hourly C–W coupling reflect the C–W relationship from different perspectives. Both showed the strongest response to VPD but with different sensitivity.

ET Partitioning Assessment Using the TSEB Model and sUAS Information across California Central Valley Vineyards

Evapotranspiration (ET) is a crucial part of commercial grapevine production in California, and the partitioning of this quantity allows the separate assessment of soil and vine water and energy fluxes. This partitioning has an important role in agriculture since it is related to grapevine stress, yield quality, irrigation efficiency, and growth. Satellite remote sensing-based methods provide an opportunity for ET partitioning at a subfield scale. However, medium-resolution satellite imagery from platforms such as Landsat is often insufficient for precision agricultural management at the plant scale. Small, unmanned aerial systems (sUAS) such as the AggieAir platform from Utah State University enable ET estimation and its partitioning over vineyards via the two-source energy balance (TSEB) model. This study explores the assessment of ET and ET partitioning (i.e., soil water evaporation and plant transpiration), considering three different resistance models using ground-based information and aerial high-resolution imagery from the Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX). We developed a new method for temperature partitioning that incorporated a quantile technique separation (QTS) and high-resolution sUAS information. This new method, coupled with the TSEB model (called TSEB-2TQ), improved sensible heat flux (H) estimation, regarding the bias, with around 61% and 35% compared with the H from the TSEB-PT and TSEB-2T, respectively. Comparisons among ET partitioning estimates from three different methods (Modified Relaxed Eddy Accumulation—MREA; Flux Variance Similarity—FVS; and Conditional Eddy Covariance—CEC) based on EC flux tower data show that the transpiration estimates obtained from the FVS method are statistically different from the estimates from the MREA and the CEC methods, but the transpiration from the MREA and CEC methods are statistically the same. By using the transpiration from the CEC method to compare with the transpiration modeled by different TSEB models, the TSEB-2TQ shows better agreement with the transpiration obtained via the CEC method. Additionally, the transpiration estimation from TSEB-2TQ coupled with different resistance models resulted in insignificant differences. This comparison is one of the first for evaluating ET partitioning estimation from sUAS imagery based on eddy covariance-based partitioning methods.

High-resolution crop yield and water productivity dataset generated using random forest and remote sensing

Accurate and high-resolution crop yield and crop water productivity (CWP) datasets are required to understand and predict spatiotemporal variation in agricultural production capacity; however, datasets for maize and wheat, two key staple dryland crops in China, are currently lacking. In this study, we generated and evaluated a long-term data series, at 1-km resolution of crop yield and CWP for maize and wheat across China, based on the multiple remotely sensed indicators and random forest algorithm. Results showed that MOD16 products are an accurate alternative to eddy covariance flux tower data to describe crop evapotranspiration (maize and wheat RMSE: 4.42 and 3.81 mm/8d, respectively) and the proposed yield estimation model showed accuracy at local (maize and wheat rRMSE: 26.81 and 21.80%, respectively) and regional (maize and wheat rRMSE: 15.36 and 17.17%, respectively) scales. Our analyses, which showed spatiotemporal patterns of maize and wheat yields and CWP across China, can be used to optimize agricultural production strategies in the context of maintaining food security.

Remote sensing-based evapotranspiration modeling using geeSEBAL for sugarcane irrigation management in Brazil

Irrigated agriculture requires the implementation of innovative tools to improve irrigation water management and accurate estimation of actual evapotranspiration (ETa) such as remote sensing-based methodology. This study aimed to evaluate the irrigation management and estimating evapotranspiration through the geeSEBAL, a new tool for automated estimation of ETa based on the Surface Energy Balance Algorithm for Land (SEBAL) and a simplified version of the Calibration using Inverse Modeling at Extreme Conditions (CIMEC) process for the endmembers selection, implemented into the Google Earth Engine (GEE) environment. GeeSEBAL has not been used yet in Brazil for irrigation proposes, and in this research, it was applied to estimate ETa using Landsat images and ERA5-Land as meteorological inputs in the largest sugarcane producing region of the world in Brazil for two ratoon seasons by comparing daily ETa with values obtained from eddy covariance (EC) data, Energy balance components using geeSEBAL were consistent with the measured data and daily ETa presenting RMSE of 0.46 mm with R2 = 0.97. Modeled ETa and Kc were similar for the two seasons, although somewhat overestimated for the fifth ratoon when compared to the EC data, mainly during high atmospheric demand (crop mid-stage). Still, the Kc values were similar to the standard values available in the literature and measured flux tower data for the two ratoons seasons. With ETa from geeSEBAL it was possible to identify water stress over the growing seasons using the remote sensing-based soil water balance, which occurred mainly during the phase after the crop reached the peak Kc (full cover stage) when the irrigation depth required was very high. This analysis showed that geeSEBAL has a significant potential for assessment of ETa for irrigation monitoring and management, even in missing climate data areas, allowing important advances in water resources management for sugarcane and other irrigated crops at field or regional scales.

Cropland Net Ecosystem Exchange Estimation for the Inland Pampas (Argentina) Using EVI, Land Cover Maps, and Eddy Covariance Fluxes

Estimations of Net Ecosystem Exchange (NEE) are crucial to assess the carbon sequestration/carbon source capacity of agricultural systems. Although several global models have been built to describe carbon flux patterns based on flux tower data, South American ecosystems (and croplands in particular) are underrepresented in the databases used to calibrate these models, leading to large uncertainties in regional and global NEE estimation. Despite the fact that almost half of the land surface is used worldwide for agricultural activities, these models still do not include variables related to cropland management. Using enhanced vegetation index (EVI) derived from MODIS imagery (250 m) and monthly CO2 exchange from a 9-year record of an eddy covariance (EC) flux tower in a crop field in the Inland Pampas region, we developed regression models to predict monthly NEE. We tested whether including a term for crop identity/land cover as a categorical variable (maize, soybean, wheat, and fallow) could improve model capability in capturing monthly NEE dynamics. NEE measured at the flux tower site was scaled to croplands across the Inland Pampa using crop-type maps, from which annual NEE maps were generated for the 2018–2019, 2019–2020, and 2020–2021 agricultural campaigns. The model based solely on EVI showed to be a good predictor of monthly NEE for the study region (r2 = 0.78), but model adjustment was improved by including a term for crop identity (r2 = 0.83). A second set of maps was generated taking into account carbon exports during harvest to estimate Net Biome Productivity (NBP) at the county level. Crops across the region as a whole acted as a carbon sink during the three studied campaigns, although with highly heterogeneous spatial and temporal patterns. Between 60% and 80% of the carbon sequestered was exported during harvest, a large decrease from the carbon sequestration capacity estimated using just NEE, which further decreased if fossil carbon emissions from agricultural supplies are taken into account. Estimates presented in this study are a first step towards upscaling carbon fluxes at the regional scale in a South American cropland area, and could help to improve regional to global estimations of carbon fluxes and refine national greenhouse gas (GHG) inventories.

Downwelling longwave radiation and sensible heat flux observations are critical for surface temperature and emissivity estimation from flux tower data

Land surface temperature (LST) is a preeminent state variable that controls the energy and water exchange between the Earth’s surface and the atmosphere. At the landscape-scale, LST is derived from thermal infrared radiance measured using space-borne radiometers. In contrast, plot-scale LST estimation at flux tower sites is commonly based on the inversion of upwelling longwave radiation captured by tower-mounted radiometers, whereas the role of the downwelling longwave radiation component is often ignored. We found that neglecting the reflected downwelling longwave radiation leads not only to substantial bias in plot-scale LST estimation, but also have important implications for the estimation of surface emissivity on which LST is co-dependent. The present study proposes a novel method for simultaneous estimation of LST and emissivity at the plot-scale and addresses in detail the consequences of omitting down-welling longwave radiation as frequently done in the literature. Our analysis uses ten eddy covariance sites with different land cover types and found that the LST values obtained using both upwelling and downwelling longwave radiation components are 0.5–1.5 K lower than estimates using only upwelling longwave radiation. Furthermore, the proposed method helps identify inconsistencies between plot-scale radiometric and aerodynamic measurements, likely due to footprint mismatch between measurement approaches. We also found that such inconsistencies can be removed by slight corrections to the upwelling longwave component and subsequent energy balance closure, resulting in realistic estimates of surface emissivity and consistent relationships between energy fluxes and surface-air temperature differences. The correspondence between plot-scale LST and landscape-scale LST depends on site-specific characteristics, such as canopy density, sensor locations and viewing angles. Here we also quantify the uncertainty in plot-scale LST estimates due to uncertainty in tower-based measurements using the different methods. The results of this work have significant implications for the combined use of aerodynamic and radiometric measurements to understand the interactions and feedbacks between LST and surface-atmosphere exchange processes.