Next Article in Journal
Contourlet-CNN for SAR Image Despeckling
Previous Article in Journal
Evaluation of OMI NO2 Vertical Columns Using MAX-DOAS Observations over Mexico City
Previous Article in Special Issue
Corn Nitrogen Status Diagnosis with an Innovative Multi-Parameter Crop Circle Phenom Sensing System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Coupled Elimination Method of Soil Moisture and Particle Size Interferences on Predicting Soil Total Nitrogen Concentration through Discrete NIR Spectral Band Data

Key Laboratory of Modern Precision Agriculture System Integration Research, Ministry of Education, China Agricultural University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(4), 762; https://doi.org/10.3390/rs13040762
Submission received: 28 December 2020 / Revised: 14 February 2021 / Accepted: 17 February 2021 / Published: 19 February 2021
(This article belongs to the Special Issue Remote Sensing for Precision Nitrogen Management)

Abstract

:
Rapid and accurate measurement of high-resolution soil total nitrogen (TN) information can promote variable rate fertilization, protect the environment, and ensure crop yields. Many scholars focus on exploring the rapid TN detection methods and corresponding soil sensors based on spectral technology. However, soil spectra are easily disturbed by many factors, especially soil moisture and particle size. Real-time elimination of the interferences of these factors is necessary to improve the accuracy and efficiency of measuring TN concentration in farmlands. Although, many methods can be used to eliminate soil moisture and particle size effects on the estimation of soil parameters using continuum spectra. However, the discrete NIR spectral band data can be completely different in the band attribution with continuum spectra, that is, it does not have continuity in the sense of spectra. Thus, relevant elimination methods of soil moisture and particle size effects on continuum spectra do not apply to the discrete NIR spectral band data. To solve this problem, in this study, moisture absorption correction index (MACI) and particle size correction index (PSCI) methods were proposed to eliminate the interferences of soil moisture and particle size, respectively. Soil moisture interference was decreased by normalizing the original spectral band data into standard spectral band data, on the basis of the strong soil moisture absorption band at 1450 nm. For the PSCI method, characteristic bands of soil particle size were identified to be 1361 and 1870 nm firstly. Next, normalized index Np, which calculated wavelengths of 1631 and 1870 nm, was proposed to eliminate soil particle size interference on discrete NIR spectral band data. Finally, a new coupled elimination method of soil moisture and particle size interferences on predicting TN concentration through discrete NIR spectral band data was proposed and evaluated. The six discrete spectral bands (1070, 1130, 1245, 1375, 1550, and 1680 nm) used in the on-the-go detector of TN concentration were selected to verify the new method. Field tests showed that the new coupled method had good effects on eliminating interferences of soil moisture and soil particle size.

1. Introduction

Excessive application of commercial fertilizers has been cited as a source of contamination of surface and groundwater [1,2,3,4]. As one of the important parts of precision agriculture (PA), based on the needs of the plant and the soil nutrients state, variable rate (VR) fertilization can reduce the overuse of manure in the farmlands to protect the soil and environment [5,6,7,8]. Rapid, real-time, and accurate acquisition of high-resolution soil parameter information is the basis of applying VR fertilization operation [9].
As the main nutrient for plant and soil, soil total nitrogen (TN) is one of the key factors that determine plant nutrient levels and soil fertility [10,11]. The estimation of TN concentration status is crucial from agricultural and environmental points of view [12,13,14]. Traditional wet chemical methods of measuring TN concentration are time-consuming, costly, have low efficiency, and also cause environmental pollution due to their consumption of harmful chemicals [15,16,17]. Spectral analysis technology is a fast, nondestructive, on-line, and free pollution method that shows an increasing development potential in TN concentration measurement [18,19,20,21]. Many researchers focused on exploring rapid TN detection methods and corresponding soil sensors based on spectral technology in the past four decades. In the 1980s, Dalal and Henry [22] showed it is possible to measure TN concentration through near-infrared (NIR) reflectance spectroscopy. Since then, many scholars have succeeded to predict TN concentration based on NIR reflectance spectroscopy in the laboratory [23,24,25].
However, collecting soil samples in the field requires a lot of human labor. It will also take a relatively long time to obtain the measurement results [17]. Thus, estimation of TN concentration through NIR reflectance spectroscopy in the laboratory does not satisfy the requirement of VR fertilization operation. Simultaneously, many scholars focus on developing soil parameter sensors to supply rapid and high-resolution acquisition of soil parameters. Sudduth and Hummel [26] designed and tested a portable, NIR spectrophotometer intended for in-field use. Mouazen et al. [27,28] developed a soil sensing device with the measuring range of the spectrometer from 306.5 to 1710.9 nm. Decreasing the capital is significant for the development of soil sensors. Thus, cheaper light-emitting diodes (LEDs), laser diodes, and tungsten halogen light sources are used as the alternative component of the spectrometer to develop soil sensors [29,30,31]. Tang et al. [32] developed a portable SOC analyzer using 850 nm as a sensitive wavelength. An et al. [31] developed a portable TN sensor based on six discrete NIR bands (940, 1050, 1100, 1200, 1300, and 1550 nm) to detect TN content, which used the LEDs as the detection light source. Zhou et al. [30] developed an in-situ TN detector with seven discrete NIR bands (1260, 1330, 1360, 1430, 1530, 1580, and 1660 nm) based on the laser diodes, which realized in-situ detection of TN concentration in the farmlands. Li et al. [33] developed the portable TN concentration detector with six discrete NIR bands (1108, 1248, 1336, 1450, 1537, and 1696 nm), which used the tungsten halogen light sources as the detection light source.
However, several problems remain through spectral technology to predict soil parameters. First, the soil spectra are easily affected by interactions with various soil constituents (e.g., minerals and soil nutrients) [34,35]. Especially, field soil samples have more variations in soil moisture and other physical conditions (i.e., soil particle size, and water). Water and hydroxyls (O-H) have strong influences in the NIR region. Overtones of O-H stretching absorb near 1450 nm, and it conceals the absorption information of soil N-H bonds in the NIR spectroscopy, resulting in the spectral curve only reflecting the changing trend of soil moisture. However, the absorption information of soil N-H bonds reflects the concentration of TN. Reflectance from a soil sample varies with its soil particle size in the NIR regions of the electromagnetic spectrum. As the soil particle size decreased, the reflectance also gradually increased. It also conceals the absorption information of soil N-H bonds in the NIR spectroscopy. Thus, the soil moisture and particle size effects on the NIR spectroscopy must be eliminated before performing the estimation of TN by the NIR spectroscopy [36,37,38].
Many scholars have analyzed the interferences of soil particle size and moisture on soil parameter prediction, and several studies have proposed methods to eliminate these disturbances. For the soil moisture, Tekin et al. [39] used NIR reflectance spectroscopy to study the effects of different soil moisture contents on soil organic matter (SOM) prediction. The root mean square error of validation (RMSEP) and the residual prediction deviation (RPD) for dry soil samples reached 1.26% and 2.83, respectively, and the soil samples with the highest moisture content had values of 1.55% and 2.29, respectively. Zhang et al. [23] studied the real-time NIR spectra of soil samples without artificial drying and sieving, and six sensitive bands were determined for predicting TN concentration through wavelet analysis and continuum removal techniques. Prediction models of TN concentration were established through partial least squares regression (PLSR) and support vector machine (SVM), respectively. For the PLSR model, the coefficient of determination of calibration ( R c 2 ) was 0.602, and the coefficient of determination of prediction ( R v 2 ) was 0.634. For the SVM model, R c 2 reached 0.823, and R v 2 reached 0.810. It revealed that the soil moisture interference was removed from the real-time soil NIR reflectance spectroscopy during TN prediction using the proposed method of wavelet analysis and continuum removal. External parameter orthogonalization, direct standardization, global moisture modeling, slope bias correction, and selective wavelength modeling have also been applied to eliminate the effect of soil moisture on soil NIR reflectance spectroscopy [40]. For the soil particle size, Barthès et al. [41] analyzed the NIR reflectance spectroscopy of sandy soils with different particle sizes in Burkina Faso and Congo and established the prediction models of TN concentration under different soil particle sizes. The results showed that the model with less than 2 mm soil particle size was better than the model established using other soil particle sizes. Bogrekci and Lee [42] investigated the effects of soil particle size on the prediction of soil rapid available phosphorus (P) using ultraviolet, visible, NIR spectroscopy. They found that for the sandy soil, with the decrease in soil particle size, the soil reflectivity increased and the absorbance decreased. The individual prediction method for each soil particle size produced better soil P predictions with lower standard errors of prediction for coarse, medium, and fine sand particles than the method that removed the effect of particle size on absorbance spectra and predicted P concentration of the soil samples.
Summarizing the above research, there was no relevant research on the elimination of soil moisture and particle size effects on predicting soil parameters through the discrete NIR spectral band data. Soil parameter detectors developed based on discrete NIR spectral band data also face serious interferences caused by soil moisture and particle size. The methods used in the laboratory to eliminate the interference of soil particle size and moisture by sieving after grinding and drying cannot be applied to the discrete NIR spectral band data in the field sensing. The above elimination methods of soil particle size and moisture are not suitable for discrete NIR spectral band data acquisition in the field sensing. An et al. [36] studied the effects of soil moisture and particle size on the prediction of TN concentration using a portable TN detector developed based on six discrete NIR bands (940, 1050, 1100, 1200, 1300, and 1550 nm). The moisture absorption index and mixed calibration set were proposed to eliminate the soil moisture and particle size interferences, respectively. However, the mixed calibration set method entailed a complicated process and could not perform real-time elimination of the interference of soil particle size.
The combination effects and real-time elimination methods of soil moisture and particle size on predicting TN concentration through the discrete NIR spectral band data have rarely been explored. The discrete NIR spectral band data can be completely different in the band attribution with continuum spectra, it does not have continuity in the sense of spectra. Therefore, the study of methods for eliminating the effects of soil particle size and moisture on the discrete NIR spectral band data has important theoretical significance and practical application value for improving the performance of detectors based on the discrete NIR spectral band data. The main goal of the paper is to establish an on-line elimination method of soil moisture and particle size interferences, which could be used in the discrete NIR spectral band data to improve the measurement accuracy of soil parameters. Thus, in this study, six discrete NIR spectral bands (1070, 1130, 1245, 1375, 1550, and 1680 nm) used in the on-the-go detector of TN concentration were chosen as the object to perform the research [29]. A new coupled elimination method of soil moisture and particle size interferences on predicting TN concentration through discrete NIR spectral band data was proposed and evaluated.

2. Materials and Methods

2.1. Major Steps of the Study

Figure 1 shows the major procedures in this research, which included three stages. During stage 1, 60 soil samples were collected from Field 1. In order to research the effect of soil moisture on discrete NIR spectral band data, the soil particle size was eliminated by artificial sieving. Additionally, the moisture absorption correction index (MACI) method to eliminate the disturbance of soil moisture on discrete NIR spectral band data was proposed and validated. During stage 2, 24 soil samples were collected from Field 2. Before performing analysis, soil moisture disturbance was eliminated through drying by the oven. Next, four groups of soil samples with different particle sizes were obtained, each group contained six TN concentration levels, and each TN concentration level contained four soil samples. A total of 96 soil samples with different particle sizes and TN concentrations were generated from the 24 soil samples. Next, the effect of soil particle size on discrete NIR spectral band data was explored and the characteristic wavebands of soil particle size were determined. The particle size correction index (PSCI) method was proposed and validated to eliminate the effect of soil particle size on discrete NIR spectral band data.
During stage 3, a field experiment was performed in Field 3 to evaluate the coupled elimination method through the on-the-go detector, which included five processing methods. Finally, five models of TN concentration were calibrated and compared.

2.2. On-The-Go Detector of TN Concentration

An on-the-go detector of TN concentration was previously developed through discrete NIR spectral band data. This detector was reported in our early research [29]. Figure 2 shows the structure of the on-the-go detector, which consisted of mechanical, optical, and control units. The mechanical unit was the basis of the entire detector and supplied a downward force for the work of the subsoiler. Thus, the electrical control cabinet of the mechanical unit was designed as a three-tier structure to integrate the optical and control units. The optical unit is used for sensing soil in the soil detection trench and a control unit embedded into an SVM model for obtaining real-time readings of TN concentration. A tungsten halogen light source was a better choice for the vehicle equipment than LED light and laser sources with respect to light intensity and stability, which provided a broadband source over the 300–2500 nm range. The selected performance parameters of the tungsten halogen light were as follows: <25% of radiation uniformity, <1% of unstable time, 20 × 20 mm radiation area, adjustable radiation intensity, 150 W of maximum power, and 2000 h of average life. The light had to be transmitted to the soil surface with minimal light loss. Therefore, an NIR optical fiber with a transmission range of 800–1900 nm was used. The diameter of the selected filter was 5 mm, and the average full width at half maximum (FWHM) bandwidth was 20 nm. The selected photoelectric sensor was InGaAs with a responding range of 800–1800 nm and a spectral responsiveness of 0–0.9. The sapphire glass attached to the bottom of the sensor mounting plate was used to ensure that the sensors and the filters were free from soil pollution. When the on-the-go detector performed detection of TN concentration in farmland, it was mounted on the tractor through a three-point suspension structure. The depth of the subsoiler that pushed into the soil exceeded 200 mm, which was used to loosen the soil and smash large soil particles. The flattener made the soil more flat and suitable for spectrum detection. The NIR light emitted by the tungsten halogen light source was transmitted through an NIR fiber to the flat detection soil surface. The diffused reflection light of the detected soil was changed into a single-band light by each optical filter, and each photoelectric detector converted the single-band light into an electrical signal. All detected electrical signals were amplified, filtered, and digitized into absorbance data, and then all data were collected and displayed in real-time by the data receiving software embedded with the whiteboard data, by which and detected spectral data the absorbance data were calculated at each sensitive wavelength. For this on-the-go detector, the group of the discrete NIR spectral bands of 1070, 1130, 1245, 1375, 1550, and 1680 nm was selected as the sensitive wavebands to predict TN concentration. Furthermore, the on-the-go detector was used to embed the new coupled elimination method of soil moisture and particle size interferences on predicting TN concentration through discrete NIR spectral band data.

2.3. Experimental Materials and Methods

During stages 1 and 2 (Figure 1), soil samples were collected from the experimental farm of China Agricultural University on the outskirts of Beijing to perform the research of soil moisture and particle size interference elimination method. The soil type was common brown soil in northern China. The experimental area is in the range of north latitude 40.15395721–40.1534362 and east longitude 116.2155061–116.2150130. Since the fertilization usually occurs from the soil surface to the soil depth of 30 cm, soil samples were collected from 5 to 30 cm after removing 5 cm of the topsoil [43]. A total of 60 soil samples used to study the effect of soil moisture were collected at Field 1 (F1, 10 ha). Corn and wheat were cropped in F1 with tillage. In total, 24 soil samples were collected from the Field 2 (F2, 6 ha), which was a standard field without fertilization (the soil nutrient concentration is almost zero). The soil textural fractions (clay, silt, and sand) of F2 were the same with F1, both the two fields were sandy loam. In total, 2–2.5 kg of soil at each sampling point was collected and put it into double-layer cowhide bags for sealing to prevent the loss of soil moisture. A total of 84 soil samples were collected.
The first set of soil samples (N = 60) was used to study the effect of soil moisture. All soil samples were sieved under laboratory conditions using a 20-mesh sifter (0.9 mm) to produce pretreated samples with consistent particle sizes [36].
Another set of soil samples (N = 24) was used to study the effect of soil particle size. In order to obtain soil samples with different TN concentration, the standard soil samples were first dried in a drying box, the temperature was set to 85 ℃, and the time for drying was 24 h [44,45]. The air-dried soil samples were randomly divided into six groups, four in each group. The soil nitrogen solution was configured with urea [46,47,48]. The concentration gradient of the soil nitrogen solution is divided into 1–6 grades, grade 1 nitrogen concentration is 0 g·kg−1, grade 2 nitrogen concentration is 0.04 g·kg−1, grade 3 nitrogen concentration is 0.08 g·kg−1, grade 4 nitrogen concentration is 0.12 g·kg−1, grade 5 nitrogen concentration is 0.16 g·kg−1, and grade 6 nitrogen concentration is 0.2 g·kg−1. When configuring the soil samples, we simulated the average summer soil moisture content (7%) of the cultivated layer of the Shangzhuang Experimental Station to configure the soil samples. After the soil configuration was completed, the drying process was carried out, and all the soil samples were sieved, respectively, 10-mesh sieves (2.0 mm), 20-mesh sieves (0.9 mm), 40-mesh sieves (0.45 mm), and 80-mesh sieves (0.2 mm). Finally, four groups of soil samples with different particle sizes were obtained, each group contained six TN concentration levels, and each TN concentration level contained four soil samples. A total of 96 soil samples with different particle sizes and TN concentrations were generated from the 24 soil samples.
During stage 3 (Figure 1), soil samples for verification were used to evaluate the new coupled elimination method. A verification experiment was conducted on the experimental farm of China Agricultural University. The soil type was common brown soil in northern China Field 3 (F3, 5.5 ha) which was one long-term research site located within 1.0 km of F1. Corn and wheat were cropped in F3 with tillage. The soil textural fractions (clay, silt, and sand) of F3 were the same as that of F1, which were sandy loam. Figure 3 shows the location and soil samples collection point in Field 3.
The experiment consisted of two parts which were data of discrete NIR spectral band data acquisition and soil samples collection. The experimental plot was 240 × 190 m, which was divided into 108 sampling cells. Each cell had a dimension of 5 × 20 m. Data and soil samples were collected every 10 m in each cell, blue cells in Figure 3 show the location of acquiring data and soil samples. All samples were placed in a double-thickness sampling bag to prevent moisture dispersion and then sent to the laboratory for optical and chemical measurements.

2.4. Laboratory Measurement

For each soil sample, it was divided into two parts, one was used for TN concentration measurement, and another was used for NIR reflectance spectroscopy. The TN concentration was tested with the FOSS Kjeltec™2300 Nitrogen Analyzer produced by FOSS, Sweden [49]. After grinding each soil sample, weighing 2.0 g, was put into a long test tube. The long test tube was also filled with 6.2 g catalyst. K2SO4:CuSO4.5H2O was mixed and ground at 30:1. Finally, 20 mL of concentrated sulfuric acid was added to the long test tube for nitrification. The temperature of the nitration furnace was set to 420 °C for 2.5 h [50]. After nitrification, it was cooled down and the total nitrogen concentration was measured with the FOSS Kjeltec™2300 nitrogen analyzer. The NIR reflectance spectroscopy of soil samples was collected by the MATRIX-I Fourier Transform NIR Spectroscopy Analyzer produced by Bruker, Germany. Before soil spectra collection, the analyzer parameters needed to be set. The spectrum collection range, spectrum sampling interval, and scan times were set to 780–2550 nm, 3 nm, and 64 times, respectively. Approximately 20 g of soil samples were put into a quartz cuvette with a diameter of 50 mm, and put into the rotating sample cell of the analyzer for spectrum detection [17,51].

2.5. Model Accuracy and Methodology

In addition to the coefficient of determination (R2), the root mean square error of validation (RMSEP), and the residual prediction deviation (RPD) were used to evaluate the prediction model of TN concentration [52]. Equations (1) and (2) present their calculation Equations.
RMSEP   =   j   =   1 n p y j y ^ p j 2 n p 1 ,
RPD   =   i   =   1 n p y j y ¯ j 2 / n p 1 / i   =   1 n p y j y ^ p j 2 / n p 1 ,
where j is the number of soil samples, yj is the TN concentration measured with the FOSS Kjeltec™2300 Nitrogen Analyzer, y ^ c j is the TN concentration predicted by the modeling set, y ^ p j is the TN concentration predicted by the validation set, y ¯ j is the average value of the TN concentration measured with the FOSS Kjeltec™2300 Nitrogen Analyzer, nc is the model set sample number, and np is the validation set sample number.

3. Results

3.1. Research on Eliminating the Interference of Soil Moisture

3.1.1. Interference of Soil Moisture on Discrete Spectral Band Data

The first set of 60 calibration samples with consistent particle sizes were measured through the developed on-the-go detector. To analyze the effect of soil moisture on soil discrete NIR spectral band data quantitatively, the spectral data of a soil sample with TN concentration of 0.1135 g·kg−1 was used as the research object. Figure 4 shows the absorbance values at seven sensitive wavelengths under different soil moisture contents. The absorbance curve obtained by the on-the-go detector from 1070 to 1680 nm showed a consistent trend, indicating that the curve was a good reflection of the soil spectral information [53]. Taking the soil absorbance at 1450 nm as an example, when the soil moisture content increased from 3.18% to 14.16%, the soil absorbance increased from 0.79 to 1.11, and the soil spectral absorbance increased by 40.51%. Absorbance data analysis indicated that along with the soil moisture content increase, the absorbance of the soil spectrum increased rapidly. Soil moisture had a serious influence on soil discrete NIR spectral band data [54,55,56]. The effects of soil moisture must be decreased when predicting TN concentration.
Our goal is to eliminate the influence of soil moisture on the prediction of TN concentration through the on-the-go detector. Soil NIR spectroscopy indicates that three distinct moisture absorption peaks are present in the soil NIR spectroscopy at 1450, 1940, and 2210 nm [57,58]. However, in our research, we only focused on the absorption peak of O-H in H2O at 1450 nm. O-H have strong influences in the NIR spectroscopy. Overtones of O-H stretching absorb near 1450 nm, and it conceals the absorption information of soil N-H bonds, resulting in the spectral curve only reflecting the change trend of soil moisture. However, the absorption information of soil N-H bonds reflects the concentration of TN. If the raw spectrum of the measured soil can be normalized to a standard spectrum with the same soil moisture content, the absorption information of soil N-H bonds on the NIR spectroscopy could be effectively extracted. Thus, the interference of soil moisture on predicting TN concentration could be decreased. However, when using the on-the-go detector to measure soil discrete NIR spectral band data, we were unable to directly obtain the soil moisture content in the fields. Thus, an alternative method must be explored to solve this problem.

3.1.2. Eliminating the Interference of Soil Moisture on Discrete NIR Spectral Band Data

To eliminate the interference of soil moisture, we need to classify different soil samples based on soil moisture content. The previous analysis revealed the difficulty of obtaining soil moisture content in real-time. Therefore, the MACI method was proposed to eliminate the effect of soil moisture on predicting TN concentration through the on-the-go detector. It was used to convert raw soil spectral data into standard spectral data. The MACI value of the on-the-go detector was calculated using Equation (3).
MACI   =   1 A 1450 1 + A 1450 × 100 ,
where A1450 is the spectral absorbance value of the soil samples at 1450 nm and MACI is the soil moisture absorption index correlation. Soil samples with different moisture contents can be divided into groups based on MACI values; thus, we divided the detected soil samples into five types based on MACI values. Table 1 shows the five MACI values when the MACI value is used to represent the soil moisture level. Their corresponding ranges of soil moisture content and accuracy are also shown. When the soil moisture was between 0% and 3.0%, the MACI value was less than 8.66. When the soil moisture was 3.0–6.0%, 6.0–10.0%, and 10.0–13.0%, the MACI values were 8.66–12.66, 12.66–14.66, and 14.66–16.66, respectively. When the soil moisture content was higher than 13%, the MACI value was greater than 16.6.
The classification of soil samples was the best according to the classification of MACI values. Although the soil moisture content could not be directly obtained in real-time when the on-the-go detector was used to measure soil in the fields, the soil samples could be classified by MACI.
Each classification group required a correction factor to normalize the spectral data, and Equation (5) was used to calculate the correction factor. In this equation, j (j = 1, 2,..., 5) is the serial number of the taxonomic group, Wj is the correction factor for the taxonomic group j, Aavgj is the soil sample average absorbance of group j, and Aavgs is the soil sample average absorbance of group 1. The soil sample average absorbance value of group 1 was used as the reference absorbance value. Table 1 also shows all five correction factors.
W j   =   A a v g j A a v g s
After determining the correction factor Wj (j = 1, 2,…, 5), the raw absorbance data were corrected according to
A j c   =   A j × W j ,
where Aj is the soil sample original absorbance value of group j, Ajc is the soil sample corrected normalized absorbance value, and Wj is the correction factor.
Absorbance data at 1070, 1130, 1245, 1375, 1450, 1550, and 1680 nm were detected using the on-the-go detector. Using the proposed MACI method to eliminate the interference of soil moisture, the wavelength absorbance at 1450 nm was selected to calculate the MACI value, and the raw spectral data were normalized to eliminate the interference caused by soil moisture content. Figure 5 shows the different absorbance curves for the change in soil moisture content after normalization. When the soil moisture content was 3.18% and 14.16%, the absorbance values at 1450 nm were 0.72 and 0.75, respectively; the absorbance value only increased by 4.16%. Compared with the soil spectrum without correction by the MACI method, the effect of soil moisture on soil discrete NIR spectral band data was greatly reduced, and the MACI method proposed in this study was proven to be effective in reducing soil moisture on the discrete NIR spectral band data.

3.2. Research on Eliminating the Interference of Soil Particle Size

3.2.1. Eliminating the Interference of Soil Moisture on Discrete NIR Spectral Band Data

Table 2 statistically analyzes the 96 soil samples prepared. The TN concentration in the soil was detected by the FOSS Kjeltec™2300 Nitrogen Analyzer [49]. The TN concentration value shows that the TN concentration of the prepared soil sample is distributed in the range of 0.003–0.206 g·kg−1, with approximately 0.04 g·kg−1 as the TN concentration interval, divided into six TN concentration grades, the gradient distribution of TN concentration grade is reasonable.
According to the analysis in Table 2, although the soil samples are collected from standard fields that are not fertilized all year, there is still TN in the soil samples, which will affect the TN concentration of the prepared soil samples.
Using the drying oven at 80 °C for 24 h, the second set of 96 soil samples was dried to eliminate the influence of soil moisture [17]. Subsequently, all soil samples in the second set were sieved by sifters of different meshes, and discrete NIR spectral band data were obtained by the on-the-go detector. Figure 6 shows the absorbance value of a single soil sample at six discrete NIR spectral bands with the TN concentration of 0.086 g·kg−1, including five curves, which are the original soil spectrum, the soil spectra with soil particle sizes of 2.0, 0.9, 0.45, and 0.2 mm. When the soil particle size decreased from 2 to 0.2 mm, the corresponding soil samples’ average absorbance also decreased from 0.74 to 0.63, and the overall absorbance changed by 14.86%. The results showed that soil particle size had a considerable effect on the soil discrete NIR spectral band data. When the soil particle size was between 0.45 and 0.9 mm, the absorbance of soil discrete NIR spectral band data was relatively stable. Larger and smaller particle sizes will cause a large difference on soil absorbance. Figure 6 also shows that the air-dried method could significantly decrease the interference of soil moisture on the spectral band of 1450 nm. Therefore, no absorption peak was observed at 1450 nm.

3.2.2. Identification of Characteristic Wavebands of Soil Particle Size

Figure 7 shows the spectral standard deviation curves of four different soil particle sizes. It can be concluded from Figure 7 that the standard deviation value in the range of 850–1000 nm was smaller than in the spectral range of 1000–2500 nm. As each spectral standard deviation curve is under the same particle size, the interference of soil particle size is eliminated, and only the TN concentration is different. Thus, the 1000–2500 nm range is the sensitive spectral region of TN concentration. Simultaneously, the black curve in Figure 7 is the spectral standard deviation curve of all 96 soil samples, by the comparative analysis between the five standard deviation curves in the spectral range of 1000–2500 nm. After adding the difference in soil particle size, the black standard deviation curve in the spectral range between 1340–1420 nm and 1795–1955 nm changed up and down. In order to show the changes between different spectral standard deviation curves, spectral ranges 1340–1420 and 1795–1955 nm of the spectral standard deviation curve were extracted in Figure 7. Obviously, this change indicates that the above two spectral regions are closely related to the change of soil particle size and are sensitive spectral regions of soil particle size. Therefore, in this study, the extreme points 1361 and 1870 nm of the two spectral regions were selected as the characteristic wavebands of soil particle size.

3.2.3. Research on the Classification of Soil Particle Size

The 96 soil samples were divided into two groups, with 48 soil samples in each group. Each group of soil samples includes four soil particle size categories, each soil particle size category includes six TN concentration grades, and each soil concentration grade includes two soil samples. The two characteristic wavebands 1361 and 1870 nm were used to classify the soil particle size. The normalized index Np was used as a single variable to predict the soil particle size, and the Np is calculated by Equation (6). The soil classification model is constructed by the support vector machine (SVM) algorithm. After many attempts, the penalty parameter c of the SVM classification model is 2 and the kernel function is 1 to obtain the optimal classification result. Figure 8 shows the prediction results of soil particle size based on SVM. The overall classification accuracy rate of the four soil categories is 97.92%. Among the four soil categories, the classification accuracy rate when the soil particle sizes are 2.0, 0.45, and 0.2 mm is 100%, and the classification accuracy rate when the soil particle size is 0.90 mm is 91.67%.
N p   =   A 1870 A 1361 A 1870 + A 1361 ,
In Equation (6): A1870 and A1361 are the absorbance of the soil at 1870 and 1361 nm.

3.2.4. Eliminating the Interference of Soil Particle Size

Based on the soil particle size classification model established by SVM, the normalized index Np was used as a single input to classify the soil particle size. The results show that the method is feasible. Based on the above results, a soil particle size correction coefficient Pc was proposed to correct the origin discrete NIR spectral band data. The soil particle size correction coefficient Pc is obtained by Equation (7),
P c   =   A 1870 A 1361 A 1870 + A 1361   N ¯ p ,
In Equation (7): A1870 and A1361 are the absorbance of the corrected soil spectrum at 1870 and 1361 nm, and N ¯ p   is the normalized index of the absorbance of the 0.20 mm soil at 1870 and 1361 nm as the reference value. It is generally believed that the smaller soil particle size can eliminate the influence of soil particle size to the greatest extent, so the spectrum data with a particle size of 0.20 mm was determined as the reference spectrum.
The absorbance Ac of discrete NIR spectral band data corrected by the Np was obtained by Equation (8),
A c   =   A × P c ,
In Equation (8): A is the soil absorbance of discrete NIR spectral band data before correction, and Ac is the soil absorbance after correction. In this study, the absorbance at six discrete NIR bands (1070, 1130, 1245, 1375, 1550, 1680 nm) of the on-the-go detector were used to verify the PSCI method. Figure 9 shows the absorbance of discrete NIR spectral band data after soil particle size correction. A comparative analysis of Figure 6 and Figure 9 was computed. For example, at 1245 nm, the difference between the absorbance value with a soil particle size of 0.2 mm and the absorbance value with soil particle sizes 2.0, 0.9, 0.45 mm, and original spectrum were 0.0737, 0.0506, 0.0361, and 0.0607, respectively. After correction with the PSCI method, the difference was 0.0263, 0.0111, −0.0011, and 0.0228, and the difference was reduced by 64.31%, 78.06%, 103.04%, and 62.43% after correction. The results show that the soil particle size correction coefficient can significantly reduce the interference of soil particle size to a greater extent.

3.3. Coupled Elimination Method of Soil Moisture and Particle Size Interferences

In order to achieve on-line coupled elimination of soil moisture and particle size interferences on predicting TN concentration through discrete NIR bands data, we performed detection of TN concentration through the soil sensors developed based on discrete spectral bands in the farmlands. MACI was firstly used to decrease the interference of soil moisture on discrete spectral band data. Then, PSCI was calculated to decrease the effect of soil particle size on the discrete spectral band data after corrected with the MACI method. Finally, TN concentration was obtained through the discrete spectral band data with coupled elimination of soil moisture and particle size interferences. Figure 10 shows the flowchart of the coupled elimination method.

3.4. Evaluation of the New Coupled Elimination Method

Discrete NIR spectral band data was obtained used the on-the-go detector by in-field and laboratory measurement. The spectral data after different moisture and particle size removal processes were used to establish the TN prediction model. The elimination effects of different soil particle size and moisture removal methods on the prediction of TN concentration by discrete NIR spectral band data were evaluated. Figure 11 showed the comparison of elimination effects. The data used for prediction models included (a) raw discrete NIR spectral band data, model 1; (b) discrete NIR spectral band data using the proposed MACI correlation method to eliminate soil moisture disturbance, model 2; (c) discrete NIR spectral band data using the proposed PSCI method to eliminate soil particle size disturbance, model 3; (d) discrete NIR spectral band data using the proposed new coupled elimination method to eliminate soil moisture and particle size disturbances, model 4; (e) discrete NIR spectral band data on dried and ground soil samples, model 5. SVM was used to establish a prediction model for TN concentration based on data treated by different methods. In the established SVM model, the kernel function is the radial basis kernel function, three parameters affecting model regression results: penalty parameter C, radial basis kernel function γ parameter λ, and loss parameter ε, and they are 1.3, 0.25, and 0.01, respectively. Table 3 and Figure 10 show the calibration results. For the raw discrete NIR spectral band data, the prediction R2 was 0.65. However, when the model was established by discrete NIR spectral band data after eliminating soil moisture and soil particle size interference, the RPD value of prediction model 4 reached 2.59, and the prediction R2 reached 0.84, respectively. For models 2 and 3, there is only soil moisture or particle size disturbance eliminated by the MACI or PSCI method, respectively. However, both of them have an improvement of prediction R2 for 12.3% and 9.2% compared with model 1. The above analysis shows that the effects of soil moisture and soil particle size on the discrete NIR spectral band were serious. Further, the TN concentration prediction model 5 was established under laboratory conditions through the discrete NIR spectral band data obtained from the soil samples after drying and grinding. The prediction accuracy of model 4 is still below that of model 5. Although, the proposed MACI and PSCI methods were quantitatively proven to have good effects on eliminating soil moisture and soil particle size interferences on predicting TN concentration through discrete NIR spectral band data. However, the 8.5% decrease of prediction was R2 between models 4 and 5, which means that the effects of soil moisture and particle size were not eliminated using the new coupled elimination method. Thus, it is necessary to continue to explore other elimination methods of soil moisture and particle size effects on predicting TN concentration through discrete NIR spectral band data.

4. Discussion

4.1. Role of New Coupled Elimination Method in Predicting TN through Discrete NIR Spectral Band Data

From the perspective of cost, the development of real-time rapid soil parameter detectors based on discrete NIR spectral band has become a current research hotspot [29,30,31,32,33]. However, for discrete NIR band data, adjacent spectral data can be completely different in the attribution of the bands, and they do not have continuity in the sense of spectrum. Thus, there is no theoretical basis for direct differentiation of a set of discrete NIR spectral band data that has no connection in the spectral sense [59]. In addition, due to the lack of discrete NIR spectral data, it does not meet the requirements of the Savitzky-Golay convolution derivation method for the number of window points, and its operation cannot be performed [60,61]. Soil particle size and moisture will cause serious interference with discrete NIR band data and affect its use [36]. Therefore, the study of effective correction methods to eliminate the interference of soil particle size and moisture on discrete NIR spectral band data has important theoretical significance and practical application value for improving the performance of instruments developed based on discrete NIR spectral band data.
There is no relevant research on the elimination of soil moisture and particle size effects on predicting soil parameters through the discrete NIR spectral band data. This study analyzes the effects of soil moisture and soil particle size on the discrete NIR spectral band data using qualitative and quantitative methods. The absorbance of soil samples increased with the increase in soil moisture content at the same TN concentration level. As the soil particle size decreased, the absorbance of soil samples decreased gradually at the same TN concentration level [36]. Correspondingly, we proposed a method to eliminate the effects of soil moisture and soil particle size on discrete NIR spectral band data. A field experiment also was conducted to evaluate the elimination effects of the proposed MACI and PSCI methods. When the model was established by discrete NIR spectral band data after eliminating soil moisture and soil particle size interference, the RPD value of prediction model 4 reached 2.59, and the prediction R2 reached 0.84, respectively. The results showed that the new coupled elimination method could effectively reduce the interferences of soil particle size and moisture on the prediction of TN concentration.
From the perspective of actual use, the new coupled elimination method can be applied to other on-the-go and portable soil parameter testing equipment developed based on discrete NIR spectral band data. This study facilitates the development of soil sensors based on discrete NIR spectral band data.

4.2. Comparison of the New Coupled Elimination Method to the Similar

Table 4 shows the comparison results of the new coupled elimination method to the similar. According to the evaluation parameters ( R v 2 and RMSEP), our study obtains slightly high estimation accuracies compared to the similar in Table 4. Even though the new coupled elimination method was not applied to the on-the-go detector, the estimation accuracies of TN were higher than the portable detector and the Veris P4000. Applying the new coupled elimination method to the on-the-go detector, the estimation accuracy was improved significantly. The R v 2 increased from 0.64 to 0.82, and the RMSEP decreased from 0.278 to 0.166. For the portable detector, the R v 2 was 0.76 with the moisture absorbance correction method and mixed calibration set used in the portable detector. Compared to the similar, the new coupled elimination method achieved better elimination effects of soil particle size and moisture on NIR spectroscopy. These comparison results further verify the validity of the new coupled elimination method in eliminating the effects of soil particle size and moisture on discrete NIR spectral band data.

4.3. Uncertainty in Current Work and Future Work

Although the new coupled elimination method was quantitatively proven to have good effects on eliminating soil moisture and soil particle size interferences on predicting TN concentration using discrete NIR spectral band data, the 8.5% decrease of prediction R2 was between models 4 and 5, which means that the effects of soil moisture and particle size were not eliminated using the new coupled elimination method. Thus, it is necessary to continue to explore other elimination methods of soil moisture and particle size effects on predicting TN concentration through discrete NIR spectral band data.
Additionally, the experimental conditions in fields are complex, and different regions and soil types also will affect the measurement of discrete NIR spectral band data of soil.
Moreover, field vehicle-mounted detectors are still limited in terms of coving large fields in a short time. The unmanned aerial vehicles (UAVs) have provided the potential alternative to the quick estimation of soil properties coving large fields in a short time. Thus, the hyperspectral cameras or the soil detectors based on the discrete NIR spectral band data mounted on the UAVs will be explored in future research. Therefore, further experiments are needed to verify the effectiveness and practicability of the proposed methods.

5. Conclusions

Commercialized soil nutrients detectors developed based on discrete NIR spectral band data cannot perform on-line elimination of soil moisture and particle size disturbances on predicting TN concentration. Bringing the detected data back to the laboratory for particle size and moisture interferences’ elimination is necessary. This process severely restricts the efficiency of field testing and cannot achieve a quick estimation of TN concentration in the farmlands. This paper reports on the online coupled elimination method of soil moisture and particle size interferences on predicting TN concentration through the discrete NIR spectral band data. The field test also was performed to evaluate the proposed coupled elimination method through the on-the-go detector.

Author Contributions

P.Z. and W.W. designed experiments, conducted field data and sample collection, and supervised soil processing and laboratory analyses; P.Z. and W.Y. implemented spectral data processing and calibration modeling under the supervision of W.Y. and M.L.; the original draft of the manuscript was written by P.Z. with editorial contributions from W.Y., and M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Key Research Projects (2017YFD0201500–2017YFD0201501, 2016YFD0700300–2016YFD0700304) and the National Natural Science Foundation of China (31801265).

Data Availability Statement

Not applicable.

Acknowledgments

The authors wish to thank teachers and students in the Key Laboratory of Modern Precision Agriculture System Integration Research for the help of Sample processing and chemical analysis. We would like to thank Weichao Wang for his help with field data collection. Additionally, great thanks to Sudduth from the USDA-ARS Cropping Systems and Water Quality Research Unit for his suggestions to improve this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bordoli, J.M.; Mallarino, A.P. Deep and shallow banding of phosphorus and potassium as alternatives to broadcast fertilization for no-till Corn. Agron. J. 1998, 90, 27–33. [Google Scholar] [CrossRef] [Green Version]
  2. Domangue, R.J.; Mortazavi, B. Nitrate reduction pathways in the presence of excess nitrogen in a shallow eutrophic estuary. Environ. Pollut. 2018, 238, 599–606. [Google Scholar] [CrossRef]
  3. Shi, P.; Zhang, Y.; Song, J.; Li, P.; Wang, Y.; Zhang, X. Response of nitrogen pollution in surface water to land use and social-economic factors in the Weihe River watershed, northwest China. Sustain. Cities Soc. 2019, 50, 101658. [Google Scholar] [CrossRef]
  4. Staver, K.W.; Brinsfield, R.B. Patterns of soil nitrate availability in corn production systems: Implications for reducing groundwater contamination. J. Soil Water Conserv. 1990, 45, 318–323. [Google Scholar]
  5. Martins, R.N.; Sárvio, D.; Valente, M.; Tadeu, J.; Rosas, F.; Santos, F.S.; Ferreira, F.; Dos, L.; Carolina, A.; Nascimento, C. Communications in Soil Science and Plant Analysis Site-specific Nutrient Management Zones in Soybean Field Using Multivariate Analysis: An Approach Based on Variable Rate Fertilization Site-specific Nutrient Management Zones in Soybean Field Using Multivariate Analysis: An Approach Based on Variable Rate. Commun. Soil Sci. Plant. Anal. 2020, 51, 687–700. [Google Scholar]
  6. Yu, J.; Yin, X.; Raper, T.B.; Jagadamma, S.; Chi, D. Nitrogen Consumption and Productivity of Cotton under Sensor-based Variable-rate Nitrogen Fertilization. Agron. J. 2019, 111, 3320–3328. [Google Scholar] [CrossRef]
  7. Qi, J.; Tian, X.; Li, Y.; Fan, X.; Yuan, H.; Zhao, J.; Jia, H. Design and experiment of a subsoiling variable rate fertilization machine. Int. J. Agric. Biol. Eng. 2020, 13, 118–124. [Google Scholar] [CrossRef]
  8. Mouazen, A.M.; Kuang, B. Soil & Tillage Research On-line visible and near infrared spectroscopy for in-field phosphorous management. Soil Tillage Res. 2016, 155, 471–477. [Google Scholar]
  9. Nawar, S.; Corstanje, R.; Halcro, G.; Mulla, D.; Mouazen, A.M. Delineation of Soil Management Zones for Variable-Rate Fertilization: A Review. Adv. Agron. 2017, 143, 175–245. [Google Scholar]
  10. Batjes, N.H. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 1996, 47, 151–163. [Google Scholar] [CrossRef]
  11. Wang, S.; Zhuang, Q.; Jin, X.; Yang, Z.; Liu, H. Predicting Soil organic carbon and soil nitrogen stocks in topsoil of forest ecosystems in northeastern china using remote sensing data. Remote Sens. 2020, 12, 1115. [Google Scholar] [CrossRef] [Green Version]
  12. Maeda, Y.; Tashiro, N.; Enoki, T.; Urakawa, R.; Hishi, T. Effects of species replacement on the relationship between net primary production and soil nitrogen availability along a topographical gradient: Comparison of belowground allocation and nitrogen use efficiency between natural forests and plantations. For. Ecol. Manag. 2018, 422, 214–222. [Google Scholar] [CrossRef]
  13. Jin, Z.; Chen, C.; Chen, X.; Hopkins, I.; Zhang, X.; Han, Z.; Jiang, F.; Billy, G. The crucial factors of soil fertility and rapeseed yield–A five year field trial with biochar addition in upland red soil, China. Sci. Total Environ. 2019, 649, 1467–1480. [Google Scholar] [CrossRef]
  14. Yoshida, H.; Takehisa, K.; Kojima, T.; Ohno, H.; Nakagawa, H. Modeling the effects of N application on growth, yield and plant properties associated with the occurrence of chalky grains of rice. Plant. Prod. Sci. 2016, 1008, 30–42. [Google Scholar] [CrossRef]
  15. Kuang, B.; Mouazen, A.M. Non-biased prediction of soil organic carbon and total nitrogen with vis-NIR spectroscopy, as affected by soil moisture content and texture Keywords. Biosyst. Eng. 2013, 114, 249–258. [Google Scholar] [CrossRef] [Green Version]
  16. Nocita, M.; Stevens, A.; Toth, G.; Panagos, P.; Montanarella, L. Prediction of soil organic carbon content by diffuse reflectance spectroscopy using a local partial least square regression approach. Soil Biol. Biochem. 2014, 68, 337–347. [Google Scholar] [CrossRef]
  17. Zhang, Y.; Li, M.; Zheng, L.; Qin, Q.; Suk, W. Spectral features extraction for estimation of soil total nitrogen content based on modified ant colony optimization algorithm. Geoderma 2019, 333, 23–34. [Google Scholar] [CrossRef]
  18. Bao, Y.; Meng, X.; Ustin, S.; Wang, X.; Zhang, X.; Liu, H. Vis-SWIR spectral prediction model for soil organic matter with diff erent grouping strategies. Catena 2020, 195, 104703. [Google Scholar] [CrossRef]
  19. Debaene, G.; Nied, J.; Pecio, A.; Anna, Ż. Effect of the number of calibration samples on the prediction of several soil properties at the farm-scale. Geoderma 2014, 215, 114–125. [Google Scholar] [CrossRef]
  20. Lin, L.; Gao, Z.; Liu, X. Estimation of soil total nitrogen using the synthetic color learning machine (SCLM) method and hyperspectral data. Geoderma 2020, 380, 114664. [Google Scholar] [CrossRef]
  21. Johnson, J.; Vandamme, E.; Senthilkumar, K.; Sila, A.; Shepherd, K.D.; Saito, K. Near-infrared, mid-infrared or combined diff use reflectance spectroscopy for assessing soil fertility in rice fields in sub-Saharan Africa. Geoderma 2019, 354, 113840. [Google Scholar] [CrossRef]
  22. Dalal, R.C.; Henry, R.J. Simultaneous determination of moisture, organic carbon, and total nitrogen by near infrared reflectance spectrophotometry. Soil Sci. Soc. Am. J. 1986, 50, 120–123. [Google Scholar] [CrossRef]
  23. Zhang, Y.; Li, M.Z.; Zheng, L.H.; Zhao, Y.; Pei, X. Soil nitrogen content forecasting based on real-time NIR spectroscopy. Comput. Electron. Agric. 2016, 124, 29–36. [Google Scholar] [CrossRef]
  24. Jiang, Q.; Li, Q.; Wang, X.; Wu, Y.; Yang, X.; Liu, F. Estimation of soil organic carbon and total nitrogen in different soil layers using VNIR spectroscopy: Effects of spiking on model applicability. Geoderma 2017, 293, 54–63. [Google Scholar] [CrossRef]
  25. Conforti, M.; Matteucci, G.; Buttafuoco, G. Using laboratory Vis-NIR spectroscopy for monitoring some forest soil properties. J. Soils Sediments 2018, 18, 1009–1019. [Google Scholar] [CrossRef]
  26. Sudduth, K.A.; Hummel, J.W. Portable, Near-infrared spectrophotometer for rapid soil analysis. Trans. ASAE 1993, 36, 185–193. [Google Scholar] [CrossRef]
  27. Mouazen, A.M.; Alhwaimel, S.A.; Kuang, B.; Waine, T. Multiple on-line soil sensors and data fusion approach for delineation of water holding capacity zones for site specific irrigation. Soil Tillage Res. 2014, 143, 95–105. [Google Scholar] [CrossRef]
  28. Mouazen, A.M.; Maleki, M.R.; De Baerdemaeker, J.; Ramon, H. On-line measurement of some selected soil properties using a VIS–NIR sensor. Soil Tillage Res. 2007, 93, 13–27. [Google Scholar] [CrossRef]
  29. Zhou, P.; Li, M.; Yang, W.; Ji, R.; Meng, C. Development of vehicle-mounted in-situ soil parameters detector based on NIR diffuse reflection. Spectrosc. Spect. Anal. 2020, 40, 2856–2861. [Google Scholar]
  30. Zhou, P.; Zhang, Y.; Yang, W.; Li, M.; Liu, Z.; Liu, X. Development and performance test of an in-situ soil total nitrogen-soil moisture detector based on near-infrared spectroscopy. Comput. Electron. Agric. 2019, 160, 51–58. [Google Scholar] [CrossRef]
  31. An, X.; Li, M.; Zheng, L.; Liu, Y.; Sun, H. A portable soil nitrogen detector based on NIRS. Precis. Agric. 2014, 15, 3–16. [Google Scholar] [CrossRef]
  32. Tang, N.; Li, M.Z.; Sun, J.Y.; Zheng, L.H.; Pan, L. Development of soil-organic-matter fast-determination instrument based on spectroscopy. Spectrosc. Spectr. Anal. 2007, 27, 2139. [Google Scholar]
  33. Li, M.; Yao, X.; Yang, W.; Zhou, P.; Hao, Z.; Zheng, L. Design of New Portable Detector for Soil Total Nitrogen Content Based on High-power Tungsten Halogen Lamp and “One-Six” Special Optical Fiber. Trans. Chinese Soc. Agric. Mach. 2019, 50, 169–174. [Google Scholar]
  34. Wang, Y.P.; Lee, C.K.; Dai, Y.H.; Shen, Y. Effect of wetting on the determination of soil organic matter content using visible and near-infrared spectrometer. Geoderma 2020, 376, 114528. [Google Scholar] [CrossRef]
  35. Stenberg, B. Effects of soil sample pretreatments and standardised rewetting as interacted with sand classes on Vis-NIR predictions of clay and soil organic carbon. Geoderma 2010, 158, 15–22. [Google Scholar] [CrossRef] [Green Version]
  36. An, X.; Li, M.; Zheng, L.; Sun, H. Eliminating the interference of soil moisture and particle size on predicting soil total nitrogen content using a NIRS-based portable detector. Comput. Electron. Agric. 2015, 112, 47–53. [Google Scholar] [CrossRef]
  37. Rienzi, E.A.; Mijatovic, B.; Mueller, T.G.; Matocha, C.J.; Sikora, F.J.; Castrignanò, A. Prediction of Soil Organic Carbon under Varying Moisture Levels Using Reflectance Spectroscopy. Soil Sci. Soc. Am. J. 2014, 78, 958–967. [Google Scholar] [CrossRef]
  38. Kaiser, M.; Kleber, M.; Berhe, A.A. How air-drying and rewetting modify soil organic matter characteristics: An assessment to improve data interpretation and inference. Soil Biol. Biochem. 2015, 80, 324–340. [Google Scholar] [CrossRef]
  39. Tekin, Y.; Tumsavas, Z.; Mouazen, A.M. Effect of moisture content on prediction of organic carbon and pH using visible and near-infrared spectroscopy. Soil Sci. Soc. Am. J. 2012, 76, 188–198. [Google Scholar] [CrossRef]
  40. Wijewardane, N.K.; Ge, Y.; Morgan, C.L. Morgan Moisture insensitive prediction of soil properties from VNIR reflectance spectra based on external parameter orthogonalization. Geoderma 2016, 267, 92–101. [Google Scholar] [CrossRef] [Green Version]
  41. Barthès, B.G.; Brunet, D.; Hien, E.; Enjalric, F.; Conche, S.; Freschet, G.T.; d’Annunzio, R.; Toucet-Louri, J. Determining the distributions of soil carbon and nitrogen in particle size fractions using near-infrared reflectance spectrum of bulk soil samples. Soil Biol. Biochem. 2008, 40, 1533–1537. [Google Scholar] [CrossRef]
  42. Bogrekci, I.; Lee, W.S. Improving phosphorus sensing by eliminating soil particle size effect in spectral measurement. Trans. ASAE. 2005, 48, 1971–1978. [Google Scholar] [CrossRef]
  43. Lee, K.S.; Sudduth, K.A.; Drummond, S.T.; Lee, D.H.; Kitchen, N.R.; Chung, S.O. Calibration methods for soil property estimation using reflectance spectroscopy. Trans. ASABE 2010, 53, 675–684. [Google Scholar] [CrossRef]
  44. Yao, X.; Yang, W.; Li, M.; Zhou, P.; Chen, Y.; Hao, Z.; Liu, Z. Prediction of Total Nitrogen in Soil Based on Random Frog Leaping Wavelet Neural Network. IFAC-PapersOnLine 2018, 51, 660–665. [Google Scholar] [CrossRef]
  45. Yao, X.; Yang, W.; Li, M.; Zhou, P.; Liu, Z. Prediction of Total Nitrogen Content in Different Soil Types Based on Spectroscopy. IFAC-PapersOnLine 2019, 52, 270–276. [Google Scholar] [CrossRef]
  46. Mcdowell, W.H.; Magill, A.H.; Aitkenhead-Peterson, J.A.; Aber, J.D.; Merriam, J.L.; Kaushal, S.S. Effects of Chronic Nitrogen Amendment on Dissolved Organic Matter and Inorganic Nitrogen in Soil Solution. For. Ecol. Manag. 2004, 196, 29–41. [Google Scholar] [CrossRef]
  47. Jämtgård, S.; Näsholm, T.; Huss-danell, K. Soil Biology and Biochemistry Nitrogen Compounds in Soil Solutions of Agricultural Land. Soil Biol. Biochem. 2010, 42, 2325–2330. [Google Scholar] [CrossRef]
  48. Stark, J.M.; Hart, S.C. Diffusion Technique for Preparing Salt Solutions, Kjeldahl Digests, and Persulfate Digests for Nitrogen-15 Analysis. Soil Sci. Soc. Am. J. 1996, 60, 1846–1855. [Google Scholar] [CrossRef]
  49. Jia, S.; Zhang, J.; Li, G.; Yang, X. Predicting Soil Nitrogen and Organic Carbon Using Near Infrared SpectroscopyCoupled with Variable Selection. Appl. Eng. Agric. 2014, 30, 641–647. [Google Scholar]
  50. Zheng, L.; Li, M.; Pan, L.; Sun, J.; Tang, N. Application of wavelet packet analysis in estimating soil parameters based on NIR spectra. Spectrosc. Spectr. Anal. 2009, 29, 1549–1552. [Google Scholar]
  51. Zheng, L.; Li, M.; Pan, L.; Sun, J.; Tang, N. Estimation of soil organic matter and soil total nitrogen based on NIR spectroscopy and BP neural network. Spectrosc. Spectr. Anal. 2008, 28, 1160–1164. [Google Scholar]
  52. Zhou, P.; Sudduth, K.A.; Veum, K.S.; Li, M. Selection of characteristic wavebands to minimize soil moisture effects with in-situ soil spectroscopy. In 2020 ASABE Annual International Virtual Meeting; ASABE: Saint Joseph, MI, USA, 2020. [Google Scholar]
  53. Ogen, Y.; Faigenbaum-golovin, S.; Granot, A.; Shkolnisky, Y. Removing Moisture Effect on Soil Reflectance Properties: A Case Study of Clay Content Prediction. Pedosph. An. Int. J. 2019, 29, 421–431. [Google Scholar] [CrossRef]
  54. Jiang, Q.; Chen, Y.; Guo, L.; Fei, T.; Qi, K. Estimating Soil Organic Carbon of Cropland Soil at Different Levels of Soil Moisture Using VIS-NIR Spectroscopy. Remote Sens. 2016, 8, 755. [Google Scholar] [CrossRef] [Green Version]
  55. Castaldi, F.; Palombo, A.; Pascucci, S.; Pignatti, S.; Santini, F. Reducing the Influence of Soil Moisture on the Estimation of Clay from Hyperspectral Data: A Case Study Using Simulated PRISMA Data. Remote Sens. 2015, 7, 15561–15582. [Google Scholar] [CrossRef] [Green Version]
  56. Lobell, D.B.; Asner, G.P. Moisture Effects on Soil Reflectance. Soil Sci. Soc. Am. J. 2002, 66, 722–727. [Google Scholar] [CrossRef]
  57. Li, M.Z.; Han, D.H.; Wang, X. Spectral Analysis and Application; Science Press: Beijing, China, 2006. [Google Scholar]
  58. Lu, W.Z. Modern Near Infrared Spectroscopy Analytical Technology; China Petrochemical Press: Beijing, China, 2010. [Google Scholar]
  59. Hu, A.Q.; Yuan, H.F.; Song, C.F.; Li, X.Y. A correction method of baseline drift of discrete spectrum of NIR. Spectrosc. Spectr. Anal. 2014, 34, 2606–2611. [Google Scholar]
  60. Carle, C. Comments on Smoothing and Differentiation of Data by Simplified Least Square Procedure. Anal. Chem. 1972, 44, 1906–1909. [Google Scholar]
  61. Savitzky, A.; Golay, M.J. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  62. Morellos, A.; Pantazi, X.; Moshou, D.; Alexandridis, T.; Whetton, R.; Tziotzios, G.; Wiebensohn, J.; Bill, R.; Mouazen, A.M. Direct Special Issue: Proximal Soil Sensing Machine learning based prediction of soil total nitrogen, organic carbon and moisture content by using VIS-NIR spectroscopy. Biosyst. Eng. 2016, 152, 104–116. [Google Scholar] [CrossRef] [Green Version]
  63. Veum, K.S.; Parker, P.A.; Sudduth, K.A.; Holan, S.H. Predicting Profile Soil Properties with Reflectance Spectra via Bayesian Covariate-Assisted External Parameter Orthogonalization. Sensor 2018, 18, 3869. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flow chart of major steps adopted in this study. In the diagram, MACI refers to the moisture absorption correction index and PSCI is the particle size correction index.
Figure 1. Flow chart of major steps adopted in this study. In the diagram, MACI refers to the moisture absorption correction index and PSCI is the particle size correction index.
Remotesensing 13 00762 g001
Figure 2. Structure of the on-the-go detector. (a) Working scene of the on-the-go detector; (b) subsoiler and optical path of the detection unit.
Figure 2. Structure of the on-the-go detector. (a) Working scene of the on-the-go detector; (b) subsoiler and optical path of the detection unit.
Remotesensing 13 00762 g002
Figure 3. Soil sampling locations and treatments of the experiments at Field 3 (a) in Beijing (b) of China (c).
Figure 3. Soil sampling locations and treatments of the experiments at Field 3 (a) in Beijing (b) of China (c).
Remotesensing 13 00762 g003
Figure 4. Absorbance curves of soil samples with different soil moisture contents.
Figure 4. Absorbance curves of soil samples with different soil moisture contents.
Remotesensing 13 00762 g004
Figure 5. Corrected absorbance at seven discrete NIR spectral bands.
Figure 5. Corrected absorbance at seven discrete NIR spectral bands.
Remotesensing 13 00762 g005
Figure 6. Absorbance curves of soil samples with different soil particle sizes.
Figure 6. Absorbance curves of soil samples with different soil particle sizes.
Remotesensing 13 00762 g006
Figure 7. Standard deviation spectra. (a) Standard deviation spectra curve in the range 850–2500 nm; (b) standard deviation spectra curves in the ranges 1340–1420 nm and 1795–1955 nm.
Figure 7. Standard deviation spectra. (a) Standard deviation spectra curve in the range 850–2500 nm; (b) standard deviation spectra curves in the ranges 1340–1420 nm and 1795–1955 nm.
Remotesensing 13 00762 g007
Figure 8. Classification of soil particle size based on support vector machine (SVM).
Figure 8. Classification of soil particle size based on support vector machine (SVM).
Remotesensing 13 00762 g008
Figure 9. Corrected absorbance at six discrete NIR spectral bands.
Figure 9. Corrected absorbance at six discrete NIR spectral bands.
Remotesensing 13 00762 g009
Figure 10. Flowchart of the coupled elimination method.
Figure 10. Flowchart of the coupled elimination method.
Remotesensing 13 00762 g010
Figure 11. Comparison of elimination effects. (a) Estimation of TN concentration with raw discrete NIR spectral band data; (b) estimation of TN concentration with discrete NIR spectral band data through the proposed MACI correlation method to eliminate soil moisture disturbance; (c) estimation of TN concentration with discrete NIR spectral band data through the proposed PSCI method to eliminate soil particle size disturbance; (d) estimation of TN concentration with discrete NIR spectral band data through the proposed MACI and PSCI methods to eliminate soil moisture and particle size disturbances; (e) estimation of TN concentration with discrete NIR spectral band data on dried and ground soil samples.
Figure 11. Comparison of elimination effects. (a) Estimation of TN concentration with raw discrete NIR spectral band data; (b) estimation of TN concentration with discrete NIR spectral band data through the proposed MACI correlation method to eliminate soil moisture disturbance; (c) estimation of TN concentration with discrete NIR spectral band data through the proposed PSCI method to eliminate soil particle size disturbance; (d) estimation of TN concentration with discrete NIR spectral band data through the proposed MACI and PSCI methods to eliminate soil moisture and particle size disturbances; (e) estimation of TN concentration with discrete NIR spectral band data on dried and ground soil samples.
Remotesensing 13 00762 g011aRemotesensing 13 00762 g011b
Table 1. Accuracy of MACI classification and correction factors.
Table 1. Accuracy of MACI classification and correction factors.
Serial NumberSoil Moisture Content GradeMACIAccuracy Rate (%)Wj
1Low level (0–3.0)>9.111001
2Low-medium level (3.0–6.0)9.11–6.05880.92
3Middle level (6.0–10.0)6.05–2.03860.86
4High-medium level (10.0–13.0)2.03–1.12810.76
5High level (>13.0)<1.121000.68
Note: MACI = moisture absorption correction index.
Table 2. Statistics of 96 soil samples prepared. The range of soil total nitrogen (TN) represents the measured TN of soil samples prepared.
Table 2. Statistics of 96 soil samples prepared. The range of soil total nitrogen (TN) represents the measured TN of soil samples prepared.
Soil Particle SizeTN Level
g·kg−1
Range of TN
g·kg−1
Particle SizeTN Level
g·kg−1
Range of TN
g·kg−1
2.0 mm0 0.005–0.012 0.9 mm0 0.003–0.012
0.04 0.036–0.047 0.04 0.029–0.037
0.08 0.071–0.077 0.08 0.076–0.091
0.12 0.109–0.123 0.12 0.119–0.133
0.16 0.162–0.171 0.16 0.153–0.176
0.2 0.186–0.194 0.2 0.191–0.206
0.45 mm0 0.006–0.02 0.2 mm0 0.009–0.018
0.04 0.031–0.042 0.04 0.026–0.04
0.08 0.063–0.071 0.08 0.060–0.067
0.12 0.119–0.131 0.12 0.128–0.138
0.16 0.142–0.153 0.16 0.153–0.169
0.2 0.178–0.190 0.2 0.187–0.199
Note: TN = soil total nitrogen.
Table 3. Comparison of elimination efforts.
Table 3. Comparison of elimination efforts.
Elimination Method R v 2 RMSEP(g·kg−1)RPD
Model 10.640.2781.59
Model 20.710.2212.06
Model 30.730.2022.12
Model 40.820.1662.53
Model 50.890.1212.72
Note: R v 2 = coefficient of determination of validation; RMSEP = root mean square error of validation; RPD = residual prediction deviation; model 1 = raw discrete NIR spectral band data; model 2 = discrete NIR spectral band data using the proposed MACI correlation method to eliminate soil moisture disturbance; model 3 = discrete NIR spectral band data using the proposed PSCI method to eliminate soil particle size disturbance; model 4 = discrete NIR spectral band data using the proposed new coupled elimination method to eliminate soil moisture and particle size disturbances; model 5 = discrete NIR spectral band data on dried and ground soil samples.
Table 4. Validation results of soil total nitrogen (TN) concentration.
Table 4. Validation results of soil total nitrogen (TN) concentration.
NSpectral Measuring EquipmentSpectral RangeElimination FactorElimination MethodCalibration Algorithm R v 2 RMSEPReferences
48Portable detector940, 1050, 1100, 1200, 1300, and 1550 nm//BPNN0.450.215[36]
48Portable detector940, 1050, 1100, 1200, 1300, 1450, and 1550 nmSoil moisture and soil particle sizePMAI and mixed calibration setBPNN0.760.030[36]
90FT-NIR analyzer (MATRIX-I, Bruker corp., Germany)800–2564 nmSoil moistureWavelet decompositionsSVM0.810.053[23]
140AgroSpec portable VIS-NIR spectrophotometer (Tec5 Technology for Spectroscopy, Germany)305–2200 nmSoil moisture and soil particle sizeFD transformation with 31 smoothing points and SNVCubist method0.730.071[62]
708Veris P4000 (Veris Technologies, Inc., Salina, KS, USA)343–2202 nmSoil moistureEPOPLS0.630.024[63]
108On-the-go detector1070, 1130, 1245, 1375, 1550, and 1680 nm//SVM0.640.278This study
108On-the-go detector1070, 1130, 1245, 1361, 1375, 1450, 1550, 1680, and 1870 nmSoil moisture and soil particle sizeNew coupled elimination methodSVM0.820.166This study
Note: N = number of samples; R v 2 = coefficient of determination of validation; RMSEP = root mean scheme.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, P.; Yang, W.; Li, M.; Wang, W. A New Coupled Elimination Method of Soil Moisture and Particle Size Interferences on Predicting Soil Total Nitrogen Concentration through Discrete NIR Spectral Band Data. Remote Sens. 2021, 13, 762. https://doi.org/10.3390/rs13040762

AMA Style

Zhou P, Yang W, Li M, Wang W. A New Coupled Elimination Method of Soil Moisture and Particle Size Interferences on Predicting Soil Total Nitrogen Concentration through Discrete NIR Spectral Band Data. Remote Sensing. 2021; 13(4):762. https://doi.org/10.3390/rs13040762

Chicago/Turabian Style

Zhou, Peng, Wei Yang, Minzan Li, and Weichao Wang. 2021. "A New Coupled Elimination Method of Soil Moisture and Particle Size Interferences on Predicting Soil Total Nitrogen Concentration through Discrete NIR Spectral Band Data" Remote Sensing 13, no. 4: 762. https://doi.org/10.3390/rs13040762

APA Style

Zhou, P., Yang, W., Li, M., & Wang, W. (2021). A New Coupled Elimination Method of Soil Moisture and Particle Size Interferences on Predicting Soil Total Nitrogen Concentration through Discrete NIR Spectral Band Data. Remote Sensing, 13(4), 762. https://doi.org/10.3390/rs13040762

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop