Combining Remote Sensing and Meteorological Data for Improved Rice Plant Potassium Content Estimation

: Potassium (K) plays a significant role in the formation of crop quality and yield. Accurate estimation of plant potassium content using remote sensing (RS) techniques is therefore of great interest to better manage crop K nutrition. To improve RS of crop K, meteorological information might prove useful, as it is well established that weather conditions affect crop K uptake. We aimed to determine whether including meteorological data into RS-based models can improve K estimation accuracy in rice ( Oryza sativa L.). We conducted field experiments throughout three growing seasons (2017–2019). During each year, different treatments (i.e., nitrogen, potassium levels and plant varieties) were applied and spectra were taken at different growth stages throughout the growing season. Firstly, we conducted a correlation analysis between rice plant potassium content and transformed spectra (reflectance spectra (R), first derivative spectra (FD) and reciprocal loga-rithm-transformed spectra (log [1/R])) to select correlation bands. Then, we performed the genetic algorithms partial least-squares and linear mixed effects model to select important bands (IBs) and important meteorological factors (IFs) from correlation bands and meteorological data (daily average temperature, humidity, etc.), respectively. Finally, we used the spectral index and machine learning methods (partial least-squares regression (PLSR) and random forest (RF)) to construct rice plant potassium content estimation models based on transformed spectra, transformed spectra + IFs and IBs, and IBs + IFs, respectively. Results showed that normalized difference spectral index (NDSI (R 1210 , R 1105 )) had a moderate estimation accuracy for rice plant potassium content (R 2 = 0.51; RMSE = 0.49%) and PLSR (FD-IBs) (R 2 = 0.69; RMSE = 0.37%) and RF (FD-IBs) (R 2 = 0.71; RMSE = 0.40%) models based on FD could improve the prediction accuracy. Among the meteorological factors, daily average temperature contributed the most to estimating rice plant potassium content, followed by daily average humidity. The estimation accuracy of the optimal rice plant potassium content models was improved by adding meteorological factors into the three RS models, with model R 2 increasing to 0.65, 0.74, and 0.76, and RMSEs decreasing to 0.42%, 0.35%, and 0.37%, respectively, suggesting that including meteorological data can improve our ability to remotely sense plant potassium content in rice.


Introduction
Potassium (K) is an essential nutrient for crop growth. Key plant functions involving K include regulation of stomatal function, enzyme activation, and the formation and transport of photosynthate [1]. Thus, K content can affect the yield, quality, and water use efficiency of crops [2][3][4][5]. Meanwhile, the price and utilization rate of K fertilizer is high, and overfertilization may pollute water resources [6]. Therefore, accurate estimation of K nutritional status of crops and improving the K fertilizer utilization rate has great potential for increasing agricultural productivity and environmental sustainability [7,8].
Indeed, changes in meteorological factors can affect K nutrition in crops. For example, Hao et al. [39] found that rising temperatures increased K accumulation of soybean (Glycine max L.). Similarly, Wang et al. [40] showed that the increase in temperature improved the plant K nutrition accumulation of rice (Oryza sativa L.). Yin et al. [41] showed that higher daily average temperature, total precipitation and potential evaporation increased the demand for potash fertilizer in rice. Nieves-Cordones et al. [42] summarized that drought stress would enhance K cation uptake in crop leaves to improve osmotic adjustment and retain water within cells. In addition, precipitation controls over soil moisture may affect passive transport of soil K ions as they transfer from root to leaf to regulate the osmotic potential of guard cells. This process helps to regulate stomatal opening and ultimately affects transpiration rates to retain water [43][44][45].
To consider the important effect of meteorological factors on plant growth, meteorological variables have been incorporated into RS-based models to improve the estimation accuracy of crop yield and biochemical parameters [46][47][48][49]. For example, Ramoelo et al. [49] found that using non-linear PLSR to combine RS and meteorological data could effectively improve the estimation accuracy of nitrogen (N) and phosphorus (P) content in graminoid leaves. Similarly, Gao et al. [47] used machine learning to combine RS and meteorological data, which improved the estimation accuracy of foliar P content in graminoids. Furthermore, using machine learning to combine satellite RS and meteorological data helped to significantly improve the estimation accuracy of wheat (Triticum aestivum L.) and maize (Zea mays L.) yield [33,46,48]. However, to date little work has been focused on using remote sensing for monitoring K in rice. As part of this work, we address the following specific questions: Question 1: which spectral index (existing and newly proposed spectral index) provides the most reliable plant potassium content estimates in rice; Question 2: which type of remote sensing approach (e.g., machine learning, spectral index) results in the most reliable crop K estimates; and Question 3: Can the combination of remote sensing and meteorological data improve K estimates?
The overarching objective of this study was to explore whether the inclusion of meteorological data could improve the estimation accuracy of RS-based K nutrition models. To reach this overarching goal, we addressed the following three objectives: (1) to compare the ability to detect variations in rice plant potassium content (PKC) between simple spectral indices and machine learning methods; (2) to build rice PKC models that combine RS and meteorological data; and (3) to test the stability and temporal transferability of RS models for predicting rice PKC. For objective (2), we hypothesized that the inclusion of meteorological data will significantly improve the estimation accuracy of K models using remotely sensed data.
To test this, a field experiment with different K, nitrogen and plant variety treatments was designed and data were collected throughout the growing season during three consecutive years (2017-2019). Data collected included canopy spectra, meteorological information, and rice plant potassium content.

Experimental Design
We conducted this study at the experimental station of the National Engineering and Technology Center for Information Agriculture (NETCIA), Rugao, China (120°05′E, 32°06′N), where the annual average temperature, average humidity, average wind speed, number of precipitation days, precipitation, and annual average frost-free period are 14.6 °C, 79.37%, 1.91 m/s, 121.3 days, 1055.5 mm, and 215.6 days, respectively. The predominant soil type was a loam (before the beginning of experiment in the year 2017), with an organic carbon concentration of 12.95 g/kg, a total nitrogen content of 1.20 g/kg, a total potassium content of 12.26 g/kg, an available nitrogen content of 75.73 mg/kg, and an available potassium content of 65.40 mg/kg (general potassium deficiency soil: 50-70 mg/kg).

Acquisition of Canopy Reflectance and Agronomic Parameters
We collected canopy hyperspectral reflectance of rice using an ASD FieldSpec Pro spectrometer (Analytical Spectral Devices, Boulder, CO, USA) during key growth stages of rice development (tillering, productive tiller critical stage (N-n), jointing, booting, heading, and filling stage). We conducted measurements between 10:00 and 14:00 under clear skies to minimize confounding effects caused by variability in illumination geometry. Before data acquisition for each plot, we collected a white reference panel measurement using a Spectralon panel (Spectralon ® , Labsphere, Inc., North Sutton, NH, USA) (approximate reflectance values of 100%). We used a one-meter-long plastic pipe to ensure the fiber optic of the spectrometer was always one meter above the top of the canopy during each sampling campaign. The field of view of bare fiber was 25°. At each plot, three reflectance readings were taken and averaged. The spectrometer collects 2151 bands from 350 to 2500 nm at a sampling interval 1.4 nm and 3 nm spectral resolution in the 350-1050 nm region, and 2 and 10 nm in the 1050-2500 nm region. We removed very short wavelength bands (350-399 nm) and the spectral reflectance values from 1351-1440 nm, 1801-1960 nm, and 2401-2500 nm due to the influence of atmospheric moisture absorption characteristics.
After collecting the spectral data, we selected three whole representative rice plants from each plot to be uprooted, placed in a bucket of water, and brought back to the laboratory for PKC analysis. The rice plants were cut at the roots and put in paper bags, ovendried at 105 °C for 30 min, and kept at 80 °C until the weight was constant. After drying, the rice plants were ground, and chemical analysis was conducted. The PKC was measured using an atomic absorption spectrometer (ICE 3000; ThermoFisher, Sunnyvale, CA, USA).

Meteorological Data Collection and Analysis
We collected meteorological data using the automated weather station Dynameta-1K (Dynamax Inc., Houston, TX, USA) from May to October (Figure 1), which included daily average temperature (Temdaily) (°C), average relative humidity (%), accumulative precipitation (mm) and average wind speed (m/s). The temperature data was divided into Temdaily and accumulated growing degree day (AGDD). AGDD represents the sum of growing degree days (GDD) from transplanting to test day, GDD = Temdaily −12.5 °C [50].
To determine the effect of meteorological factors on rice PKC, we used a linear mixed effects model and set rice PKC as a dependent variable, meteorological factors as fixed effects, and day as a random effect with a random intercept. To avoid collinearity and overfitting, meteorological factors were removed using a stepwise multiple linear regression (SMLR) until the variance inflation factor (VIF) of all meteorological factors was less than 3 [51]. The linear mixed effects analysis was performed with the "lme4" package [52] in the R language (R Core Team, 2018). We considered meteorological factors at p < 0.001 to have a statistically significant effect on rice PKC and selected them as important meteorological factors (Ifs). The standardized slope (b*) was used to determine the effect size of each meteorological factor.

Modeling and Validation of Potassium Nutrition
The workflow associated with our approach to identify a reliable remote sensing model for predicting rice PKC is shown in Figure 2. First, using normalized difference spectral index (NDSI) and correlation analysis, we selected two-band indices and correlation bands (CBs) based on three transformed spectra (reflectance spectra (R), first derivative spectra (FD) and reciprocal logarithm-transformed (log [1/R])). Next, we used the genetic algorithms partial least-squares (GA-PLS) (more details see below) and linear mixed effects model to select important bands (IBs) and important meteorological factors (IFs), respectively. Finally, we used spectral indices and machine learning (partial least-squares regression (PLSR) and random forest (RF)) to build K nutrition RS models and combined RS models and IFs to build new K nutrition models. Furthermore, we used a validation dataset to evaluate the performance of K nutrition models and divided all data into different datasets to test the stability and temporal transferability of K models and selected optimal K models according to simplicity and estimation accuracy.
To select important bands that contribute the most to the variability of rice PKC, we performed variable selection with the commonly used GA-PLS [53]. GA-PLS was implemented through the "plsVarSel" package [54] in the R language. To compare the ability of simple spectral indices and machine learning methods to detect variability in rice PKC, we did the following. First, we built simple linear regression models between rice PKC and spectral indices proposed by previous studies and this paper. The optimal spectral index was selected by root mean square error (RMSE), bias and relative root mean square error (RE). Meanwhile, a SMLR model was developed between rice PKC and the optimal spectral index and IFs. To avoid overfitting, two statistical parameters of the SMLR model (VIF and tolerance (TOL)) were set less than 3 and greater than 0.33, respectively [51]. We also employed machine learning (i.e., PLSR and RF), which selected IBs and IFs as input variables. PLSR was implemented through the "pls" package [55] in the R language. RF included two important tuning parameters: (1) the number of regression trees (ntree) and (2) the number of predictors tested at each node (mtry). The ntree was set to 1000 according to previous studies [37,56] and mtry was determined by the "randomForest" package [57] in the R language. Figure 2. Workflow of the methodology used here to identify optimal remote sensing models for predicting rice PKC. Note: R = reflectance, FD = first derivative spectra, LOG = reciprocal logarithmtransformed spectra, NDSI = normalized difference spectral index, PLSR = partial least-squares regression, RF = random forest, GDD = growing degree day, and GA-PLS = genetic algorithms partial least-squares.
To build the rice PKC estimation models and test their performance, a total of 864 samples from three growing seasons (2017-2019) were randomly split into a calibration (2/3 of data, 576 samples) and validation (1/3 of data, 288 samples) dataset ( Table 2). The coefficient of determination (R 2 ), RMSE, bias, RE, and Akaike information criterion (AIC) [58] of all model types were calculated in R (R Core Team, 2018) and used to evaluate the prediction accuracy of these models. Note: STD and CV are standard deviation and coefficient of variation (STD/Mean), respectively. Table 2 showed the descriptive statistics of PKC for different datasets. The 2018 dataset exhibited a higher mean PKC, and the 2019 dataset exhibited a higher degree of dispersion. The PKC decreased with the growth of rice plant; the no K fertilization (K0) had the lowest PKC and K2 had the highest PKC (Figure 3a,b). We selected the optimal treatment (V1N2K2) to show the mean reflectance with six growth stages ( Figure 3c); before the booting stage, the spectral reflectance decreased in the visible region and increased in the near-infrared region, after the booting stage, the spectral reflectance decreased in the near-infrared region due to the influence of ears. Then, we showed the variation pattern of mean NDSI for R, FD, and LOG with different AGDD (calculation of NDSI see Section 3.4), in summary, NDSI values of R and FD showed a similar variation, they decreased with the growth of rice plant, no K fertilization treatment (K0) had lowest NDSI value, and K2 had the high NDSI value (Figure 3d,e). For the NDSI value of LOG, it showed an opposite variation to R (Figure 3f).

Correlation of Rice PKC with Spectral and Meteorological Data
The reflectance spectra (R) were positively correlated with rice PKC at 400-723, 1441-1800, and 1961-2400 nm (maximum R 2 value of 0.38 at 1969 nm) and negatively correlated at 750-955 and 1008-1130 nm (Figure 4a). The log (1/R) spectra had an opposite trend of the R spectra. The FD spectra had a negative correlation with rice PKC at 759-810, 985-1066, and 1268-1341 nm (maximum R 2 value of 0.38 at 1045 nm) and a positive correlation at 920-973 and 1077-1264 nm. Correlation bands (CBs) were selected by significant correlation analysis and used for selecting important bands (IBs) (see Section 3.3 below).
After excluding highly collinear meteorological factors based on TOL and VIF values (Table 3), we fit a linear mixed effects model with daily average temperature, humidity, and wind speed as fixed effects and day as a random effect with a random intercept on rice PKC. The results from the linear mixed effects model revealed that daily average temperature had the strongest (b* = 0.233, p < 0.001) significant effect on the rice PKC, followed by average humidity (b* = −0.065, p < 0.001) (Figure 4b). Thus, daily average temperature and humidity were selected as important meteorological factors (Ifs).

Figure 4.
The correlation between rice plant potassium content (PKC) (a), reflectance spectra (R), first derivative spectra (FD), and reciprocal logarithm-transformed spectra (log (1/R)). Dotted horizontal lines illustrate significance level (p value < 0.001). Standardized slope (b*) estimates for meteorological factors of the linear mixed effects model with PKC as the dependent variable and the following fixed effects: daily average temperature, average wind, and humidity. The symbol of *** represents statistical significance levels (p < 0.001) (b). (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article). Note: VIF = variance inflation factor, TOL = tolerance, and R 2 adj = adjusted coefficient of determination.

Important Bands Selection
The GA-PLS was used to select IBs that were sensitive to variability in rice PKC from the CBs. The results showed that the number of IBs were 40, 57, and 30 for the R, FD, and Log (1/R), respectively (Table 4). Some IBs for transformed spectra were located in the visible region that was influenced by pigment content. Most IBs for transformed spectra were located in near-infrared and shortwave infrared regions that were influenced by water, protein, starch, lignin, cellulose, and sugar (Figure 5d).
For our machine learning models (PLSR and RF), we used IBs identified with GA-PLS to predict rice PKC. The PLSR model based on FD-IBs showed to be the strongest predictor of rice PKC (R 2 = 0.69, RMSE = 0.38%) (Figure 8b), followed by the PLSR models based on R-IBs and LOG-IBs (R 2 = 0.67, RMSE = 0.39%) (Figure 8a,c). RF modelling showed similar results to PLSR modelling and the RF model based on FD spectra showed to be the strongest predictor of rice PKC, followed by RF models based on R-IBs and LOG-IBs (R 2 = 0.71, 0.62 and 0.61, RMSE = 0.37%, 0.42% and 0.43%, respectively) (Figure 9a-c). Compared to the ability of spectral indices to detect the variations in rice PKC, the machine learning methods showed better estimation accuracy (R 2 and RMSE increase and decrease of more than 0.1 and 0.06%, respectively).

Estimation of Rice PKC with a Combination of Remote Sensing and Meteorological Data
Here, we added meteorological data identified in our correlation analyses (see Section 3.2) into RS models. The optimal spectral indices among different transformed spectra and IFs (daily average temperature and humidity) were selected to develop SMLR models ( Table 6). We found that SMLR models that combined spectral indices and IFs improved the estimation accuracy of rice PKC, and the model based on the FD index performed better than R and log (1/R) indices for estimating rice PKC. The SMLR model of NDSI (FD1505, FD805) + AGDD showed the best performance for estimating rice PKC, followed by NDSI (R1210, R1105) + AGDD and NDSI (LOG1210, LOG1180) + AGDD. Compared to the models using spectral index alone, the R 2 of these SMLR models increased by 17.24%, 24.45%, and 45.45% to 0.68, 0.65, and 0.64, respectively. Note: R = reflectance, FD = first derivative spectra, LOG = reciprocal logarithm-transformed spectra, AGDD = accumulated growing degree day, VIF = variance inflation factor, TOL = tolerance, and R 2 adj = adjusted coefficient of determination.
To test the stability and temporal transferability of rice PKC models further, all the data was divided into eleven different datasets according to N level, K level, variety, and experiment year (Table 8). We selected NDSI (R1210, R1105) and NDSI (FD1505, FD805) using fewer bands and PLSR (FD-IBs) and RF (FD-IBs) using many bands to test their performance. Our results showed that the rice PKC models of spectra combined with meteorological factors could better detect the variability of rice PKC than using spectra alone (with higher R 2 and lower RMSE). Considering the ability of RS models for detecting the variability in rice PKC, they performed better when used in japonica rice than indica rice; Low N level (N1) could be better detected by them than high N level (N2); K stress levels could be better detected by them than normal K level: K1 > K3 > K0 > K2. Meanwhile, the K RS models have good performances under different years and performed better when used in Exp. 3 (2019) than Exp. 2 (2018) and Exp. 1 (2017). Finally, we selected the optimal PKC models: one was NDSI (FD1505, FD805) + AGDD; another was PLSR (FD-IBs + IFs) and RF (FD-IBs + IFs). 0.039 −432 Note: R = reflectance, FD = first derivative spectra, LOG = reciprocal logarithm-transformed spectra, IBs = important bands, IFs = important meteorological factors, NDSI = normalized difference spectral index, LR = linear regression, SMLR = stepwise multiple linear regression, PLSR = partial least-squares regression, RF = random forest, AGDD = accumulated growing degree day, IFs-Temdaily represents daily average temperature and humidity, IFs-AGDD represents AGDD and daily average humidity, and the italic font represents the best performances of various methods. Table 8. Correlations between rice plant potassium content (PKC) and different remote sensing methods under the different datasets. Note: R = reflectance, FD = first derivative spectra, IBs = important bands, IFs = important meteorological factors, NDSI = normalized difference spectral index, PLSR = partial least-squares regression, RF = random forest, AGDD = accumulated growing degree day, and IFs-AGDD represents AGDD and daily average humidity.

Canopy Sensitive Bands and Transformed Spectra for Plant K Estimation
Our results showed that the spectral bands of rice canopy that were sensitive to PKC were mainly concentrated in the near-infrared and shortwave infrared regions (1000-1340 nm, 1450-1800 nm, and 1961-2400 nm) (Figure 6a). This result was consistent with previous studies that the shortwave infrared spectra could be used to identify crop K nutrition levels in other crops [9,15,62,63]. The K spectral index proposed here included bands at 1210 nm and 1505 nm (Table 5), which have previously been used to construct a spectral index to monitor vegetation moisture status [64]. Meanwhile, the K spectral indices NDSI (R1645, R1715) and NDSI (R870, R1450) proposed by Pimstein et al. [15] and NDSI (R1705, R1320) proposed by Lu et al. [9] performed well in this paper, included shortwave infrared bands and these bands were also closely related to vegetation moisture status [65]. In addition, most of the important bands selected by GA-PLS were also located in the near-infrared and shortwave infrared regions, and only a few important bands were located in visible regions ( Figure 5). The near-infrared and shortwave infrared bands contributed the most to the K estimation ( Figure 10). For K element, it does not have strong absorption features in spectra, and their predictions are due in part to correlations with other traits such as pigments, water, starch, protein, cellulose, and sugar [66,67]. One physiological explanation of why important spectral regions for the remote sensing of K appear to be mainly located in the near-infrared and shortwave infrared spectral regions is that K exists in the form of K ions in the plant, which regulates stomatal conductance and osmotic potential, and enhances the synthesis of many enzymes, proteins, starches, and cellulose [1]. Finally, the variations in the plant K element will influence the spectral changes of corresponding absorption features.
The rice PKC models using a novel spectral index and machine learning based on the FD spectra showed better estimation accuracy than these rice PKC models based on the R spectra (Table 7). This was consistent with the conclusion by Lu et al. [9] who found that the FD spectral index was better than the R spectral index on estimating the rice leaf K content, and Gao et al. [47] found that the performance of machine learning based on the FD spectra on estimating leaf phosphorus content in forage was better than the R spectra. The K RS models using FD spectra had a better performance than using the original reflectance spectra in that FD spectra could reduce noise signals generated by atmospheric water absorption, soil background, and data redundancy, and enhance the absorption characteristics of plant leaf biochemical parameters [23].  (Table 4) for reflectance spectra (a), first derivative spectra (b), and reciprocal logarithm-transformed spectra (c). We used the importance function in the "Random-Forest" package to perform this analysis. The %IncMSE represents the percent increase in mean squared error.

Performance of Spectral Index and Machine Learning for Estimating Plant K
The results of this paper showed that the optimal estimation model of rice PKC could be divided into two types: one was the spectral index NDSI (FD1505, FD805) using two bands, and the R 2 and RMSE of the model were 0.58 and 0.47%, respectively. Another was the PLSR (FD-IBs) and RF (FD-IBs) of the machine learning method using multi-band spectral information, the R 2 of the models were 0.69 and 0.71, and the RMSEs were 0.37% and 0.40%, respectively. The spectral index-based approach requires only two bands, hence it seems to be feasible to design a relatively low-cost sensor that could be installed on a tractor or Unmanned Aerial Vehicle to inform K fertilizer decisions. Even though the estimation accuracy was higher for machine learning approaches, they rely on data from hyperspectral sensors that are expensive, fragile, and require expert knowledge to be processed. Hence, for K nutrition estimation, machine learning approaches might be most useful for scientists, but might be less practical for operational decision making.

Incorporating Remote Sensing and Meteorological Data for Plant K Estimation
Of the meteorological variables in our study, the daily average temperature contributed the most to estimating rice PKC (b* = 0.233), followed by the daily average humidity (b* = −0.065). This was consistent with previous studies [39,40,42] that found rising temperature would increase crop K uptake. Similarly, Gisleröd et al. [68] showed that reducing air humidity would increase plant K uptake. Meanwhile, Yin et al. [41] also found an increase in average temperature and evaporation would increase the demand for potash fertilizer. The increase or decrease in average temperature and humidity will affect stomata opening and closing, thereby influencing plant transpiration and water transfer. K is closely related to water [9], which mainly regulated the opening and closing of stomata and the penetration potential to affect plant moisture status [69]. Thus, the daily average temperature and humidity will affect crop K uptake. Our study suggested that combining RS and meteorological data could effectively improve the estimation accuracy of crop K nutrition. Specifically, combining spectral and IFs-AGDD in models outperformed models with spectral and IFs-Temdaily (Table 7). AGDD represents the effective accumulated temperature that crops require to complete their different growth stages [70], which suggests the accumulated temperature is more important than daily conditions for our results. Our finding that meteorological information could improve the RS of rice PKC in all model types was consistent with previous works which found that combining spectra and meteorological factors could improve the estimation accuracy of forage P content [47,49]. Similarly, other studies found that combining spectra and meteorological factors could improve the estimation accuracy of crop yield [33,46,48,71]. Meteorological factors mainly represent plant environmental information that plants are dynamically responding to and the spectra represent the current status of the plant. This might explain why combining spectra and meteorological factors can improve the accuracy of crop K nutrition estimation.

Evaluation of K Models and Future Directions
Compared to K estimation models using remote sensing data alone, the optimal K estimation models proposed in this study which incorporated meteorological data showed better stability under different datasets and better temporal transferability between different years ( Table 8). The reason why RS models performed better when used in japonica rice than in indica rice is that japonica rice is more sensitive to K, which leads to more abundant K accumulation [9,19]. Low N levels can be better detected using RS than high N levels due to ammonium ions in plant roots under high N concentrations, which have an inhibitory effect on K uptake and thus reduce plant K accumulation [72]. K stress levels can be better detected using RS than normal K levels, indicating that K stress could affect the stomatal opening and closing of plant leaves and osmotic potential regulation, thus affecting leaf water status [73]. Therefore, vegetation signals under K stress are easier to remotely sense. Though beyond the scope of this study, phosphorus levels may also influence the transfer of soil potassium ions. Therefore, it would be valuable if future work could focus on the interaction between phosphorus and potassium.
When we tested the performance of RS models in different years, our result showed that RS models performed better when used in Exp. 3 (2019) than Exp. 2 (2018) and Exp. 1 (2017); this is as our study carried out the same K treatment experiment of a rice-wheat rotation system in the same experiment location. With the progress of the experiment, the difference of soil K gradient is more obvious, which can be seen from the coefficient of variation of the dataset (Table 2). Hence, the K gradient difference of Exp. 3 is the largest and the estimation accuracy of the K nutrition model is the highest. These experiments were conducted at only one crop and ecological site and the band used in the novel K spectral indices had been used for monitoring vegetation water status. Additionally, for other crops at different field sites, variability in topography will likely influence water movement and hence the availability of K. Thus, future studies should consider the topography attributes as predictor variables and test their importance in different crops, field sites (including different soil types), and drought stress conditions. Meanwhile, the ease of access to meteorological data continues to increase thanks to the establishment of environmental monitoring networks in many agricultural regions [74], fusing NDSI from a multi-band remote sensing device with air temperature and humidity information from the environmental monitoring networks could help to more effectively manage PKC in rice.

Conclusions
This study systematically compared the performance of different transformed spectra, spectral indices, machine learning, and spectra combined with meteorological factors for estimating rice K nutrition. Results showed that the K sensitive bands for rice canopy were mainly located in the near-infrared and shortwave infrared regions of the electromagnetic spectrum. The K nutrition models of spectral index and machine learning based on the FD spectra showed good estimation accuracy. Meanwhile, among the meteorological factors, daily average temperature explained most of the variability of observed rice PKC, followed by daily average humidity. The estimation accuracy of rice PKC could be significantly improved by combining spectra and meteorological data in all model types. The rice PKC estimation models of NDSI (FD1505, FD805) + AGDD, PLSR (FD-IBs + IFs-AGDD) and RF (FD-IBs + IFs-AGDD) showed the best performance, with R 2 of 0.68, 0.74, and 0.76, and RMSEs of 0.41, 0.35, and 0.37, respectively. Furthermore, the optimal models showed better stability and temporal transferability under different years and cultivation factors (N and K levels and varieties). Our findings will provide novel evidence for the improved estimation of crop K nutrition based on RS combined with meteorological data.