Land-surface evapotranspiration derived from a first-principles primary production model

Estimation of crop gross primary production (GPP): I. impact of MODIS observation footprint and impact of vegetation BRDF characteristics

Thermally derived evapotranspiration from the Surface Temperature Initiated Closure (STIC) model improves cropland GPP estimates under dry conditions

Global estimates of evapotranspiration and gross primary production based on MODIS and global meteorology data

Gross Primary Production (GPP) and evapotranspiration (ET) play a critical role of the global carbon, water and energy cycle. Accurate quantification of the global GPP and ET could improve our ability to understand global climate change and energy budget. However, most of the GPP and ET remote sensing models fail to take the coupled relationship between vegetation transpiration (Et) and photosynthesis into consideration. More importantly, these models might ignore the difference of transpiration and photosynthesis rate in different groups of leaves (sunlit and shaded). Here, we coupled the estimates of daily GPP and ET at 84 global forest sites based on the Two-Leaf Light Use Efficiency model and the Penman-Monteith equation that were linked by the Ball-Berry conductance model. The developed model was well calibrated with the root mean square error (RMSE) and the coefficient of determination (R2) being 1.97 gC/m2 day and 0.77 for GPP respectively, and being 21.91 W/m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> and 0.65 for ET, respectively. In the meantime, the validation results demonstrated the good performance of the coupled model, with the RMSE and R2 being 1.95 gC/m2 day and 0.77 for GPP, respectively, and being 21.34 W/m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> and 0.66 for ET, respectively.

Maize is one of the most important crops cultivated on the global scale. Accurate estimation of maize Gross Primary Production (GPP) can provide valuable information for regional and global carbon budget studies. From site level to regional/global scales, GPP estimation depends on remote sensing or eddy covariance flux data. In this research, the 8-day composite GPP of maize was estimated by Moderate Resolution Imaging Spectroradiometer (MODIS) and flux tower data at eight study sites using a Regional Production Efficiency Model (REG-PEM). The performance of the model was assessed by analyzing the linearly regression of GPP estimated from the REG-PEM model (GPPEST) with the GPP predicted from the eddy covariance data (GPPEC). The coefficient of determination, root mean squared error and mean absolute error of the regression model were calculated. The uncertainties of the model are also discussed in this research. The seasonal dynamics (phases and magnitudes) of the GPPEST reasonably agreed with those of GPPEC, indicating the potential of the satellite-driven REG-PEM model for up-scaling the GPP in maize croplands. Furthermore, the maize GPP estimated by this model is more accurate than the MODIS GPP products (MOD17A2). In particular, MOD17A2 significantly underestimated the GPP of maize croplands. The uncertainties in the REG-PEM model are mostly contributed by the maximum light use efficiency and the fraction of photosynthetically active radiation. Keywords: gross primary production, maize, light use efficiency, MODIS, remote sensing, flux data DOI: 10.3965/j.ijabe.20160902.2245 Citation: Li S H, He P, Liu B S, Ni P, Han X. Modeling of maize gross primary production using MODIS imagery and flux tower data. Int J Agric & Biol Eng, 2016; 9(2): 110－118.

Accurate high-resolution actual evapotranspiration (ET) and gross primary production (GPP) information is essential for understanding the large-scale water and carbon dynamics. However, substantial uncertainties exist in the current ET and GPP datasets in China because of insufficient local ground measurements used for model constrain. This study utilizes a water-carbon coupled model, Penman-Monteith-Leuning Version 2 (PML-V2), to estimate 500 m ET and GPP at a daily scale. The parameters of PML-V2 (China) were well-calibrated against observations of 26 eddy covariance flux towers across nine plant functional types in China, indicated by Nash–Sutcliffe Efficiency (NSE) of 0.75 and Root Mean Square Error (RMSE) of 0.69 mm d−1 for daily ET respectively, and NSE of 0.82 and RMSE of 1.71 g C m-2 d−1 for daily GPP. The model estimates get a small bias of 6.28 % and a high NSE of 0.82 against water‐balance annual ET estimates across 10 major river basins in China. Further evaluations suggest that the newly developed product outperforms its global version and other typical products (MOD16A2, SEBAL, GLEAM, MOD17A2H, VPM, and EC-LUE) in estimating both ET and GPP. Moreover, PML-V2 (China) accurately monitors the intra-annual variations in ET and GPP in the croplands with a dual-cropping system. Using the new data showed that, over the last 20 years, the annual GPP and water use efficiency experienced a significant (p < 0.001) increase (8.51 g C m-2 yr-1 and 0.02 g C mm-1 H2O yr-1, respectively), but annual ET showed a non-significant (p > 0.05) increase (0.65 mm yr-1). This indicates that vegetation in China exhibits a huge potential for carbon sequestration with little cost in water resources. The PML-V2 (China) product provides a great opportunity for academic communities and various agencies for scientific studies and applications, freely available at http://dx.doi.org/10.11888/Terre.tpdc.272389.

According to 230 flux site observations, intra-day and diurnal ground heat flux (G) accounted for 19.2 % and 28.8 % of the corresponding net radiation, respectively. This indicates that G plays an important role in remote sensing (RS) evapotranspiration (ET) models. The accuracy of the five G simulation methods in the surface energy balance-based RS ET models was evaluated using half-hourly observations. The linear coefficient (LC) method and the methods embedded with the normalized difference vegetation index (NDVI) were able to accurately simulate a half-hourly G series at most sites. The mean and median Nash-Sutcliffe efficiency (NSE) values of all sites were generally higher than 0.50 in each half-hour period. However, the two methods embedded by fractional vegetation coverage showed poor performance in most half-hour periods, except the sunrise and sunset periods, with the mean and median NSE values of all sites below 0.20. The accuracy of each method varied significantly at different sites and at half-hour intervals. The highest accuracy was exhibited during sunrise periods (6:00–7:00), followed by sunset periods (17:00–18:00). There were 92 % (211/230) sites with an NSE of the LC method greater than 0.50 at 6:30. It showed a slightly higher accuracy during night periods than during daytime periods. The lowest accuracy was observed at noon periods (10:00–15:30). The sites with an NSE exceeding 0.50 only accounted for 51 % (118/230) and 43 % (100/230) at 10:30 and 13:30, respectively. The accuracy of Northern Hemisphere sites was generally higher than that of Southern Hemisphere sites. In general, the highest and lowest accuracies were observed at the high- and low-latitude sites, respectively. The performance of the LC method and the methods embedded with NDVI were generally satisfactory at the Eurasian and North American sites, with the NSE values of most sites exceeding 0.70. Conversely, it exhibited relatively poor performance at the African, South American, and Oceanian sites, especially the African sites. Both the temporal and spatial distributions of the accuracy of the G simulation were positively correlated with the correlation between G and the net radiation. Although the G simulation methods accurately simulated the G series at most sites and time periods, their performance was poor at some sites and time periods. The application of RS ET datasets covering these sites requires caution. The optimal parameter value for each method varied greatly at different sites. This indicates that the fixed parameter values in the G simulation methods are not appropriate. Further improvement of G simulations at these sites and time periods is recommended for the RS ET modelers.

An attempt to introduce atmospheric CO2 concentration data to estimate the gross primary production by the terrestrial biosphere and analyze its effects

An improved Biome-BGC model for estimating net primary productivity of alpine meadow on the Qinghai-Tibet Plateau

<strong class="journal-contentHeaderColor">Abstract.</strong> Earth's drylands are home to more than two billion people, provide key ecosystem services, and exert a large influence on the trends and variability in Earth's carbon cycle. However, modeling dryland carbon and water fluxes with remote sensing suffers from unique challenges not typically encountered in mesic systems, particularly in capturing soil moisture stress. Here, we develop and evaluate an approach for the joint modeling of dryland gross primary production (GPP), net ecosystem exchange (NEE), and evapotranspiration (ET) in the western United States (US) using a suite of AmeriFlux eddy covariance sites spanning major functional types and aridity regimes. We use artificial neural networks (ANNs) to predict dryland ecosystem fluxes by fusing optical vegetation indices, multitemporal thermal observations, and microwave soil moisture and temperature retrievals from the Soil Moisture Active Passive (SMAP) sensor. Our new dryland ANN (DrylANNd) carbon and water flux model explains more than 70 % of monthly variance in GPP and ET, improving upon existing MODIS GPP and ET estimates at most dryland eddy covariance sites. DrylANNd predictions of NEE were considerably worse than its predictions of GPP and ET likely because soil and plant respiratory processes are largely invisible to satellite sensors. Optical vegetation indices, particularly the normalized difference vegetation index (NDVI) and near-infrared reflectance of vegetation (NIR<span class="inline-formula">_v</span>), were generally the most important variables contributing to model skill. However, daytime and nighttime land surface temperatures and SMAP soil moisture and soil temperature also contributed to model skill, with SMAP especially improving model predictions of shrubland, grassland, and savanna fluxes and land surface temperatures improving predictions in evergreen needleleaf forests. Our results show that a combination of optical vegetation indices and thermal infrared and microwave observations can substantially improve estimates of carbon and water fluxes in drylands, potentially providing the means to better monitor vegetation function and ecosystem services in these important regions that are undergoing rapid hydroclimatic change.

The measurement of carbon exchange between vegetation and the atmosphere is vital to quantify the impact of environmental variables on the carbon sequestration capacity of forests, and to predict how they will respond to future climate. In this study we use proximal remote sensing, defined as observations made from non-contact radiometric or imaging sensors in close proximity to the forest canopy (10–20 m), as an intermediate upscaling tool between direct measurements of carbon fluxes and satellite-derived estimations of primary productivity in a tropical dry forest (TDF) in Jalisco, Mexico. Two broad-band vegetation indices (VIs), the normalized difference VI and the enhanced vegetation index 2 (EVI2), were calculated from proximally sensed canopy properties, validated with field estimates of the fraction of absorbed photosynthetically active radiation by photosynthetic tissue (fAPARgreen), and compared to estimates of gross primary productivity (GPP) and net ecosystem exchange of CO2, measured from a flux tower. The VIs captured the phenology of the TDF, both under typical summer rainfall and during an atypically-dry wet season in El Niño of 2009. The VIs also tracked a secondary leaf-flushing in the dry season of 2010. Our study suggests that (1) VIs are the best predictors of gross carbon uptake, able to explain up to 86% of variations in GPP; (2) VIs are accurate predictors of the photosynthetic capacity of green tissue, able to explain up to 99% of fAPARgreen variation; and (3) VIs and soil water content can be used to develop an empirical model that captures the seasonal trajectory of GPP from high respiration after the rain pulses, to rapid leaf development, and finally to slow senescence as the soil dries out. Proximal remote sensing constitutes a useful tool to link field-base measurements of carbon fluxes to satellite- or airborne-derived estimates of carbon exchange.

Accurate estimation of reference evapotranspiration (ETo) plays a vital role in irrigation and water resource planning. The Penman–Monteith method recommended by the Food and Agriculture Organization (FAO PM56) is widely used and considered a standard to calculate ETo. However, FAO PM56 cannot be used with limited meteorological variables, so it is compulsory to choose an alternative model for ETo estimation, which requires fewer variables. This study built ten machine learning (ML) models based on multi-function, neural network, and tree-based structure against the FAO PM56 method. For this purpose, time series temperature data on a monthly scale are only used to train ML models. The developed ML models were applied to estimate ETo at different test stations and the obtained results were compared with the FAO PM56 method to verify and validate their performance in ETo estimation for the selected stations. In addition, multiple statistical indicators, including root-mean-square error (RMSE), coefficient of determination (R2), mean absolute error (MAE), Nash–Sutcliffe efficiency (NSE), and correlation coefficient (r) were calculated to compare the performance of each ML model on ETo estimation. Among the applied ML models, the ETo tree boost (TB) ML model outperformed the other ML models in estimating ETo in diverse climatic conditions based on statistical indicators (R2, NSE, r, RMSE, and MAE). Moreover, the observed R2, NSE, and r were the highest for the TB ML model, while RMSE and MAE were found to be the lowest at the study sites compared to other applied ML models. Lastly, ETo point data yielded from the TB ML model was used in an interpolation process to create monthly and annual ETo maps. Based on the ETo maps, this study suggests mainly a focus on areas with high ETo values and proper irrigation scheduling of crops to ensure water sustainability.

Ecosystem water use efficiency in a warm-temperate mixed plantation in the North China

<strong class="journal-contentHeaderColor">Abstract.</strong> We develop high resolution (1 km) estimates of Gross Primary Productivity (GPP), Ecosystem Respiration (ER) and Net Ecosystem Exchange (NEE) over the Australian continent for the period January 2003 to June 2022 by empirical upscaling of flux tower measurements. We compare our estimates with nine other products that cover the three broad categories that define current methods for estimating the terrestrial carbon cycle and assess if consiliences between datasets can point to the correct dynamics of Australia&rsquo;s carbon cycle. Our results indicate that regional empirical upscaling greatly improves upon the existing global empirical upscaling efforts, outperforms process-based models, and agrees much better with the dynamics of CO₂ flux over Australia as estimated by two regional atmospheric inversions. Our nearly 20-year estimates of terrestrial carbon fluxes revealed Australia is a strong net carbon sink of -0.44 (IQR=0.42) PgC/year on-average, with an inter-annual variability of 0.18 PgC/year and an average seasonal amplitude of 0.85 PgC/yr. Annual mean carbon uptake estimated from other methods ranged considerably, while carbon flux anomalies showed much better agreement between methods. NEE anomalies were predominately driven by cumulative rainfall deficits and surpluses, resulting in larger anomalous responses from GPP over ER. In contrast, we show that the long-term average seasonal cycle is dictated more by the variability in ER than GPP, resulting in peak carbon uptake typically occurring during the cooler, drier Austral autumn, and winter months. This new estimate of Australia&rsquo;s terrestrial carbon cycle provides a benchmark for assessment against Land Surface Model simulations, and a means for monitoring of Australia&rsquo;s terrestrial carbon cycle at an unprecedented high-resolution. We call this new estimate of Australia&rsquo;s terrestrial carbon cycle, &ldquo;AusEFlux&rdquo; (Australian Empirical Fluxes).

Reliable modeling of carbon and water fluxes is essential for understanding the terrestrial carbon and water cycles and informing policy strategies aimed at constraining carbon emissions and improving water use efficiency. We used an assimilation framework (LPJ-Vegetation and soil moisture Joint Assimilation, or LPJ-VSJA) to improve gross primary production (GPP) and evapotranspiration (ET) estimates globally. The terrestrial biosphere model that we used is the integrated model – LPJ-PM coupled from the Lund-Potsdam-Jena Dynamic Global Vegetation Model (LPJ-DGVM) and a hydrology module (i.e., the updated Priestley–Taylor Jet Propulsion Laboratory model, PT-JPLSM). Satellite-based soil moisture products derived from the Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active and Passive (SMAP) and leaf area index (LAI) from the global Land and Ground satellite (GLASS) product were assimilated into LPJ-PM to improve GPP and ET simulations using a Proper Orthogonal Decomposition-based ensemble four-dimensional variational assimilation method (PODEn4DVar). The joint assimilation framework LPJ-VSJA achieved the best model performance (with an R2 of 0.91 and 0.81 and an RMSD reduced by 50.4 % and 38.4 % for GPP and ET, respectively, compared with those of LPJ-DGVM at the monthly scale). The assimilated GPP and ET demonstrated a better performance in the arid and semiarid regions (GPP: R2 = 0.73, ubRMSD = 1.05 g C m−2 d−1; ET: R2 = 0.73, ubRMSD = 0.61 mm d−1) than in the humid and sub-dry humid regions (GPP: R2 = 0.61, ubRMSD = 1.23 g C m−2 d−1; ET: R2 = 0.66; ubRMSD = 0.67 mm d−1). The ET simulated by LPJ-PM that assimilated SMAP or SMOS had a slight difference, and the ET that assimilated SMAP soil moisture data was more improved than that assimilated SMOS data. Our global simulation modeled by LPJ-VSJA was compared with several global GPP and ET products (e.g., GLASS GPP, GOSIF GPP, GLDAS ET, GLEAM ET) using the triple collocation (TC) method. Our products, especially ET, exhibited advantages in the overall error distribution (estimated error (μ): 3.4 mm month−1; estimated standard deviation of μ: 1.91 mm month−1). Our research showed that the assimilation of multiple datasets could reduce model uncertainties, while the model performance differed across regions and plant functional types. Our assimilation framework (LPJ-VSJA) can improve the model simulation performance of daily GPP and ET globally, especially in water-limited regions.