Potassium Simulation Using HYDRUS-1D with Satellite-Derived Meteorological Data under Boro Rice Cultivation

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.

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.

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

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.

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.

Information about the state of the forest is of vital importance in forest management planning. To enable high-precision modelling, many forest planning systems demand input data at the single-tree level. The conventional strategy for collecting such data is a plot-wise field inventory. This is expensive and, thus, cost-efficient alternatives are of interest. During recent years, the focus has been on remote sensing techniques. The k nearest neighbour (kNN) estimation method is a way to assign plot-wise data to all stands in a forest area, using remotely sensed data in connection with a sparse sample of field reference plots. Plot-wise aerial photograph interpretations combined with information from a stand register were used in this study. Nearness to a reference plot was decided upon using a regression transform distance. Standing stem volume was estimated with a relative root mean square error (RMSE) equal to 20% at the stand level, while age could be estimated with a RMSE equal to 15%. A cost-efficient data-capturing strategy could be to assign plot data with the presented kNN method to some types of forest, while using traditional field inventories in other, more valuable, stands.

수공구조물의 위험도에 관한 불확실성을 검토하기 위하여 본 연구에서는 빈도해석을 통하여 추정되는 설계홍수량의 분산량을 고려한 불확실성 해석을 실시하였다. Gumbel 분포형을 기본 분포형으로 가정하였으며, 모멘트법, 최우도법, 확률가중모멘트법을 이용하여 각 매개변수 추정방법별로 추정된 설계홍수량에 대한 이론적인 분산량을 산정하였다. 이론적으로 유도한 분산량의 특성을 규명하기 위하며 다양한 표본크기와 설계연한, 비초과확률 및 변동계수조건에 대하여 Monte-Carlo 모의를 실시하고 각 매개변수 추정방법별 비교를 실시하였다. 그 결과 확률 가중 모멘트법을 사용한 경우 위험도에 대하여 상대적으로 가장 작은 상대편의 및 상대제곱근오차를 발생시키는 것으로 나타났으며, 최우도법의 경우에는 상대적으로 큰 표본자료에 대해서는 설계연한 및 비초과확률에 관계없이 작은 상대편의 및 상대제곱근오차를 발생시키는 것으로 나타났다. 또한 다양한 변동계수 조건은 상대편의 및 상대제곱근오차의 측면에서 고려하여 볼 때 거의 영향을 주지 않는 것으로 나타났다. The uncertainty of the risk of failure of hydraulic structures can be determined by estimating the variance of the risk of failure based on the methods of moments, probability weighted moments, and maximum likelihood assuming that the underlying model is the Gumbel distribution. In this paper, the variance of the risk of failure was derived. Monte Carlo simulation was peformed to verify the characteristics of the derived formulas for various sample size, design life, nonexceedance probability, and variation coefficient. As the results, PWM showed the smallest relative bias and root mean square error than the others while ML showed the smallest ones for relatively large sample siBes regardless of design life and nonexceedance probability. Also, it was found that variation coefficient does not effect on the relative bias and relative root mean square error.