Leaf Area Index Data Research Articles

A Physically Based Soil Moisture Index From Passive Microwave Brightness Temperatures for Soil Moisture Variation Monitoring

Soil moisture is a pivotal hydrological variable that links the terrestrial water, energy, and carbon cycles. In this article, a new soil moisture (SM) index (SMI), which aims to capture the temporal variability of SM, irrespective of cloud cover and solar illumination, was developed by using the L-band SM active passive (SMAP) radiometer observations. The SMI was proposed on the basis of two key foundations: 1) vegetation and roughness have similar effects on “depolarization” of microwave emission, while SM enhances polarization differences and 2) vegetation and roughness generally impose positive effects on surface emissivity, while SM and emissivity are negatively correlated. Based on the two physical principles, it is possible to decouple the effects of SM and those of vegetation and surface roughness in a 2-D space independent of vegetation type and roughness condition. The proposed SMI was then validated by in situ measurements from five dense SM networks covering different vegetation and climatic conditions and also compared with SMAP passive and European space agency climate change initiative (ESA CCI) SM products at a coarse resolution of 36 km, and SMAP-enhanced passive and Japan Aerospace Exploration Agency (JAXA) advanced microwave scanning radiometer (AMSR2) SM products at a medium resolution of 9 km. The results show that the new SMI is able to well reproduce the temporal dynamic of SM with a favorable averaged correlation coefficient value of 0.87 and 0.84 at 36 and 9 km, respectively, higher than that of SMAP passive (0.80), SMAP-enhanced passive (0.77), ESA CCI (0.69), and JAXA AMSR2 (0.53). After removing the systematic differences between satellite and site-specific SM data by using the cumulative distribution function (CDF) matching technique, the SMI can achieve an average root mean squared error (RMSE) of 0.031 and 0.036 m3m−3 at 36 and 9 km during the validation period, respectively, lower than that of the satellite SM products. In addition to surface temperature, the SMI does not need any further information from other sensors [e.g., the optical normalized difference vegetation index (NDVI) or leaf area index (LAI) data] to guarantee an all-weather monitoring. Therefore, it has great potential to estimate SM variability on a global scale.

Outage prediction models for snow and ice storms

Abstract This paper describes the development of two Outage Prediction Models (OPMs) for power outages caused by snow and ice storms on electric distribution networks, and their performance evaluation in the Northeastern United States. The first is a Machine Learning (ML) based model that is set up to predict power outages on a regular grid (4-km grid spacing). The second is a Generalized Linear Model (GLM) that is set up to predict total power outages at the town level. Both models are fed with Numerical Weather Prediction (NWP) outputs, satellite-derived leaf area index (LAI) and land cover data, and utility specified infrastructure and historical outage data. For both models the most important variables are the amount of assets on the territory, LAI, snow density, and – for the ice model – the amount of freezing rain. Results from cross-validation experiments based on 54 outage events spanning three orders of magnitude show median absolute percentage errors around 70% for both models. The GLM is able to predict extreme events (such as the devastating October 2011 nor’easter) better than ML models do, while ML models have better performance for lower impact events and present better characterization of the spatial distribution of power outages.

Adjusting Ceptometer Data to Improve Leaf Area Index Measurements

Leaf Area Index (LAI) is an important plant parameter for both farmers and plant scientists to monitor and/or model the growth and the well-being of plants. Since direct LAI measurement techniques are relatively laborious and time-consuming, various indirect methods have been developed and widely used since the early 1990s. The LP-80 ceptometer uses a linear array of PAR (photosynthetically active radiation) sensors for non-destructive LAI measurements that is backed by 15 years of research. Despite this, considerable discrepancy can be found between the expert opinions regarding the optimal illumination conditions recommended for the measurement. The sensitivity of ceptometer-based LAI values to PAR was investigated, and a simple method was devised to correct raw ceptometer data collected under non-ideal light conditions. Inadequate light conditions (PAR < 1700 µmol m−2 s−1) could cause an underestimation of LAI. Using the corrected LAI values, the ceptometer data showed a significantly better fit (higher R2, smaller mean average error and closer to zero mean signed error values) to the destructive LAI data for both wheat and maize. With the help of the correction equations, the use of the LP-80 ceptometer could be extended to days when light conditions are not ideal.

ESTIMATION METHODOLOGY FOR FOREST BIOMASS IN MONGOLIA USING REMOTE SENSING

Abstract. The forest biomass is one of the most important parameters for the global carbon stock. Information on the forest volume, coverage and biomass are important to develop the global perspective on the CO2 concentration changes. Objective of this research is to estimate forest biomass in the study area. The study area is Hangal sum, Bulgan province, Mongolia. Backscatter coefficients for vertical transmit and vertical receive (VV), for vertical transmit and horizontal receive (VH) from Sentinel data and Leaf Area Index (LAI) from Landsat data were used in the study area. We developed biomass estimation approach using ground truth data which is DBH, height and soil moisture. The coefficient α, β, δ, γ were found from the approach. The output map from the approach was compared with VV and VH, LAI data. The relationship between output map and VH data shows a positive result R2 = 0.61. This study suggests that the biomass estimation using Remote sensing data can be applied in forest region in the North.

Multi-Scale Validation of MODIS LAI Products Based on Crop Growth Period

Leaf area index (LAI) is one of the most important canopy structure parameters utilized in process-based models of climate, hydrology, and biogeochemistry. In order to determine the reliability and applicability of satellite LAI products, it is critical to validate satellite LAI products. Due to surface heterogeneity and scale effects, it is difficult to validate the accuracy of LAI products. In order to improve the spatio-temporal accuracy of satellite LAI products, we propose a new multi-scale LAI product validation method based on a crop growth cycle. In this method, we used the PROSAIL model to derive Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) LAI data and Gaofen-1 (GF-1) for the study area. The Empirical Bayes Kriging (EBK) interpolation method was used to perform a spatial multi-scale transformation of Moderate Resolution Imaging Spectroradiometer (MODIS) LAI products, GF-1 LAI data, and ASTER LAI data. Finally, MODIS LAI satellite products were compared with field measured LAI data, GF-1 LAI data, and ASTER LAI data during the growing season of crop field. This study was conducted in the agricultural oasis area of the middle reaches of the Heihe River Basin in northwestern China and the Conghua District of Guangzhou in Guangdong Province. The results suggest that the validation accuracy of the multi-scale MODIS LAI products validated by ASTER LAI data were higher than those of the GF-1 LAI data and the reference field measured LAI data, showing a R2 of 0.758 and relative mean square error (RRMSE) of 28.73% for 15 m ASTER LAI and a R2 of 0.703 and RRMSE of 30.80% for 500 m ASTER LAI, which imply that the 15 m MODIS LAI product generated by the EBK method was more accurate than the 500 m and 8 m products. This study provides a new validation method for satellite remotely sensed products.

High temporal resolution of leaf area data improves empirical estimation of grain yield

Empirical yield estimation from satellite data has long lacked suitable combinations of spatial and temporal resolutions. Consequently, the selection of metrics, i.e., temporal descriptors that predict grain yield, has likely been driven by practicality and data availability rather than by systematic targetting of critically sensitive periods as suggested by knowledge of crop physiology. The current trend towards hyper-temporal data raises two questions: How does temporality affect the accuracy of empirical models? Which metrics achieve optimal performance? We followed an in silico approach based on crop modelling which can generate any observation frequency, explore a range of growing conditions and reduce the cost of measuring yields in situ. We simulated wheat crops across Australia and regressed six types of metrics derived from the resulting time series of Leaf Area Index (LAI) against wheat yields. Empirical models using advanced LAI metrics achieved national relevance and, contrary to simple metrics, did not benefit from the addition of weather information. This suggests that they already integrate most climatic effects on yield. Simple metrics remained the best choice when LAI data are sparse. As we progress into a data-rich era, our results support a shift towards metrics that truly harness the temporal dimension of LAI data.

Spatially and Temporally Continuous Leaf Area Index Mapping for Crops through Assimilation of Multi-resolution Satellite Data

As a key parameter that represents the structural characteristics and biophysical changes of crop canopy, the leaf area index (LAI) plays a significant role in monitoring crop growth and mapping yield. A considerable amount of farmland is dispersed with strong spatial heterogeneity. The existing time series satellite LAI products fail to capture spatial distributions and growth changes of crops due to coarse spatial resolutions and spatio-temporal discontinuities. Therefore, it becomes crucial for fine resolution LAI mapping in time series over crop areas. A two-stage data assimilation scheme was developed for dense time series LAI mapping in this study. A LAI dynamic model was first constructed using multi-year MODIS LAI data. This model coupled with the PROSAIL radiative transfer model, and MOD09A1 reflectance data were used to retrieve temporal LAI profiles at the 500 m resolution with the assistance of the very fast simulated annealing (VFSA) algorithm. Then, the LAI dynamics at the 500 m scale were incorporated as prior information into the Landsat 8 OLI reflectance data for time series LAI mapping at the 30 m resolution. Finally, the spatio-temporal continuities and retrieval accuracies of assimilated LAI values were assessed at the 500 m and 30 m resolutions respectively, using the MODIS LAI product, fine resolution LAI reference map and field measurements. The results indicated that the assimilated the LAI estimations at the 500 m scale effectively eliminated the spatio-temporal discontinuities of the MODIS LAI product and displayed reasonable temporal profiles and spatial integrity of LAI. Moreover, the 30 m resolution LAI retrievals showed more abundant spatial details and reasonable temporal profiles than the counterparts at the 500 m scale. The determination coefficient R2 between the estimated and field LAI values was 0.76 with a root mean square error (RMSE) value of 0.71 at the 30 m scale. The developed method not only improves the spatio-temporal continuities of the LAI at the 500 m scale, but also obtains 30 m resolution LAI maps with fine spatial and temporal consistencies, which can be expected to meet the needs of analysis on crop dynamic changes and yield mapping in fragmented and highly heterogeneous areas.

Estimating Rainfall Interception of Vegetation Canopy from MODIS Imageries in Southern China

The interception of rainfall by vegetation canopies plays an important role in the hydrologic process of ecosystems. Most estimates of canopy rainfall interception in present studies are mainly through field observations at the plot region. However, it is difficult, yet important, to map the regional rainfall interception by vegetation canopy at a larger scale, especially in the southern rainy areas of China. To obtain a better understanding of the spatiotemporal variation of vegetation canopy rainfall interception with regard to the basin scale in this region, we extended a rainfall interception model by combining the observed rainfall data and moderate resolution imaging spectroradiometer leaf area index (MODIS_LAI) data to quantitatively estimate the vegetation canopy rainfall interception rate (CRIR) at small/medium basin scales in Guangdong Province, which is undergoing large changes in vegetation cover due to rapid urban expansion in the area. The results showed that the CRIR in Guangdong declined continuously during 2004–2012, but increased slightly in 2016, and the spatial variability of CRIR showed a diminishing yearly trend. The CRIR also exhibited a distinctive spatial pattern, with a higher rate to the east and west of the mountainous areas and a lower rate in the central mountainous and coastal areas. This pattern was more closely related to the spatial variation of the LAI than that of rainfall due to frequent extreme rainfall events saturating vegetation leaves. Further analysis demonstrated that forest coverage, instead of background climate, has a certain impact on the canopy rainfall interception, especially the proportion of broad-leaved forests in the basin, but more in-depth study is warranted in the future. In conclusion, the results of this study provide insights into the spatiotemporal variation of canopy rainfall interception at the basin scale of the Guangdong Province, and suggest that forest cover should be increased by adjusting the species composition to increase the proportion of native broad-leaved species based on the local condition within the basin. In addition, these results would be helpful in accurately assessing the impacts of forest ecosystems on regional water cycling, and provide scientific and practical implications for water resources management.

Red-Edge Band Vegetation Indices for Leaf Area Index Estimation From Sentinel-2/MSI Imagery

The estimation of leaf area index (LAI) from optical remotely sensed data based on vegetation indices (VIs) is a quick and practical approach to acquire LAI over vast areas. Reflectance in the red-edge bands is sensitive to vegetation status, and its information is thought to be useful in agricultural applications. Based on three red-edge band observations (represented as RE1, RE2, and RE3 for bands 5–7) from the Multispectral Instrument (MSI) onboard the Sentinel-2 satellite, this article aims to investigate the feasibility and performance of using red-edge bands for LAI estimates with the VI method and ground-measured LAI data sets. Sensitivity analysis from PROSAIL simulations revealed that RE1 is mainly affected by the influence of the leaf chlorophyll content, and this uncertainty should not be ignored during LAI estimation. For the normalized difference vegetation index (NDVI), modified simple ratio (MSR), chlorophyll index (CI), and wide dynamic range vegetation index (WDRVI), the optimal combination of Sentinel-2 bands for LAI estimation was RE2 and RE3, with a minimum root-mean-square error (RMSE) of 0.75. Four 3-band red-edge VIs were proposed to exploit the full content of the red-edge bands of Sentinel-2, and their performance in LAI estimation improved slightly. However, both 2-band red-edge VIs and 3-band red-edge VIs remained slightly saturated at high LAI levels; therefore, a segmental estimation with a threshold was suggested for large LAIs. The results indicate that the optimal 2-band red-edge VIs and proposed 3-band red-edge VIs are effective tools for crop LAI estimation in multiple-growth stages with Sentinel-2 MSI images.

Toward High-Resolution Soil Moisture Monitoring by Combining Active-Passive Microwave and Optical Vegetation Remote Sensing Products with Land Surface Model.

The assimilation of radiometer and synthetic aperture radar (SAR) data is a promising recent technique to downscale soil moisture products, yet it requires land surface parameters and meteorological forcing data at a high spatial resolution. In this study, we propose a new downscaling approach, named integrated passive and active downscaling (I-PAD), to achieve high spatial and temporal resolution soil moisture datasets over regions without detailed soil data. The Advanced Microwave Scanning Radiometer (AMSR-E) and Phased Array-type L-band SAR (PALSAR) data are combined through a dual-pass land data assimilation system to obtain soil moisture at 1 km resolution. In the first step, fine resolution model parameters are optimized based on fine resolution PALSAR soil moisture and moderate-resolution imaging spectroradiometer (MODIS) leaf area index data, and coarse resolution AMSR-E brightness temperature data. Then, the 25 km AMSR-E observations are assimilated into a land surface model at 1 km resolution with a simple but computationally low-cost algorithm that considers the spatial resolution difference. Precipitation data are used as the only inputs from ground measurements. The evaluations at the two lightly vegetated sites in Mongolia and the Little Washita basin show that the time series of soil moisture are improved at most of the observation by the assimilation scheme. The analyses reveal that I-PAD can capture overall spatial trends of soil moisture within the coarse resolution radiometer footprints, demonstrating the potential of the algorithm to be applied over data-sparse regions. The capability and limitation are discussed based on the simple optimization and assimilation schemes used in the algorithm.

Agricultural Policy Environmental eXtender (APEX) Simulation of Spring Peanut Management in the North China Plain

Spring peanut is a valuable alternative crop to mitigate water scarcity caused by excessive water use in conventional cropping systems in the North China Plain (NCP). In the present study, we evaluated the capability of the Agricultural Policy Environmental eXtender (APEX) model to predict spring peanut response to sowing dates and seeding rates in order to optimize sowing dates, seeding rates, and irrigation regimes. Data used for calibration and validation of the model included leaf area index (LAI), aboveground biomass (ABIOM), and pod yield data collected from a field experiment of nine sowing dates and seeding rate combinations conducted from 2017 to 2018. The calibrated model was then used to simulate peanut yield responses to extended sowing dates (5 April to 4 June with a 5-day interval) and seeding rates (15 plants m−2 to 50 plants m−2 with a 5 plants m−2 interval) using 38 years of weather data as well as yield, evapotranspiration (ET), and water stress days under different irrigation regimes (rainfed, one irrigation before planting (60 mm) or at flowering (60 mm), and two irrigation with one time before planting and one time at flowering (60 mm each time) or at pod set (60 mm each time)). Results show that the model satisfactorily simulates pod yield of peanut based on R2 = 0.70, index of agreement (d value) being 0.80 and percent bias (PBIAS) values ≤4%. Moreover, the model performed reasonably well in predicting the emergence, LAI and ABIOM, with a R2 = 0.86, d = 0.95 and PBIAS = 8% for LAI and R2 = 0.90, d = 0.97 and PBIAS = 1% for ABIOM, respectively. Simulation results indicate that the best combination of sowing dates and seeding rates is a density of 35–40 plants m−2 and dates during early-May to mid-May due to the influence of local climate and canopy structure to the growth and yield of peanut. Under the optimal sowing date and plant density, an irrigation depth of 60 mm during flowering gave a pod yield (5.6 t ha−1) and ET (464 mm), which resulted in the highest water use efficiency (12.1 kg ha−1 mm−1). The APEX model is capable of assessing the effects of management practices on the growth and yield of peanut. Sowing 35–40 plants m−2 during early-May to mid-May with 60 mm irrigation depth is the recommended agronomic practice for peanut production in the water-constrained NCP.

Assimilating spatiotemporal MODIS LAI data with a particle filter algorithm for improving carbon cycle simulations for bamboo forest ecosystems

Bamboo forests are an important part of the forest ecosystem, which has strong carbon sequestration potential and plays an important role in the global carbon cycle. As a key parameter for simulating the carbon cycle using forest ecosystem models, the quality of leaf area index (LAI) data has a direct influence on the accuracy of modelling results. Here, we used the particle filter (PF) algorithm and PROSAIL model to assimilate MODIS LAI products, which were then used to drive a boreal ecosystem productivity simulator model to simulate the bamboo forest carbon cycle. The results showed that the relationship between the assimilated and observed LAI values was very significant, with an R2 of 0.95 and an RMSE of 0.28, greatly improving the precision of MODIS LAI products. The R2 values for the gross primary productivity (GPP), net ecosystem exchange (NEE), and total ecosystem respiration (TER) simulated by the assimilated LAI values and observed carbon fluxes were 0.65, 0.45 and 0.70, respectively, and the RMSE values were 1.10 g C m−2 day−1, 1.00 g C m−2 day−1 and 0.35 g C m−2 day−1, respectively. Compared with the results of the carbon cycle simulated by non-assimilated LAI, the R2 values of the GPP, NEE and TER values that were simulated by assimilated LAI increased by 27.5%, 45.2% and 6.1%, and the RMSE values decreased by 29.9%, 23.7% and 22.2%, respectively. Therefore, coupling the PF and PROSAIL models can greatly improve the simulation precision for the large-scale bamboo forest carbon cycle. This study laid the foundation for simulating the carbon cycle over a large-scale bamboo forest based on low-resolution data in the future.

Assimilation of Remotely Sensed LAI Into CLM4CN Using DART

AbstractPlant leaves play an important role in water, carbon, and energy exchanges between terrestrial ecosystems and atmosphere. Assimilating remotely sensed leaf area index (LAI) into land surface models is a promising approach to improve our understanding of those processes. Toward this goal, this study uses the Community Land Model with carbon and nitrogen components (CLM4CN) coupled with the Data Assimilation Research Testbed (DART). Global Land Surface Satellite (GLASS) LAI data are assimilated via the Ensemble Adjustment Kalman Filter. A random 40‐member atmospheric forcing ensemble is used to drive the CLM4CN to provide background error covariance. The results show that assimilating GLASS LAI and updating both LAI and leaf C/N is an effective method to provide a high‐accuracy estimate of LAI. The simulations always systematically overestimate LAI, especially in low‐latitude regions, with the largest bias up to 5 m2/m2, which are effectively corrected in the analyzed LAI, with the bias reduced to ±1 m2/m2. Significantly improved regions are located in central Africa, Amazonia, southern Eurasia, northeastern China, and western Europe, where evergreen/deciduous forests and mixed forests are dominant. Except for the temperate zone in the Southern Hemisphere, the analyzed LAI can well represent seasonal variations. The most pronounced assimilation impact in low‐latitude regions is attributed to large initial forecast error covariance and sufficient background errors. The MOD 16 evapotranspiration estimates and upscaled gross primary production have been used to evaluate the assimilation impact, which highlight neutral to highly positive improvement.

A Density-Based Approach for Leaf Area Index Assessment in a Complex Forest Environment Using a Terrestrial Laser Scanner

Forests are an important part natural ecosystems, by for example providing food, fiber, habitat, and biodiversity, all of which contribute to stable natural systems. Assessing and modeling the structure and characteristics of forests, e.g., Leaf Area Index (LAI), volume, biomass, etc., can lead to a better understanding and management of these resources. In recent years, Terrestrial Laser Scanning (TLS) has been recognized as a tool that addresses many of the limitations of manual and traditional forest data collection methods. In this study, we propose a density-based approach for estimating the LAI in a structurally-complex forest environment, which contains variable and diverse structural attributes, e.g., non-circular stem forms, dense canopy and below-canopy vegetation cover, and a diverse species composition. In addition, 242 TLS scans were collected using a portable low-cost scanner, the Compact Biomass Lidar (CBL), in the Hawaii Volcanoes National Park (HAVO), Hawaii Island, USA. LAI also was measured for 242 plots in the site, using an AccuPAR LP-80 ceptometer. The first step after cleaning the point cloud involved detecting the higher forest canopy in the light detection and ranging (lidar) point clouds, using normal change rate assessment. We then estimated Leaf Area Density (LAD), using a voxel-based approach, and divided the canopy point cloud into five layers in the Z (vertical) direction. These five layers subsequently were divided into voxels in the X direction, where the size of these voxels were obtained based on inter-quartile analysis and the number of points in each voxel. We hypothesized that the intensity returned to the lidar system from woody materials, like branches, would be higher than from leaves, due to the liquid water absorption feature of the leaves and higher reflectance for woody material at the 905 nm laser wavelength. We also differentiated between foliar and woody materials using edge detection in the images from projected point clouds and evaluated the density of these regions to support our hypothesis. Density of points, or the number of points divided by the volume of a grid, in a 3D grid size of 0.1 m, was calculated for each of the voxels. The grid size was determined by investigating the size of the branches in the lower portion of the canopy. Subsequently, we fitted a Kernel Density Estimator (KDE) to these values, with the threshold set based on half of the area under the curve in each of the density distributions. All the grids with a density below the threshold were labeled as leaves, while those grids above the threshold were identified as non-leaves. Finally, we modeled LAI using the point densities derived from the TLS point clouds and the listed analysis steps. This model resulted in an R 2 value of 0.88. We also estimated the LAI directly from lidar data using the point densities and calculating LAD, which is defined as the total one-sided leaf area per unit volume. LAI can be obtained as the sum of the LAD values in all the voxels. The accuracy of LAI estimation was 90%, with an RMSE value of 0.31, and an average overestimation of 9 % in TLS-derived LAI, when compared to field-measured LAI. Algorithm performance mainly was affected by the vegetation density and complexity of the canopy structures. It is worth noting that, since the LAI values cannot be considered spatially independent throughout all the plots in this site, we performed semivariogram analysis on the field-measured LAI data. This analysis showed that the LAI values can be assumed to be independent in plots that are at least 30 m apart. As a result, we divided the data into six subsets in which the plots were 30 m spaced. The R 2 values for these subsets, based on modeling of the field-measured LAI using leaf point density values, ranged between 0.84–0.96. The results bode well for using this method for efficient, automatic, and accurate/precise estimation of LAI values in complex forest environments, using a low-cost, rapid-scan TLS.

Watershed-scale-based forecast method for leaf area index data based on the integration of time series MODIS products and meteorological data

The leaf area index (LAI) is an important parameter used to characterize land vegetation growth patterns, and LAI forecasts are highly beneficial to natural resource management decisions. However, the LAI can be affected by biophysical, climate, and anthropogenic factors, thus, its future watershed-scale patterns remain difficult to assess. Herein, an object-oriented watershed-scale LAI forecasting (OWLF) method is proposed to solve the forecasting difficulties encountered when using only past meteorological data. The model was constructed with filtered multi-year MODIS LAI time series and climate index variables. Then, the dependence of the LAI on four mean monthly climate index variables was evaluated, namely, latent evapotranspiration, rainfall capacity, average air temperature, and sunlight hours, in the Poyang Lake basin. By combining the climate indexes derived from the past 1–3 months, the model forecasts the subsequent months LAI for different land cover types. The LAI forecasting results derived with a 13-year time series (2000–2012) suggest that the OWLF method can effectively recognize the expected spatial patterns, and the data agreed reasonably well with LAI dynamics and phenological periods. This work offers a promising way to exploit combined satellite and climate index data in novel and more accurate watershed-scale forecasting studies.

Comparison of different sequential assimilation algorithms for satellite-derived leaf area index using the Data Assimilation Research Testbed (version Lanai)

Abstract. The leaf area index (LAI) is a crucial parameter for understanding the exchanges of mass and energy between terrestrial ecosystems and the atmosphere. In this study, the Data Assimilation Research Testbed (DART) has been successfully coupled to the Community Land Model with explicit carbon and nitrogen components (CLM4CN) by assimilating Global Land Surface Satellite (GLASS) LAI data. Within this framework, four sequential assimilation algorithms, including the kernel filter (KF), the ensemble Kalman filter (EnKF), the ensemble adjust Kalman filter (EAKF), and the particle filter (PF), are thoroughly analyzed and compared. The results show that assimilating GLASS LAI into the CLM4CN is an effective method for improving model performance. In detail, the assimilation accuracies of the EnKF and EAKF algorithms are better than those of the KF and PF algorithm. From the perspective of the average and RMSD, the PF algorithm performs worse than the EAKF and EnKF algorithms because of the gradually reduced acceptance of observations with assimilation steps. In other words, the contribution of the observations to the posterior probability during the assimilation process is reduced. The EAKF algorithm is the best method because the matrix is adjusted at each time step during the assimilation procedure. If all the observations are accepted, the analyzed LAI seem to be better than that when some observations are rejected, especially in low-latitude regions.

Joint Assimilation of Leaf Area Index and Soil Moisture from Sentinel-1 and Sentinel-2 Data into the WOFOST Model for Winter Wheat Yield Estimation.

It is well known that timely crop growth monitoring and accurate crop yield estimation at a fine scale is of vital importance for agricultural monitoring and crop management. Crop growth models have been widely used for crop growth process description and yield prediction. In particular, the accurate simulation of important state variables, such as leaf area index (LAI) and root zone soil moisture (SM), is of great importance for yield estimation. Data assimilation is a useful tool that combines a crop model and external observations (often derived from remote sensing data) to improve the simulated crop state variables and consequently model outputs like crop total biomass, water use and grain yield. In spite of its effectiveness, applying data assimilation for monitoring crop growth at the regional scale in China remains challenging, due to the lack of high spatiotemporal resolution satellite data that can match the small field sizes which are typical for agriculture in China. With the accessibility of freely available images acquired by Sentinel satellites, it becomes possible to acquire data at high spatiotemporal resolution (10–30 m, 5–6 days), which offers attractive opportunities to characterize crop growth. In this study, we assimilated remotely sensed LAI and SM into the Word Food Studies (WOFOST) model to estimate winter wheat yield using an ensemble Kalman filter (EnKF) algorithm. The LAI was calculated from Sentinel-2 using a lookup table method, and the SM was calculated from Sentinel-1 and Sentinel-2 based on a change detection approach. Through validation with field data, the inverse error was 10% and 35% for LAI and SM, respectively. The open-loop wheat yield estimation, independent assimilations of LAI and SM, and a joint assimilation of LAI + SM were tested and validated using field measurement observation in the city of Hengshui, China, during the 2016–2017 winter wheat growing season. The results indicated that the accuracy of wheat yield simulated by WOFOST was significantly improved after joint assimilation at the field scale. Compared to the open-loop estimation, the yield root mean square error (RMSE) with field observations was decreased by 69 kg/ha for the LAI assimilation, 39 kg/ha for the SM assimilation and 167 kg/ha for the joint LAI + SM assimilation. Yield coefficients of determination (R2) of 0.41, 0.65, 0.50, and 0.76 and mean relative errors (MRE) of 4.87%, 4.32%, 4.45% and 3.17% were obtained for open-loop, LAI assimilation alone, SM assimilation alone and joint LAI + SM assimilation, respectively. The results suggest that LAI was the first-choice variable for crop data assimilation over SM, and when both LAI and SM satellite data are available, the joint data assimilation has a better performance because LAI and SM have interacting effects. Hence, joint assimilation of LAI and SM from Sentinel-1 and Sentinel-2 at a 20 m resolution into the WOFOST provides a robust method to improve crop yield estimations. However, there is still bias between the key soil moisture in the root zone and the Sentinel-1 C band retrieved SM, especially when the vegetation cover is high. By active and passive microwave data fusion, it may be possible to offer a higher accuracy SM for crop yield prediction.

Erratum: Novelli, F., et al. Assimilation of Sentinel-2 Leaf Area Index Data into a Physically-Based Crop Growth Model for Yield Estimation. Agronomy 2019, 9, 255

The authors wish to correct the following erratum in this paper [...]

Response of surface air temperature to the change of leaf area index in the source region of the Yellow River by the WRF model

Leaf area index (LAI) is a crucial land-atmosphere exchange parameter. LAI change can cause a variation of other land surface parameter. In this research, three experiments were conducted to investigate the impact of LAI and albedo change on surface air temperature in the source region of the Yellow River (SRYR) using the Weather Research and Forecasting (WRF) model. Three experiments used the same settings, initial and boundary conditions except for the LAI and albedo data. The control simulation (CTL) used the WRF model climatological LAI data, the second simulation (LAI) used the Moderate Resolution Imaging Spectroradiometer (MODIS) LAI data, and the third simulation (LAIALB) used the MODIS LAI and MODIS albedo. The results show MODIS LAI is greater by 34.5% than WRF climatological LAI, and the MODIS albedo is lower by 24.3% than WRF climatological albedo over the whole growing season of 2006 in SRYR. All the experiments can simulate the surface air temperature (Ta) spatial distribution characteristics, but underestimate the values of 1.3 °C in CTL experiment and 0.6 °C in LAIALB experiment in SRYR. The simulated Ta by LAI experiment is lower 0.1 °C than by the CTL experiment, but the simulated Ta by the LAIALB experiment is obviously higher 0.6 °C than the CTL experiment. The LAI experiment shows a cooling effect because the higher MODIS LAI decreases the canopy resistance, which induces a positive average 2.1 Wm−2 latent heat flux (LH) and a negative average − 2.1 Wm−2 sensible heat flux (Hs). The LAIALB experiment presents a warming effect because of low MODIS albedo comparing WRF albedo, which changes the radiation components and results in an obvious negative − 14.0 Wm−2 upward short wave radiation, a positive 11.7 Wm−2 net radiation, and a positive 10.9 Wm−2 heat flux. In fact, more precipitation produces more snow and high surface albedo in the WRF model, which results in a cold temperature bias in SRYR.

Development and evaluation of a simple hydrologically based model for terrestrial evapotranspiration simulations

Following the scheme of the groundwater-soil-plant-atmosphere continuum (GSPAC), a simple process-based model (SiTH; Simple Terrestrial Hydrosphere) was developed to estimate the dynamics of terrestrial evapotranspiration (ET) at the daily step. The input data include net radiation, air temperature, precipitation, leaf area index, vegetation type and soil data, most of which are readily available. Locally, the model performed well in simulating the dynamics of ET and soil moisture over selected FLUXNET sites. Globally, the daily 0.25° ET and groundwater table depth estimations in year 2005 were determined using available gridded datasets. The spatial pattern was reasonable and the range of values corresponded well with other global ET and groundwater table depth products. In future studies, we will produce a long-term daily 0.25° global ET and groundwater table depth products that spans from 1984 to present by using different input datasets.