Next Article in Journal
A Deep Learning-Based Robust Change Detection Approach for Very High Resolution Remotely Sensed Images with Multiple Features
Next Article in Special Issue
Analysis of Chlorophyll Concentration in Potato Crop by Coupling Continuous Wavelet Transform and Spectral Variable Optimization
Previous Article in Journal
Low-cost Fiberoptic Probe for Ammonia Early Detection in Fish Farms
Previous Article in Special Issue
Prediction of Early Season Nitrogen Uptake in Maize Using High-Resolution Aerial Hyperspectral Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Developing a Proximal Active Canopy Sensor-based Precision Nitrogen Management Strategy for High-Yielding Rice

1
International Center for Agro-Informatics and Sustainable Development (ICASD), College of Resources and Environmental Sciences, China Agricultural University, Beijing 100193, China
2
Agro-Geoinformatics Research Center (ARC), School of Surveying and Land Information Engineering, Henan Polytechnic University, Jiaozuo 454000, China
3
Precision Agriculture Center, Department of Soil, Water and Climate, University of Minnesota, St. Paul, MN 55108, USA
4
Beijing Agriculture Technology Extension Station, Beijing 100029, China
5
Center for Precision Agriculture, Norwegian Institute of Bioeconomy Research (NIBIO), Nylinna 226, 2849 Kapp, Norway
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(9), 1440; https://doi.org/10.3390/rs12091440
Submission received: 16 February 2020 / Revised: 18 April 2020 / Accepted: 29 April 2020 / Published: 2 May 2020
(This article belongs to the Special Issue Remote Sensing for Precision Nitrogen Management)

Abstract

:
RapidSCAN is a portable active canopy sensor with red, red-edge, and near infrared spectral bands. The objective of this study is to develop and evaluate a RapidSCAN sensor-based precision nitrogen (N) management (PNM) strategy for high-yielding rice in Northeast China. Six rice N rate experiments were conducted from 2014 to 2016 at Jiansanjiang Experiment Station of China Agricultural University in Northeast China. The results indicated that the sensor performed well for estimating rice yield potential (YP0) and yield response to additional N application (RIHarvest) at the stem elongation stage using normalized difference vegetation index (NDVI) (R2 = 0.60–0.77 and relative error (REr) = 6.2–8.0%) and at the heading stage using normalized difference red edge (NDRE) (R2 = 0.70–0.82 and REr = 7.3–8.7%). A new RapidSCAN sensor-based PNM strategy was developed that would make N recommendations at both stem elongation and heading growth stages, in contrast to previously developed strategy making N recommendation only at the stem elongation stage. This new PNM strategy could save 24% N fertilizers, and increase N use efficiencies by 29–35% as compared to Farmer N Management, without significantly affecting the rice grain yield and economic returns. Compared with regional optimum N management, the new PNM strategy increased 4% grain yield, 3–10% N use efficiencies and 148 $ ha−1 economic returns across years and varieties. It is concluded that the new RapidSCAN sensor-based PNM strategy with two in-season N recommendations using NDVI and NDRE is suitable for guiding in-season N management in high-yield rice management systems. Future studies are needed to evaluate this RapidSCAN sensor-based PNM strategy under diverse on-farm conditions, as well as to integrate it into high-yield rice management systems for food security and sustainable development.

Graphical Abstract

1. Introduction

As one of the major cereal crops in the world, more than half of the world’s population takes rice (Oryza sativa L.) as the staple food [1]. The area under rice cultivation in Asia accounts for 90% of the world’s total rice area [2]. At the same time, inappropriate nitrogen (N) fertilizer application rates and timing result in low N use efficiency (NUE) in this area [3]. Northeast China is a major rice production region in China and the abovementioned management problems are common [4,5,6]. Facing these challenges, Chinese agricultural scientists have developed regional optimum N management (RONM) systems, aiming to obtain higher yields with less resources and N losses suitable for different regions [4,6,7]. The RONM system using fixed N rates and timing optimum for a region may not be optimal for a specific site, year, and variety in that region [8,9]. Precision N management (PNM) strategies consider both spatial and temporal variability in soil N supply and crop N demand. They have the potential to further improve NUE over the RONM strategy [8].
Active crop canopy sensors have been increasingly used to develop in-season site-specific N management strategies, allowing non-destructive real-time diagnosis of crop N status and N recommendations. They have their own light sources and are not affected by environmental light conditions [8]. The GreenSeeker active canopy sensor (Trimble Navigation Limited, Sunnyvale, CA, USA) is a commonly used sensor for guiding in-season N management [8]. It has red (R) and near-infrared (NIR) spectral wavebands and two default vegetation indices (VI): normalized difference vegetation index (NDVI) and ratio vegetation index (RVI) [8]. A GreenSeeker sensor-based PNM strategy has been previously developed to improve NUE while maintaining rice yield in Northeast China [10]. In this strategy, basal and tillering N rates were the same as RONM, while panicle fertilizer rates at the stem elongation stage were adjusted based on N fertilizer optimization algorithm using the GreenSeeker sensor [10]. The key components of this strategy include in-season estimation of yield potential without additional topdressing N application (YP0) and N response index (RIHarvest) [10,11,12]. The potential yield with sufficient topdressing N application (YPN) can be estimated by multiplying YP0 and RIHarvest [11] and then the N topdressing requirement is estimated by dividing the yield increase (YPN−YP0) by the average NUE [10]. Because of the saturation problem of NDVI at moderate to high biomass conditions, the estimation of YP0 and RIHarvest using GreenSeeker NDVI was not very satisfactory across all stages [10] or at later growth stages (e.g., the heading stage) [13]. A previous study using the Crop Circle ACS 470 sensor (Holland Scientific, Inc., Lincoln, NE, USA) indicated that red edge-based VIs had the potential to overcome the NDVI saturation problem and improve the estimation of YP0 and RIHarvest, especially at later rice growth stages [13].
The RapidSCAN CS-45 sensor (Holland Scientific Inc., Lincoln, Nebraska, USA) is a relatively new alternative active crop canopy sensor available in the market. It is a lightweight and convenient portable sensor with built-in global positioning system and red edge (RE) band in addition to red and near infrared bands. Another advantage of the RapidSCAN sensor is that the sensor data collection is not influenced by measurement height in the range of 0.3 to 3 m [14]. It provides NDVI and normalized difference red edge (NDRE) as two default VIs, in addition to the R, RE, and NIR waveband reflectance. Besides NDVI and NDRE, many different VIs can be calculated. This sensor was found to perform well for estimating rice N status indicators at different growth stages [14]. Zhang et al. [15] reported that NDRE had a better rice yield prediction accuracy than NDVI from stem elongation to booting stage using the RapidSCAN sensor. More studies are needed to develop RapidSCAN sensor-based PNM strategies for rice.
To ensure both food security and agricultural sustainable development, integrated precision rice management systems have been developed to increase rice yield and NUE simultaneously [16]. In such systems, in addition to N applications before transplanting, at the tillering and stem elongation stages, grain N fertilizer was also applied at the heading stage to better meet the N demand of high-yielding rice. The previously developed GreenSeeker-based PNM strategies did not perform well to guide grain N fertilizer application at the heading stage because of the NDVI saturation [13]. The RapidSCAN sensor has the potential to overcome the saturation problem of NDVI and research is needed to develop a RapidSCAN sensor-based PNM strategy for high-yielding rice management systems that will guide topdressing N applications at both stem elongation and heading stages.
Therefore, the objectives of this study are to (i) evaluate the potential of in-season estimation of the rice yield potential and the response to N application at different growth stages using the RapidSCAN CS-45 sensor, (ii) develop a RapidSCAN sensor-based PNM strategy for high-yielding rice, and (iii) evaluate the RapidSCAN sensor-based PNM strategy for different varieties, N status, and years in Northeast China.

2. Materials and Methods

2.1. Study Site

The study was conducted in Sanjiang Plain, Heilongjiang Province, Northeast China (47.2°N, 132.6°E). The main soil type in this area is Albic soil, classified as Mollic Planosols in the FAO-UNESCO system, and typical Argialbolls in the soil taxonomy [17]. The study site is located in a cool-temperate sub-humid continental monsoon climate zone. The temperature ranges from −41 °C in the winter to 38 °C in the summer, with a mean annual temperature of 1.9 °C. About 72% of its annual precipitation (500–600 mm) occurs from June to September. The annual frost-free period is about 120–140 days long [16].

2.2. Calibration and Validation Experiments

Six plot experiments (Exp.) were conducted from 2014 to 2016 at Jiansanjiang Experiment Station of China Agricultural University, involving two different varieties, N rates and sensor-based N management strategies (Table 1). Each experiment had the same five N rates (0, 40, 80, 120, and 160 kg N ha−1). In addition, the experiments in 2015–2016 consisted of a sensor-based PNM treatment using the RapidSCAN sensor. The N fertilizer was applied in five N rate treatments (except the control treatment without N application) as three splits: 40% as basal N before transplanting, 30% at tillering stage, and the remaining 30% N at the stem elongation stage. According to previous studies [6,9], the N rate treatment for 120 kg N ha−1 was used as RONM system in this region. The sensor-based PNM treatments were also based on the RONM system, with the same basal (48 kg ha−1) and tillering (36 kg ha−1) N rates. The panicle and grain fertilizer rates were determined according to active canopy sensor-based N recommendation algorithm and applied at the stem elongation and heading stages, respectively.
Exp. 1 to 3 used Longjing 31, which is an 11-leaf variety requiring about 130 days to reach maturity. Exp. 4 to 6 used Longjing 21, which is a 12-leaf variety that needs about 133 days to maturity. All plot experiments were replicated three times in a randomized complete block design. The N source was granular urea. To evaluate the potential of the crop canopy sensors to estimate rice YP0 and RIHarvest at the stem elongation and heading stages, each plot of all experiments (except the 0 kg N ha−1 treatment) was divided into two parts: 4.5 × 9 m as the main plot and 2.5 × 9 m as the subplot without receiving the third N application. For all the treatments, 50 kg P2O5 ha−1 in the form of Ca(H2PO4)2 was applied before transplanting and 105 kg K2O ha−1 in the form of KCl was applied as two splits: 50% before transplanting and 50% at the stem elongation stage. Rice seedlings were prepared in a greenhouse and transplanted into the experimental fields in mid-May. The field and crop management in these experiments followed the regional recommendations.

2.3. Proximal Sensing Data Collection

The RapidSCAN CS-45 sensor was used to collect reflectance data in this study. The sensor with modulated light emitting diodes irradiates the crop canopy and determines a portion of the radiation reflected from the crop canopy, without being affected by ambient illumination. The internal polychromatic light source includes three spectral bands centered at R (670 nm), RE (730 nm), and NIR (780 nm) wavelengths. According to the manufacturer, the sensor has the unique feature of Pseudo Solar Reflectance measurements that are independent of height in the range of 0.3 m to 3 m. Considering the potential influence of viewing angle and measurement area on sensor readings, the sensor footprint was parallel to the plant rows with the beam of light being perpendicular to rice canopy about 0.7–0.9 m above the canopy. The sensor was carried at a consistent speed to collect sensor readings from four different rows (3 m per row) in the middle of each plot. The reflectance values were then averaged to represent the reflectance for each plot.
Following the methodology established in the previous study on the RapidSCAN sensor [14], fifty-one VIs were evaluated in this study for estimating YP0 and RIHarvest and the best performing VIs for calibration and validation are listed in Table 2. NDVI and NDRE were provided as two default indices for this sensor (see Table 2). Reflectance data were collected at stem elongation and heading stages, which were the key stages for panicle and grain fertilizer applications.

2.4. Plant Sampling and Measurements

At the stem elongation, heading, and maturity stages, 3 hills with tillers representative of each plot were randomly selected for assessing the aboveground biomass. After cleaning with water, all roots were removed. The plant samples were then oven dried for 30 min at 105 °C and then at 70 °C until constant weight, and weighed to determine their biomass. They were later ground to pass a 0.5 mm sieve. Plant N concentration was determined using the Kjeldahl-N method.
Rice was harvested at the end of September or early October. Grain yield was determined by hand harvesting three 1 m2 areas in each plot where spectral reflectance data were collected. Grains were separated from straw using a small grain thresher and then weighed. Grain moisture was determined immediately after weighing. The rice grain weight was adjusted to a moisture content of 140 g kg−1.
Agronomic efficiency of N (AEN) and partial factor productivity of N (PFPN) were calculated using the following equations:
AE N   ( kg   kg 1 ) = Grain   yield     Grain   yield   at   control N   rate × 100
PFP N (   kg   kg 1 ) = Grain   yield N   rate

2.5. Development and Evaluation of RapidSCAN-Based Precision Nitrogen Management Strategies

Based on Yao et al. [10], the RapidSCAN-based PNM strategy in this study was developed by first establishing the models to estimate YP0 and RIHarvest using in-season estimate of yield (INSEY) and in-season N response index based on VI (RI-VI), respectively. INSEY can be regarded as an estimate of average daily biomass production from the time of transplanting to the day of sensing [11]. It was calculated as NDVI divided by the number of growing degree days > 0 [10]. In this study, however, the number of days from transplanting to sensing was used instead of growing degree days to calculate INSEY, similarly to the method of Cao et al. [13]. With respect to their study, the selected RapidSCAN VIs were used here to replace the GreenSeeker NDVI or RVI. RIHarvest indicates the actual crop yield response to additional N within a given year [28,29] and was calculated as follows [13]:
RI Harvest = Yield _ N rich Yield _ CK
where Yield_Nrich is the average yield of plots receiving sufficient N application (the 160 kg N ha−1 treatment in this study), and Yield_CK is the average yield of plot without receiving the third N application at the stem elongation stage or the fourth N application at the heading stage.
RI-VI was calculated in the same way as RIHarvest, with the exception that VIs derived from RapidSCAN sensor were used instead of yield. YPN was calculated by multiplying YP0 and RIHarvest. Finally, the N topdressing requirement is estimated by dividing the yield gap (YPN − YP0) by the average AE of topdressing N (AEtopdressing) [10]. The AEtopdressing should be higher than the one for the whole season, and will be predicted during the growing season using the predicted RIHarvest [30].
To ensure sufficient N supply for grain filling and higher NUE in high-yield rice management systems, a strategy for in-season site-specific N management of rice using RapidSCAN at stem elongation and heading stages was developed in this study (Figure 1). First, the topdressing N application rate (Nrate) at stem elongation stage was determined as mentioned above, and then this rate was split in two doses, 2/3 as panicle fertilizer at stem elongation stage (SE_Nrate) and 1/3 for grain fertilizer at the heading stage. Second, the RapidSCAN sensor was used to estimate the potential yield with added N application at the heading stage (HD_YPN). The difference between estimated YPN at stem elongation and heading stages (HD_YPN - SE_YPN) was used to adjust the remaining 1/3 Nrate to match the crop N demand at the heading stage. Therefore, the recommended N topdressing application rate at the heading stage (HD_Nrate) can be determined as follows (Figure 1):
HD _ N rate = HD _ YP N     SE _ YP N AE topdressing + 1 3   N rate
where HD_Nrate is the topdressing N application rate at heading stage, HD_YPN is the predicted yield potential with topdressing N application at heading stage, SE_YPN is the predicted yield potential with topdressing N application at the stem elongation stage, and AEtopdressing is the topdressing N agronomic efficiency.
In addition, two restrictions were applied after considering the rice production situation in Northeast China:
YP N     YP max   ( 12   t   ha 1 )
N min ( 0   kg   ha 1 )     N rate     N max   ( 48   kg   ha 1 )
where YPmax is the maximum obtainable yield, Nmin and Nmax are the minimum and maximum topdressing N rates.
According to the definition and the methods of Raun et al. [11] and [12], the YPmax in the study region was set to 12 t ha−1 based on previous studies and farmer survey data in this region [10,13,16,31]. The topdressing N application rates at the stem elongation stage was set to 0 to 48 kg ha−1 based on farmer surveys and previous studies in this region [10,16,31].
In order to determine whether the restrictions applied to the rice PNM strategies were suitable, three N rate treatments (80, 120, and 160 kg ha−1) from Exp. 1–6 were chosen to represent three rice N status before topdressing (deficient, optimal, and surplus, respectively). They were used to evaluate the RONM and the developed PNM strategies by calculating the differences between economically optimum N rate (EONR) and N rates recommended by RONM or PNM strategies. In order to evaluate the potential of the developed PNM strategy, the RapidSCAN-based PNM treatment in Exp. 2–3 and 5–6 was compared for yield, N rate, and NUE with the control treatment (0 kg N ha−1), the 160 kg N ha−1 treatment reflecting the farmer N management (FNM) and the 120 kg N ha−1 identical with the RONM. For the RapidSCAN sensor-based PNM treatment, the topdressing N rate was estimated based on the PNM strategy developed in this study using NDVI at the stem elongation stage and NDRE at the heading stage (Exp. 2 and 5 using data up to 2015; Exp. 3 and 6 using data up to 2016).
Economic return to N (E, $ ha−1) was used to evaluate the profitability of different N management systems, and was calculated as follows:
E = ( Y N Y 0 ) × P Y     N total × P N
where YN (kg ha−1) is the rice grain yield with N application, Y0 (kg ha−1) is the rice grain yield of the check treatment without any N application, PY is rice grain price (0.44 $ kg−1). Ntotal is the total N fertilizer application rate (kg ha−1). PN is the N fertilizer price (0.54 $ kg−1).

2.6. Statistical Analysis

Data collected from the three-year experiments were pooled together and then randomly divided into calibration dataset (67% of the observations) and validation dataset (33% of the observations) for the estimation of YP0 and RIHarvest using RapidSCAN sensor. The coefficients of determination (R2) for the relationships between VIs and agronomic parameters were calculated using SPSS 18.0 (SPSS Inc., Chicago, Illinois, USA), and the models with the highest R2 were selected. In addition to R2, the performance of the models for predicting YP0 and RIHarvest was also evaluated using the root mean square error (RMSE) and relative error (REr). Analysis of variance were conducted using the SAS software package Version 9.0 (SAS Institute Inc., Cary, NC, USA). The means for treatments were compared with least significant difference (LSD) test at the 0.05 probability level (at p < 0.05).

3. Results

3.1. Changes in NDVI vs. NDRE among Different N Rates, Varieties, Stages, and Years

Rice grain yield was significantly affected by the factors of N rates, varieties, and years, and RapidSCAN-based NDVI and NDRE also showed similar results (Table 3). Except for the variation of year, NDVI was significantly affected across different N rates, varieties, and growth stages. NDRE was significantly affected by these factors, except for variety. The changes in NDVI vs. NDRE were also shown in Figure 2. YP0, NDVI, and NDRE all increased with N rates (Figure 3). The average NDVI showed significant difference between Longjing 31 (0.68) and Longjing 21 (0.72) across years, growth stages, and N levels.

3.2. Yield Without Additional Topdressing N Application

The performance of the INSEY calculated with NDVI, NDRE, and best performing VIs to estimate rice YP0 varied with different growth stages across N rate treatments, sites, and years (Table 4 and Figure 2). At the stem elongation growth stage, two varieties performed similarly. INSEY-NDVI explained 72–76% of YP0 variability, which was better than INSEY-NDRE (60–66%). The INSEY calculated with best performing VIs (nonlinear index (NLI) (INSEY_NLI), modified simple ratio (MSR) (INSEY_MSR), and modified soil-adjusted vegetation index (MSAVI) (INSEY_MSAVI)) performed similarly (R2 = 0.74–0.78) to INSEY_NDVI.
At the heading stage, however, INSEY_NDVI did not perform very well, explaining only 28-59% of the YP0 variability. Moreover, INSEY_NDRE performed consistently better than INSEY_NDVI, explaining 73–89% of the YP0 variability (Table 4). The INSEY calculated with red edge optimal soil adjusted vegetation index (REOSAVI) (INSEY_REOSAVI), modified red edge soil adjusted vegetation index (MRESAVI) (INSEY_MRESAVI), and NDVI*RVI (INSEY_NDVI*RVI), did not perform significantly better than INSEY_NDRE (Table 4). Furthermore, the YP0 was better estimated using INSEY for Longjing 31 than for Longjing 21.
The validation results were similar to the calibration results (Table 4). At the stem elongation stage, the INSEY_NDVI and INSEY calculated with best performing VI had similar performance for predicting YP0, with R2, RMSE, and REr of 0.65–0.66, 0.47, and 7.4% across both varieties. At the heading stage, INSEY_NDRE and INSEY calculated with best performing VI performed similarly for predicting YP0, with R2, RMSE, and REr of 0.70, 0.65–0.66, and 8.6–8.7%, respectively.

3.3. The Responsiveness to Topdressing N Application

The performance of the RI_VI calculated with NDVI (RI_NDVI), NDRE (RI_NDRE), and best performing VIs to estimate rice RIHarvest varied with growth stages across N rate treatments, sites, and years (Table 5 and Figure 4). At the stem elongation stage, RI_NDRE and RI calculated with best performing VIs (optimal vegetation index (VIopt), red edge wide dynamic range vegetation index (REWDRVI) and RVI) did not perform significantly better than RI_NDVI (R2 = 0.67–0.78). The RIHarvest was estimated better for Longjing 31 (R2 = 0.71–0.79) than for Longjing 21 (R2 = 0.68–0.71) or across varieties (R2 = 0.64–0.68). The validation results showed similar pattern to the calibration results (Table 5). The prediction of RIHarvest for Longjing 31 (R2 = 0.72–0.78, RMSE = 0.08–0.09 and REr = 6.1–6.8) was better than for Longjing 21 (R2 = 0.60–0.67, RMSE = 0.10–0.11 and REr = 7.0–7.4) or across varieties (R2 = 0.61–0.63, RMSE = 0.11 and REr = 7.8–8.0).
At the heading stage, RI_NDRE and RI calculated with best performing VIs (normalized NIR index (RI_NNIRI), red edge optimal VI (RI_REVIopt), REOSAVI (RI_REOSAVI)) performed similarly for the prediction of RIHarvest (Table 5). They worked better than RI_NDVI (R2 = 0.79–0.85 vs. R2 = 0.67–0.79) for Longjing 21 or across varieties based on the calibration results. However, they performed similarly for Longjing 31 (R2 = 0.92–0.95). With validation, RI_NDRE (R2 = 0.75–0.82, RMSE = 0.10 and REr = 8.3–8.7) and RI calculated with best performing VIs (R2 = 0.73–0.82, RMSE = 0.10 and REr = 8.3–8.7) were better than RI_NDVI (R2 = 0.60–0.65, RMSE = 0.12–0.16 and REr = 10.2–13.4) for a specific variety or across varieties (Table 5). The RI calculated with best performing VIs performed better at the heading stage than the stem elongation stage for either calibration or validation.

3.4. In-Season Prediction of Nitrogen Use Efficiency

AEtopdressing varied significantly between Longjing 31 (39.5 kg kg−1) and Longjing 21 (59.0 kg kg−1) (Figure 5a). RIHarvest had a strong positive relationship with AEtopdressing, with the correlation being stronger for Longjing 21 (R2 = 0.76) than for Longjing 31 (R2 = 0.61) (Figure 5b). Therefore, in-season predicted RIHarvest was used to predict AEtopdressing for side-dress N fertilizer recommendations in this study.

3.5. Evaluating Different Precision Nitrogen Management Strategies Under Variable Nitrogen Status Using Scenario Analysis

Based on the abovementioned results (Table 4 and Table 5), different PNM strategies were developed, as explained in Methods and Figure 3. To evaluate the performance of these strategies under different N status, the N rate treatments (80, 120, and 160 kg ha−1) in Exp.1–6 were selected for scenario analysis to determine the difference between recommended N rates (ΔNrate) based on PNM strategies and EONR calculated using the N responses in each variety-year (Figure 6, Figure 7 and Figure 8). The RONM strategy did not consider the variation of years, varieties, and rice N status, and used a fixed N topdressing rate (36 kg N ha−1). Four RapidSCAN sensor-based PNM strategies were evaluated. The tested PNM strategies all had the same basal and tillering N application rates as RONM, but panicle fertilizer at the stem elongation stage was recommended using NDVI or the best performing VIs or panicle and grain fertilizer rates at the stem elongation and heading stages were recommended using NDVI and NDRE or the best performing VIs.
There were significant differences between RONM and N rates recommended by different RapidSCAN sensor-based PNM strategies under different years, varieties, and rice N status. The N topdressing rate of RONM was consistently lower than EONR under deficient N status (Figure 6), but higher than EONR under surplus N conditions (Figure 8). Under relatively optimum N conditions, RONR was slightly higher than EONR for Longing 31, but lower for Longjing 21 (Figure 7).
All the four PNM strategies recommended higher N rates than RONM under deficient N conditions and lower N rates under surplus N conditions, while under relatively optimal N conditions, the PNM strategies would recommend slightly higher, lower or similar N rates depending on the year, variety, and PNM strategy. In a specific year across different N conditions, the ΔNrate for the PNM strategies was in a range of 0–15 kg N ha−1. Among all the four PNM strategies, the PNM strategies based on in-season N recommendations at both stem elongation and heading stages (with 3-year cumulative ΔNrate of 40-50 kg N ha−1) performed better than the PNM strategies making N recommendation only at the stem elongation stage for Longjing 21 under different N conditions (with 3-year cumulative ΔNrate of 68-69 kg N ha−1). For Longjing 31, the PNM strategies making N recommendations at both stem elongation and heading stages performed much better under deficient N conditions, with 3-year cumulative ΔNrate (46-48 kg N ha−1) being lower than the PNM strategies making N recommendations only at the stem elongation stage (63-64 kg N ha−1). Irrespective of the varieties, using NDVI and NDRE for panicle and grain N fertilizer recommendations had better or similar performance as the best performing VIs.

3.6. Evaluation Experiments

Based on the above results, we chose NDVI and NDRE for panicle and grain fertilizer recommendation to evaluate the potential of RapidSCAN sensor-based PNM strategy in the evaluation experiments. These experiments included the 0 kg N ha−1 treatment as control (CK), the 160 kg N ha−1 treatment as FNM treatment, the 120 kg N ha−1 treatment as RONM treatment, and RapidSCAN-based PNM treatment in Exp. 2–3 and 5–6 (Table 6). The RapidSCAN sensor-based PNM treatment in evaluation experiments had the same N rate as RONM (N120) before topdressing and made two in-season adjustments at the stem elongation and heading growth stages.
Making two in-season N recommendations to better meet rice N requirements, the RapidSCAN-based PNM strategy recommended different N rates for different years and varieties. For Longjing 31, the total recommended N rates by the PNM strategy were 118–122 kg N ha−1, which were similar to the rates of the RONM strategy but 24–26% lower than the N rates of the FNM strategy. For Longjing 21, the total N rates recommended by the PNM strategy were similar to or higher than the rates given by the RONM strategy but 21–26% lower than the N rate of the FNM strategy (Table 6).
The RapidSCAN-based PNM strategy resulted in higher yield than RONM strategy across years and varieties. For Longjing 31, the yield of the PNM strategy was not significantly different from FNM and RONM strategies in a specific year, but it was significantly higher (5%) than the RONM strategy across two years. For Longjing 21, the PNM strategy significantly increased grain yield by 4% with respect to the RONM strategy in 2016 or across two years, but the increase was not significantly different from the FNM strategy.
The PNM strategy led to the highest AEN and PFPN across years and varieties. For Longjing 31, the PNM strategy increased AEN and PFPN by 14% and 5% over the RONM strategy and both by 33% over the FNM strategy across two years, respectively. For Longjing 21, the PNM strategy increased AEN by 33% and PFPN by 32% over the FNM strategy across two years. Compared with RONM, the PNM strategy increased AEN by 7%, but did not improve PFPN significantly. Across varieties and years, the PNM strategy increased AEN and PFPN by 32–33% over FNM and by 3–10% over RONM.
The PNM strategy consistently increased economic returns compared with RONM by 101–168 $ ha−1 for specific year-variety combination or across years and varieties. The PNM strategy increased economic returns by 0-91 $ ha−1 over the FNM strategy, but none of the increase was statistically significant.

4. Discussion

4.1. Selection of Vegetation Indices for Precision N Management Strategy Development

NDVI and NDRE are two widely used VIs in precision agriculture and have good relationships with aboveground biomass and plant N uptake [14,32]. They are the default VIs of the RapidSCAN sensor, and increased with N application rates in this study, which was in agreement with Aranguren et al. [33]. Besides NDVI and NDRE, many other VIs can be calculated from the three wavebands of RapidSCAN [13,34]. It is important to know if any other VI will perform consistently better than the default VIs and should be selected for PNM strategy development.
In this study, we found that the two default indices (NDVI and NDRE) were sufficient for estimating key parameters to develop the PNM strategy at different key growth stages and other VIs would not be needed.
The stem elongation stage is the key stage to apply panicle fertilizer for rice in Northeast China. The rice canopy was not closed at this stage, and the biomass was not very high yet. The NDVI values were not saturated and achieved similar performance for predicting YP0 (R2 = 0.66-0.70; REr = 6.3-7.7%) and RIHarvest (R2 = 0.60–0.77; REr = 6.2–8.0%) as best performing VIs (R2 = 0.63–0.70 and 0.63–0.78; REr = 6.5–7.6% and 6.1–7.8%, respectively). The recommended N rates using NDVI were also similar to the best performing VIs across both varieties based on the scenario analysis results. This agreed with many studies achieving good results with GreenSeeker NDVI-based PNM systems at early growth stages [10,12,34,35,36,37,38].
The heading stage is the key stage to apply grain fertilizer for rice in Northeast China in high-yield rice management systems [13,16]. GreenSeeker NDVI become saturated at this stage and cannot be used to make N recommendations [10,13,16,34,39]. As a result, the PNM strategy making two in-season N recommendations at stem elongation and heading stages only using NDVI had the worst performance among all the PNM strategies evaluated in the scenario analysis. One approach to overcome the NDVI saturation problem is to combine NDVI with relative plant height data [30]. Another approach is to replace NDVI with other VIs that can overcome the saturation effect [32]. The RapidSCAN sensor has a RE band in addition to the R and NIR bands, which can be used to calculate more RE-based VIs to overcome the saturation problems of NDVI [13,14,40]. In this study, good relationships between RE-based VIs and YP0 or RIHarvest were found at the heading stage, with R2 being 0.75–0.95 for calibration, and R2 and REr being 0.70–0.82 and 7.3–9.0% for validation, respectively. Encouragingly, the NDRE index achieved similar performance for predicting YP0 (R2 = 0.70–0.76; REr = 7.3–8.6%) and RIHarvest (R2 = 0.75–0.82; REr = 8.3–8.7%) as best performing VIs (R2 = 0.70–0.76 and 0.73–0.82; REr = 7.3–9.0% and 8.3–8.7%, respectively). Furthermore, the N rates recommended by the PNM strategy using NDVI at the stem elongation stage and NDRE at heading stage were very similar to the N rates recommended by the PNM strategy using best performing VIs at the stem elongation and heading stages.

4.2. In-Season Prediction of Agronomic Efficiency of Topdressing Nitrogen

The NUE is an important parameter for developing crop sensing-based PNM strategies, and can directly influence the final recommended N rates. Previous crop sensor-based PNM strategies generally used a constant NUE value. Raun et al. [12] used the expected recovery efficiency to calculate the N topdressing rate for wheat with the range of 0.5 to 0.7 in America. Cao et al. [34] used 0.4 as the expected recovery efficiency in developing PNM strategy for winter wheat in North China Plain. According to the multi-site-year data in N plot experiments, Yao et al. [10] used fixed AEtopdressing (26.79 kg kg−1) for the development of rice PNM strategy.
Studies indicated that NUE is influenced by environmental factors (e.g., soil fertility, water supply, sunshine, accumulated temperature, etc.) and the plant genotypes (e.g., ability to absorb N, photosynthetic rate, stress resistance, root system structure, etc.) [41,42,43,44]. Wang et al. [30] considered the influence of early N application and soil types on AEN and used in-season predicted RIHarvest with GreenSeeker sensor to predict AEN for developing crop sensor-based PNM strategy for spring maize in Northeast China. Following their idea, this study found that early N application and rice variety were two key factors influencing the AEtopdressing of rice in the study region. If early season N application rate is low, N deficiency will lead to a higher efficiency of N use, while a high N application rate during the early growth stage will result in N surplus conditions before topdressing, which will lead to a lower use efficiency of the topdressing N fertilizer. RIHarvest is an indicator of crop N response to additional N application. If early season N status is deficient, the crop will be more responsive to additional topdressing or side-dressing N application. On the other hand, if the early season N status is surplus, the crop will not be responsive to additional N application, and the NUE for the topdressing or side-dressing N application will be lower. Therefore, the RIHarvest and AEtopdressing should be positively correlated, as found in this study (R2 = 0.61 ~ 0.76). This finding was consistent with the results of Wang et al. [30] for spring maize (R2 = 0.72 ~ 0.74). On the other hand, the big difference of AEtopdressing between Longjing 31 (39.5 kg kg−1) and Longjing 21 (59.0 kg kg−1) is shown in Figure 5a. The variety with longer growing period will require more N supply and the photosynthetic utilization will be higher. Therefore, in-season prediction of AEtopdressing can be based on crop N status using in-season active sensor predicted response index and variety- or variety group-specific models should be developed.

4.3. Evaluation of Different Precision N Management Strategies

In scenario analysis, the four PNM strategies for N rate recommendations performed much better than the RONM strategy, especially under deficient or surplus N status. The PNM strategies making two in-season N recommendations at both stem elongation and heading stages performed better than PNM strategies making only one in-season N recommendation at the stem elongation stage, and their recommended N rates were 15-28 kg N ha−1 closer to EONR under deficient N status across years and varieties. For the longer-growing variety (Longjing 21), the PNM strategy making two N recommendations performed consistently better under different N conditions, being 2-28 kg N ha−1 closer to EONR across years. However, no significant difference was found under optimal or surplus N status for Longjing 31, the shorter-growing variety. One of the reasons might be that Longjing 21 typically required a higher N rate and more time (3-7 days) to booting and grain filling than Longjing 31 [9]. The PNM strategy with two in-season N recommendations can give us two opportunities to adjust N application rates according to in-season crop growth and weather conditions, and can better meet rice N needs to increase grain yield. This is more important for the long-season variety Longjing 21 than the shorter-season variety of Longjing 31 [45].
In the evaluation experiments, the RapidSCAN-based PNM strategy could achieve similar rice yield as the FNM strategy, save 24% N fertilizers, and increase NUE by an average of 32–33% across years and varieties (Table 6). Because of the low N fertilizer prices under agricultural subsidies by the government in China, the difference of economic return between PNM and FNM was small (average 36 $ ha−1 annually) across varieties. These results agreed with the results of Yao et al. [10], who reported that PNM strategies based on GreenSeeker or chlorophyll meter increased PFPN by 48-65% without significantly increasing the yield in comparison to FNM.
Although the total N rates of PNM and RONM strategies in the evaluation experiments were not significantly different, the PNM strategy allowed us to better adjust N distribution for different varieties and plant N status, resulting in the increased grain yield, AEN, PFPN, and economic return by an average of 4%, 10%, 3%, and 148 $ ha−1,respectively (Table 6). Zhao et al. [16] found that a chlorophyll meter-based PNM with two in-season N recommendations at the same growth stages as the PNM strategy developed in this study had higher N accumulation at later growth stages to achieve higher yield, and larger panicle size and grain fill percentage. Wang et al. [45] integrated a GreenSeeker sensor-based PNM strategy making one in-season N recommendation at stem elongation stage into a high-yield management system (with optimized transplanting density and water management) increased rice grain yield and NUE by 10% and 1–33% over regional optimum rice management, respectively. Future studies are needed to integrate the RapidSCAN sensor-based PNM strategy developed in this study into similar high-yield management systems to further evaluate its potential to simultaneously improve rice yield, NUE, and economic returns over the regional optimum management system under diverse on-farm conditions.

5. Conclusions

The RapidSCAN sensor has the potential to predict rice grain yield and response to additional N application in Northeast China at the stem elongation and heading stages during the growing season, with R2 and REr being 0.63–0.82 and 6.1–9.0%, respectively. At the stem elongation stage, NDVI could be used to predict YP0 and RIHarvest (R2 = 0.66–0.70 and 0.60–0.77; REr = 6.3–7.7% and 6.2–8.0%) with similar performance as the best performing indices (R2 = 0.63–0.70 and 0.63–0.78; REr = 6.5–7.6% and 6.1–7.8%, respectively). At the heading stage, NDRE could be used to predict YP0 and RIHarvest (R2 = 0.70–0.76 and 0.75–0.82; REr = 7.3–8.6% and 8.3–8.7%, respectively) as accurately as the best performing indices (R2 = 0.70–0.76 and 0.73–0.82; REr = 7.3–9.0% and 8.3–8.7%, respectively). The combination of the two default indices (NDVI and NDRE) of the RapdiSCAN sensor was sufficient for rice PNM strategy development without the need of incorporating other indices. Across the two tested varieties under deficient N conditions or for Longjing 21 under different N conditions, the PNM strategy making two in-season N recommendations performed better than the PNM strategies making only one in-season N recommendation at the stem elongation stage. The RapidSCAN sensor-based PNM strategy with panicle and grain fertilizer recommendations could lead to similar yield as FNM, save 24% N application rate, and increase AEN and PFPN by an average of 33% and 32% across years and varieties, respectively. Compared with RONM, the PNM strategy increased grain yield, AEN, PFPN, and economic return by an average of 4%, 10%, 3%, and 148 $ ha−1 across years and varieties, respectively. Future studies are needed to further evaluate the RapidSCAN sensor-based PNM strategy under diverse on-farm conditions, and integrate it into high-yield rice management systems.

Author Contributions

Conceptualization, J.L. (Junjun Lu) and Y.M.; methodology and validation, Y.M., J.L. (Junjun Lu), and X.W.; investigation, J.L. (Junjun Lu), W.S., J.L. (Jingxin Li), and X.H.; resources: Y.M., and J.L. (Junjun Lu); writing—original draft preparation, J.L. (Junjun Lu); writing—review and editing, Y.M., K.K. and Z.C.; supervision, Y.M.; project administration and funding acquisition, Y.M., J.L. (Junjun Lu), and K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key National Research and Development Program (2016YFD0200600; 2016YFD0200602), Norwegian Ministry of Foreign Affairs (SINOGRAIN II, CHN-17/0019), Doctoral Foundation of Henan Polytechnic University (B2019-5), Key Scientific Research Projects of Universities Supported by Foundation of Henan Educational Committee (20A210013), and Key Scientific and Technological Project of Henan Province (192102310038).

Acknowledgments

The kind assistance and support provided for this research by Yuan Gao, Wen Yang, Yong Yang, Haisheng He, Honglin Li, Fengyan Liu, and Xiaolong Wang at Jiansanjiang Institute of Agricultural Research and Mr. Guojun Li from Jiansanjiang Department of Agriculture, Jiansanjiang Branch Bureau of Agricultural Reclamation are highly appreciated.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AbbreviationEnglish Full NameUnit
AEtopdressingThe topdressing nitrogen agronomic efficiency%
AENAgronomic efficiency of nitrogen%
EONREconomically optimum nitrogen ratekg N ha−1
FNMFarmer nitrogen management-
HD_NrateThe recommended nitrogen topdressing application rate at the heading stagekg N ha−1
HD_YPNThe potential yield with added nitrogen topdressing application at the heading stagekg ha−1
INSEYIn-season estimate of yield-
NNitrogen-
NDRENormalized difference red edge-
NDVINormalized difference vegetation index-
NIRNear infrared-
NrateThe recommended N topdressing application ratekg N ha−1
NUENitrogen use efficiency-
PFPNPartial factor productivity of nitrogenkg kg−1
PNMPrecision nitrogen management-
Q, E or PThe quadratic, exponential, or power fit-
RRed-
R2The coefficients of determination-
RERed edge-
RErRelative error%
RIHarvestYield responsiveness to additional nitrogen fertilizer applications-
RI-VIIn-season nitrogen response index based on vegetation index-
RMSEThe root mean square errorSame units as statistical data used
RONMRegional optimum nitrogen management-
RVIRatio vegetation index-
SE_NrateThe recommended nitrogen topdressing application rate at the stem elongation stagekg N ha−1
SE_YPNThe potential yield with added nitrogen topdressing application at the stem elongation stagekg ha−1
SE-NDVIThe recommended nitrogen rate using normalized difference vegetation index at the stem elongation stage kg N ha−1
VI0Vegetation index at plots without additional nitrogen topdressing application-
VINVegetation index at plots with sufficient nitrogen fertilization plots-
VIsVegetation indices-
YP0The yield potential without additional nitrogen applicationkg ha−1
YPmaxThe maximum obtainable yieldkg ha−1
YPNThe potential yield with added nitrogen fertilizationkg ha−1
ΔNrateThe difference between recommended nitrogen rate and economically optimum nitrogen ratekg N ha−1

References

  1. Cantrell, R.P.; Reeves, T.G. The cereal of the world’s poor takes center stage. Science 2002, 296, 53. [Google Scholar] [CrossRef] [PubMed]
  2. Normile, D. Reinventing rice to feed the world. Science 2008, 321, 330–333. [Google Scholar] [CrossRef] [PubMed]
  3. Foley, J.A.; Ramankutty, N.; Brauman, K.A.; Cassidy, E.S.; Gerber, J.S.; Johnston, M.; Mueller, N.D.; O’Connell, C.; Ray, D.K.; West, P.C.; et al. Solutions for a cultivated planet. Nature 2011, 478, 337–342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zhang, F.; Chen, X.; Vitousek, P. Chinese agriculture: An experiment for the world. Nature 2013, 497, 33–35. [Google Scholar] [CrossRef] [PubMed]
  5. Peng, X.; Yang, Y.; Yu, C.; Chen, L.; Zhang, M.; Liu, Z.; Sun, Y.; Luo, S.; Liu, Y. Crop management for increasing rice yield and nitrogen use efficiency in Northeast China. Agron. J. 2015, 107, 1682–1690. [Google Scholar] [CrossRef]
  6. Wu, L.; Chen, X.; Cui, Z.; Wang, G.; Zhang, W. Improving nitrogen management via a regional management plan for Chinese rice production. Environ. Res. Lett. 2015, 10, 095011. [Google Scholar] [CrossRef]
  7. Chen, X.; Cui, Z.; Fan, M.; Vitousek, P.; Zhao, M.; Ma, W.; Wang, Z.; Zhang, W.; Yan, X.; Yang, J. Producing more grain with lower environmental costs. Nature 2014, 514, 486–489. [Google Scholar] [CrossRef]
  8. Cao, Q.; Miao, Y.; Feng, G.; Gao, X.; Liu, B.; Liu, Y.; Li, F.; Khosla, R.; Mulla, D.J.; Zhang, F. Improving nitrogen use efficiency with minimal environmental risks using an active canopy sensor in a wheat-maize cropping system. Field Crop. Res. 2017, 214, 365–372. [Google Scholar] [CrossRef]
  9. Zhang, J.; Miao, Y.; Batchelor, W.; Lu, J.; Wang, H.; Kang, S. Improving high-latitude rice nitrogen management with the CERES-rice crop model. Agronomy 2018, 8, 263. [Google Scholar] [CrossRef] [Green Version]
  10. Yao, Y.; Miao, Y.; Huang, S.; Gao, L.; Ma, X.; Zhao, G.; Jiang, R.; Chen, X.; Zhang, F.; Yu, K.; et al. Active canopy sensor-based precision N management strategy for rice. Agron. Sustain. Dev. 2012, 32, 925–933. [Google Scholar] [CrossRef] [Green Version]
  11. Raun, W.R.; Solie, J.B.; Johnson, G.V.; Stone, M.L.; Mullen, R.W.; Freeman, K.W.; Thomason, W.E.; Lukina, E.V. Improving nitrogen use efficiency in cereal grain production with optical sensing and variable rate application. Agron. J. 2002, 94, 815–820. [Google Scholar] [CrossRef] [Green Version]
  12. Raun, W.R.; Solie, J.B.; Stone, M.L.; Martin, K.L.; Freeman, K.W.; Mullen, R.W.; Zhang, H.; Schepers, J.S.; Johnson, G.V. Optical sensor-based algorithm for crop nitrogen fertilization. Commun. Soil Sci. Plan. 2005, 36, 2759–2781. [Google Scholar] [CrossRef] [Green Version]
  13. Cao, Q.; Miao, Y.; Shen, J.; Yu, W.; Yuan, F.; Cheng, S.; Huang, S.; Wang, H.; Yang, W.; Liu, F. Improving in-season estimation of rice yield potential and responsiveness to topdressing nitrogen application with Crop Circle active crop canopy sensor. Precis. Agric. 2016, 17, 136–154. [Google Scholar] [CrossRef]
  14. Lu, J.; Miao, Y.; Shi, W.; Li, J.; Yuan, F. Evaluating different approaches to non-destructive nitrogen status diagnosis of rice using portable RapidSCAN active canopy sensor. Sci. Rep. 2017, 7, 14073. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, K.; Ge, X.; Shen, P.; Li, W.; Liu, X.; Cao, Q.; Zhu, Y.; Cao, Q.; Tian, Y. Predicting rice grain yield based on dynamic changes in vegetation indexes during early to mid-growth stages. Remote Sens. 2019, 11, 387. [Google Scholar] [CrossRef] [Green Version]
  16. Zhao, G.; Miao, Y.; Wang, H.; Su, M.; Fan, M.; Zhang, F.; Jiang, R.; Zhang, Z.; Liu, C.; Liu, P.; et al. A preliminary precision rice management system for increasing both grain yield and nitrogen use efficiency. Field Crop. Res. 2013, 154, 23–30. [Google Scholar] [CrossRef]
  17. Xing, B.; Dudas, M.J.; Zhang, Z.; Xu, Q. Pedogenetic characteristics of albic soils in the Three River Plain, Heilongjiang province. Acta. Pedolog. Sin. 1994, 31, 95–104, (In Chinese with English Abstract). [Google Scholar]
  18. Rouse, J.W.; Haas, J.R.H.; Schell, J.A.; Deering, D.W. Monitoring vegetation systems in the Great Plains with ERTS. In Third Earth Resources Technology Satellite-1 Symposium; NASA Special Publication: Washington, DC, USA, 1974; pp. 309–317. [Google Scholar]
  19. Barnes, E.M.; Clarke, T.R.; Richards, S.E.; Colaizzi, P.D.; Haberland, J.; Kostrzewski, M.; Waller, P.; Choi, C.; Riley, E.; Thompson, T. Coincident detection of crop water stress, nitrogen status and canopy density using ground-based multispectral data. In Proceedings of the International Conference on Precision Agriculture and Other Resource Management, Bloomington, IN, USA, 16–19 July 2000; pp. 16–19. [Google Scholar]
  20. Jordan, C.F. Derivation of leaf-area index from quality of light on the forest floor. Ecology 1969, 50, 663–666. [Google Scholar] [CrossRef]
  21. Qi, J.; Chehbouni, A.; Huete, A.; Kerr, Y.; Sorooshian, S. A modified soil adjusted vegetation index. Remote Sens. Environ. 1994, 48, 119–126. [Google Scholar] [CrossRef]
  22. Chen, J.M. Evaluation of vegetation indices and a modified simple ratio for boreal applications. Can. J. Remote Sens. 1996, 22, 229–242. [Google Scholar] [CrossRef]
  23. Reyniers, M.; Walvoort, D.J.J.; Baardemaaker, J.D. A linear model to predict with a multi-spectral radiometer the amount of nitrogen in winter wheat. Int. J. Remote Sens. 2006, 27, 4159–4179. [Google Scholar] [CrossRef]
  24. Goel, N.S.; Qin, W. Influences of canopy architecture on relationships between various vegetation indices and LAI and FPAR. Remote Sens. Rev. 1994, 10, 309–347. [Google Scholar] [CrossRef]
  25. Gong, P.; Pu, R.; Biging, G.S.; Larrieu, M.R. Estimation of forest leaf area index using vegetation indices derived from Hyperion hyperspectral data. IEEE Trans. Geosci. Remote 2003, 41, 1355–1362. [Google Scholar] [CrossRef] [Green Version]
  26. Cao, Q.; Miao, Y.; Wang, H.; Huang, S.; Cheng, S.; Khosla, R.; Jiang, R. Non-destructive estimation of rice plant nitrogen status with Crop Circle multispectral active canopy sensor. Field Crop. Res. 2013, 154, 133–144. [Google Scholar] [CrossRef]
  27. Jasper, J.; Reusch, S.; Link, A. Active sensing of the N status of wheat using optimized wave-length combination: Impact of seed rate, variety and growth stage. In Precision Agriculture’09; Henten, E.J.V., Goense, D., Lokhorst, C., Eds.; Wageningen Academic Publishers: Wageningen, The Netherlands, 2009; pp. 23–30. [Google Scholar]
  28. Johnson, G.V.; Raun, W.R. Nitrogen response index as a guide to fertilizer management. J. Plant Nutr. 2003, 26, 249–262. [Google Scholar] [CrossRef] [Green Version]
  29. Mullen, R.W.; Freeman, K.W.; Raun, W.R.; Johnson, G.V.; Stone, M.L.; Solie, J.B. Identifying an in-season response index and then potential to increase wheat yield with nitrogen. Agron. J. 2003, 95, 347–351. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, X.; Miao, Y.; Dong, R.; Chen, Z.; Guan, Y.; Yue, X.; Fang, Z.; Mulla, D.J. Developing active canopy sensor-based precision nitrogen management strategies for maize in Northeast China. Sustainability 2019, 11, 706. [Google Scholar] [CrossRef] [Green Version]
  31. Su, M.; Huang, S.; Zhao, G.; Zhao, Q.; Yao, Y.; Miao, Y. Current rice management practices of farmers in Heilongjiang land reclamation area and improvement strategies. North. Rice 2012, 42, 28–33, (In Chinese with English Abstract). [Google Scholar]
  32. Li, S.; Ding, X.; Kuang, Q.; Ata-Ul-Karim, S.T.; Cheng, T.; Liu, X.; Tian, Y.; Zhu, Y.; Cao, W.; Cao, Q. Potential of UAV-based active sensing for monitoring rice leaf nitrogen status. Front. Plant Sci. 2018, 9, 1834. [Google Scholar] [CrossRef] [Green Version]
  33. Aranguren, M.; Castellón, A.; Aizpurua, A. Crop sensor-based in-season nitrogen management of wheat with manure application. Remote Sens. 2019, 11, 1094. [Google Scholar] [CrossRef] [Green Version]
  34. Cao, Q.; Miao, Y.; Li, F.; Gao, X.; Liu, B.; Lu, D.; Chen, X. Developing a new Crop Circle active canopy sensor-based precision nitrogen management strategy for winter wheat in North China Plain. Precis. Agric. 2017, 18, 2–18. [Google Scholar] [CrossRef]
  35. Tubaña, B.S.; Arnall, D.B.; Walsh, O.; Chung, B.; Solie, J.B.; Girma, K.; Raun, W.R. Adjusting midseason nitrogen rate using a sensor-based optimization algorithm to increase use efficiency in corn. J. Plant Nutr. 2008, 31, 1393–1419. [Google Scholar] [CrossRef]
  36. Macnack, N.; Khim, B.C.; Mullock, J.; Raun, W. In-season prediction of nitrogen use efficiency and grain protein in winter wheat (Triticum aestivum L.). Commun. Soil Sci. Plan. 2014, 45, 2480–2494. [Google Scholar] [CrossRef] [Green Version]
  37. Purba, J.; Sharma, R.K.; Jat, M.L.; Thind, H.S.; Gupta, R.K.; Chaudhary, O.P.; Chandna, P.; Khurana, H.S.; Kumar, A.; Uppal, H.S.; et al. Site-specific fertilizer nitrogen management in irrigated transplanted rice (Oryza sativa) using an optical sensor. Precis. Agric. 2015, 16, 455–475. [Google Scholar]
  38. Li, F.; Miao, Y.; Zhang, F.; Cui, Z.; Li, R.; Chen, X.; Zhang, H.; Schroder, J.; Raun, W.R.; Jia, L. In-season optical sensing improves nitrogen-use efficiency for winter wheat. Soil Sci. Soc. Am. J. 2009, 73, 1566–1574. [Google Scholar] [CrossRef] [Green Version]
  39. Bonfil, D.J. Wheat phenomics in the field by RapidScan: NDVI vs. NDRE. Isr. J. Plant Sci. 2016, 1–14. [Google Scholar] [CrossRef]
  40. Lu, J.; Miao, Y.; Shen, J.; Cao, Q.; Huang, S.; Wang, H.; Wu, H.; Hu, S.; Hu, X. Improving estimation of rice yield potential using active canopy sensor Crop Circle ACS 430 in Northeast China. In Precision Agriculture’15; Stafford, J.V., Ed.; Wageningen Academic Publishers: Wageningen, The Netherlands, 2015; pp. 102–115. [Google Scholar] [CrossRef]
  41. Fan, M.; Shen, J.; Yuan, L.; Jiang, R.; Chen, X.; Davies, W.J.; Zhang, F. Improving crop productivity and resource use efficiency to ensure food security and environmental quality in China. J. Exp. Bot. 2011, 63, 13–24. [Google Scholar] [CrossRef]
  42. Ju, C.; Buresh, R.J.; Wang, Z.; Zhang, H.; Liu, L.; Yang, J.; Zhang, J. Root and shoot traits for rice varieties with higher grain yield and higher nitrogen use efficiency at lower nitrogen rates application. Field Crop. Res. 2015, 175, 47–55. [Google Scholar] [CrossRef]
  43. Ladha, J.K.; Pathak, H.; Krupnik, T.J.; Six, J.; Kessel, C.V. Efficiency of fertilizer nitrogen in cereal production: Retrospects and prospects. Adv. Agron. 2005, 87, 85–156. [Google Scholar]
  44. Li, S.; Wang, Z.; Hu, T.; Gao, Y.; Stewart, B.A. Nitrogen in dryland soils of China and its management. Adv. Agron. 2009, 101, 123–181. [Google Scholar]
  45. Wang, H.; Miao, Y.; Zhao, G.; Yao, Y.; Khosla, R. Evaluating different integrated precision rice management strategies in Northeast China. In Proceedings of the Second International Conference on Agro-Geoinformatics (Agro-Geoinformatics 2013), Fairfax, VA, USA, 12–16 August 2013; pp. 419–422. [Google Scholar] [CrossRef]
Figure 1. RapidSCAN sensor-based in-season N recommendation algorithm developed for determining topdressing N rates at stem elongation stage and heading stage of high-yielding rice in this study. VIN: vegetation index at plots with sufficient N fertilization plots; VI0: vegetation index at plots without additional N topdressing application; RI-VI: in-season N response index based on vegetation index; RIHarvest: N response index based on yield; DAT: the number of days from planting to the date of sensing; INSEY: in-season estimate of yield; YP0: the potential yield without additional topdressing N application; Nrate: recommended topdressing N application rate at stem elongation stage; AEtopdressing: agronomic efficiency of topdressing N; SE_YPN or HD_YPN: the potential yield with added topdressing N application at the stem elongation stage or heading stage, respectively; SE_Nrate or HD_Nrate: recommended topdressing N application rate at the stem elongation stage or heading stage, respectively.
Figure 1. RapidSCAN sensor-based in-season N recommendation algorithm developed for determining topdressing N rates at stem elongation stage and heading stage of high-yielding rice in this study. VIN: vegetation index at plots with sufficient N fertilization plots; VI0: vegetation index at plots without additional N topdressing application; RI-VI: in-season N response index based on vegetation index; RIHarvest: N response index based on yield; DAT: the number of days from planting to the date of sensing; INSEY: in-season estimate of yield; YP0: the potential yield without additional topdressing N application; Nrate: recommended topdressing N application rate at stem elongation stage; AEtopdressing: agronomic efficiency of topdressing N; SE_YPN or HD_YPN: the potential yield with added topdressing N application at the stem elongation stage or heading stage, respectively; SE_Nrate or HD_Nrate: recommended topdressing N application rate at the stem elongation stage or heading stage, respectively.
Remotesensing 12 01440 g001
Figure 2. The relationships between yield without additional topdressing N application (YP0) and in-season estimate of yield (INSEY) calculated with NDVI, NDRE, MSAVI, and NDVI*RVI across all varieties at stem elongation stage and heading stage.
Figure 2. The relationships between yield without additional topdressing N application (YP0) and in-season estimate of yield (INSEY) calculated with NDVI, NDRE, MSAVI, and NDVI*RVI across all varieties at stem elongation stage and heading stage.
Remotesensing 12 01440 g002
Figure 3. Rice yield without additional topdressing N application (YP0) or NDVI (a; c; e; g) and YP0 or NDRE (b; d; f; h) obtained by RapidSCAN at the stem elongation (ad) and heading (eh) stages as affected by different N rates for Longjing 31 (a; b; e; f) and Longjing 21 (c; d; g; h) in 2014–2016, respectively. Different color bars represent the value of NDVI or NDRE in different years. Different dots represent YP0 in different years. The red curves were the curves of rice yield potential without additional topdressing N application (YP0). Vertical bars represent the LSD value (p = 0.05) among different N levels.
Figure 3. Rice yield without additional topdressing N application (YP0) or NDVI (a; c; e; g) and YP0 or NDRE (b; d; f; h) obtained by RapidSCAN at the stem elongation (ad) and heading (eh) stages as affected by different N rates for Longjing 31 (a; b; e; f) and Longjing 21 (c; d; g; h) in 2014–2016, respectively. Different color bars represent the value of NDVI or NDRE in different years. Different dots represent YP0 in different years. The red curves were the curves of rice yield potential without additional topdressing N application (YP0). Vertical bars represent the LSD value (p = 0.05) among different N levels.
Remotesensing 12 01440 g003
Figure 4. The relationships between response index calculated with yield (RIHarvest) and response index calculated with NDVI, NDRE, RVI, and REOSAVI across all varieties at stem elongation stage and heading stage.
Figure 4. The relationships between response index calculated with yield (RIHarvest) and response index calculated with NDVI, NDRE, RVI, and REOSAVI across all varieties at stem elongation stage and heading stage.
Remotesensing 12 01440 g004
Figure 5. The difference in AEtopdressing between Longjing 31 and Longjing 21 (a) and its relationship with RIHarvest (b). The red ✳ indicates the average of AEtopdressing in Figure 5a. The red lines are the different regression models for Longjing 31 and Longjing 21 in Figure 5b; *** indicates significance at the level of p < 0.001.
Figure 5. The difference in AEtopdressing between Longjing 31 and Longjing 21 (a) and its relationship with RIHarvest (b). The red ✳ indicates the average of AEtopdressing in Figure 5a. The red lines are the different regression models for Longjing 31 and Longjing 21 in Figure 5b; *** indicates significance at the level of p < 0.001.
Remotesensing 12 01440 g005
Figure 6. The difference of recommended topdressing N application rate (ΔNrate) from economically optimal N rates (EONR) between RONM strategy and RapidSCAN sensor-based PNM strategies for the rice varieties of Longjing 31 (a) and Longjing 21 (b) under deficient N status before topdressing (N80) in 2014–2016. RONM: the regional optimum N management topdressing N rate (36 kg N ha−1 in this study). SE-NDVI and SE-VItop indicated recommended topdressing N application rate calculated by the models based on NDVI and the top performing VIs (VItop) for different varieties at stem elongation stage, respectively. SE&HD-NDVI&NDRE and SE&HD-VItop indicated recommended topdressing N application rates calculated by the models with NDVI (at stem elongation stage) and NDRE (at heading stage), and VItop at stem elongation and heading stage, respectively. Vertical bars represent the LSD value (p = 0.05). Different letters indicate significant difference at p < 0.05 level within the same year.
Figure 6. The difference of recommended topdressing N application rate (ΔNrate) from economically optimal N rates (EONR) between RONM strategy and RapidSCAN sensor-based PNM strategies for the rice varieties of Longjing 31 (a) and Longjing 21 (b) under deficient N status before topdressing (N80) in 2014–2016. RONM: the regional optimum N management topdressing N rate (36 kg N ha−1 in this study). SE-NDVI and SE-VItop indicated recommended topdressing N application rate calculated by the models based on NDVI and the top performing VIs (VItop) for different varieties at stem elongation stage, respectively. SE&HD-NDVI&NDRE and SE&HD-VItop indicated recommended topdressing N application rates calculated by the models with NDVI (at stem elongation stage) and NDRE (at heading stage), and VItop at stem elongation and heading stage, respectively. Vertical bars represent the LSD value (p = 0.05). Different letters indicate significant difference at p < 0.05 level within the same year.
Remotesensing 12 01440 g006
Figure 7. The difference of recommended topdressing N application rate (ΔNrate) from economically optimal N rates (EONR) between RONM strategy and RapidSCAN sensor-based PNM strategies for the rice varieties of Longjing 31 (a) and Longjing 21 (b) under optimal N status before topdressing (N120) in 2014–2016. RONM: the regional optimum N management topdressing N rate (36 kg N ha−1 in this study). SE-NDVI and SE-VItop indicated recommended topdressing N application rates calculated by the models using NDVI and the top performing VIs (VItop) for different varieties at stem elongation stage, respectively. SE&HD-NDVI&NDRE and SE&HD-VItop indicated recommended topdressing N application rates calculated by the models using NDVI (at stem elongation stage) and NDRE (at heading stage), and VItop at stem elongation and heading stage, respectively. Vertical bars represent the LSD value (p = 0.05). Different letters indicate significant difference at p < 0.05 level within the same year.
Figure 7. The difference of recommended topdressing N application rate (ΔNrate) from economically optimal N rates (EONR) between RONM strategy and RapidSCAN sensor-based PNM strategies for the rice varieties of Longjing 31 (a) and Longjing 21 (b) under optimal N status before topdressing (N120) in 2014–2016. RONM: the regional optimum N management topdressing N rate (36 kg N ha−1 in this study). SE-NDVI and SE-VItop indicated recommended topdressing N application rates calculated by the models using NDVI and the top performing VIs (VItop) for different varieties at stem elongation stage, respectively. SE&HD-NDVI&NDRE and SE&HD-VItop indicated recommended topdressing N application rates calculated by the models using NDVI (at stem elongation stage) and NDRE (at heading stage), and VItop at stem elongation and heading stage, respectively. Vertical bars represent the LSD value (p = 0.05). Different letters indicate significant difference at p < 0.05 level within the same year.
Remotesensing 12 01440 g007
Figure 8. The difference of recommended N topdressing application rate (ΔNrate) from economically optimal N rates (EONR) between RONM strategy and RapidSCAN sensor-based PNM strategies for the rice varieties of Longjing 31 (a) and Longjing 21 (b) under surplus N status before topdressing (N160) in 2014–2016. RONM: the regional optimum N management topdressing N rate (36 kg N ha−1 in this study). SE-NDVI and SE-VItop indicated recommended topdressing N application rates calculated by the models using NDVI and the top performing VIs (VItop) for different varieties at stem elongation stage, respectively. SE&HD-NDVI&NDRE and SE&HD-VItop indicated recommended topdressing N application rates calculated by the models using NDVI (at stem elongation stage) and NDRE (at heading stage), and VItop at stem elongation and heading stage, respectively. Vertical bars represent the LSD value (p = 0.05). Different letters indicate significant difference at p < 0.05 level within the same year.
Figure 8. The difference of recommended N topdressing application rate (ΔNrate) from economically optimal N rates (EONR) between RONM strategy and RapidSCAN sensor-based PNM strategies for the rice varieties of Longjing 31 (a) and Longjing 21 (b) under surplus N status before topdressing (N160) in 2014–2016. RONM: the regional optimum N management topdressing N rate (36 kg N ha−1 in this study). SE-NDVI and SE-VItop indicated recommended topdressing N application rates calculated by the models using NDVI and the top performing VIs (VItop) for different varieties at stem elongation stage, respectively. SE&HD-NDVI&NDRE and SE&HD-VItop indicated recommended topdressing N application rates calculated by the models using NDVI (at stem elongation stage) and NDRE (at heading stage), and VItop at stem elongation and heading stage, respectively. Vertical bars represent the LSD value (p = 0.05). Different letters indicate significant difference at p < 0.05 level within the same year.
Remotesensing 12 01440 g008
Table 1. The details about the crop growth and crop sensing within the experiments performed in this study.
Table 1. The details about the crop growth and crop sensing within the experiments performed in this study.
ExperimentYearVarietyTransplanting DateHarvest DateStem Elongation StageHeading Stage
Sensing DateDATSensing DateDAT
Exp. 12014Longjing 3119 May29 September3 July4526 July68
Exp. 2201520 May4 October 6 July4730 July71
Exp. 3201619 May25 September5 July4725 July67
Exp. 42014Longjing 2119 May29 September7 July4926 July68
Exp. 5201520 May4 October7 July482 August.74
Exp. 6201619 May25 September5 July4725 July67
Note: DAT: the number of days from transplanting to sensing.
Table 2. NDVI, NDRE, and the best performing vegetation indices for calibration and validation from RapidSCAN used in this study.
Table 2. NDVI, NDRE, and the best performing vegetation indices for calibration and validation from RapidSCAN used in this study.
IndexFormulaReference
Normalized difference vegetation index (NDVI) ( NIR     R ) / ( NIR + R ) Rouse et al. [18]
Normalized difference red edge (NDRE) ( NIR     RE ) / ( NIR + RE ) Barnes et al. [19]
Ratio vegetation index (RVI) NIR / R Jordan [20]
Modified soil-adjusted vegetation index (MSAVI) [ ( 2   ×   NIR + 1 )     ( 2   ×   NIR + 1 ) 2     8 ( NIR     R ) ] / 2 Qi et al. [21]
Modified simple ratio (MSR) ( NIR / R     1 ) / NIR / R + 1 Chen [22]
Optimal vegetation index (VIopt) 1.45   ×   ( NIR 2   + 1 ) / ( R + 0.45 ) Reyniers et al. [23]
Nonlinear index (NLI) ( NIR 2     R ) / ( NIR 2 + R ) Goel and Qin [24]
NDVI*RVI ( NIR 2     R ) / ( NIR + R 2 ) Gong et al. [25]
Red edge wide dynamic range vegetation index (REWDRVI) ( 0.12 INIR     RE ) / ( 0.12 RNIR + RE ) Cao et al. [26]
Red edge optimal soil adjusted vegetation index (REOSAVI) ( 1 + 0.16 ) ( NIR     RE ) / ( NIR + RE + 0.16 ) Cao et al. [26]
Modified red edge soil adjusted vegetation index (MRESAVI) [ ( 2   ×   NIR + 1 )     ( 2   ×   NIR + 1 ) 2     8 ( NIR     RE ) ] / 2 Cao et al. [26]
Optimized red edge vegetation index (REVIopt) 100   ×   ( lnNIR     lnRE ) Jasper et al. [27]
Normalized near infrared index (NNIRI) NIR / ( NIR + RE + R ) Lu et al. [14]
Table 3. Significance of mean squares in the analysis of yield without additional topdressing N application (YP0), RapidSCAN NDVI, and NDRE under five N rates (0, 40, 80, 120, 180 kg ha−1) combined at two stages (stem elongation and heading stage), across three years (2014-2016) for two varieties (Longjing 31 and Longjing 21).
Table 3. Significance of mean squares in the analysis of yield without additional topdressing N application (YP0), RapidSCAN NDVI, and NDRE under five N rates (0, 40, 80, 120, 180 kg ha−1) combined at two stages (stem elongation and heading stage), across three years (2014-2016) for two varieties (Longjing 31 and Longjing 21).
Source of VariationDegree of FreedomSignificance of Mean Square
NDVINDREYP0
N level4*********
Variety1**ns***
Stage1*********
Year2ns******
Note: *, **, and *** indicate significance at 0.05, 0.01, and 0.001 probability levels, respectively. ns = non-significant.
Table 4. Calibration and validation results for predicting yield without additional topdressing N application (YP0) using the in-season estimate of yield (INSEY) calculated with the RapidSCAN’s default indices (NDVI and NDRE) and the best performing vegetation indices for two varieties at the stem elongation (SE) and heading (HD) stages in 2014–2016.
Table 4. Calibration and validation results for predicting yield without additional topdressing N application (YP0) using the in-season estimate of yield (INSEY) calculated with the RapidSCAN’s default indices (NDVI and NDRE) and the best performing vegetation indices for two varieties at the stem elongation (SE) and heading (HD) stages in 2014–2016.
VarietyStageCalibration ResultsValidation Results
IndexModelR2R2RMSEREr
Longjing 31SEINSEY_NDVIE0.760.700.497.7
INSEY_NDREP0.660.620.558.5
INSEY_NLIE0.780.700.497.6
HDINSEY_NDVIQ0.590.540.749.9
INSEY_NDREQ0.890.760.547.3
INSEY_REOSAVIQ0.890.760.557.3
Longjing 21SEINSEY_NDVIE0.720.660.396.3
INSEY_NDREP0.620.340.548.7
INSEY_MSRE0.740.630.416.5
HDINSEY_NDVIQ0.280.640.9312.1
INSEY_NDREQ0.770.750.668.5
INSEY_ MRESAVIQ0.780.720.699.0
Across varietiesSEINSEY_NDVIE0.730.660.477.4
INSEY_NDREP0.600.490.568.8
INSEY_MSAVIE0.740.650.477.4
HDINSEY_NDVIQ0.410.590.8110.7
INSEY_NDREQ0.730.700.658.6
INSEY_NDVI*RVIP0.750.700.668.7
Note: Q, E, and P: stand for quadratic, exponential, and power models. RMSE: root mean square error. REr: relative error (%).
Table 5. Calibration and validation results for the response index calculated with yield (RIHarvest) predicted by response index calculated with NDVI, NDRE, and the best performing vegetation indices for two varieties at stem elongation stage (SE) and heading stage (HD) in 2014–2016.
Table 5. Calibration and validation results for the response index calculated with yield (RIHarvest) predicted by response index calculated with NDVI, NDRE, and the best performing vegetation indices for two varieties at stem elongation stage (SE) and heading stage (HD) in 2014–2016.
VarietyStageCalibration ResultsValidation Results
IndexModelR2R2RMSEREr
Longjing 31SERI_NDVIQ0.780.770.086.2
RI_NDREQ0.710.720.096.8
RI_VIoptQ0.790.780.086.1
HDRI_NDVIQ0.920.650.1210.2
RI_NDREQ0.920.750.108.4
RI_NNIRIQ0.950.730.108.7
Longjing 21SERI_NDVIQ0.680.600.117.3
RI_NDREQ0.690.630.117.4
RI_REWDRVIQ0.710.670.107.0
HDRI_NDVIE0.670.600.1613.4
RI_NDREQ0.790.820.108.7
RI_REVIoptQ0.790.820.108.7
Across varietiesSERI_NDVIQ0.670.610.118.0
RI_NDREQ0.640.620.118.0
RI_RVIQ0.680.630.117.8
HDRI_NDVIQ0.790.620.1411.5
RI_NDREQ0.850.780.108.3
RI_REOSAVIQ0.850.780.108.3
Note: Q, E, and P stand for the quadratic, exponential, and power models. RMSE: root mean square error. REr: relative error (%).
Table 6. N application rate (Nrate), grain yield, agronomy efficiency of N (AEN), partial factor productivity of N (PFPN), and economic returns of different rice N management strategies in the evaluation experiments conducted in 2015 and 2016.
Table 6. N application rate (Nrate), grain yield, agronomy efficiency of N (AEN), partial factor productivity of N (PFPN), and economic returns of different rice N management strategies in the evaluation experiments conducted in 2015 and 2016.
Treatment20152016Across Two Years
Longjing 31Longjing 21Longjing 31Longjing 21Longjing 31Longjing 21Average
Nrate (kg ha−1)CK0000000
FNM160160160160160160160
RONM120120120120120120120
PNM118119122127120123121
Yield (t ha−1)CK5.26 b5.46 b5.25 b5.17 c5.25 c5.31 c5.28 c
FNM8.19 a8.60 a8.40 a8.90 ab8.30 a8.75 ab8.52 a
RONM7.78 a8.36 a8.05 a8.68 b7.91 b8.52 b8.22 b
PNM8.14 a8.59 a8.43 a9.06 a8.29 a8.83 a8.56 a
AEN (kg kg−1)CK-------
FNM18.3 b19.6 c19.7 c23.3 b19.0 c21.5 c20.3 c
RONM21.0 ab24.2 b23.3 b29.3 a22.2 b26.8 b24.5 b
PNM24.4 a26.4 a26.1 a30.8 a25.3 a28.6 a26.9 a
PFPN (kg kg−1)CK-------
FNM51.2 c53.7 c52.5 c55.6 b51.9 c54.7 b53.3 c
RONM64.8 b69.7 b67.1 b72.3 a65.9 b71.0 a68.5 b
PNM68.9 a72.4 a69.1 a71.6 a69.0 a72.0 a70.5 a
Economic return ($ ha−1)CK-------
FNM1205 a1297 a1300 a1556 a1253 a1426 a1339 a
RONM1045 a1214 a1166 a1481 a1106 b1348 b1227 b
PNM1205 a1315 a1334 a1647 a1269 a1481 a1375 a
Note: CK: the control treatment with no N application; RONM: regional optimum N management; FNM: farmer N management; PNM: RapidSCAN-based precision N management strategy. Within a column for parameter, values followed by different letters are significantly different (p < 0.05).

Share and Cite

MDPI and ACS Style

Lu, J.; Miao, Y.; Shi, W.; Li, J.; Hu, X.; Chen, Z.; Wang, X.; Kusnierek, K. Developing a Proximal Active Canopy Sensor-based Precision Nitrogen Management Strategy for High-Yielding Rice. Remote Sens. 2020, 12, 1440. https://doi.org/10.3390/rs12091440

AMA Style

Lu J, Miao Y, Shi W, Li J, Hu X, Chen Z, Wang X, Kusnierek K. Developing a Proximal Active Canopy Sensor-based Precision Nitrogen Management Strategy for High-Yielding Rice. Remote Sensing. 2020; 12(9):1440. https://doi.org/10.3390/rs12091440

Chicago/Turabian Style

Lu, Junjun, Yuxin Miao, Wei Shi, Jingxin Li, Xiaoyi Hu, Zhichao Chen, Xinbing Wang, and Krzysztof Kusnierek. 2020. "Developing a Proximal Active Canopy Sensor-based Precision Nitrogen Management Strategy for High-Yielding Rice" Remote Sensing 12, no. 9: 1440. https://doi.org/10.3390/rs12091440

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