Evaluation of Coffee Plants Transplanted to an Area with Surface and Deep Liming Based on Multispectral Indices Acquired Using Unmanned Aerial Vehicles

Abstract

The use of new technologies to monitor and evaluate the management of coffee crops allowed for a significant increase in productivity. Precision coffee farming has leveraged the development of this commodity by using remote sensing and Unmanned Aerial Vehicles (UAVs). However, the success of coffee farming in the country also resulted from management practices, including liming management in the soils. This study aimed to evaluate the response of coffee seedlings transplanted to areas subjected to deep liming in comparison to conventional (surface) liming, using vegetation indices (VIs) generated by multispectral images acquired using UAVs. The study area was overflown bimonthly by UAVs to measure the plant height, crown diameter, and chlorophyll content in the field. The VIs were generated and compared with the data measured in the field using linear time graphs and a correlation analysis. Linear regression was performed to predict the biophysical parameters as a function of the VIs. A significant difference was found only in the chlorophyll content. Most indices were correlated with the biophysical parameters, particularly the green chlorophyll index (GCI) and the canopy area calculated via vectorization. Therefore, UAVs proved to be effective coffee monitoring tools and can be recommended for coffee producers.

1. Introduction

The coffee crop has a great impact on the general scope of Brazilian exportation and has a great influence on the national agricultural chain. Despite recent advances in agriculture, the rate of increase in the crop yield has slowed down so that progress in specific areas has become insufficient to meet the challenges of the increasing demands for food, energy, and environmental resources. In this context, applications from different areas must be integrated into a multidisciplinary approach to obtain the most appropriate solutions to meet ecological and economic needs for coffee crop management [1,2]. Considering this, the generation of information through the use of images obtained with Unmanned Aerial Vehicles has become a fundamental tool in the monitoring and evaluation of coffee crops. In this context, precision coffee farming, defined in [3] as a set of techniques and technologies that help coffee growers manage their crops by considering the spatial soil parameters and individual plant characteristics, has been crucial. In addition, constant management practices and the farmers’ knowledge over the years have been essential for the success of coffee farming in the country.

In farming management, liming is a highly important soil preparation practice for coffee growing. Limestone contains calcium and magnesium. The former plays a key role in deepening the root system and increasing plant sustainability and resistance to drought, in addition to helping leaf retention, fruit ripening, and budding. The latter is a chlorophyll component and is therefore necessary for photosynthesis. Thus, a deficiency in these nutrients may compromise plant development and even lead to plant death [4,5,6].

In addition to providing nutrients, liming can increase the cation exchange capacity (CTC) and the soil pH, thereby neutralizing its acidity, which is a recurring problem in most Brazilian soils. However, calcium has limited mobility in soils. Consequently, the effect of conventional liming is restricted to the topsoil layer, that is, from 0 to 20 cm. For this reason, the root system only benefits from liming in this depth range [7,8]. Several studies have tested gypsum as an alternative to neutralize aluminum and transport calcium to deeper soil layers [9], but gypsum consists of a neutral salt and therefore is not a soil corrective [10]. Furthermore, the presence of sulfate ions may leach nutrients and the bases Ca²⁺, Mg²⁺, and K⁺, rendering them inaccessible to the roots [9,11]. Accordingly, liming deeper soil layers may be a good alternative to correct the soil pH and, consequently, increase the CTC and nutrient availability along the root system of the coffee tree, in addition to neutralizing the aluminum content.

Studies have assessed different liming depths for crops such as corn [12,13], sugarcane [14], and olive trees [15], generally finding good results for both soil pH correction and aluminum neutralization, or even for physical attributes of the soil. However, little research has been conducted on coffee growing at greater depths.

In recent years, remote sensing has been widely used in monitoring crops and decision making for better management practices. Spectral remote sensing allows for the early, efficient, objective, and non-destructive assessment of plant responses to environmental stress factors [16]. Plants react to biotic and abiotic stresses through biophysical and biochemical changes, such as reductions in the biomass and chlorophyll content and changes in the internal structures of the leaves, which can be easily detected through differences in the energy reflected in the visible spectral regions and near-infrared regions [17]. However, coffee culture still has some gaps in coffee farming because this perennial crop is highly complex. Furthermore, the emergence of Unmanned Aerial Vehicles (UAVs), which can be equipped with multispectral cameras, opens up even more possibilities since their images are superior to satellite images in spatial and temporal resolution [18].

Remote sensing has been addressed in coffee studies for disease monitoring, water stress, and the correlation of other biotic factors with abiotic factors through vegetation indices (VIs) from satellite images [19,20,21,22,23]. Recently, some studies have been conducted with UAVs equipped with multispectral cameras for crop monitoring, monitoring errors of planting, and the evaluation of fruit ripening, planting failures, and correlations between field measurements and images [24,25,26,27,28,29,30]. However, few studies have used remote sensing and UAVs to monitor agricultural practices, including the response of transplanted seedlings in areas with surface and deep liming.

Considering the above, the aim of this study was to evaluate the response of coffee plants transplanted to an area with surface and deep liming based on multispectral indices acquired using UAVs. For this purpose, a comparative study was performed between indices derived from aerial images and field data, such as the plant height, crown diameter, the leaf area index (LAI) derived from these measurements, and the chlorophyll content, obtained via temporal and correlation analyses.

2. Materials and Methods

2.1. Study Area

The experiment was conducted at Samambaia Farm, located in the municipality of Santo Antônio do Amparo in the region of Campos das Vertentes, southeastern Minas Gerais [31]. The regional climate, according to the Köppen climate classification modified in [32], is Cwb, subtropical humid, with dry winters and temperate summers, mean temperatures ranging from 18 °C to 22 °C, and a mean annual rainfall of 1650 mm. The dry season stretches from May to September, and the rainy season normally lasts until late March or early April.

The study area (Figure 1) has 0.31 ha and an average altitude of 1026 m (a.s.l.), and it is located at the geographic coordinates 507,868.215 W and 7,690,442 S (longitude and latitude in meters) in the UTM zone, 23S projection and SIRGAS 2000 geodetic reference. Coffee seedlings of the cultivar Catucaí 2SL, species Coffea arabica L., were transplanted to this area. The soil of the region is red-yellow latosol [33].

2.2. Crop Treatments and Transplanting

For coffee pre-planting, initially, a dose of 4 tons of dolomitic limestone per hectare was applied to the soil surface throughout the study area. This limestone was incorporated via disc harrowing, first with a heavy-duty disc harrow and then a lighter secondary disc harrow. Then, subsoiling was performed using a subsoiler type (Ikeda company, Marília, Brazil) with a 60 cm deep rod.

Subsequently, 35 cm deep furrows were opened, with 3.8 m spacing between rows, which were treated with gypsum, limestone, and magnesite doses (150 g/m, 600 g/m, and 100 g/m, respectively). Furrowing was repeated to incorporate the aforementioned inputs and to apply phosphate fertilizer (monoammonium phosphate (MAP)—200 g/m) and chicken manure (2000 g/m). Then, the furrows were closed using a 3-shank ripper, which also broke up the remaining clods. Last, 400 g/m dolomitic limestone was applied to the soil at depths ranging from 0.60 to 0.80 m using the implement inter-row weeder, leaving 3 rows without management, adopted as controls, to monitor this experiment.

Seedlings, previously grown in a nursery with 0.50 m × 3.8 m spacing (between plants and between rows, respectively), in perforated bags, were transplanted in November and December 2018 using a Mafes transplanter (Gralha model). Ground control points (GCPs) were used to mark the sample units to more easily identify them in the UAV images.

2.3. UAV, Flight Plans, and Multispectral Sensor

Aerial images were acquired with a UAV DJI Matrice 100, a quadcopter drone with rotating propellers, Wi-Fi connection to the base station, autonomous operation from 22 to 28 min, a 5 km range, and a 17 m/s maximum speed in Global Positioning System (GPS) mode [34].

The Precision Flight software (Version 1.3.2, Precision Hawk, Raleigh, NC, USA) [35] was used to plan the flight missions, in which the flight directions (across the rows), image overlay (80 × 80%), velocity (8 m/s), and height (50 m) were configured.

The images were acquired using the Parrot Sequoia sensor, which has an RGB sensor with a resolution of 16 megapixels (4608 × 3456) and four extra sensors with a resolution of 1.5 megapixels (1280 × 960) in the spectral bands of green (550 nm, BP 40), red (660 nm, BP 40), red-edge (735 nm, BP 10), and near-infrared (790 nm, BP 40). The camera was equipped with a solar radiation sensor coupled to a global navigation satellite system (GNSS) and a radiometric calibration panel to calibrate EMR and to standardize the data in reflectance values. The spatial resolution of the mosaics and the reflectance maps was 4.85 cm pixels.

2.4. Image Processing and Spectral Data

The images were processed in the PIX4D Mapper software, version 4.4.10 [36], according to the process flow diagram shown in Figure 2. The QGIS 3.4.14 Madeira software [37] and eCognition Developer 9.0 [38] were used to generate Vis and segment images and produce and store geospatial vector data (shapefiles) on the coffee plant samples. The values of the parameters measured in the field—the height, crown diameter, and chlorophyll—were entered into the tables outlining the attributes of the vector data.

The parameters adopted in each image processing step are outlined in

Table 1. Parameters used in the image processing flow of the Pix4D software (Mapper version 4.4.12).

Parameter	Configuration
Initial processing (alignment)
Keypoint image scale	Full
Matching image pairs	Aerial grid or corridor
Number of keypoints	Automatic
Calibration method	Alternative
Point cloud and mesh texture
Image scale	Original image size
Point density	High
Minimum number of matches	3
3D textured mesh	High resolution
DSM and orthomosaic
Resolution	Automatic
Orthomosaic	GeoTIFF
Reflectance maps and index calculation
Radiometric calibration	Irradiance sensor and calibration grid
Bands	Green, red, red edge, near infrared
Resolution	Automatic
Reflectance map	GeoTIFF

.

2.5. Vegetation Indices

The following VIs were calculated to study the spectral responses of coffee plants subjected to different treatments: normalized difference vegetation index (NDVI), green NDVI (GNDVI), normalized difference red edge (NDRE), modified soil-adjusted vegetation index 2 (MSAVI2), modified chlorophyll absorption in reflectance index 1 (MCARI1), and green chlorophyll index (GCI) (Equations (1)–(6), respectively). These selected spectral indicators are commonly used to assess productivity based on the chlorophyll content, water stress, nitrogen, and vigor, among other factors [39,40,41,42,43].

$$\mathrm{NDVI}=\frac{{\mathsf{\rho }}_{\mathrm{nir}}-{\mathsf{\rho }}_{\mathrm{red}}}{{\mathsf{\rho }}_{\mathrm{nir}}+{\mathsf{\rho }}_{\mathrm{red}}}$$

(1)

$$\mathrm{GNDVI}=\frac{{\mathsf{\rho }}_{\mathrm{nir}}-{\mathsf{\rho }}_{\mathrm{green}}}{{\mathsf{\rho }}_{\mathrm{nir}}+{\mathsf{\rho }}_{\mathrm{green}}}$$

(2)

$$\mathrm{NDRE}=\frac{{\mathsf{\rho }}_{\mathrm{nir}}-{\mathsf{\rho }}_{\mathrm{red}-\mathrm{edge}}}{{\mathsf{\rho }}_{\mathrm{nir}}+{\mathsf{\rho }}_{\mathrm{red}-\mathrm{edge}}}$$

(3)

$$\mathrm{MSAVI}2=0.5\times \left[{(2\mathsf{\rho }}_{\mathrm{nir}}+1)-\sqrt{{(2\mathsf{\rho }}_{\mathrm{nir}}+1)\times 2-8\times {(\mathsf{\rho }}_{\mathrm{nir}}-{\mathsf{\rho }}_{\mathrm{red}})}\right]$$

(4)

$$\mathrm{MCARI}1=1.5\times \left[2.5{(2\mathsf{\rho }}_{\mathrm{nir}}-{\mathsf{\rho }}_{\mathrm{red}})-1.3{(2\mathsf{\rho }}_{\mathrm{nir}}-{\mathsf{\rho }}_{\mathrm{red}})\right]$$

(5)

$$\mathrm{CIg}=\frac{{\mathsf{\rho }}_{\mathrm{nir}}}{{\mathsf{\rho }}_{\mathrm{green}}}- 1$$

(6)

where: ρ_green—the reflectance in the green spectral range; ρ_red—the reflectance in the red spectral range; ρ_(red-edge)—the reflectance in the red-edge spectral range; ρ_nir—the reflectance in the near-infrared spectral range.

The indices were calculated in QGIS using the “Raster Calculator” tool and extracted using the “Zonal Statistics” tool. The vectors were generated in the eCognition Developer software (version 9.0.1) using a process tree, according to [44], with the following sequence of tasks: image segmentation including “multi-resolution segmentation” and “spectral difference segmentation”; coffee plant sampling including “sample selection”; training and classification including “support vector machine” (parameters: brightness, shape, compactness, spectral bands, and VIs); and export of vectors in shapefile format.

2.6. Field Sampling

The plant height and crown diameter (average values in cm) measurements were taken with a graduated measuring tape, and the chlorophyll content was measured with a chlorophyll sensor (Figure 3) every two months on the same day of the flights.

The chlorophyll content was assessed using an atLEAF+ Chl meter (FT Green LLC, Wilmington, DE, USA), which returns dimensionless values. These values were converted into chlorophyll a and b contents using Equations (7) and (8), respectively [45].

$${\mathrm{Chl}}_{\mathrm{a}}=-5.774 + 0.430 \times \mathrm{atLEAF} + 0.0045 \times {\mathrm{atLEAF}}^2$$

(7)

$${\mathrm{Chl}}_{\mathrm{b}}=0.040 \times {\mathrm{atLEAF}}^{1.57}$$

(8)

where atLEAF is the chlorophyll content recorded using the chlorophyll meter (dimensionless), ${\mathrm{C}\mathrm{h}\mathrm{l}}_{\mathrm{a}}$ is chlorophyll content a (μg/cm²), and ${\mathrm{C}\mathrm{h}\mathrm{l}}_{\mathrm{b}}$ is chlorophyll content b (μg/cm²).

The LAI was calculated according to the method in [46], using the plant height and crown diameter to relate the leaf area to the crop area (surface area occupied by the crop) according to Equation (9).

$$\mathrm{LAI}=0.0134+0.7276\times {\mathrm{D}}^2\times \mathrm{h}$$

(9)

where D is the plant diameter (m) and h is the plant height (m).

Samples were collected as shown in Figure 4. To evaluate the response of the coffee plants in an area with deep liming, in the treatment area, 3 rows were marked by dividing each row into 5 sampling units, each of which consisted of 5 sampled plants, on average, spaced apart by 20 plants, totaling 15 samples. In the control area, that is, with surface liming, the same procedure was followed with one row of spacing from the previously sampled area to prevent border row effects.

2.7. Meteorological Variables

Monthly cumulative rainfall and monthly mean temperature data were retrieved from the National Aeronautics and Space Administration (NASA) initiative, Prediction Of Worldwide Energy Resources (POWER) [47], to assess the effects of meteorological conditions on coffee growth.

2.8. Statistical Analysis

Data collected in the field and extracted from images were tabulated in the Microsoft Excel software (2013) [48], where graphs were also plotted to evaluate variations in coffee plants over time.

The Shapiro–Wilk test [49] was performed in the SISVAR software to assess data normality, resulting in non-parametric values. The data were submitted to the Mann–Whitney test [50], which is recommended for non-parametric variables, in the software environment R, with a significance level ρ = 0.05.

Spearman’s rank order correlation [51] (${\widehat{\rho }}_s$) was calculated to evaluate the correlation levels between the parameters measured in the field and those extracted from the UAV images. Based on the most relevant correlations, a linear regression analysis was performed to predict the physical variables sampled in the field based on digital measurements from the UAV images. The models were developed in the software environment R [52].

3. Results and Discussion

3.1. Temporal Analysis of Meteorological Conditions and Coffee Growth Response

A temporal analysis of the meteorological conditions during the coffee cycle can be seen in Figure 5. It is possible to observe variations in the cumulative rainfall and mean air temperature during the experiment. The highest accumulation of rain was observed in the summer, and the air temperature remained constant. On the other hand, in the winter, the volume of rain and the air temperature were reduced. Climatic conditions are of great importance in the period of establishment of seedlings in the field. According to the results presented, it is possible to highlight that the meteorological conditions were favorable for the establishment and growth of the coffee plants.

For a clearer understanding of the coffee plants’ temporal dynamics, images of the same plot of coffee plants during the study months can be seen in Figure 6, showing the VIs evaluated in this study and the false-color composite RGB 4, 2, and 1 (near infrared, red, and green, respectively), which were widely used to evaluate the leaf response in terms of vegetation.

The coffee growth variables and vegetation indices’ responses for the areas with deep liming and the areas with surface liming can be seen in

Table 2. Coffee growth variables including height, crown diameter, chlorophyll content, LAI, chlorophyll a, chlorophyll b, total chlorophyll (a + b), and canopy area and the VIs, including the GCI, MSAVI2, MCARI1, GNDVI, NDRE, and NDVI, of the coffee seedlings transplanted to areas treated with deep liming and with surface liming.

	Sampled Months
Variables	May	July	September	November	January	March
	Treatment with deep liming
Height (cm)	31.01	35.72	39.87	50.35	55.07	70.45
Crown diameter (cm)	29.49	30.45	31.65	49.24	58.69	81.01
Chlorophyll content (IRC)	69.95	73.36	59.80	64.55	59.81	65.41
LAI	0.0356	0.0408	0.0460	0.1122	0.1639	0.3715
Chlorophyll a (μg/cm2)	46.7594	50.1935	36.4486	40.8418	36.1390	41.7047
chlorophyll b (μg/cm2)	31.7791	34.0795	24.9167	27.8441	24.7043	28.4203
Total chlorophyll (μg/cm2)	80.0600	86.0390	62.2095	69.7472	61.6010	71.2432
Canopy area (cm2)	631.56	777.16	721.66	1583.75	1729.22	3725.61
GCI	3.9045	2.9215	2.9382	4.3389	4.3003	6.7518
MSAVI2	0.5526	0.3604	0.2927	0.6516	0.4991	0.6086
MCARI1	0.5648	0.3193	0.2425	0.6830	0.4511	0.5375
GNDVI	0.6541	0.5824	0.5873	0.6787	0.6739	0.7660
NDRE	0.3158	0.1451	0.1143	0.1902	0.2230	0.3244
NDVI	0.7032	0.5132	0.5647	0.7865	0.7799	0.8462
	Treatment with surface liming
Height (cm)	33.16	35.53	40.25	50.11	54.33	69.11
Crown diameter (cm)	31.87	30.57	32.98	50.17	58.44	78.22
Chlorophyll content (IRC)	60.20	72.20	59.20	62.97	56.98	67.28
LAI	0.0414	0.0414	0.0518	0.1215	0.1659	0.3500
Chlorophyll a (μg/cm2)	36.6896	48.9462	35.8778	39.2296	33.3926	43.6815
chlorophyll b (μg/cm2)	25.0791	33.2493	24.5370	26.7668	22.8703	29.7396
Total chlorophyll (μg/cm2)	62.6016	83.8670	61.2292	66.9435	56.8532	74.6893
Canopy area (cm2)	614.84	824.29	818.97	1647.60	1925.71	3825.22
GCI	3.9389	2.7637	2.5180	4.3164	4.3529	6.6507
MSAVI2	0.5491	0.3496	0.2579	0.6811	0.5188	0.5990
MCARI1	0.5561	0.3095	0.2097	0.7371	0.4721	0.5320
GNDVI	0.6551	0.5697	0.5510	0.6752	0.6749	0.7607
NDRE	0.3095	0.1398	0.1140	0.1957	0.2178	0.3244
NDVI	0.7047	0.5021	0.5106	0.8014	0.7916	0.8311

. The results show the temporal variation in the studied variables, including the parameters measured via remote sensing, such as the VIs and canopy area calculated via vectorization, and the growth variables, such as the LAI and chlorophyll content, based on field sampling.

According to the results, after establishing the seedlings in the field, it was possible to observe linear increases in the values of the growth variables of height, crown diameter, LAI, and canopy area for the treatments with limestone on the surface and in depth. It was not possible to observe differences in the growth response between the studied treatments. This response may be associated with the effect of limestone on the soil, since the interaction and effect occur gradually and slowly, resulting in a more evident response throughout the crop cycle.

Figure 7 shows the damage on the plants caused by unfavorable conditions during dry and cool months. Since the photographs were recorded in July, the plants’ conditions can be associated with precipitation/temperature (Figure 5) and the spectral response (Table 2) over the sampled months, considering that plants normally take 30 to 60 days to respond to rainfall effects [53]. Therefore, in this study, with bimonthly measurements, the coffee crop year was divided into the dry season, from May to September, and the rainy season, from November to March.

In general, the physical measurements increased bimonthly, as shown by the variation in the LAI (derived from height and diameter), the canopy area calculated via vectorization, chlorophyll contents, and VIs. In September, the values of the LAI and canopy area calculated via vectorization stagnated, possibly due to winter defoliation. The LAI findings corroborate those of [46], which assessed a mean of 0.27 for the coffee plants after 15 months of planting, and [54], which showed a variation in the LAI over time similar to that found in the present study when studying evapotranspiration in young coffee plants.

The graphs also show that, in July and January, which are critical winter and summer months, the chlorophyll values were the highest and lowest, respectively. This response is related to the fact that under water stress, chloroplasts act to increase CO₂ uptake to maintain the physiological activities of plants, thereby increasing the concentration of photosynthetic pigments whilst decreasing the leaf water content, which results in higher values in chlorophyll meter readings.

The VIs decreased during the dry season (from May to September) and increased again during the rainy season. The dry season is also when coffee plants grow the least because most of this period overlaps with winter, which is a period of stress during which photoassimilates were used to grow and maintain the root system to the detriment of vegetative growth, which also explains the drop in the VI values. In addition, the drop in temperature may damage the leaves, as shown in Figure 7, which is ultimately reflected in the values of the VIs. These occurrences corroborate the studies of [55,56,57], which also reported lower VI values in the satellite images of coffee plants during the dry season, especially in the winter, and in the monitoring of coffee trees that were recently transplanted via UAVs, as evidenced in [30].

During the rainy season, the values of the VIs increased again. Some of them, such as MCARI1 and MSAVI2, were higher in November than in January. This finding suggests that these indices were more sensitive to flowering in this experimental area. Flowering usually occurs between September and November, and despite its minimum intensity in new plants, such as the coffee seedlings sampled in this study (1 year after transplanting), it can affect the responses of the coffee plants, as shown in [58], which observed an increase in the VI values in colza crops, and in [59], which reported that the NDVI and enhanced vegetation index (EVI) of Coffea arabica peaked during flowering.

In general, the variation in the values of the VIs evaluated in this study is similar to the dynamics found in [59]. The author used the moderate resolution imaging spectroradiometer (MODIS) NDVI to propose an agrometeorological–spectral monitoring model based on the responses of coffee seedlings and mature and pruned coffee plants over time. Other studies reported a similar temporal variation when evaluating the VIs of coffee plants [55,56,60]. Therefore, although this experiment involves young coffee plants, this pattern can be associated with the phenological phases of the coffee tree reported in [61]. Nevertheless, the values of this VI are lower in young coffee plants than in adult coffee plants, primarily due to the lower biomass concentration [59,62].

In the evaluation of the difference between coffee plants subjected to different treatments, that is, deep and surface liming, in Table 2, no clear differences in variations over time are identified in most variables. The lack of differences related to growth and the LAI may be explained by the failure of the seedlings’ roots to reach depths of 70 to 80 cm (liming horizon) during the study period and thus benefit from the input. Only the chlorophyll parameter showed higher values in May and January, especially in the former, in the plantation treated with deep liming than in the plantation treated without deep liming.

The VIs of GCI, NDVI, and GNDVI, which are frequently associated with plant vigor and chlorophyll a, are slightly higher in the plantation with deep liming in the dry season. The ability to maintain plant vigor in periods of stress indicates limestone action. Liming decreases soil acidity, which increases nutrient availability, thereby providing nutrition to coffee plants during critical periods. Furthermore, the inter-row weeder, an implement used during liming, may have increased the soil porosity, favoring water infiltration and storage [9,63] whilst ensuring longevity and reducing losses.

3.2. Statistical Analysis for Coffee Growth Variables

To assess the differences between the transplanted coffee seedling growth of the treatments evaluated in this study, that is, deep liming vs. surface liming, the Mann–Whitney U test was performed at a significance level of p = 0.05.

Table 3. Results of Mann–Whitney U test for non-parametric distributions assessing the existence of significant differences between plantations with and without deep liming.

Study Variable	Null Hypothesis (H0)	Mean Rank	U Value (Sum of Ranks)	Significance Value p	Decision
Height (cm)	The distribution of the variable is the same for the different managements	SL = 91.33 DL = 89.67	8219.50 8070.50	0.831	Accepts the null hypothesis (H0)
Diameter (cm)	The distribution of the variable is the same for the different managements	SL = 92.88 DL = 88.12	8359.00 7931.00	0.540	Accepts the null hypothesis (H0)
LAI	The distribution of the variable is the same for the different managements	SL = 93.98 DL = 87.02	8458.00 7832.00	0.371	Accepts the null hypothesis (H0)
Chlorophyll a + b (μg/cm2)	The distribution of the variable is the same for the different managements	SL = 82.18b DL = 98.82a	7396.50b 8893.50a	0.032 *	Rejects the null hypothesis (H0); the managements differ
Estimated crown area (cm2)	The distribution of the variable is the same for the different managements	SL = 93.43 DL = 87.57	8409.00 7881.00	0.450	Accepts the null hypothesis (H0)
GCI	The distribution of the variable is the same for the different managements	SL = 88.91 DL = 92.09	8002.00 8288.00	0.682	Accepts the null hypothesis (H0)
MSAVI2	The distribution of the variable is the same for the different managements	SL = 90.03 DL = 90.97	8103.00 8187.00	0.904	Accepts the null hypothesis (H0)
MCARI1	The distribution of the variable is the same for the different managements	SL = 90.11 DL = 90.89	8110.00 8180.00	0.920	Accepts the null hypothesis (H0)
GNDVI	The distribution of the height variable is the same for the different managements	SL = 88.32 DL = 92.68	7949.00 8341.00	0.575	Accepts the null hypothesis (H0)
NDRE	The distribution of the variable is the same for the different managements	SL = 89.78 DL = 91.22	8080.00 8210.00	0.852	Accepts the null hypothesis (H0)
NDVI	The distribution of the variable is the same for the different managements	SL = 88.58 DL = 92.42	7972.00 8318.00	0.621	Accepts the null hypothesis (H0)

* Significant differences at the significance level ρ = 0.05; SL= surface liming; DL= deep liming.

shows no significant differences between the treatments evaluated in this study, except for the chlorophyll parameter. The plantation with deep liming shows a higher chlorophyll concentration than that without deep liming. Therefore, although the other variables do not show any clear differences between managements, the chlorophyll meter can be an effective tool in identifying these differences in this type of evaluation. The U value (Mann–Whitney sum of ranks) of the plantation treated with deep liming is 8893.50 in contrast to 7396.50 in the plantation without deep liming, with a mean total chlorophyll (a + b) content of 86.04 μg/cm² versus 83.87 μg/cm, respectively.

Few studies on deep liming in coffee plantations were found in the literature, especially long-term studies or research relating deep liming to VIs, thus requiring comparisons to other cultures. According to the authors of [15], who evaluated deep liming in olive trees, which is a perennial crop, as well as in coffee trees, the use of European soil, which has naturally lower pH values, is a requirement for the success of these plants. The authors observed an improvement in the soil’s physical–chemical attributes, but the only parameter with a positive and limiting correlation with productivity was the concentration of boron; that is, applying lime did not result in any significant differences compared to conventional planting. The authors [64] observed a 66% increase in root growth up to the 60 cm horizon by increasing the limestone dose in this profile and observed a 140% increase in production.

The authors of [65] evaluated different subsurface limestone doses for the Catuaí and Icatu varieties (up to 34 cm below the surface) in soils with a high aluminum content and observed increases in the total root length/root dry matter values (cm/g) and in the total root surface/root dry matter (cm²/g), in addition to a decrease in the aluminum content at the treatment site, mainly in the Catuaí cultivar. However, biophysical parameters such as the leaf area, stem length, and shoot dry matter production remained unaffected by management.

In [66], the authors evaluated coffee plants for one year after transplanting them to a cambisol subjected to furrowing or subsoiling and to different deep liming doses at 40, 60, and 80 cm and compared the results to those of coffee seedlings transplanted to soil without these preparations. All of the treatments resulted in physical and water improvements up to the 20 cm layer. The authors observed a better distribution of water and pore volume up to 40 cm in depth. The 80 cm deep liming stood out for its good water availability in this horizon, thus proving to be an alternative for mitigating possible damage caused by water deficits during the winter. The height and NDVI measurements were performed using a Green Seeker sensor, and the readings correlated well with the physical parameters measured in the field.

The study conducted in [66] showed that adopting treatments was effective, although no significant differences were found between the treatments regarding residual amounts of calcium and magnesium. These findings may therefore be associated with the soil factor because this natural cambisol had a high clay content, mainly below 20 cm of depth, in contrast to the latosol of the present study. Thus, the management adopted in comparison to the control represented a significant increase in the soil physical attributes rather than a chemical improvement related to liming.

3.3. Spearman’s Rank Correlation

To determine the degree of correlation between the physical variables measured in the field and those derived from the UAV images in this experiment, a Spearman’s rank correlation (${\widehat{\rho }}_s$) analysis was performed, and the results are presented in

Table 4. Spearman’s rank correlation (${\widehat{\rho }}_s$) between the variables measured in the field and the variables assessed in UAV images of seedling plantations with deep liming and without deep liming.

	Field Parameters
VIs (UAV)	Height (cm)	Crown Diameter (cm)	Chlorophyll Content (IRC)	LAI	Chlorophyll a (μg/cm2)	Chlorophyll b (μg/cm2)	Chlorophyll a + b (μg/cm2)
Treatment with deep liming
GCI	0.8735	0.9281	−0.1337	0.9525	−0.1625	−0.1624	−0.1654
MSAVI	0.5117	0.5911	0.0440	0.5290	0.0139	0.0137	0.0106
MCARI1	0.3320	0.4088	0.0975	0.3351	0.0707	0.0705	0.0678
GNDVI	0.8393	0.9033	−0.1780	0.8943	−0.2076	−0.2077	−0.2107
NDRE	0.3858	0.5263	0.1873	0.5756	0.1705	0.1702	0.1685
NDVI	0.7627	0.8199	−0.3415	0.7462	−0.3701	−0.3703	−0.3732
Crown area	0.9583	0.9753	−0.2097	0.9953	−0.2385	−0.2384	−0.2413
Treatment without deep liming
GCI	0.8688	0.9224	0.0850	0.9324	0.0619	0.0618	0.0594
MSAVI	0.5167	0.5953	−0.0640	0.5128	−0.0932	−0.0936	−0.0966
MCARI1	0.3630	0.4362	−0.0903	0.3419	−0.1187	−0.1190	−0.1220
GNDVI	0.8076	0.8787	−0.0238	0.8541	−0.0492	−0.0494	−0.0521
NDRE	0.4474	0.5621	−0.0499	0.5927	−0.0673	−0.0671	−0.0689
NDVI	0.7314	0.7993	−0.2605	0.7102	−0.2877	−0.2881	−0.2909
Crown area	0.9740	0.9776	0.1797	0.9974	0.1610	0.1607	0.1588

.

All of the variables are positively correlated, except for chlorophyll, which has low or negative values across all correlated pairs. Some parameters stand out, such as the Vis of GCI, GNDVI, and NDVI. The latter has a 0.7627 correlation with the height of coffee plants growing in soil treated with deep liming, thus corroborating the findings of the authors of [66], who, when studying deep liming in coffee plantations, found that the Greenseeker NDVI was best correlated with height for liming at 60 cm of depth. The first two indices stand out, with correlation values close to 1 for parameters such as height, diameter, and LAI, which also have a correlation of 0.99 with the canopy area estimated in the software for both managements. The results indicate that these variables are good parameters for estimating or modeling some physical variables. Indices that use the green band have been applied in studies to estimate the biomass and other biophysical parameters in different crops [67,68]. The authors of [43] highlighted that the GNDVI, for example, is more sensitive to the chlorophyll concentration than the NDVI.

Other VIs do not obtain such promising correlation values, such as the MSAVI2, with moderate values (below 0.60), and the MCARI1 and NDRE, with low correlations (${\widehat{\rho }}_s$) between 0.30 and 0.40, indicating that they are not good indices for monitoring coffee plants during development.

Following the non-parametric tests performed in the previous section, no major differences were found in the correlations between the study treatments, except for chlorophyll. Although low, the correlation between chlorophyll and extra-field parameters is more positive in the plantation with deep liming than in that without deep liming. The VI of NDRE has the highest value of correlation with this parameter (${\widehat{\rho }}_s$ = 0.1705). The low correlations, even negative, with most of the VIs calculated in this study, may be related to the leaf sampling, which was performed randomly. According to [4], even in the same coffee tree, more exposed leaves (the outer section of plagiotropic branches or leaves located in the crown) are more susceptible to burning and chlorosis during periods of lower temperatures or may suffer more during intense summers. Conversely, the more internal leaves are more protected from these factors or may receive less solar radiation, and these factors affect the response of photosynthetic pigments.

The correlation values between the VIs and chlorophyll range from 0.55 to 0.78 in wheat [69]. For coffee plants, the authors of [20] developed models for predicting chlorophyll using Sentinel-2 images. The model developed by these authors explained 77% of the chlorophyll variation by using the 10 m bands of the sensor, unlike this study. The authors reported that mature coffee plants, older than 4 years, improved the model considerably compared to younger coffee plants because the near-infrared band of adult coffee plants is more reflective, and thus, the canopy area is more representative in the images. Young coffee plants, such as those sampled in this study, reduced the performance of the entire data set when modeling the total chlorophyll. Furthermore, the authors used, in addition to an SPAD chlorophyll meter, leaf samples of chlorophyll from the laboratory, which increased the accuracy of the readings.

Considering the VIs, UAV geometric measurements, and field measurements, both treatments performed similarly in this experiment. This result may be associated with the shallow root system of growing coffee plants, which therefore cannot benefit from deep liming. In turn, the positive correlation between most of the variables measured using UAVs—such as the VIs and geometric measurements—and the field measurements highlights the potential of using UAVs to estimate biophysical parameters, which shows the potential of this technology for monitoring coffee trees.

The authors suggest that further studies be carried out under conditions in which monitoring is carried out during the coffee plant’s adult stage, since the effect of deep liming treatment could be visible after a few seasons, which would require further monitoring. Also, deep liming could have a greater impact on coffee productivity instead of directly affecting the growth variables.

4. Conclusions

The results presented in this study contribute to Brazilian and world coffee growing. The use of remote sensing is widely used in agriculture. However, for coffee production areas, some gaps exist regarding the effectiveness of remote sensing in the evaluation of agricultural practices and the plant response after their establishment in the field. With our results, it is possible to observe that the use of VIs generated via multispectral images acquired using UAVs enabled the precise monitoring of the coffee crops after the seedlings were transplanted in the field. Also, the relationship between the evaluations performed in the field and the Vis’ results was satisfactory. Even when not observing the differences between the treatments studied—that is, liming application on the soil surface and in depth—this study showed a great applicability of remote sensing in the evaluation and monitoring of coffee crops.

Most of the variables determined using UAVs were highly correlated with the parameters measured in the field, except for chlorophyll, whose correlations were generally low or negative. The vegetation indices, particularly GCI, GNDVI, and NDVI, showed correlation values above 0.70 with the LAI and with the canopy area calculated via vectorization, whose correlation values with the LAI, crown diameter, and height were close to 1 in both treatments.