INTRODUCTION
The use of remote sensing is increasingly widespread in the analysis and management of natural resources. Its use is related to the availability of time series satellites images in visible, infrared and ultraviolet spectra. In the same way, the number of satellites placed in orbit with thermal imaging, multispectral and hyperspectral sensors have contributed in the uses of spectral images (Triantafyllou et al., 2019; Segarra et al., 2020). Through monitoring with remote sensors, it has been possible to carry out biophysical a characterization of the surface vegetation, the state of the soil and humidity. Due to the optical properties of leaves, their absorption in the red spectral band of the electromagnetic spectrum is very high, while they reflect in the near infrared band. (Xie et al., 2019).
Numerous investigations in this field have shown the existence of empirical relationships between sensor information and biophysical properties such as biomass, leaf cover and surface moisture content. (Gutierrez et al., 2018; Hatfield et al., 2019; De Grave et al., 2020; Sishodia et al., 2020).
The vegetation indices (VI) are calculated through the combination of the red and near infrared bands of the electromagnetic spectrum. Thus, the empirical relationships calculated between biophysical variables of crops and VI, reflect general patterns or trends between crops optical properties and their architecture with respect to the spatial and angular distribution of the foliage phytoelements for a given sensor-illumination geometry. (Kharuf-Gutierrez et al., 2018).
In Cuba, the cultivation of sugarcane is carried out extensively and it has a high economical value with the production of raw and refined sugar. The by-products of the industrial process are used as raw material in food, chemical, and energy industries, among others. As part of the modernization of sugarcane cultivation, improvements have been made in varieties, irrigation technologies and harvesting and transportation. (Alvarez, 2014). Precision agriculture techniques have also been applied in fertilization, irrigation and sanitary control systems. (Hernández et al., 2003; Quevedo et al., 2006; Gutierrez et al., 2018).
In this framework, the present work aims to evaluate the effectiveness of monitoring sugarcane cultivation through spectral indices, obtained from satellite images and based on the correlation with the crop development in the field.
MATERIALS AND METHODS
Location of Area of Study
The study was carried out in the Basic Unit of Cooperative Production (UBPC) ¨Juan Verdecia¨, Camajuaní Municipality, Cuba. For crop monitoring, field number 3 was selected, located at 22°26'55'' North and 79°41'48'' West coordinates, with an area of 5,5 ha planted with sugarcane variety of CP 52-43 (Figure 1) and average yield of 35 t/ha. The previous harvest took place on 09-01-2021 and a combined KTP-2M and the ZIL-130 truck were used for it.
Methodology for Field Measurements
Crop growth control was performed to determine the effectiveness of satellite measurements. For this, the double diagonal or X distribution method was used (Flores, 2006), nine experimental points were taken, which GPS georeferenced for precision of 0,2 m. The controlled variables in the plants were leaf width, stem diameter and stem height. In addition, the measurement of soil moisture and thickness of the vegetation cover was carried out. For plants dimensions and ground cover, eight samples were taken at each experimental point and four repetitions were carried out. Tape measure, graduated rule and Vernier caliper with an appreciation of 1 mm and 0,05 mm, respectively, were used in measuring (Figures 2a and 2b).
Moisture based on dry soil (hbss)%, was determined by gravimetric method. Measurements were made in the ridge next to the cane stump (Figure 2c). Four measurements were made at each experimental point at 0 to 10 cm depth. The humidity of the samples was determined in the laboratory according to NC 3437:2003 standard.
Methodology for Satellite Monitoring
For satellite image monitoring of the area, a georeferenced study plot was configured using the services available in the Earth Observed System (EOS) on the site https://eos.com. The indices monitored during the vegetative period of the sugarcane plantation were: Normalized Difference Vegetation Index (NDVI), Soil Adjusted Vegetation Index (SAVI), Enhanced Vegetation Index (EVI) and the Normalized Differential Water Index (NDWI). The spatial distribution images for indices were obtained from LandViewer Repository, calculated from spectral images taken by terrestrial reconnaissance satellites LandSat and Sentinel-2. With data of field measurements and the spectral indices of crop, the descripted statistics of variables were determined, multiple correlation and regression analyses were carried out for a confidence interval of 95% using Centurions Plus software.
RESULTS AND DISCUSSION
Result of Measurements in the Field
When monitoring the evolution of the crop physiological indicators through the field measurements, reduced increases were seen until April, as shown in Figure 3. In this sense, a gradual increase in the width of the leaf, diameter and height of the stem took place. This increase become more evident in measurements made to the stem height. The three indicators showed a similar pattern for the nine months of cultivation, after which a gradual decrease happens.
The soil cover has a maximum value after the harvesting the field, as a result of the cleaning process by the harvester during the cutting (Aguilar-Pardo et al., 2016). The leaf fragments and shoot reach an average height of 5,6 cm, forming a regular residues layer distributed on the ground. In the second month, a reduction in coverage happens, as a result of the first mechanized cultivation, later a gradual reduction of it takes place gradually as part of the decomposition process, with increase in soil moisture. A further decrease in soil cover takes place in the month seven of cultivation, caused by the second mechanized cultivation operation which reduces the superficial layer of the soil to 0,2 cm.
The measurements made in soil moisture throughout the crop growth (Figure 4), show the transition from dry to rainy period, imposing two fundamental moisture conditions. As a result, the average value of humidity in the first four months was 18,9%, because the absence of rainfall, which reduced the availability of micronutrients and minerals for new shoots (Taiz y Zeiger, 2006). In the fifth month of the crop, the first rainfall in the area took place, soil moisture increases an average of 42.9% for the rest of the season and remained stable due to new rains contribution and water retention by plantation foliage.
Result of the Spectral Indices
Figure 5 shows the spectral images of the indices evaluated at the beginning of the study. The EVI index reaches values between 0,3 and 0,6, located fundamentally in the central zone of the field and values between 0,2 and 0,3 in periphery. In both cases, these values are related to the presence of low density vegetation, without damage from pests or other cause (Wang et al., 2010). With a similar distribution pattern to the EVI, the NDVI index shows values between 0,4 and 0,5 in the central area of the crop and values between 0,3 and 0,4 in the periphery, which characterizes a vegetation of limited development for both cases (Gilabert et al., 2002). The SAVI vegetation index shows a regular distribution in the range of 0,2 to 0,3 over the entire surface of the field. The NDWI index shows values in the range of -0,6 and 0,3, which corresponds to dry soil and a low moisture content in the biomass (Senay et al., 2014).
The spectral images of the crop final stage (Figure 6), show EVI index reaches values between 0,6 and 0,8 up to 74,4% of the crop, in the same way, NDVI index shows values between 0,7 and 0,8 which represents abundant vegetation, corresponding with the foliar development of plantation according to the field. SAVI index shows values between 0,3 and 0,6 which underestimates the real density of the existing biomass over surface. Regarding the water content of biomass and soil, NDWI index shows values between -1 and -0.6. Although the soil at this stage has moisture content greater than 40%, and the water content in the grass during the maturation stage is greater than 70% (Rosell Pardo y Ramírez Rubio, 2021).
Analysis of Spectral Indices Variability
As shown in NDVI index from LanSat and Sentinel satellites (Figure 7), between the frames under study, the values obtained in both cases represent an increases in vegetation with respect to initial values. Two stages are identified, divided by the fifth month of the crop, where NDVI goes from 0,33 to 0,65. This increase coincides with the beginning of the rainfall, where biomass increase is reflected by plants biophysical variables in according to data from field. However, by Sentinel monitoring, a more stable representation of the increases in vegetation is achieved, as well as more adequate values of initial and final states of the crop. The fluctuation of LandSat data in the fifth and tenth month does not correspond to changes in vegetation, being attributable to the dispersion of electromagnetic waves due to the effect of radiation from the atmosphere (Wang et al., 2010).
In moisture monitoring by NDWI index (Figure 8), in both satellites, initial values between -0,35 and -0,4, are obtained identifying a dry soil which coincide with values from field soil moisture. From the fifth month, NDWI decreases to -0.5 in both cases, and remains so until the tenth month where it shows a brief increase, to end again at low values of -0,47. With the foliar growth, the acquisition of soil reflectance is difficult, until the fifth month NDWI represent a dry soil due to the little foliar cover of soil. With the beginning of spring season, the increase in biomass happens and NDWI begins to reflect the moisture contained on it, which interferes with soil moisture readings. Despite the impossibility of monitoring soil moisture directly, through the NDWI is possible to identify the moisture values for uncovered soil as well as the stage where changes take place.
Correlation between Indices and Field Measurements
Table 1 shows the correlation between the spectral indices and growth variables through the Pearson correlation coefficient. In this sense, the strongest correlation with NDVI is obtained with respect to the stem diameter, reaching a value of 0,93 followed by leaf width, soil moisture and stem height, all the coefficient have positive values except coverage. Soil moisture and soil cover do not have a direct link with the spectral response of plants vegetative state and correlation with NDVI is established more clearly than the other variables. Similarly, in the correlation between NDVI and NDWI, a value close to -1 is obtained, which shows the low interdependence between indices. NDWI, only has a positive correlation with land cover of 0,67; negative values are obtained for the other parameters. Although soil moisture is the parameter most closely related to NDWI, the correlation obtained is -0,817, the development of plant foliage is one of the factors that prevent establishing the adequate prediction of moisture.
Index | Stem height | Cover | Stem diameter | Leaf width | Soil moisture | NDVI | NDWI | |
---|---|---|---|---|---|---|---|---|
NDVI | Coef.Pearson | 0.782 | -0.896 | 0.936 | 0.911 | 0.873 | -0.873 | |
P-Value | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | ||
NDWI | Coef.Pearson | -0.491 | 0.675 | -0.734 | -0.686 | -0.817 | -0.873 | |
P-Value | 0.020 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 | ||
Stem diameter | Coef.Pearson | 0.9273 | -0.9585 | 0.9900 | 0.8883 | 0.9363 | -0.733 | |
P-Value | 0.0000 | 0.0000 | 0.0000 | 0.0000 | 0.0000 | 0.0001 |
The linear regression adjusted model between stem diameter and NDVI index is shown in Figure 9. There, it is possible to make an adequate forecast of the stage of development of the crop from NDVI with a determination coefficient R2 of 87,6%. Likewise, Table 1 shows a high correlation of stem diameter with height and width of the leaf, whichmakes possible to predict the general state of development of the crop at different times.
CONCLUSIONS
In the evolution of the physiological indicators of the crop, a first stage of four months is distinguished, characterized by a limited development, and soil moisture average of 18,9 %.
The EVI, NDVI and NDWI indices shows results according to the prevailing conditions in the first four months of the crop.
After starting the rainfall, an increase in biophysical variables in sugarcane plants takes place, high values of EVI and NDVI indices are obtained by mean of satellite image monitoring, while NDWI index does not agree with the field moisture conditions.
A high correlation is obtained between the stem diameter and NDVI index described by a linear regression model at 97% adjustment.
Through monitoring with Sentinel-2 satellite, a more stable representation of the increase in vegetation is achieved, as well as more adequate values of its initial and final state are obtained.