Nonlinear Hyperspectral Mixture Analysis for tree cover estimates in orchards

ABSTRACTThe normalized difference vegetation index (NDVI) is a commonly used index for monitoring crop growth status. Previous studies have shown that the leaf area index (LAI) estimation based on NDVI is limited by saturation that occurs under conditions of relatively dense canopies (LAI > 2 m2 m–2). To reduce the saturation effect, we suggested new spectral indices through the spectral indices approach. The results suggested that the two-band normalized difference spectral index (NDSI = ((ρ940 – ρ730) /(ρ940 + ρ730))) resulted from the two-band spectral indices approach and the three-band modified normalized difference spectral index (mNDSI = ((ρ940 – 0.8 × ρ950) – ρ730) /((ρ940 – 0.8 × ρ950) + ρ730)) resulted from the three-band spectral indices approach, and they were able to mitigate saturation and improve the LAI prediction with a determination coefficient (R2) of 0.77 and 0.78, respectively. In the validation based on data from independent experiments, these new indices exhibited an accuracy with relative root mean square error (RRMSE) lower than 23.38% and bias higher than –0.40. These accuracies were significantly higher than those obtained with some existing indices with good performance in LAI estimation, such as the enhanced vegetation index (EVI) (RRMSE = 30.19%, bias = –0.34) and the modified triangular vegetation index 2 (MTVI2) (RRMSE = 29.30%, bias = –0.28), and the indices with the ability to mitigate the saturation, such as the wide dynamic range vegetation index (WDRVI) (RRMSE = 31.37%, bias = –0.54), the red-edge wide dynamic range vegetation index (red-edge WDRVI) (RRMSE = 26.34%, bias = –0.54), and the normalized difference red-edge index (NDRE) (RRMSE = 28.41%, bias = –0.56). Additionally, these new indices were more sensitive under moderate to high LAI conditions (between 2 and 8 m2 m–2). Between these two new developed spectral indices, there was no significant difference in the accuracy and sensitivity assessments. Considering the index structure and convenience in application, we demonstrated that the two-band spectral index NDSI((ρ940 – ρ730) /(ρ940 + ρ730)) is efficient in mitigating saturation and has considerable potential for estimating the LAI of canopies throughout the entire growing season of wheat (Triticum aestivum L.), whereas the three-band spectral index contributes lesser in the saturation mitigation provided the red-edge band has been contained.

Abstract. This paper examines various geometric shapes to determine the most suitable one for calculating tree crown area in field measurements. The study was conducted in an Eldarica pine plantation forest, which was digitally mapped using RGB images captured by a Phantom 4 Pro drone. Tree crowns were manually digitized (MD) from these images to serve as reference data. Field measurements, including the large and small crown diameters, were collected to evaluate crown areas derived from different geometric shapes. The geometric shapes considered were: Oval with Both Diameters (OBD), Circle with Small Diameter (CSD), Circle with Large Diameter (CLD), and Circle with Mean Diameters (CMD). Three analyses were performed to assess the results: correlation analysis (R2), relative root mean square error (RRMSE), and shape analysis, which included overestimation (OverID) and underestimation (UnderID) indices. The results revealed that the choice of geometric shape significantly impacts the accuracy of crown area calculations. The OBD model based on the outer boundary diameter yielded the best results with RRMSE = 0.29, R2 = 0.84, OverID = 0.18, and UnderID = 0.23, followed closely by the CMD method. In contrast, the CSD and CLD models performed less effectively, with RRMSE = 0.52, R2 = 0.42, OverID = 0.11, UnderID = 0.35 (CSD), and RRMSE = 0.59, R2 = 0.37, OverID = 0.46, UnderID = 0.22 (CLD). These differences in performance are likely due to the inclusion of empty spaces within the crown area in some models. However, the findings of this study are not universally applicable to all tree crown area calculations. The geometric shape used for crown area estimation must align with the structural characteristics of the forest and the specific geometry of the tree species under consideration.

The error and bias of supplementing a short, arid climate, rainfall record with regional vs. global frequency analysis

Abstract. A dual-lidar system comprising two WindCube v2.1 Doppler lidars – one oriented horizontally at a 45° angle relative to the other – was deployed to estimate atmospheric turbulence using the variance method. This approach derives second-order velocity statistics directly from line-of-sight (LOS) velocity variances, enabling the reconstruction of the three-dimensional velocity variance components without retrieving instantaneous wind vectors. Its performance is evaluated against the traditional method commonly used in the wind energy sector, which reconstructs instantaneous velocity components from single-lidar LOS measurements prior to computing turbulence statistics. Both methods are assessed at a single measurement height using a 30 d collocated dataset, with sonic anemometer measurements as reference and classification by atmospheric stability. The analysis focuses on turbulence intensity (TI), and performance is quantified using the mean relative bias error (MRBE) and relative root mean square error (RRMSE), following Det Norske Veritas (DNV) load-validation acceptance criteria (MRBE within ±5 % and RRMSE ≤ 15 %). Results show that the variance method demonstrates improved agreement with the reference measurements across nearly all wind-speed bins. For along-wind TI, MRBE values range between 0 % and −5 %, with RRMSE remaining below 15 % over the full wind-speed range, fully satisfying the DNV criteria. For cross-wind TI, the MRBE criterion is met, while the RRMSE threshold is exceeded. In contrast, the traditional method meets both DNV criteria for along-wind TI only at wind speeds above approximately 7 m s−1 and fails to satisfy either criterion for cross-wind TI across the investigated range. Overall, the results demonstrate that the variance method provides more robust and DNV-compliant TI estimates than the traditional reconstruction approach, particularly for along-wind turbulence relevant to wind turbine load validation.

China has been facing serious land degradation and desertification in its north and northwest arid and semi-arid areas. Monitoring the dynamics of percentage vegetation cover (PVC) using remote sensing imagery in these areas has become critical. However, because these areas are large, remote, and sparsely populated, and also because of the existence of mixed pixels, there have been no accurate and cost-effective methods available for this purpose. Spectral unmixing methods are a good alternative as they do not need field data and are low cost. However, traditional linear spectral unmixing (LSU) methods lack the ability to capture the characteristics of spectral reflectance and scattering from endmembers and their interactions within mixed pixels. Moreover, existing nonlinear spectral unmixing methods, such as random forest (RF) and radial basis function neural network (RBFNN), are often costly because they require field measurements of PVC from a large number of training samples. In this study, a cost-effective approach to mapping PVC in arid and semi-arid areas was proposed. A method for selection and purification of endmembers mainly based on Landsat imagery was first presented. A probability-based spectral unmixing analysis (PBSUA) and a probability-based optimized k nearest-neighbors (PBOkNN) approach were then developed to improve the mapping of PVC in Duolun County in Inner Mongolia, China, using Landsat 8 images and field data from 920 sample plots. The proposed PBSUA and PBOkNN methods were further validated in terms of accuracy and cost-effectiveness by comparison with two LSU methods, with and without purification of endmembers, and two nonlinear approaches, RF and RBFNN. The cost-effectiveness was defined as the reciprocal of cost timing relative root mean square error (RRMSE). The results showed that (1) Probability-based spectral unmixing analysis (PBSUA) was most cost-effective and increased the cost-effectiveness by 29.3% 29.3%, 33.5%, 50.8%, and 53.0% compared with two LSU methods, PBOkNN, RF, and RBFNN, respectively; (2) PBSUA, RF, and RBFNN gave RRMSE values of 22.9%, 21.8%, and 22.8%, respectively, which were not significantly different from each other at the significance level of 0.05. Compatibly, PBOkNN and LSU methods with and without purification of endmembers resulted in significantly greater RRMSE values of 27.5%, 32.4%, and 43.3%, respectively; (3) the average estimates of the sample plots and predicted maps from PBSUA, PBOkNN, RF, and RBFNN fell in the confidence interval of the test plot data, but those from two LSU methods did not, although the LSU with purification of endmembers improved the PVC estimation accuracy by 25.2% compared with the LSU without purification of endmembers. Thus, this study indicated that the proposed PBSUA had great potential for cost-effectively mapping PVC in arid and semi-arid areas.

Soil pH is a key factor that controls soil nutrient availability, soil microbial activities, and crop growth and development. However, studies on the soil pH variations of cultivated lands in different horizons at the regional scale remain limited. In this work, 348 soil samples were collected from three soil horizons (A, B, and C) at 120 sites over the hilly region of Chongqing, southwestern China. Six topographic indicators, four climate parameters, and parent material were considered. Classification and regression trees (CARTs) were applied to investigate the relationships between soil pH and the variables in the A, B, and C horizons. Model performances were evaluated by root mean square error (RMSE), relative root mean square error (RRMSE), and coefficient of determination (R2). Results showed that soil pH increased obviously from the A to C horizons. Soil pH was predicted well by the forcing factors with the CART models in all horizons. RMSE, RRMSE, and R2 varied between 0.37 and 0.435, between 5.93 and 7.23%, and between 0.71 and 0.80, respectively. The relative importance of the studied variables to soil pH differed with the horizons. Annual temperature range (ATR), terrain wetness index (TWI), and Melton ruggedness number were critical factors that controlled soil pH variability in the A horizon. Parent material, precipitation of warmest quarter (PWQ), ATR, and TWI were important variables in the B horizon. Parent material, PWQ, ATR, and precipitation were key factors in the C horizon. The results are expected to provide valuable information for designing appropriate measurements for agricultural practices and preventing soil acidification.

Unsaturated zone N fate and transport were evaluated at four sites to identify the predominant pathways of N cycling: an almond [Prunus dulcis (Mill.) D.A. Webb] orchard and cornfield (Zea mays L.) in the lower Merced River study basin, California; and corn-soybean [Glycine max (L.) Merr.] rotations in study basins at Maple Creek, Nebraska, and at Morgan Creek, Maryland. We used inverse modeling with a new version of the Root Zone Water Quality Model (RZWQM2) to estimate soil hydraulic and nitrogen transformation parameters throughout the unsaturated zone; previous versions were limited to 3-m depth and relied on manual calibration. The overall goal of the modeling was to derive unsaturated zone N mass balances for the four sites. RZWQM2 showed promise for deeper simulation profiles. Relative root mean square error (RRMSE) values for predicted and observed nitrate concentrations in lysimeters were 0.40 and 0.52 for California (6.5 m depth) and Nebraska (10 m), respectively, and index of agreement (d) values were 0.60 and 0.71 (d varies between 0 and 1, with higher values indicating better agreement). For the shallow simulation profile (1 m) in Maryland, RRMSE and d for nitrate were 0.22 and 0.86, respectively. Except for Nebraska, predictions of average nitrate concentration at the bottom of the simulation profile agreed reasonably well with measured concentrations in monitoring wells. The largest additions of N were predicted to come from inorganic fertilizer (153-195 kg N ha(-1) yr(-1) in California) and N fixation (99 and 131 kg N ha(-1) yr(-1) in Maryland and Nebraska, respectively). Predicted N losses occurred primarily through plant uptake (144-237 kg N ha(-1) yr(-1)) and deep seepage out of the profile (56-102 kg N ha(-1) yr(-1)). Large reservoirs of organic N (up to 17,500 kg N ha(-1) m(-1) at Nebraska) were predicted to reside in the unsaturated zone, which has implications for potential future transfer of nitrate to groundwater.

Soil water content (SWC) is an important factor restricting crop growth and yield in cropland ecosystems. The observation and simulation of soil moisture contribute greatly to improving water-use efficiency and crop yield. This study was conducted at the Shandong Yucheng Agro-ecosystem National Observation and Research Station in the North China Plain. The study period was across the winter wheat (Triticum aestivum L.) growth stages from 2017 to 2019. A cosmic-ray neutron probe was used to monitor the continuous daily SWC. Furthermore, the crop leaf area index (LAI), yield, and aboveground biomass of winter wheat were determined. The root zone quality model 2 (RZWQM2) was used to simulate and validate the SWC, crop LAI, yield, and aboveground biomass. The results showed that the simulation errors of SWC were minute across the wheat growth stages and mature stages in 2017–2019. The root mean square error (RMSE) and relative root mean square error (RRMSE) of the SWC simulation at the jointing stage of winter wheat were 0.0296 and 0.1605 in 2017–2018, and 0.0265 and 0.1480 in 2018–2019, respectively. During the rain-affected days, the RMSE (0.0253) and RRMSE (0.0980) for 2017–2018 were significantly lower than those of 2018–2019 (0.0301 and 0.1458, respectively), indicating that rain events decreased the model accuracy in the dry years compared to the wet years. The simulated LAIs were significantly higher than the measured values. The simulated yield value of winter wheat was 5.61% lower and 3.92% higher than the measured yield in 2017–2018 and in 2018–2019, respectively. The simulated value of aboveground biomass was significantly (45.48%) lower than the measured value in 2017–2018. This study showed that, compared with the dry and cold wheat growth period of 2018–2019, the higher precipitation and temperature in 2017–2018 led to a poorer simulation of SWC and crop-growth components. This study indicated that annual abnormal rainfall and temperature had a significant influence on the simulation of SWC and wheat growth, especially under intensive climate-change stress conditions.

Reliable information concerning stand volume is fundamental to making strategic decisions in sustainable forest management. A variety of remotely sensed data and different inventory methods have been used for the estimation of forest biometric parameters. Particularly, airborne laser scanning (ALS) point clouds are widely used for the estimation of stand volume and forest biomass using an area-based approach (ABA) framework. This method relies on the reference measurements of field plots with the necessary prerequisite of a precise co-registration between ground reference plots and the corresponding ALS samples. In this research, the allometric area-based approach (AABA) is proposed in the context of stand volume estimation of Scots pine (Pinus sylvestris L.) stands. The proposed method does not require detailed information about the coordinates of the field plots. We applied Polish National Forest Inventory data from 9400 circular field plots (400 m2) to develop a plot level stand volume allometric model using two independent variables: top height (TH) and relative spacing index (RSI). The model was developed using the multiple linear regression method with a log–log transformation of variables. The hypothesis was that, the field measurements of TH and RSI could be replaced with corresponding ALS-derived metrics. It was assumed that TH could be represented by the maximum height of the ALS point cloud, while RSI can be calculated based on the number of tree crowns delineated within the ALS-derived canopy height model. Performance of the developed AABA model was compared with the semi-empirical ABASE (with two predictors: TH and RSI) and empirical ABAE (several point cloud metrics as predictors). The models were validated at the plot level using 315 forest management inventory plots (400 m2) and at the stand level using the complete field measurements from 42 Scots pine dominated forest stands in the Milicz forest district (Poland). The AABA model showed a comparable accuracy to the traditional ABA models with relatively high accuracy at the plot (relative root mean square error (RMSE) = 22.8 per cent; R2 = 0.63) and stand levels (RMSE = 17.8 per cent, R2 = 0.65). The proposed novel approach reduces time- and cost-consuming field work required for the classic ABA method, without a significant reduction in the accuracy of stand volume estimations. The AABA is potentially applicable in the context of forest management inventory without the necessity for field measurements at local scale. The transportability of the approach to other species and more complex stands needs to be explored in future studies.

In this study, we developed a model for photosynthetic production in double cropping rice by integrating the advantages in current crop models (including the models of canopy structure, canopy light distribution, canopy photosynthesis and dry matter production). The canopy light distribution and dry matter accumulation were preliminarily validated with independent field experiment datasets. The distribution of direct radiation both on a level surface and on the leaf surface within canopy, the canopy daily photosynthate and its characteristics with varying leaf area index for three typical plant types (erect both upper and lower, upper erect and lower flat, and flat both upper and lower) were quantitatively analyzed by the model. The results indicated that there was a good agreement between the simulated and observed values. The root mean square error (RMSE), relative root mean square error (RRMSE) and correlation coefficient (r) of prediction of canopy light distribution in double cropping rice were 12.01 J ·m-2·s-1, 8.2% and 0.9929, respectively. Meanwhile, the RMSE, RRMSE and r of prediction of dry matter accumulation were 0.83 t·hm-2, 14.6% and 0.9772, respectively. It was indicated that the model had a performance. The upper erect and lower flat plant type had highest canopy daily photosynthate due to higher incident sun light received on the leaf surface, leaf photosynthetic efficiency and leaf area index.

Chillers operation energy efficiency is the main factors affecting the efficiency of air conditioning systems. A chillers operation energy efficiency model and its method for a real chiller plant were proposed in this paper, whose parameters were optimized by GA optimization algorithm. Moreover, as R_RMSE (Relative Root Mean Square Error) was adopted as the evaluation of the prediction accuracy, the results showed that the prediction accuracy of SVR model based on GA optimization algorithm was better than other models, which average R_RMSE was 4.52%. The method proposed in this paper can predict the operation energy efficiency of the chiller accurately to provide basic for chiller energy efficiency analysis, fault detection and diagnosis, and optimizing control.

Enhancement of bioenergy estimations within forests using airborne laser scanning and multispectral line scanner data

For the planning and sustainable management of forest resources, well-managed plantations are of great significance to mitigate the decrease of forested areas. Monitoring these planted forests is essential for forest resource inventories. In this study, two ALOS PALSAR-2 quad-polarimetric synthetic aperture radar (SAR) images and ground measurements were employed to estimate growing stem volume (GSV) of Chinese fir plantations located in a hilly area of southern China. To investigate the relationships between forest GSV and polarization characteristics, single and fused variables were derived by the Yamaguchi decomposition and the saturation value of GSV was estimated using a semi-exponential empirical model as a base model. Based on the estimated saturation values and relative root mean square error (RRMSE), the single and fused characteristics and corresponding models were selected and integrated, which led to a novel saturation-based multivariate method used to improve the GSV estimation and mapping of Chinese fir plantations. The new findings included: (1) All the original polarimetric characteristics, statistically, were not significantly correlated with the forest GSV, and their logarithm and ratio transformation fused variables greatly improved the correlations, thus the estimation accuracy of the forest GSV. (2) The logarithm transformation of surface scattering resulted in the greatest saturation, value but the logarithm transformation of double-bounce scattering resulted in the smallest RRMSE of the GSV estimates. (3) Compared with the single transformations, the fused variables led to more reasonable saturation values and obviously reduced the values of RRMSE. (4) The saturation-based multivariate method led to more accurate estimates of the forest GSV than the univariate method, with the smallest value (29.64%) of RRMSE achieved using the set of six variables. This implied that the novel saturation-based multivariate method provided greater potential to improve the estimation and mapping of the forest GSV.

Resistance value (R-value) is one of the basic subgrade stiffness characterizations that express a material’s resistance to deformation. In this paper, artificial intelligence (AI)-based models—especially M5P, support vector machine (SVM), and Gaussian process regression (GPR) algorithms—are built for R-value evaluation that meets the high precision and rapidity requirements in highway engineering. The dataset of this study comprises seven parameters: hydrated lime-activated rice husk ash, liquid limit, plastic limit, plasticity index, optimum moisture content, and maximum dry density. The available data are divided into three parts: training set (70%), test set (15%), and validation set (15%). The output (i.e., R-value) of the developed models is evaluated using the performance measures coefficient of determination (R2), mean absolute error (MAE), relative squared error (RSE), root mean square error (RMSE), relative root mean square error (RRMSE), performance indicator (ρ), and visual framework (Taylor diagram). GPR is concluded to be the best performing model (R2, MAE, RSE, RMSE, RRMSE, and ρ equal to 0.9996, 0.0258, 0.0032, 0.0012, 0.0012, and 0.0006, respectively, in the validation phase), very closely followed by SVM, and M5P. The application used for the aforementioned approaches for predicting the R-value is also compared with the recently developed artificial neural network model in the literature. The analysis of performance measures for the R-value dataset demonstrates that all the AI-based models achieved comparatively better and reliable results and thus should be encouraged in further research. Sensitivity analysis suggests that all the input parameters have a significant influence on the output, with maximum dry density being the highest.

Conservation agriculture, especially no‐tillage, has proven to become sustainable farming in many agricultural environments globally. In spite of advantages of no‐till systems, this practice may result in excess infiltration into the soil and can enhance the movement of mobile nutrients and some pesticides to subsurface drains and groundwater along preferential pathways. The goal of this study was to evaluate the capacity of DRAINMOD‐N II to simulate subsurface nitrate‐N leaching in no‐till fields in Truro, Nova Scotia, Canada, from 2003 to 2006. The model performance was first evaluated by comparing observed and simulated drain outflow data that is an essential prerequisite for the model to obtain a proper prediction of NO3‐N movement, and then by comparing observed and simulated NO3‐N concentration in no‐till fields using three statistical indices, relative root mean square error (RRMSE), average absolute deviation (AAD) and the correlation coefficient (R2). The RRMSE, AAD and R2 for the validation period were determined to be 1.09, 1.85 and 0.83 mm for drain outflow, and 1.43, 0.51 and 0.79 mg l−1 for NO3‐N concentration respectively. The results showed that DRAINMOD‐N II predicted NO3‐N leaching reasonably well in drainage outflow of no‐till fields over the whole period. Copyright © 2018 John Wiley & Sons, Ltd.