Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 18 July 2025

Sec. Technical Advances in Plant Science

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1621080

This article is part of the Research TopicIntegration of Advanced Technologies in Orchard ManagementView all 6 articles

Orchard variable-rate spraying method integrating GNSS and wind-excited audio-conducted leaf area density

Hangxing Zhao,,&#x;Hangxing Zhao1,2,3†Shenghui Yang,,&#x;Shenghui Yang1,2,3†Wenwei Li,Wenwei Li2,3Han Feng,Han Feng2,3Shijie Jiang,Shijie Jiang3,4Weihong Liu,Weihong Liu2,3Jingbin LiJingbin Li5Yongjun Zheng,*Yongjun Zheng2,3*Songchao Zhang*Songchao Zhang1*
  • 1Key Laboratory of Modern Agricultural Equipment, Ministry of Agriculture and Rural Affairs, Nanjing, China
  • 2College of Engineering, China Agricultural University, Beijing, China
  • 3State Key Laboratory of Intelligent Agricultural Power Equipment, Beijing, China
  • 4College of Agricultural Equipment Engineering, Henan University of Science and Technology, Luoyang, China
  • 5College of Mechanical and Electrical Engineering, Shihezi University, Shihezi, China

Introduction: Conventional air-assisted sprayers used in orchards often suffer from excessive pesticide waste, high residue levels, and uneven droplet distribution on fruit tree canopies. Precision spraying technologies have emerged to address these limitations by enabling dynamic regulation of spray parameters according to canopy characteristics. Among these, leaf area density is a key indicator for describing canopy sparseness. However, accurate and automated measurement of canopy leaf area density remains challenging due to leaf shading effects. As a result, few fully functional variable-rate spraying systems have been developed based on this parameter.

Methods: This study presents a variable-rate spraying method that integrates global navigation satellite system (GNSS) positioning with wind-excited audio-conducted estimation of canopy leaf area density. A self-propelled orchard spraying platform was developed to acquire real-time GNSS positioning and audio-conducted canopy leaf area density data. Based on this, a method was established for generating prescription maps that integrate spatial positioning and canopy density information. A variable-rate spray control model and algorithm were then constructed to regulate spray flow according to the spatial distribution of leaf area density across the orchard.

Results: Field experiments demonstrated that the system achieved a mean relative error of only 5.52% in spray flow rate regulation. Compared with conventional constant-rate spraying, the variable-rate mode reduced the longitudinal coefficient of variation (CV) of droplet deposition by 55.75% on adaxial leaf surfaces and by 33.22% on abaxial surfaces, with a maximum reduction of 62.32% in transverse CV. Ground runoff of spray solution was reduced by 62.29%, and droplet deposition density on leaf surfaces exceeded 25 droplets/cm², meeting the standard for low-volume insecticide application.

Discussion: The integration of GNSS and wind-excited audio sensing for real-time canopy density assessment enables more precise and efficient pesticide application in orchards. This system significantly improves droplet deposition uniformity while reducing environmental losses, offering a promising technical solution for the development of intelligent and sustainable plant protection equipment.

1 Introduction

Efficient control of orchard pests and diseases is essential for ensuring high fruit yield and quality (Jiang et al., 2022a; Hu et al., 2025). Chemical spraying remains the predominant method for plant protection due to its high efficiency (Rincón et al., 2020; Zheng et al., 2020; Jiang et al., 2022b). Air-assisted spraying, recognized by the Food and Agriculture Organization (FAO) of the United Nations as an advanced plant protection technology, significantly enhances pesticide utilization efficiency and reduces spray drift (Gu et al., 2022a; Yang et al., 2022; Feng et al., 2024). However, conventional uniform spraying methods apply the same spray dosage across all fruit trees based solely on preset parameters. These approaches fail to consider spatial variations in canopy structure and leaf area density (Planas et al., 2013; Colaço et al., 2019a, b).

As a result, they often lead to localized over-application and untreated zones, which not only waste pesticides and increase environmental pollution but also cause inconsistent pest control outcomes (Zhai et al., 2018; Li et al., 2024a). With the growing adoption of precision agriculture, variable-rate spraying (VRS) has emerged as a preferred solution for efficient and environmentally friendly orchard operations (Song et al., 2020). VRS systems integrate sensing, decision-making, and actuation to enable spatially differentiated spraying based on canopy variability. Such systems significantly improve pesticide use efficiency, reduce environmental impact, and optimize the balance between operational efficiency and pest control uniformity (Xue et al., 2023).

Accurate sensing of fruit tree canopies is fundamental to the effectiveness of variable-rate spraying systems (He X., 2020). Among various canopy characteristics, leaf area density is a key indicator of canopy sparsity and directly influences the spatial allocation of spray volume (Wei et al., 2022). Existing methods for leaf area density detection mainly include ultrasonic sensing (Gil et al., 2007; Palleja and Landers, 2015), light detection and ranging (LiDAR) (Hosoi and Omasa, 2007; Li et al., 2017), and multispectral or RGB imaging (Tu et al., 2019; Raj et al., 2021). Ultrasonic sensors estimate canopy volume or LEAF AREA DENSITY by emitting acoustic pulses and measuring their echoes (Palleja and Landers, 2017). However, they are highly susceptible to environmental noise. Even minor distance fluctuations can result in volume estimation errors up to 50 mm, with relative deviations reaching 30% (Palleja et al., 2010). LiDAR acquires three-dimensional point clouds of the canopy by emitting laser pulses, enabling precise reconstruction of canopy structure and estimation of leaf area index (LAI) and volume distribution (Gu et al., 2022b). Nevertheless, in regions where multiple leaf layers overlap and the interlayer spacing is smaller than the divergence angle of the laser beam, LiDAR often fails to penetrate dense canopies. It can only capture the surface geometry, lacking the ability to provide insights into internal leaf area density (Wang and Fang, 2020). Multispectral and RGB imaging estimate LEAF AREA DENSITY using vegetation indices (e.g., NDVI, GNDVI) or reconstruct canopy models through deep learning combined with stereo vision (Liu et al., 2021a). However, their accuracy is significantly compromised by fluctuations in lighting, background interference, and occlusion from overlapping foliage (Wang et al., 2024).Conventional sensors, such as LiDAR and multispectral cameras, therefore fall short of meeting the requirements for internal canopy sensing under field conditions (Li et al., 2024b). To address these limitations, novel sensing strategies must be explored to improve the accuracy and robustness of LEAF AREA DENSITY estimation. Audio-conducted sensing has recently demonstrated promising potential in agricultural applications (Cao et al., 2025). Our research group has observed that the amplitude of wind-induced canopy audio signals correlates with LEAF AREA DENSITY (Li et al., 2024b). This suggests that canopy vibration sounds induced by wind can indirectly reflect internal LEAF AREA DENSITY. (Li et al., 2024c). Unlike optical or laser-based methods, audio signals are inherently immune to lighting fluctuations and leaf occlusion. By extracting time–frequency features from wind-excited canopy audio, it becomes feasible to achieve leaf area density classification even under complex field noise conditions.

The decision-making and actuation mechanisms of variable-rate spraying systems are generally categorized into two types: sensor-based real-time variable-rate spraying systems and prescription map-based variable-rate spraying systems (Zhou et al., 2017). Sensor-based variable-rate spraying systems typically utilize a single type of sensor to detect canopy characteristics such as position, volume, and structural features in real time (Zhang et al., 2018; Gu et al., 2020; Luo et al., 2024). Spray valves are regulated at the millisecond level through pulse-width modulation (PWM) technology to achieve immediate response. However, in complex orchard environments, sensor-based systems are highly susceptible to environmental factors such as lighting conditions, wind speed, and leaf occlusion, which constrain both detection accuracy and decision-making speed (Whitfield et al., 2022). Moreover, supporting high-frequency data decoding and dynamic canopy feature computation requires each sprayer to be equipped with high-performance onboard computing and sensing devices. This significantly increases system investment costs and maintenance expenses (Han et al., 2023). In contrast, prescription map-based variable-rate spraying systems pre-generate spatially referenced prescription maps using multi-source sensor data (Román et al., 2021). During operation, RTK-GPS technology is employed to resolve the corresponding prescription units in real time (Tayari et al., 2015; Liu et al., 2024; Taseer and Han, 2024). In spraying execution, the system only needs to adjust the duty cycle of the solenoid valves via PWM according to the spray rates specified in the prescription map (Partel et al., 2021; Zhang et al., 2024). This approach significantly reduces the hardware requirements for online computation and sensing. It also enables the integration of multiple decision-making factors—such as terrain, canopy structure, and historical growth information—through multi-source data fusion (Tewari et al., 2020). As a result, prescription map-based systems achieve both high spray response speed and high crop information accuracy. Consequently, they exhibit greater application potential and broader promotional value in orchard scenarios (Yu et al., 2022).

In summary, to address the challenges in sensing costs, system integration, and adaptability in variable-rate spraying research, this study proposes an orchard variable-rate spraying system based on wind-induced audio recognition. The main contributions are as follows:

1. Design of a compact self-propelled variable-rate spraying platform: A small-scale, self-propelled spraying platform was developed. It is capable of acquiring orchard-wide positioning data via GNSS and capturing canopy leaf area density using wind-excited audio signals under field conditions.

2. Development of a prescription map generation method integrating canopy leaf area density and spatial positioning information: A mapping model was established to link leaf area density, required spray volume, and PWM duty cycle. Based on this model, a grid-based prescription map was generated by integrating GNSS-based spatial positioning and classified leaf area density levels, providing spatial guidance for precision spraying.

3. Establishment and integration of a variable-rate spraying control model: A control model and algorithm utilizing orchard-wide leaf area density distribution were formulated, enabling dynamic adjustment of spraying rates according to canopy structural variability. These developments led to the full integration of a complete variable-rate spraying system.

Through these innovations and implementations, this study aims to provide a cost-effective and dynamically adaptive LiDAR solution for smart orchards.

2 Materials and methods

2.1 Design of a compact self-propelled variable-rate spraying platform

The compact self-propelled variable-rate sprayer developed in this study (Figure 1a; Table 1) employs a modular integrated architecture driven by the Robot Operating System (ROS). The platform integrates an air-assisted spraying system, a variable-rate spray module, a wind-excited audio-conducted sensing module, an RTK-GNSS-based autonomous navigation system, and an electric tracked chassis. The air-assisted system consists of a high-flow, high-pressure axial-flow fan, an air chamber, and internal guide vanes. Based on the principle of flow displacement, the system generates a uniform dynamic pressure field within the canopy, ensuring efficient droplet penetration and uniform deposition. The spraying module incorporates a three-cylinder plunger pump with a recirculation function, an electric ball valve, and Teejet TP6503 flat-fan nozzles. Spray flow rates are precisely regulated by adjusting the PWM duty cycle applied to the valve. The wind-excited excitation unit is composed of a DPT30-80M centrifugal fan and multilayer vibration isolation and silencing structures. Adjustable wind speeds are used to stimulate branch and leaf vibrations. The audio-conducted sensing module is equipped with a PM420 microphone fitted with a windscreen and an anti-vibration mount. It enables real-time transmission of audio signals for canopy leaf area density estimation. The navigation module utilizes a Huace P3-DU dual-antenna RTK-GNSS receiver. After universal transverse mercator (UTM) projection, operation trajectories are smoothed using the least squares method (Liu et al., 2021b). Steering commands are generated in real time based on the pure pursuit algorithm (Hu et al., 2024), achieving centimeter-level path tracking and closed-loop control for variable-rate spraying.

Figure 1
Illustrates the design and control system of a self-propelled small-scale variable-rate orchard sprayer. Subfigure (a) presents the overall structural diagram, showing components such as the triplex plunger pump, spray nozzle, variable valve, axial fan, GNSS antenna, PM420 microphone, wind-excited device, and chemical tank, all mounted on an electric tracked chassis. Subfigure (b) displays the hardware block diagram of the control system, including audio signal acquisition via microphone and industrial computer, GNSS-based vehicle pose detection, wind excitation and relay modules, and PWM-controlled spray and drive units. The STM32F103C8T6 microcontroller coordinates all modules through TTL and CAN communication.

Figure 1. Self-propelled small-scale variable-rate orchard sprayer. 1. Triplex plunger pump 2. Spray nozzle 3. Relay 4. Variable valve 5. Axial fan 6. Air chamber 7. Muffler duct 8. Anti-vibration frame 9. PM420 microphone 10. Windproof sponge 11. GNSS antenna 12. P3-DU receiver 13. wind-excited device 14. Industrial compute 15. Inverter 16. Chemical tank 17. Electric tracked chassis. (a) Overall structural diagram (b) Hardware block diagram of the control system.

Table 1
www.frontiersin.org

Table 1. Main technical specifications of the compact self-propelled orchard variable-rate sprayer.

The hardware block diagram of the control system is shown in Figure 1b. The upper computer, a Zhanmei GK7000 industrial computer running under the ROS framework, integrates RTK-GNSS positioning, vehicle pose estimation, and audio data acquisition. These functions are used for prescription map generation and path planning. The lower computer, based on an STM32F103C8T6 microcontroller, performs closed-loop control of the excitation device, variable-rate valve, fan, and steering actuators via the CAN bus. The system supports both remote-control and autonomous operation modes.

2.2 Prescription map generation method integrating canopy leaf area density and spatial positioning information

2.2.1 Wind-excited audio signal acquisition and leaf area density classification

A preliminary field survey was carried out by randomly selecting 100 apple trees to assess canopy leaf area density. In the preliminary experiments, the leaf area density measurement procedure is illustrated in Figures 2a–c and includes the following steps:

Figure 2
Illustrates the distribution and estimation of canopy leaf area density. Subfigure (a) shows the canopy sampling area divided into left and right zones. Subfigure (b) presents the leaf area calculation process using a reference marker and pixel analysis. Subfigure (c) depicts a leaf count measurement setup using a cubic sampling frame in the field. Subfigure (d) provides a histogram of canopy counts across six classes of leaf area density ranging from 2.0 to 8.0 m²·m⁻³. Subfigure (e) displays a pie chart showing the proportional distribution of canopy counts by leaf area density class.

Figure 2. Canopy leaf area density distribution. (a) Canopy sampling area (b) Leaf area calculation (c) Leaf count measurement (d) Distribution of canopy counts by leaf area density (e) Proportion of canopy counts by leaf area density.

2.2.1.1 Selection and labeling of target canopy sections

A total of 100 apple trees were randomly selected from the demonstration orchard. Each tree canopy was divided into left and right sections (Figure 2a), resulting in 200 canopy regions in total. Each region was labeled numerically (1–200) to facilitate subsequent acquisition of leaf area density and wind-induced audio signals.

2.2.1.2 Leaf area density measurement of the tree canopy

As shown in Figure 2b, 50 leaves were randomly sampled from each selected tree. The average leaf area was calculated using Equation 1. Starting from region 1, a PVC square frame was used to delimit the canopy area to be measured (Figure 2c). The number of leaves within the frame was manually counted. This procedure was repeated three times at different locations within the same canopy region to obtain an average leaf count. Finally, the leaf area density of the region was calculated using Equation 2.

Sleaf=1nli=1nlNleaf,iSref,iNref,i(1)

where Sleaf is the average area of a single leaf (m²); nl is the number of randomly selected leaves; Sref,i is the area of the reference black square (4 × 10⁻4 m²); Nleaf,i is the number of pixels corresponding to the leaf; and Nref,i is the number of pixels corresponding to the reference black square.

ρleaf=NcSleafVc(2)

where ρleaf is the canopy leaf area density (m²·m⁻³); Nc is the number of leaves within the canopy; and Vc is the canopy volume (m³).

The results showed that the canopy leaf area density ranged from 2.0 m²·m⁻³ to 8.0 m²·m⁻³, with 85.5% of the samples falling within the range of 3.0 m²·m⁻³ to 6.0 m²·m⁻³ (Figures 2d, e). Due to the minimal differences in the required spray volume between adjacent density intervals, the leaf area density was classified into six levels (A–F) based on increments of 1.0 m²·m⁻³. This classification approach eliminates the need for precise numerical estimation and allows variable-rate spraying decisions to be made directly according to the assigned category. Consequently, it reduces model complexity and ensures the efficiency and accuracy of spraying decision-making.

Based on the leaf area density classification model previously developed by the research team (Li et al., 2024c), the leaf area density detection system in this study consists of three sequential stages: audio acquisition, audio enhancement, and audio-based recognition (Figure 3a).

Figure 3
Presents the workflow and schematic for leaf area density detection based on audio signals. Subfigure (a) shows a flowchart of the process. It begins with wind-excitation and acoustic signal collection, followed by time-domain processing and logarithmic transformation using FFT, then audio signal enhancement via separation and masking models. The enhanced signal is processed by an SCA-DCNN recognition model to classify leaf area density levels. Subfigure (b) illustrates a schematic of the physical setup, showing a tracked platform with a duct system directing wind toward a fruit tree canopy at a fixed distance for acoustic signal acquisition.

Figure 3. Workflow of leaf area density detection and schematic of audio signal acquisition scheme. (a) Workflow of leaf area density detection (b) Schematic diagram of audio signal acquisition.

Step 1: audio acquisition

Wind-induced canopy vibration was generated using a centrifugal fan operating at a wind speed of 10 m/s. A high-sensitivity PM420 microphone was positioned 0.3 m from the canopy. After the wind excitation device stabilized (approximately 5 s), time-domain excitation audio signals at each sampling unit were synchronously recorded using an audio acquisition device for a duration of 30 s (Figure 3b).

Step 2: audio enhancement

The audio signal was first processed using short-time Fourier transform (STFT) with an fast fourier transform (FFT) size of 1024, a frame shift of 512, and a Hamming window to extract the complex spectrogram. A logarithmic transformation was then applied to obtain the log-power spectrum (LPS). Each LPS frame was concatenated with three preceding and three following frames, resulting in a contextual feature vector with a total dimension of 7 × 513 = 3591. This vector was subsequently subjected to single-channel audio separation processing to derive the ideal ratio masks (IRM) corresponding to the excitation audio and background noise. Based on these masks, spectral masking was performed to enhance the original signal, which was then reconstructed into the time domain via inverse STFT (ISTFT).

Step 3: audio-based recognition

The enhanced signal underwent a second STFT and was resampled into a 256 × 256 log power spectrogram. This spectrogram served as a single-channel input image to an improved spatial attention convolutional neural network (SCA-DCNN)(Li et al., 2024c). The model consists of six convolution–pooling modules, incorporating multiple spatial attention mechanisms to enhance the weighted representation of time–frequency features. The final fully connected layer outputs a six-dimensional probability vector representing the leaf area density levels, which are mapped to the label set [‘A’, ‘B’, ‘C’, ‘D’, ‘E’, ‘F’].

This three-stage cascaded architecture combined signal denoising and feature extraction. It enabled real-time and reliable leaf area density classification under complex field conditions and provided critical data support for subsequent grid-based prescription spraying operations.

2.2.2 Prescription map generation method

The prescription map serves as the spatial reference for the precise control of variable-rate spraying. Its core elements include the grid size, the grid center coordinates, and the mapping of required spray volumes. Based on over four years of orchard data, including row spacing (3.0–5.0 m), plant spacing (2.0–3.0 m), canopy diameter (2.0–3.0 m), and the response characteristics of the electric ball valve, the grid size was ultimately set to 3 m × 3 m. This setting was selected to balance spatial resolution with actuator dynamic performance. The center coordinates of each grid are obtained in real time via RTK-GNSS and are subsequently converted into planar coordinates using UTM projection. The spray requirement level (Si) is determined according to the leaf area density classification results described in Section 2.2.1. The six classified levels (A–F) are mapped to corresponding spray requirement levels (1–6), thereby establishing a high-precision one-to-one correspondence between the leaf area density and the required spray volume.

The prescription map generation process consists of two stages (Figure 4).

Figure 4
Illustrates the methodology for prescription map generation in orchards. Subfigure (a) shows the orchard layout with rows of trees and a self-propelled sprayer moving between them, recording sampling points for canopy leaf area density using GNSS-based positioning. Red dots indicate sampling locations, and tree canopies are labeled to represent sampled areas. Subfigure (b) presents a pseudocode algorithm that generates the prescription map. The algorithm includes steps for grid partitioning, interpolation of leaf area density, and spatial averaging using weights. Final outputs include CSV files and optional GeoTIFF or Shapefile exports for use in spraying systems.

Figure 4. Methodology for prescription map generation. (a) GNSS positioning data acquisition and leaf area density estimation (b) Workflow for prescription map generation in orchards.

Step 1: global data acquisition

In remote-control mode, the spraying platform traverses the orchard rows at 3 m intervals. At each sampling point, the planar coordinates (xi, yi) and the corresponding leaf area density level (Li) are synchronously recorded. These data form a spatial distribution dataset of leaf area density covering the entire orchard.

Step 2: prescription map construction

Based on the established mapping relationship, the leaf area density level (Li) at each sampling point is converted into a spray requirement level (Si). All coordinate and spray requirement triplets (xi, yi, Si) are then sequentially structured and stored in a TXT file, generating standardized prescription map data.

This method features a compact structure and high operational efficiency. It supports dynamic updates and batch deployment, providing precise spatial decision-making references for subsequent path planning and closed-loop spraying operations.

2.3 Variable-rate spraying control method based on orchard-wide leaf area density distribution

2.3.1 Construction method of the variable control model based on leaf area density

To satisfy the low-volume spraying criterion of achieving no fewer than 25 droplets per square centimeter on the leaf surface, a quantitative relationship was formulated between the canopy leaf area density and the required spray volume. This relationship is expressed in Equation 3 (Jiang et al., 2025):

Q=5×105×kqρleafπCv36(1δbLd)Pb(3)

where Q is the required spray volume (L); kq is a correction coefficient (set to 3.3 in this study); ρleaf is the leaf area density (m²·m⁻³); Cv is a constant related to the nozzle, considering factors such as liquid surface tension and air density (set to 3.0 × 10⁻4 in this study); Ld is the distance between the nozzle and the canopy (m); δb is a correction coefficient accounting for spray losses (set to 0.1); and Pb is the spray pressure (MPa).

This equation provides the theoretical foundation for linking spatially variable canopy structure to the spray dosage needed in precision orchard operations.

2.3.2 Variable-rate spraying algorithm

A fitting function q = f(xp) between the PWM duty cycle (xp) and the spray flow rate (q) was established through calibration experiments. Combined with the prescription map grid size of 3 m and an operational speed of 0.5 m/s, a closed-loop variable-rate spraying control model was constructed. The control algorithm was designed based on a four-layer architecture of “Data–Perception–Decision–Execution” (Figure 5a). The specific workflow is described as follows.

Figure 5
Shows the variable-rate spraying control method based on orchard-wide leaf area density distribution. Subfigure (a) presents a pseudocode algorithm that loads prescription map data, calculates the distance between the sprayer’s real-time position and target spray points, and adjusts the PWM signal to control the spray valve accordingly. Subfigure (b) illustrates a dual-part workflow: the left part shows the generation of the prescription map, beginning with GNSS and audio data acquisition, followed by signal enhancement, coordinate transformation, and map construction. The right part displays the operational flow of the sprayer, including loading the map, autonomous navigation, application level reading, valve adjustment, and spray execution.

Figure 5. Variable-rate spraying control method based on orchard-wide leaf area density distribution. (a) Variable-rate spraying control algorithm (b) Workflow of the sprayer.

Step 1: data layer

The prescription map TXT file is imported, extracting the center coordinates xtarget, ytarget) of each grid along with the corresponding spray requirement level Si, which are stored as a data array(xtarget, ytarget, Si).

Step 2: perception layer

The current position of the sprayer (xcurrent, ycurrent) is obtained in real time via RTK-GNSS. The Euclidean distance between the current position and each grid center is calculated using Equation 4:

Dactual=(xcurrentxtarget)2+(ycurrentytarget)2(4)

Step 3: decision layer

When the actual distance Dactual is less than half of the predefined grid size (3 m), the system determines that the target grid has been reached and reads the corresponding spray requirement level (Si).

Step 4: execution layer

Based on the spray requirement level Si, the spray volume is calculated using the model in Equation (3) and the calibration curve q = f(xp) to determine the required duty cycle (xp). A PWM command is then issued via the CAN bus to the lower controller to drive the variable-rate valve, achieving precise flow rate control.

This control method enables real-time optimization of the spray volume for each grid, significantly improving pesticide utilization efficiency and spray uniformity.

2.3.3 Workflow of the variable-rate spraying system

As shown in Figure 5b, the operational workflow adopts an asynchronous architecture of “sensing–decision–actuation” and consists of three main stages:

1. Static multimodal sensing: After the sprayer is remotely driven into the target inter-row space and comes to a stationary position, the wind-excited device and the microphone synchronously collect canopy audio spectra and RTK-GNSS data at predefined grid intervals. The collected audio data were processed through the Audio Enhancement and Audio-Based Recognition stages, resulting in the classification of the corresponding leaf area density. The corresponding grid coordinates are determined through UTM projection. The classification results and positional data are saved as separate text files.

2. Dynamic decision-making optimization: Based on the leaf area density classification and coordinate data of each grid, the required spray level is computed. A point-based prescription map is then generated, and the optimal spraying path is fitted accordingly.

3. Precision spraying control: The system switches to autonomous navigation mode. The industrial computer (upper controller) continuously aligns the vehicle position with the corresponding prescription grid in real time. The lower controller dynamically adjusts the PWM duty cycle according to the spray volume required for each grid. This enables accurate control and ensures precision variable-rate spraying.

2.4 Orchard field testing and validation

2.4.1 Experimental site: demonstration orchard

To enhance the representativeness and generalizability of this study, a standardized demonstration orchard was selected in Hongcaohe, Fuping County, Baoding City, Hebei Province, China (Figure 6a). The orchard featured a plant spacing of 2.0 m, a row spacing of 4.0 m, and an average tree height of 3.5 m. At this site, canopy leaf area density detection, wind-induced audio signal acquisition, spray flow rate calibration, and droplet deposition tests were conducted sequentially. The layout and field environment for each experiment are illustrated in Figure 6.

Figure 6
Shows the experimental site and layout for field testing of the variable-rate sprayer. Subfigure (a) presents a satellite image with a red polygon outlining the orchard boundary and four marked GPS waypoints. Subfigure (b) shows the sprayer at the spray flow calibration site, equipped with two white containers for measuring spray output. Subfigure (c) illustrates the sprayer operating within an orchard surrounded by dense fruit trees. Subfigure (d) depicts the droplet deposition test setup, where red-and-white poles are systematically placed under the canopy to support water-sensitive papers or collectors for evaluating spraying performance.

Figure 6. Experimental sites and test layout. (a) Location of the experimental site (b) Spray flow calibration test site (c) Field operation scene of the sprayer (d) Layout of the droplet deposition test site.

2.4.2 Spray volume error analysis test

To verify the calibration accuracy and control performance of the variable-rate spraying system, spray volume calibration and precision testing were conducted in August 2024 (Figure 6b).

2.4.2.1 PWM–flow rate calibration

In the single-side variable spraying mode, only the right-side valve was calibrated. The experimental setup included the self-propelled sprayer, a water bucket, an electronic scale, and a stopwatch. The calibration procedure was conducted as follows. After the sprayer pump was activated and stabilized, the nozzle was positioned over an empty water bucket. The system was operated sequentially at duty cycle settings ranging from 0% to 100% in 10% increments, with each setting maintained for 1 minute. After closing the valve, the mass of the liquid collected in the bucket was measured. Each duty cycle condition was tested three times, and the average value was calculated. If the relationship between the flow rate and the duty cycle exhibited nonlinearity, the duty cycle intervals were refined within the relevant ranges, and additional tests were performed to improve the calibration curve. The resulting fitting function q = f(xp) provided a quantitative basis for variable-rate spray control.

2.4.2.2 Variable-rate spraying accuracy test

In a 60 m long standardized testing area, the spraying operation zone was divided into six consecutive regions according to the spray requirement levels (Levels 1–6) specified in the prescription map. A buffer zone was established before and after each operational region to minimize transitional effects (Figure 7). The theoretical spray volume Q’i for each region was calculated based on the grid area and the corresponding prescription level. During the experiment, the nozzle was positioned over an empty water bucket, and spraying operations were sequentially performed for each spray requirement region. The actual spray volume (Qi) was measured by weighing the liquid collected in the bucket. Each region was tested three times to reduce random errors, and the mean value of the three measurements was used for analysis. The relative error between the actual and theoretical spray volumes for each region was calculated using Equation 5 to quantify the accuracy of the variable-rate spraying operation.

Figure 7
Illustrates the control accuracy test layout for variable-rate spraying. A rectangular orchard area is divided into a grid of 3 m × 3 m cells, with a central 60 m spray path running vertically. Trees are positioned on both sides of the grid. Several grid cells are labeled with different spray levels, such as

Figure 7. Control accuracy test scheme for variable-rate spraying.

δs=|Qi¯Qi'|Qi'×100%(5)

This testing method effectively evaluates the stability and reliability of the variable-rate spraying algorithm and provides robust performance verification for the practical application of variable-rate spraying technology.

2.4.3 Droplet deposition and runoff test

2.4.3.1 Test location

To comprehensively evaluate the canopy coverage uniformity and ground runoff suppression performance of the variable-rate spraying system, comparative experiments between variable-rate and constant-rate spraying modes were conducted in August 2024 at the demonstration orchard.

The experiments were performed in accordance with the Chinese national standards NY/T 992–2006 Air-assisted Orchard Sprayer Operation Quality and JB/T 9782–2014 General Test Methods for Plant Protection Machinery (Figures 6c, d).

2.4.3.2 Test design

The specific implementation of the droplet deposition test is illustrated in Figure 8. Three representative apple trees, labeled Tree 1, Tree 2, and Tree 3 (Figure 8a), were selected within the test area as sample trees. In the middle and lower canopy of each tree, three vertical testing layers were established at heights of 1.8 m (upper layer), 1.2 m (middle layer), and 0.6 m (lower layer). Each testing layer was further divided horizontally into nine sampling positions, including four outer positions (A, B, C, D), four inner positions (a, b, c, d), and one central trunk position (o). At each sampling position, water-sensitive papers were affixed to both the adaxial and abaxial surfaces of leaves using paper clips (Figure 8b).

Figure 8
Illustrates the experimental layout for droplet deposition testing. (a) Shows the sprayer operation scheme in a two-row orchard setup, with red dots marking sampling positions and three fruit trees (labeled 1–3) designated for canopy deposition evaluation. The sprayer moves in opposite directions along the rows. (b) Depicts the sampling point distribution within the tree canopy, divided into upper, middle, and lower layers, and further segmented into outer and inner canopy zones. (c) Displays a top-down view of the ground-level sampling grid around the tree, with nine labeled positions (G1–G9) arranged in a 0.8 m × 0.8 m square pattern beneath the canopy for droplet collection.

Figure 8. Experimental scheme for droplet deposition testing. (a) Sprayer operation scheme (b) Sampling point layout within the canopy (c) Ground-level sampling point layout.

To assess ground runoff, nine ground sampling points, labeled G1 to G9, were evenly distributed beneath the vertical projection of the canopy (Figure 8c). Water-sensitive papers were attached to ground marker rods using bulldog clips. Both variable-rate and constant-rate spraying tests were conducted under the same operational conditions. The driving speed was maintained at 0.5 m/s, the spraying pressure was set to 0.5 MPa, and water was used as the spraying medium. In the variable-rate spraying group, valve openings were dynamically adjusted according to the prescription map. In contrast, in the constant-rate spraying group, the valves remained fully open during operation.

2.4.3.3 Data analysis method

After the experiments, all water-sensitive papers were collected and analyzed using DepositScan™ software, developed by the Chinese Ministry of Agriculture. The software was used to extract the droplet deposition volume per unit area (μL/cm²) and the droplet density (number/cm²). To minimize individual variability, the data at each sampling position were averaged across the three sampled trees. During data analysis, the droplet deposition count, droplet deposition volume, and the coefficient of variation (CV) of deposition volume were calculated for each sampling position. These indicators were used to comprehensively evaluate the uniformity and effectiveness of the spraying operation. Specifically, the droplet deposition count was used to assess whether the deposition met the minimum application requirement (≥25 droplets/cm²). The droplet deposition volume and its CV were used to quantify the uniformity of droplet distribution. Furthermore, the canopy droplet deposition characteristics were analyzed in two dimensions: longitudinally (across different canopy height layers) and transversely (across the outer, inner, and trunk regions).

The specific calculation methods are described as follows.

The droplet deposition volumes at the nine sampling positions in the upper, middle, and lower canopy layers are denoted as X1j, X2j, and X3j, respectively, where j = 1, 2,…, 9. The longitudinal coefficient of variation (CV) of droplet deposition was calculated according to Equations 69, following the steps below:

Step 1: Calculate the average droplet deposition volume for each canopy layer.

{X¯u=19j=19X1jX¯m=19j=19X2jX¯l=19j=19X3j(6)

In the equations, X¯u, X¯m, and X¯l represent the average droplet deposition volumes for the upper, middle, and lower canopy layers, respectively.

Step 2: Calculate the overall mean droplet deposition volume across the three canopy layers, denoted as μv.

μv=X¯u+X¯m+X¯l3(7)

Step 3: Calculate the standard deviation of droplet deposition volume across the three canopy layers, denoted as σv.

σv=(X¯uμv)2+(X¯mμv)2+(X¯lμv)23(8)

Step 4: Calculate the longitudinal coefficient of variation of droplet deposition volume, denoted as Cv.

Cv=σvμv×100%(9)

where Cv is the longitudinal coefficient of variation; μv is the overall mean droplet deposition volume across the three canopy layers; and σv is the standard deviation of droplet deposition volume across the three canopy layers.

The transverse coefficient of variation was calculated using Equations 1013.

X¯j=X1j+X2j+X3j9,   j=1,2,,9(10)
μh=19j=19X¯j(11)
σh=19j=19(X¯jμh)2(12)
Ch=σhμh×100%(13)

where Ch is the transverse coefficient of variation; X¯j is the mean droplet deposition volume across the three canopy layers at a single sampling point; μh is the overall mean droplet deposition volume across the nine sampling points; and σh is the standard deviation of droplet deposition volume across the nine sampling points.

This methodology systematically and quantitatively revealed the differences in canopy coverage uniformity and ground runoff between the variable-rate and constant-rate spraying modes. It provided a solid foundation for the performance validation and optimization of precision spraying technologies.

3 Results and discussion

3.1 Spray volume error analysis test results

3.1.1 Experimental calibration results of spray flow rate in response to PWM duty cycle

Figure 9a presents the calibration results between the PWM duty cycle (xp) and the spray flow rate (q), based on a total of 23 collected data sets. Analysis of the calibration data revealed the following patterns. When xp was less than 33%, the spray flow rate (q) remained approximately zero. When xp ranged from 33% to 55%, q increased approximately linearly with xp. When xp exceeded 55%, q tended to saturate. Therefore, as shown in Figure 9b, the data within the range of xp from 33% to 55% were used to fit the functional relationship between the spray flow rate (q) and the PWM duty cycle (xp), as expressed in Equation 14.

Figure 9
Displays the calibration results of the PWM (pulse-width modulation) duty cycle versus spray flow rate. (a) Shows a scatter plot of experimental data where spray flow rate (L/min) increases non-linearly with PWM duty cycle (%), eventually reaching a saturation plateau above 60%. (b) Zooms into the linear response region (approximately 33%–55%), presenting a linear regression line fitted to the data. The fitting equation is shown as q = 0.0074x − 1.9438 with R² = 0.9312, indicating a strong linear correlation in this range.

Figure 9. PWM–spray flow rate calibration results. (a) Experimental data (b) Duty cycle–spray flow rate fitting curve.

q=0.0744xp1.9438  (33x55)(14)

3.1.2 Variable-rate spraying control model

Based on the operational parameters of a driving speed of 0.5 m/s and a prescription map grid size of 3 m × 3 m, a variable-rate spraying control model was established to translate spatially distributed canopy leaf area density levels into quantitative spray control signals. The model defines a one-to-one mapping between leaf area density levels, the corresponding spray requirement levels, the required spray volume per grid, and the PWM duty cycle applied to the electric spray valve.

The leaf area density was classified into six discrete levels (A–F), based on the field survey results described in Section 2.2. Each level reflects a distinct canopy density category, with level A indicating sparse foliage and level F representing the densest canopy conditions. For each leaf area density level, the corresponding spray requirement level was assigned from 1 to 6. The required spray volume was determined using Equation 3, accounting for droplet deposition thresholds and nozzle characteristics. Subsequently, calibration experiments were conducted to establish the functional relationship between spray volume and PWM duty cycle. The PWM duty cycles corresponding to each spray requirement level were derived from the fitted flow rate curve Equation 14. The complete variable-rate spraying control model is summarized in Table 2. This model enables real-time, grid-based adjustment of spray output according to canopy structure variability, providing a practical basis for precise and efficient field operations.

Table 2
www.frontiersin.org

Table 2. Variable-rate spraying control model.

3.1.3 Prescription map integration

Following the method described in Section 2.3, a prescription map with a grid size of 3 m × 3 m was constructed for three rows of apple trees within the test area. The navigation path was generated based on the center coordinates of each grid and the corresponding spray requirement levels (Figure 10). The prescription map and the navigation path were successfully imported into the sprayer’s control system. These data provided spatial decision-making references for the subsequent closed-loop spraying operations.

Figure 10
Presents the orchard prescription map and corresponding sprayer navigation path. The chart uses X- and Y-coordinates (in meters) to depict the spatial layout. Blue dots represent sampling points arranged in two diagonal rows, with orange lines indicating the navigation path of the sprayer through these points. Each sampling point is annotated with a number from 1 to 6, signifying the required application level at that position. The path zigzags across the orchard rows, enabling variable-rate spraying based on localized prescription data.

Figure 10. Orchard prescription map and sprayer navigation path.

3.1.4 Spray volume error analysis in field tests

The accuracy test results of variable-rate spraying are summarized in Table 3. The absolute relative errors between the average actual spray volume and the theoretical spray volume across different spray requirement levels ranged from 0.00% to 14.81%. The overall mean relative error was 5.52%. This result met the field accuracy requirement for orchard variable-rate spraying.

Table 3
www.frontiersin.org

Table 3. Accuracy test results of variable-rate spraying.

3.2 Droplet deposition and runoff test results

3.2.1 Comparison of spraying distribution patterns

The canopy droplet deposition patterns under variable-rate and constant-rate spraying modes are shown in Figure 11. Both spraying modes met the minimum application requirement of 25 droplets/cm². However, the variable-rate spraying mode significantly reduced pesticide usage compared to the constant-rate spraying mode.

Figure 11
Illustrates the distribution of canopy droplet deposition under variable-rate and constant-rate spraying conditions, with both count and volume representations. Four 3D bar charts are displayed: (a) Variable-rate spraying: droplet count distribution across sampling points and canopy layers, including upper, middle, and lower, both adaxial and abaxial leaf surfaces. (b) Constant-rate spraying: droplet count distribution showing different spatial patterns compared to (a). (c) Variable-rate spraying: droplet volume distribution (µL·cm⁻²), highlighting how droplet amounts vary across canopy positions. (d) Constant-rate spraying: droplet volume distribution, showing generally higher or more uniform values across all surfaces compared to (c). Each subfigure categorizes deposition across six leaf surfaces (upper/middle/lower and adaxial/abaxial) at various sampling points, enabling comparative assessment of spraying strategies.

Figure 11. Canopy droplet deposition: count and amount distributions. (a) Variable-rate spraying - count distribution (b) Constant-rate spraying – count distribution (c) Variable-rate spraying – amount distribution (d) Constant-rate spraying – amount distribution.

3.2.2 Analysis of longitudinal droplet deposition distribution characteristics

The average droplet deposition on the adaxial and abaxial leaf surfaces was measured across the upper, middle, and lower canopy layers. The corresponding longitudinal coefficients of variation (CV) were also calculated to evaluate deposition uniformity. The results are summarized in Figure 12 and Table 4.

Figure 12
Bar chart comparing average droplet deposition counts (no./cm²) on adaxial and abaxial leaf surfaces across upper, middle, and lower canopy layers under two spraying strategies: variable-rate spraying (blue bars) and constant-rate spraying (orange bars). The x-axis represents different sampling point locations (e.g., upper adaxial, middle adaxial, lower adaxial, upper abaxial, etc.), while the y-axis indicates droplet deposition counts. Constant-rate spraying shows higher deposition on most surfaces, particularly the lower abaxial surface, whereas variable-rate spraying shows more uniform distribution across the canopy.

Figure 12. Average droplet deposition count distribution on adaxial and abaxial leaf surfaces across canopy layers.

Table 4
www.frontiersin.org

Table 4. Longitudinal coefficient of variation (CV) across canopy layers.

Analysis of Figure 12 and Table 4 indicates that the variable-rate spraying mode not only met the canopy spraying requirements but also reduced pesticide usage and improved longitudinal droplet deposition uniformity within the canopy. The detailed results are presented as follows:

3.2.2.1 Spray coverage requirement

Under both the constant-rate and variable-rate spraying modes, the average droplet deposition counts on the adaxial and abaxial surfaces of leaves in each canopy layer exceeded 25 droplets/cm². This meets the minimum application requirement. The lowest deposition count was observed on the abaxial surface of the middle canopy layer under the variable-rate spraying mode. The value was 54 droplets/cm², which still significantly exceeded the threshold.

3.2.2.2 Reduction in pesticide usage

The variable-rate spraying mode effectively reduced the amount of pesticide applied while still meeting the application requirements. Under the constant-rate spraying mode, the average droplet deposition volumes in the upper, middle, and lower canopy layers were 103.87%, 243.18%, and 370.15% higher, respectively, than those under the variable-rate spraying mode.

3.2.2.3 Improved longitudinal deposition uniformity

The variable-rate spraying mode significantly improved the longitudinal uniformity of droplet deposition. Under the constant-rate spraying mode, the longitudinal coefficients of variation (CV) for the adaxial and abaxial leaf surfaces were 62.23% and 70.32%, respectively. In contrast, under the variable-rate spraying mode, the corresponding CVs were reduced to 6.48% and 37.10%. This represents reductions of 55.75% and 33.22% for the adaxial and abaxial surfaces, respectively.

3.2.3 Analysis of transverse droplet deposition distribution characteristics

Figure 13 shows the distribution of average droplet deposition counts on the adaxial and abaxial leaf surfaces at each sampling point. The results demonstrate that, under both spraying modes, the average droplet counts at all sampling points on both leaf surfaces exceeded 25 droplets/cm², thereby meeting the minimum application requirement. The lowest droplet count was recorded at the abaxial surface of sampling point b under the variable-rate spraying mode, with a value of 44 droplets/cm². Despite being the lowest, this value still significantly surpassed the threshold.

Figure 13
Bar charts showing the average droplet deposition count (no./cm²) at each sampling point on (a) adaxial and (b) abaxial leaf surfaces under two spraying modes: variable-rate spraying (blue) and constant-rate spraying (orange). The x-axis represents sampling points labeled A–D and a–d, with

Figure 13. Distribution of average droplet deposition counts on adaxial and abaxial leaf surfaces at each sampling point. (a) Adaxial surface (b) Abaxial surface.

The average droplet deposition volumes on the adaxial and abaxial leaf surfaces at each sampling point (Figure 14) and the canopy transverse coefficients of variation (CV) (Table 5) were analyzed to further investigate the uniformity of droplet deposition distribution in the transverse direction under both spraying modes. The results are summarized as follows:

1. Under constant-rate spraying, the transverse CVs of droplet deposition on the adaxial and abaxial leaf surfaces were 52.50% and 110.13%, respectively.

2. Under variable-rate spraying, the transverse CVs on the adaxial and abaxial surfaces were 51.82% and 47.81%, respectively.

Figure 14
Four circular heatmaps illustrating the average droplet deposition amount (µL/cm²) across different sampling points on adaxial and abaxial leaf surfaces. Each circle is labeled with eight directional sampling points (A–D, a–d) and a center point (O), with a color scale from 0 (blue) to 30 (red). (a) Adaxial surface under variable-rate spraying: moderate to low deposition with spatial variability; lowest around point

Figure 14. Distribution of average droplet deposition amounts on adaxial and abaxial surfaces across sampling points. (a) Adaxial surface under variable-rate spraying (b) Adaxial surface under constant-rate spraying (c) Abaxial surface under variable-rate spraying (d) Abaxial surface under constant-rate spraying.

Table 5
www.frontiersin.org

Table 5. Transverse coefficient of variation (CV) within the canopy.

These findings indicate that the variable-rate spraying mode effectively improved the uniformity of droplet deposition across sampling points on both leaf surfaces. In particular, the transverse CV on the abaxial surface was reduced by 62.32% compared to constant-rate spraying.

3.2.4 Ground runoff analysis

The ground runoff distribution under the two spraying modes is illustrated in Figure 15. It can be visually observed that the ground pesticide runoff under the variable-rate spraying mode was significantly lower than that under the constant-rate spraying mode.

Figure 15
Two square heatmaps showing the distribution of ground-level droplet deposition counts (no./cm²) under different spraying strategies, labeled with nine sampling points (G1–G9) arranged in a 3 × 3 grid. The color scale ranges from 0 (blue) to 40 (red). (a) Variable-rate spraying: lower overall deposition with smooth spatial variation. Blue and green dominate, especially at the center (G5) and bottom edge (G8), indicating reduced ground runoff. (b) Constant-rate spraying: high and uniform deposition across most regions. Intense red at G5 and extending vertically and horizontally, indicating significant ground runoff concentration.

Figure 15. Distribution of ground runoff under variable-rate and constant-rate spraying. (a) Variable-rate spraying (b) Constant-rate spraying.

The average ground pesticide runoff was 36.91 μL/cm² for the constant-rate spraying mode and 13.92 μL/cm² for the variable-rate spraying mode. Compared to constant-rate spraying, variable-rate spraying reduced ground runoff by approximately 62.29%.

3.3 Discussion

In this study, an integrated system combining audio signal analysis, canopy leaf area density classification, and variable-rate spray control was developed. The system incorporated wind-excited audio-conducted sensing with GNSS-based prescription mapping technologies. It effectively addressed the limitations of conventional sensors in detecting leaf area density within dense canopies. Moreover, it mitigated the decision-making delays commonly encountered in real-time variable-rate spraying systems. The proposed system also reduced hardware requirements for online computation and sensing. Its feasibility was fully validated through field experiments.

Compared with previous studies, most existing research employs ultrasonic sensors (Palleja and Landers, 2015), multispectral/red-green-blue (RGB) cameras (Tu et al., 2019; Raj et al., 2021), and LiDAR (Gu et al., 2022b) to detect canopy thickness, measure canopy height, and estimate canopy volume based on known travel speed. These methods generally assume that canopy foliage is uniformly distributed, thereby neglecting the internal heterogeneity of leaf area density. In contrast, the system proposed in this study is unaffected by audio diffraction and occlusion. It enables stable acquisition of excitation signals from deep within dense canopies, thereby facilitating the estimation of internal leaf area density. Additionally, sensor-based real-time variable-rate spraying systems typically rely on a single type of sensor (Zhang et al., 2019; Luo et al., 2024). These systems require real-time decoding and dynamic analysis of complex canopy features, which often results in low detection accuracy, slow response speeds, and limited decision-making capabilities. By integrating RTK-GNSS data with orchard-wide leaf area density distribution, a predefined prescription map was generated. This prescription map allows spray decisions to be executed simply by reading grid-based text files. The approach ensures high precision while significantly reducing reliance on onboard computational power. It demonstrates notable advantages in terms of cost efficiency, system deployment, and real-time performance.

Although the proposed system exhibits excellent performance, the current leaf area density estimation model is based on a limited set of samples from specific fruit tree varieties and growth stages. Therefore, further optimization is necessary through large-scale trials across multiple regions and crop types. Considering the system’s advantages of low cost, high robustness, and ease of deployment, it holds potential for extension to precision pesticide applications in modern orchards, tea plantations, and understory economic crops. In the future, the system could be enhanced by integrating meteorological sensors and multimodal data sources, such as vision and LiDAR. Environmental adaptation could be further improved through the application of attention mechanisms or graph neural networks. Furthermore, incorporating reinforcement learning-based online optimization algorithms may enable intelligent joint scheduling of valve opening degrees and operational paths. These advancements are expected to enhance resource utilization efficiency and environmental sustainability, thereby promoting the broader industrial application of this technology.

4 Conclusion

This study presents an orchard precision spraying system that integrates wind-excited audio-conducted sensing, grid-based prescription mapping, and variable-rate spray control. A high-precision mapping relationship among leaf area density, required spray rate, and PWM duty cycle was established. The system was validated through field experiments conducted in a commercial orchard. The main conclusions are as follows:

(1) Completion of prescription map generation and system integration: A 3 m × 3 m grid-based prescription map was successfully generated using RTK-GNSS-acquired operation trajectories and classified leaf area density data. A mapping curve linking leaf area density, required spray rate, and PWM duty cycle was established through calibration experiments, achieving a goodness of fit of R² = 0.93. The system integrated an electric tracked chassis, an autonomous navigation module, variable-rate spraying actuators, and a ROS-based control platform. It enabled online loading of the prescription map, path fitting, and coordinated control between the upper and lower computers for dynamic adjustment of the spray valves.

(2) Validation of spraying performance and operational benefits: The variable-rate spraying system satisfied the canopy spraying requirements while significantly reducing pesticide usage and improving spray uniformity. Experimental results showed that, under the variable-rate spraying mode, the longitudinal coefficients of variation (CV) for droplet deposition were 6.48% on the adaxial surface and 37.10% on the abaxial surface. Compared to constant-rate spraying, these values represented reductions of 55.75% and 33.22%, respectively. The transverse CVs were 51.82% on the adaxial surface and 47.81% on the abaxial surface, with reductions of 0.68% and 62.32%, respectively. Furthermore, the ground runoff of spray solution was reduced to 13.92 μL/cm², representing a 62.29% decrease compared to constant-rate spraying.

In conclusion, this study provides a complete technical system solution for orchard plant protection. The proposed methodology can be extended to similar orchard scenarios, such as pear and peach orchards. It offers practical value in enhancing spraying precision and reducing pesticide non-point source pollution.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.

Author contributions

HZ: Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing, Conceptualization, Data curation, Formal analysis. SY: Conceptualization, Software, Validation, Writing – review & editing, Methodology, Project administration, Resources, Supervision. WWL: Conceptualization, Data curation, Formal analysis, Methodology, Supervision, Validation, Visualization, Writing – review & editing. HF: Investigation, Supervision, Validation, Writing – review & editing. SJ:. WHL: Formal analysis, Supervision, Visualization, Writing – review & editing. JL: Conceptualization, Methodology, Project administration, Supervision, Writing – review & editing. YZ: Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Resources, Supervision, Visualization, Writing – review & editing. SZ: Conceptualization, Project administration, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This paper was funded by the Open Project of Key Laboratory of Modern Agricultural Equipment, Ministry of Agriculture and Rural Affairs (grant number No. 2023005) National Research Natural Science Foundation of China (NSFC) (grant number:32171901), the Xinjiang Production and Construction Corps Agricultural Key Core Technology Research Project (Grant No. NYHXGG2023AA402), and the China Postdoctoral Science Foundation (Grant No. 2024T171010).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Cao, Y., Yin, Z., Duan, Y., Cao, R., Hu, G., Liu, Z., et al. (2025). Research on improved sound recognition model for oestrus detection in sows. Comput. Electron. Agric. 231, 109975. doi: 10.1016/j.compag.2025.109975

Crossref Full Text | Google Scholar

Colaço, A. F., Molin, J. P., Rosell-Polo, J. R., and Escolà, A. (2019a). Spatial variability in commercial orange groves. Part 1: canopy volume and height. Precis. Agric. 20, 788–804. doi: 10.1007/s11119-018-9612-3

Crossref Full Text | Google Scholar

Colaço, A. F., Molin, J. P., Rosell-Polo, J. R., and Escolà, A. (2019b). Spatial variability in commercial orange groves. Part 2: relating canopy geometry to soil attributes and historical yield. Precis. Agric. 20, 805–822. doi: 10.1007/s11119-018-9615-0

Crossref Full Text | Google Scholar

Feng, H., Xu, P., Yang, S., Zheng, Y., Li, W., Liu, W., et al. (2024). Back pressure generated by downwash and crosswind on spatial atomization characteristics during uav spraying: cfd analysis and verification. Pest Manage. Sci. 80, 1348–1360. doi: 10.1002/ps.7865

PubMed Abstract | Crossref Full Text | Google Scholar

Gil, E., Escolà, A., Rosell, J. R., Planas, S., and Val, L. (2007). Variable rate application of plant protection products in vineyard using ultrasonic sensors. Crop Prot. 26, 1287–1297. doi: 10.1016/j.cropro.2006.11.003

Crossref Full Text | Google Scholar

Gu, C., Wang, X., Yang, F., and Zhai, C. (2020). Research progress on variable-rate spraying technology in orchards. Appl. Eng. Agric. 36, 927–942. doi: 10.13031/aea.14201

Crossref Full Text | Google Scholar

Gu, C., Zhao, C., Zou, W., Yang, S., Dou, H., Zhai, C., et al. (2022b). Innovative leaf area detection models for orchard tree thick canopy based on lidar point cloud data. Agriculture 12, 1241. doi: 10.3390/agriculture12081241

Crossref Full Text | Google Scholar

Gu, C., Zou, W., Wang, X., Chen, L., and Zhai, C. (2022a). Wind loss model for the thick canopies of orchard trees based on accurate variable spraying. Front. Plant Sci. 13, 1010540. doi: 10.3389/fpls.2022.1010540

PubMed Abstract | Crossref Full Text | Google Scholar

Han, L., Wang, S., Wang, Z., Jin, L., and He, X. (2023). Method of 3d voxel prescription map construction in digital orchard management based on lidar-rtk boarded on a ugv. Drones 7, 242. doi: 10.3390/drones7040242

Crossref Full Text | Google Scholar

Hosoi, F. and Omasa, K. (2007). Factors contributing to accuracy in the estimation of the woody canopy leaf area density profile using 3d portable lidar imaging. J. Exp. Bot. 58, 3463–3473. doi: 10.1093/jxb/erm203

PubMed Abstract | Crossref Full Text | Google Scholar

Hu, C., Ru, Y., Fang, S., Rong, Z., Zhou, H., Yan, X., et al. (2025). Orchard variable rate spraying method and experimental study based on multidimensional prescription maps. Comput. Electron. Agric. 235, 110379. doi: 10.1016/j.compag.2025.110379

Crossref Full Text | Google Scholar

Hu, C., Ru, Y., Li, X., Fang, S., Zhou, H., Yan, X., et al. (2024). Path tracking control for brake-steering tracked vehicles based on an improved pure pursuit algorithm. Biosyst. Eng. 242, 1–15. doi: 10.1016/j.biosystemseng.2024.04.006

Crossref Full Text | Google Scholar

Jiang, S., Chen, B., Li, W., Yang, S., Zheng, Y., and Liu, X.. (2022a). Stereoscopic plant-protection system integrating uavs and autonomous ground sprayers for orchards. Front. Plant Sci. 13, 1040808. doi: 10.3389/fpls.2022.1040808

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, S., Li, W., and Ma, H. (2025). Design of target-variable spraying system based on favd of fruit tree canopy. Front. Plant Sci. 16, 1582664. doi: 10.3389/fpls.2025.1582664

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, S., Yang, S., Xu, J., Li, W., Zheng, Y., Liu, X., et al. (2022b). Wind field and droplet coverage characteristics of air-assisted sprayer in mango-tree canopies. Pest Manage. Sci. 78, 4892–4904. doi: 10.1002/ps.7110

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L., Zhang, R., Chen, L., Hewitt, A. J., He, X., et al. (2024a). Toward a remote sensing method based on commercial lidar sensors for the measurement of spray drift and potential drift reduction. Sci. Total Environ. 918, 170819. doi: 10.1016/j.scitotenv.2024.170819

PubMed Abstract | Crossref Full Text | Google Scholar

Li, S., Dai, L., Wang, H., Wang, Y., He, Z., and Lin, S.. (2017). Estimating leaf area density of individual trees using the point cloud segmentation of terrestrial lidar data and a voxel-based model. Remote Sens. 9, 1202. doi: 10.3390/rs9111202

Crossref Full Text | Google Scholar

Li, W., Jiang, S., Yang, S., Feng, H., Liu, W., Zheng, Y., et al. (2024b). Leaf-density estimation for fruit-tree canopy based on wind-excited audio. J. Field Robotics 41, 1469–1479. doi: 10.1002/rob.22336

Crossref Full Text | Google Scholar

Li, W., Yang, S., Zhao, H., Jiang, S., Zheng, Y., Liu, X., et al. (2024c). Deep learning method for leaf-density estimation based on wind-excited audio of fruit-tree canopies. Comput. Electron. Agric. 222, 109062. doi: 10.1016/j.compag.2024.109062

Crossref Full Text | Google Scholar

Liu, S., Jin, X., Nie, C., Wang, S., Yu, X., Cheng, M., et al. (2021a). Estimating leaf area index using unmanned aerial vehicle data: shallow vs. Deep machine learning algorithms. Plant Physiol. 187, 1551–1576. doi: 10.1093/plphys/kiab322

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, W., Li, W., Feng, H., Xu, J., Yang, S., Zheng, Y., et al. (2024). Overall integrated navigation based on satellite and lidar in the standardized tall spindle apple orchards. Comput. Electron. Agric. 216, 108489. doi: 10.1016/j.compag.2023.108489

Crossref Full Text | Google Scholar

Liu, X., Zhang, C., Zhang, H., Yang, S., Jiang, S., Zheng, Y., et al. (2021b). Inter-row automatic navigation method by combining least square and svm in forestry. Trans. Chin. Soc Agric. Eng. 37, 157–164. doi: 10.11975/j.issn.1002-6819.2021.09.018

Crossref Full Text | Google Scholar

Luo, S., Wen, S., Zhang, L., Lan, Y., and Chen, X. (2024). Extraction of crop canopy features and decision-making for variable spraying based on unmanned aerial vehicle lidar data. Comput. Electron. Agric. 224, 20. doi: 10.1016/j.compag.2024.109197

Crossref Full Text | Google Scholar

Palleja, T. and Landers, A. J. (2015). Real time canopy density estimation using ultrasonic envelope signals in the orchard and vineyard. Comput. Electron. Agric. 115, 108–117. doi: 10.1016/j.compag.2015.05.014

Crossref Full Text | Google Scholar

Palleja, T. and Landers, A. J. (2017). Real time canopy density validation using ultrasonic envelope signals and point quadrat analysis. Comput. Electron. Agric. 134, 43–50. doi: 10.1016/j.compag.2017.01.012

Crossref Full Text | Google Scholar

Palleja, T., Tresanchez, M., Teixido, M., Sanz, R., Rosell, J. R., Palacin, J., et al. (2010). Sensitivity of tree volume measurement to trajectory errors from a terrestrial lidar scanner. Agric. For. Meteorology 150, 1420–1427. doi: 10.1016/j.agrformet.2010.07.005

Crossref Full Text | Google Scholar

Partel, V., Costa, L., and Ampatzidis, Y. (2021). Smart tree crop sprayer utilizing sensor fusion and artificial intelligence. Comput. Electron. Agric. 191, 106556. doi: 10.1016/j.compag.2021.106556

Crossref Full Text | Google Scholar

Planas, S., Camp, F., Solanelles, F., Sanz, R., and Rosell-Polo, J. R. (2013). Advances in pesticide dose adjustment in tree crops (Wageningen Academic) 13, 533–539. doi: 10.3920/9789086867783_067

Crossref Full Text | Google Scholar

Raj, R., Walker, J. P., Pingale, R., Nandan, R., Naik, B., and Jagarlapudi, A.. (2021). Leaf area index estimation using top-of-canopy airborne rgb images. Int. J. Appl. Earth Observation Geoinformation 96, 102282. doi: 10.1016/j.jag.2020.102282

Crossref Full Text | Google Scholar

Rincón, V. J., Grella, M., Marucco, P., Alcatrão, L. E., Sanchez-Hermosilla, J., and Balsari, P.. (2020). Spray performance assessment of a remote-controlled vehicle prototype for pesticide application in greenhouse tomato crops. Sci. Total Environ. 726, 138509. doi: 10.1016/j.scitotenv.2020.138509

PubMed Abstract | Crossref Full Text | Google Scholar

Román, C., Arnó, J., and Planas, S. (2021). Map-based zonal dosage strategy to control yellow spider mite (eotetranychus carpini) and leafhoppers (empoasca vitis & jacobiasca lybica) in vineyards. Crop Prot. 147, 105690. doi: 10.1016/j.cropro.2021.105690

Crossref Full Text | Google Scholar

Song, L., Huang, J., Liang, X., Yang, S. X., Hu, W., and Tang, D. (2020). An intelligent multi-sensor variable spray system with chaotic optimization and adaptive fuzzy control. Sensors 20, 2954. doi: 10.3390/s20102954

PubMed Abstract | Crossref Full Text | Google Scholar

Taseer, A. and Han, X. (2024). Advancements in variable rate spraying for precise spray requirements in precision agriculture using unmanned aerial spraying systems: a review. Comput. Electron. Agric. 219, 108841. doi: 10.1016/j.compag.2024.108841

Crossref Full Text | Google Scholar

Tayari, E., Jamshid, A. R., and Goodarzi, H. R. (2015). Role of gps and gis in precision agriculture. J. Sci. Res. Dev. 2, 157–162. Available at: Available at: https://jsrad.org/

Google Scholar

Tewari, V. K., Pareek, C. M., Lal, G., Dhruw, L. K., and Singh, N. (2020). Image processing based real-time variable-rate chemical spraying system for disease control in paddy crop. Artif. Intell. Agric. 4, 21–30. doi: 10.1016/j.aiia.2020.01.002

Crossref Full Text | Google Scholar

Tu, Y., Johansen, K., Phinn, S., and Robson, A. (2019). Measuring canopy structure and condition using multi-spectral uas imagery in a horticultural environment. Remote Sens. 11, 269. doi: 10.3390/rs11030269

Crossref Full Text | Google Scholar

Wang, Y. and Fang, H. (2020). Estimation of lai with the lidar technology: a review. Remote Sens. 12, 3457. doi: 10.3390/rs12203457

Crossref Full Text | Google Scholar

Wang, Y., Kootstra, G., Yang, Z., and Khan, H. A. (2024). Uav multispectral remote sensing for agriculture: a comparative study of radiometric correction methods under varying illumination conditions. Biosyst. Eng. 248, 240–254. doi: 10.1016/j.biosystemseng.2024.11.005

Crossref Full Text | Google Scholar

Wei, Z., Xue, X., Salcedo, R., Zhang, Z., Gil, E., Sun, Y., et al. (2022). Key technologies for an orchard variable-rate sprayer: current status and future prospects. Agronomy 13, 59. doi: 10.3390/agronomy13010059

Crossref Full Text | Google Scholar

Whitfield, E., Hilbourne, O., Stamper, S., and Walklate, P (2022). Precision Orchard Dosing project: Development of a high-precision variable rate spray system in tree fruit orchards. Aspects of Applied Biology, 147, International Advances in Pesticide Application. Available at: https://www.researchgate.net/profile/E-Whitfield/publication/366016701_Precision_Orchard_Dosing_project_Development_of_a_high_precision_variable_rate_spray_system_in_tree_fruit_orchards/links/65c3868d1e1ec12eff7bb3c2/Precision-Orchard-Dosing-project-Development-of-a-high-precision-variable-rate-spray-system-in-tree-fruit-orchards.pdf

Google Scholar

He, X. (2020). Research progress and developmental recommendations on precision spraying technology and equipment in China. Smart Agric. 2, 133. doi: 10.12133/j.smartag.2020.2.1.201907-SA002

Crossref Full Text | Google Scholar

Xue, X., Luo, Q., Ji, Y., Ma, Z., Zhu, J., Li, Z., et al. (2023). Design and test of kinect-based variable spraying control system for orchards. Front. Plant Sci. 14, 1297879. doi: 10.3389/fpls.2023.1297879

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, S., Xu, P., Jiang, S., and Zheng, Y. (2022). Downwash characteristics and analysis from a six-rotor unmanned aerial vehicle configured for plant protection. Pest Manage. Sci. 78, 1707–1720. doi: 10.1002/ps.6790

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, F., Bai, J., Jin, Z., Zhang, H., Guo, Z., and Chen, C.. (2022). Research on precise fertilization method of rice tillering stage based on uav hyperspectral remote sensing prescription map. Agronomy 12, 2893. doi: 10.3390/agronomy12112893

Crossref Full Text | Google Scholar

Zhai, C., Zhao, C., Wang, N., Long, J., Wang, X., Weckler, P., et al. (2018). Research progress on precision control methods of air-assisted spraying in orchards. Trans. Chin. Soc. Agric. Eng. 34, 1–15. doi: 10.11975/j.issn.1002-6819.2018.10.001

Crossref Full Text | Google Scholar

Zhang, L., Sun, B., Zhao, D., Shan, C., Wang, G., Song, C., et al. (2024). Prediction of cotton fpar and construction of defoliation spraying prescription map based on multi-source uav images. Comput. Electron. Agric. 220, 108897. doi: 10.1016/j.compag.2024.108897

Crossref Full Text | Google Scholar

Zhang, Z., Wang, X., Lai, Q., and Zhang, Z. (2018). Review of variable-rate sprayer applications based on real-time sensor technologies provisional chapter review of variable-rate sprayer applications based on real-time sensor technologies 13. doi: 10.5772/intechopen.73622

Crossref Full Text | Google Scholar

Zheng, Y. J., Chen, B. T., Lyu, H. T., Kang, F., and Jiang, S. J. (2020). Research progress of orchard plant protection mechanization technology and equipment in China. Trans. Chin. Soc Agric. Eng. 36, 110–124. doi: 10.11975/j.issn.1002-6819.2020.20.014

Crossref Full Text | Google Scholar

Zhou, L., Xue, X., Zhou, L., Zhang, L., Ding, S., Chang, C., et al. (2017). Research situation and progress analysis on orchard variable rate spraying technology. Nongye Gongcheng Xuebao/Transactions Chin. Soc. Agric. Eng. 33, 80–92. doi: 10.11975/j.issn.1002-6819.2017.23.011

Crossref Full Text | Google Scholar

Keywords: GNSS, variable-rate spraying, leaf area density, prescription map, canopy droplet deposition

Citation: Zhao H, Yang S, Li W, Feng H, Jiang S, Liu W, Li J, Zheng Y and Zhang S (2025) Orchard variable-rate spraying method integrating GNSS and wind-excited audio-conducted leaf area density. Front. Plant Sci. 16:1621080. doi: 10.3389/fpls.2025.1621080

Received: 30 April 2025; Accepted: 30 June 2025;
Published: 18 July 2025.

Edited by:

Edney Leandro da Vitória, Federal University of Espirito Santo, Brazil

Reviewed by:

Ping Xu, Hangzhou Dianzi University, China
Daniele Sarri, University of Florence, Italy

Copyright © 2025 Zhao, Yang, Li, Feng, Jiang, Liu, Li, Zheng and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yongjun Zheng, enlqQGNhdS5lZHUuY24=; Songchao Zhang, emhhbmdzb25nY2hhb0BjYWFzLmNu

†These authors share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.