Phenotyping provides important support for corn breeding. Unfortunately, the rapid detection of phenotypes has been the major limiting factor in estimating and predicting the outcomes of breeding programs. This study was focused on the potential of phenotyping to support corn breeding using unmanned aerial vehicle (UAV) images, aiming at mining and deepening UAV techniques for comparing phenotypes and screening new corn varieties. Two geometric traits (plant height, canopy leaf area index (LAI)) and one lodging resistance trait (lodging area) were estimated in this study. It was found that stereoscopic and photogrammetric methods were promising ways to calculate a digital surface model (DSM) for estimating corn plant height from UAV images, with R² = 0.7833 (p < 0.001) and a root mean square error (RMSE) = 0.1677. In addition to a height estimation, the height variation was analyzed for depicting and validating the corn canopy uniformity stability for different varieties. For the lodging area estimation, the normalized DSM (nDSM) method was more promising than the gray-level co-occurrence matrix (GLCM) textural features method. The estimation error using the nDSM ranged from 0.8% to 5.3%, and the estimation error using the GLCM ranged from 10.0% to 16.2%. Associations between the height estimation and lodging area estimation were done to find the corn varieties with optimal plant heights and lodging resistance. For the LAI estimation, the physical radiative transfer PROSAIL model offered both an accurate and robust estimation performance both at the middle (R² = 0.7490, RMSE = 0.3443) and later growing stages (R² = 0.7450, RMSE = 0.3154). What was more exciting was that the estimated sequential time series LAIs revealed a corn variety with poor resistance to lodging in a study area of Baogaofeng Farm. Overall, UAVs appear to provide a promising method to support phenotyping for crop breeding, and the phenotyping of corn breeding in this study validated this application.

Data assimilation has been demonstrated as the potential crop yield estimation approach. Accurate quantification of model and observation errors is the key to determining the success of a data assimilation system. However, the crop growth model error is not fully taken into account in most of the previous studies. The objective of this study is to better quantify the model uncertainty in the data assimilation system. Firstly, we calibrated a crop growth model and inferred its posterior uncertainty based on the Global LAnd Surface Satellite (GLASS) 250-m LAI product, regional statistical data, station observations, and field measurements with a Markov chain Monte Carlo (MCMC) method. Secondly, the model posterior uncertainty was used in the Ensemble Kalman Filter (EnKF) algorithm to better characterize the ensemble distribution of model errors. Our results indicated the proposed Bayesian posterior-based EnKF can improve the accuracy of winter wheat yield estimation at both the point scale (the coefficient of determination R ² value increasing from 0.06 to 0.41, the mean absolute percentage error MAPE value decreasing from 12.65% to 7.82%, and the root mean square error RMSE value decreasing from 987 to 688 kg∙ha ^-1 ) and the regional scale (R ² value from 0.30 to 0.57, MAPE value from 19.67% to 10.13%, and RMSE value from 1275 to 695 kg∙ha ^-1 ) compared with the open-loop estimation. Our analysis also indicated that the Bayesian posterior-based EnKF can perform better compared to the standard Gaussian perturbation-based EnKF. The proposed framework provides an important reference for crop yield estimation at the regional scale in similar agricultural landscapes worldwide.

Light pollution, a phenomenon in which artificial nighttime light (NTL) changes the form of brightness and darkness in natural areas such as protected areas (PAs), has become a global concern due to its threat to global biodiversity. With ongoing global urbanization and climate change, the light pollution status in global PAs deserves attention for mitigation and adaptation. In this study, we developed a framework to evaluate the light pollution status in global PAs, using the global NTL time series data. First, we classified global PAs (30,624) into three pollution categories: non-polluted (5974), continuously polluted (8141), and discontinuously polluted (16,509), according to the time of occurrence of lit pixels in/around PAs from 1992 to 2018. Then, we explored the NTL intensity (e.g., digital numbers) and its trend in those polluted PAs and identified those hotspots of PAs at the global scale with consideration of global urbanization. Our study shows that global light pollution is mainly distributed within the range of 30°N and 60°N, including Europe, north America, and East Asia. Although the temporal trend of NTL intensity in global PAs is increasing, Japan and the United States of America (USA) have opposite trends due to the implementation of well-planned ecological conservation policies and declining population growth. For most polluted PAs, the lit pixels are close to their boundaries (i.e., less than 10 km), and the NTL in/around these lit areas has become stronger over the past decades. The identified hotspots of PAs (e.g., Europe, the USA, and East Asia) help support decisions on global biodiversity conservation, particularly with global urbanization and climate change.

In the Loess Plateau (LP) of China, the vegetation degradation and soil erosion problems have been shown to be curbed after the implementation of the Grain for Green program. In this study, the LP is divided into the northwestern semi-arid area and the southeastern semi-humid area using the 400 mm isohyet. The spatial–temporal evolution of the vegetation NDVI during 2000–2015 are analyzed, and the driving forces (including factors of climate, environment, and human activities) of the evolution are quantitatively identified using the geographical detector model (GDM). The results showed that the annual mean NDVI in the entire LP was 0.529, and it decreased from the semi-humid area (0.619) to the semi-arid area (0.346). The mean value of the coefficient of variation of the NDVI was 0.1406, and it increased from the semi-humid area (0.1165) to the semi-arid area (0.1926). The annual NDVI growth rate in the entire LP was 0.0079, with the NDVI growing faster in the semi-humid area (0.0093) than in the semi-arid area (0.0049). The largest increments of the NDVI were from grassland, farmland, and woodland. The GDM results revealed that changes in the spatial distribution of the NDVI could be primarily explained by the climatic and environmental factors in the semi-arid area, such as precipitation, soil type, and vegetation type, while the changes were mainly explained by the anthropogenic factors in the semi-humid area, such as the GDP density, land-use type, and population density. The interactive analysis showed that interactions between factors strengthened the impacts on the vegetation change compared with an individual factor. Furthermore, the ranges/types of factors suitable for vegetation growth were determined. The conclusions of this study have important implications for the formulation and implementation of ecological conservation and restoration strategies in different regions of the LP.

The tassel development status and its branch number in maize flowering stage are the key phenotypic traits to determine the growth process, pollen quantity of different maize varieties, and detasseling arrangement for seed maize production fields. Rapid and accurate detection of tassels is of great significance for maize breeding and seed production. However, due to the complex planting environment in the field, such as unsynchronized growth stage and tassels vary in size and shape caused by varieties, the detection of maize tassel remains challenging problem, and the existing methods also cannot distinguish the early tassels. In this study, based on the time series unmanned aerial vehicle (UAV) RGB images with maize flowering stage, we proposed an algorithm for automatic detection of maize tassels which is suitable for complex scenes by using random forest (RF) and VGG16. First, the RF was used to segment UAV images into tassel regions and non-tassel regions, and then extracted the potential tassel region proposals by morphological method; afterwards, false positives were removed through VGG16 network with the ratio of training set to validation set was 7:3. To demonstrate the performance of the proposed method, 50 plots were selected from UAV images randomly. The precision, recall rate and F1-score were 0.904, 0.979 and 0.94 respectively; 50 plots were divided into early, middle and late tasseling stages according to the proportion of tasseling plants and the morphology of tassels. The result of tassels detection was late tasseling stage > middle tasseling stage > early tasseling stage, and the corresponding F1-score were 0.962, 0.914 and 0.863, respectively. It was found that the model error mainly comes from the recognition of leaves vein and reflective leaves as tassels. Finally, to show the morphological characteristics of tassel directly, we proposed an endpoint detection method based on the tassel skeleton, and further extracted the tassel branch number. The method proposed in this paper can well detect tassels of different development stages, and support large scale tassels detection and branch number extraction.