Next Article in Journal
Edaphic Response and Behavior of Agricultural Soils to Mechanical Perturbation in Tillage
Previous Article in Journal
Vegetation Indices Applied to Suborbital Multispectral Images of Healthy Coffee and Coffee Infested with Coffee Leaf Miner
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of an Automated Linear Move Fertigation System for Cotton Using Active Remote Sensing

1
Department of Agricultural Sciences, Clemson University, Clemson, SC 29634, USA
2
Edisto Research and Education Center, Clemson University, Blackville, SC 29217, USA
*
Author to whom correspondence should be addressed.
AgriEngineering 2022, 4(1), 320-334; https://doi.org/10.3390/agriengineering4010022
Submission received: 30 January 2022 / Revised: 11 March 2022 / Accepted: 15 March 2022 / Published: 18 March 2022

Abstract

:
Optimum nitrogen (N) application is essential to the economic and environmental sustainability of cotton production. Variable-rate N fertigation could potentially help farmers optimize N applications, but current overhead irrigation systems normally lack automated site-specific variable-rate fertigation capabilities. The objective of this study was to develop an automated variable-rate N fertigation based on real-time Normalized Difference Vegetation Index (NDVI) measurements from crop sensors integrated with a lateral move irrigation system. For this purpose, NDVI crop sensors and a flow meter integrated with Arduino microcontrollers were constructed on a linear move fertigation system at the Edisto Research and Education Center in Blackville, South Carolina. A computer program was developed to automatically apply site-specific variable N rates based on real-time NDVI sensor data. The system’s ability to use the NDVI data to prescribe N rates, the flow meter to monitor the flow of N, and a rotary encoder to establish the lateral’s position were evaluated. Results from this study showed that the system could accurately use NDVI data to calculate N rates when compared to hand calculated N rates using a two-sample t-test (p > 0.05). Linear regression analysis showed a strong relationship between flow rates measured using the flow meter and hand calculations (R2 = 0.95), as well as the measured distance travelled using the encoder and the actual distance travelled (R2 = 0.99). This study concludes that N management decisions can be automated using NDVI data from on-the-go handheld GreenSeeker crop sensors. The developed system can provide an alternative N application solution for farmers and researchers.

1. Introduction

Upland cotton (Gossypium hirsutum L.) is an important fiber and oilseed crop grown in many areas worldwide [1]. The United States is the third-largest producer of cotton after India and China [2]. In 2020, the United States produced over twenty million bales of cotton representing USD 7 billion in total value [2]. Demand for cotton-based products has grown with the global population and economy. Global competition, costs associated with water, pesticides, and fertilizers used in cotton production, and diminished prices are among the major challenges of cotton producers [2,3,4,5,6].
Nitrogen fertilizer (N) is one of the major inputs of cotton production. Under- or over-application of N can adversely affect the environment, reduce lint yields, and incur monetary costs to cotton producers through wasted fertilizer and/or reduced product quality and quantity [7,8,9,10,11,12,13]. Application of N beyond cotton needs negatively impacts the environment, which contributes to the destruction of the stratospheric ozone layer, stresses ecosystems, acidifies rain, and contaminates drinking water [11,14,15,16]. The risks of environmental N loss are problematic, especially in the U.S. Southeastern Coastal Plain (SCP) region, where sandy soils with low nutrient holding capacity are prevalent [14,17].
Soil texture variability in the SCP complicates the management of N because varying soil textures differ in their ability to hold nutrients [18]. This means a uniform application rate of N may not address the needs of cotton plants in a field with varying soil types. Clemson University (CU) Extension recommends applying 112 kg/ha (100 lb/ac) of N for irrigated cotton production and 78 kg/ha (70 lb/ac) for dryland cotton production over the growing season, with approximately 70% of the N applied during side-dress application in the SCP region [19]. Multiple factors, such as crop growth stage, soil moisture content, soil texture, and yield potential, could cause sufficient N rates to vary above or below the CU Extension recommendations [12,20,21]. Nitrogen application based on the plant’s needs reduces the risks associated with N management of cotton [22].
The Normalized Difference Vegetation Index (NDVI) can be used to quantify the N requirements of cotton [23]. NDVI value is obtained from the measured magnitude of reflectance of visible red (Red) and near infrared (NIR) using Equation (1) [21]:
NDVI = ( NIR Red ) / NIR + Red
Chlorophyll is the primary pigment responsible for a plant’s leaf color. Nitrogen is a significant component of the molecules in chlorophyll; thus, N is directly correlated with leaf chlorophyll content. NDVI has been used to determine the plant relative N content [24]. A plant with higher N content absorbs more red light and reflects more infrared light than a plant with lower N content [12]. Consequently, cotton with adequate N will have higher NDVI values than N-deficient cotton [24]. In addition to the N content of the plant, other factors, such as plant stage of development, canopy structure and water deficiency, can affect the optical properties of plants. The use of reference N-rich strips helps mitigate the effects of these factors other than N in crops and increases the sensitivity of NDVI measurements [25].
Site-specific N recommendations based on cotton response to N are critical for SCP due to the inherent soil variability in the region [12]. Sensor-based N rate algorithms have been developed to increase the precision of N application [12]. The Clemson N application rate model is based on a sensor-based N rate algorithm that was initially developed for wheat (Triticum aestivum L.) and was modified for use in cotton [12,26]. The Clemson sensor-based cotton N rate algorithm is based on a prediction of the potential response, or response index (RI), to additional N. This is done by comparing NDVI data from sections of the field with high N content, called N rich strips, to NDVI data from parts of the field that have not received N fertilizer [12]. The N rich strip NDVI sensor values (NRSSV) and field NDVI sensor values (FSV) are then correlated with soil electrical conductivity to account for differing soil types [12]. Previous research has shown that the Clemson sensor-based cotton N rate algorithm reduced the total N applied by 40% compared to traditional growers’ applications of 101 kg/ha (90 lb/ac) with similar yields [12]. This side-dress N rate estimation method successfully applied liquid N using different application methods in cotton and corn [12,18,27].
The cost associated with N fertilizer application increases the necessity for reducing over- and under-application [28]. These costs could potentially be reduced by applying N using variable rate applications to meet site-specific crop requirements [29]. Commercial variable-rate applicator systems, such as the Raven RCM, John Deere Rate Controller 2000, and Trimble Field IQ, utilize both map-based and on-the-go decision making systems. These systems are available to producers and integrated with band applicators, sprayers, and spreaders; however, on-the-go decision making methods for variable rate N fertigation have not been integrated with lateral move irrigation systems to our knowledge. A map-based variable-rate N fertigation system for a lateral move irrigation system was developed at Clemson University in 2018 [14]. This map-based N fertigation system allowed different N rates to be applied across the field simultaneously by injecting liquid N fertilizer directly into the lateral move irrigation system [14]. A manually operated flow control valve was used to control the flow rate of the main N line, and pulsing solenoid values were used to apply the variable N rates in different sections of the lateral. The N rates for the map-based system were predetermined using the Clemson sensor-based cotton N rate algorithm, which required the user to collect NDVI plant data from the field before application. In addition, the manually controlled flow control valve required the user to be present as the lateral moved through the field.
In this project, the manual aspects of the system developed by Williams (2018) are replaced by (1) mounting NDVI sensors on the applicator to generate on-the-go, site-specific N rate prescriptions, and (2) automating the flow control bypass system to implement on-the-go rate variation in the absence of an operator. Producers could employ an automated variable-rate N application system to meet the unique nutrient needs of their field. Furthermore, an automated system would benefit producers by increasing efficiency and crop yields, reducing labor and input costs, and protecting the environment.
The overall goal of this study was to develop a fertigation system that automated N fertigation based on real-time NDVI measurements using a lateral move irrigation system. The specific objectives to accomplish this goal were to (a) integrate the required individual components including an NDVI crop sensor and flow meter with an Arduino microcontroller, and (b) construct the automated fertigation system and develop the software to automatically apply site-specific variable N rates based on real-time NDVI sensor data.

2. Materials and Methods

2.1. Site Description

A 76-m lateral-move irrigation system (Reinke Manufacturing, Deshler, Nebraska) located at the Clemson University Edisto Research and Education Center in Blackville, South Carolina (33° 21.530′ N, 81° 19.949′ W) was used to develop the variable-rate N fertigation system. According to the Web Soil Survey [30], the experimental field had an area of around 2.0 ha with 2 soil types, including an Orangeburg Loamy Sand (VaA) and Barnwell Loamy Sand (VaB), representing 43% and 57% of the area, respectively. The VaA soil had a 0–2% slope, a typical profile that included a Loamy Sand Ap layer (0–15 cm depth), a Sandy Clay Loam Bt layer (15–135 cm), and a Sandy Clay Loam BCt layer (135–200 cm). The VaA soil was well-drained, with a depth to the water table of more than 200 cm. The VaB soil had a 2–6% slope, a typical profile that included a Loamy Sand Ap layer (0–20 cm), a Sandy Clay Loam Bt layer (20–127 cm), and a Sandy Clay Loam BC layer (127–200 cm). The VaB soil was well-drained, with a depth to the water table of 100–200 cm. The field had a pH of 6.7 measured in 2018. Between May and October 2021, the average precipitation was 107 mm, and the average high and low temperatures were 31.0 °C and 16.7 °C in Blackville, SC [31]. The lateral-move system was used to irrigate 1.59 ha of the experimental field that was split into 110 field plots measuring 7.8 m by 15.2 m. The irrigation system was equipped with low energy precision application (LEPA) drops spaced at 1.8 m along the length of the lateral. A Veris 3100 soil electrical conductivity (EC) meter (Veris Technologies, Inc., Salina, KS, USA) was used to classify the field plots into two soil EC classification groups. Cotton was planted in May of 2021 with a row spacing of 97 cm to test the system functions under field conditions.

2.2. Hardware Integration with Microcontrollers

A microcontroller was integrated with NDVI sensors, a flow meter, a flow control valve, and a rotary encoder. These components were integrated before the full assembly of the fertigation system in the laboratory, and their functions were tested.

2.2.1. Clemson N Application Rate Model and NDVI Sensors

The automated variable-rate N fertigation system automatically determined the N application rates based on the Clemson sensor-based N recommendation algorithm shown in Equations (2)–(5). This algorithm predicts the potential response to additional N with a response index (RI) value, which is the ratio between the NDVI data from N rich strips and NDVI data obtained from the fields that have not received N fertilizer. The N rich strip NDVI sensor values (NRSSV) and field NDVI sensor values (FSV) are then correlated with soil electrical conductivity to account for differing soil types.
RI = NRSSV / FSV
INSEY = FSV / DAE
YP 0 = 431.46 × 2.72 104.98 × INSEY k
NR = YP 0 × RI YP 0 × FNCS / NUE
where RI, response index; NRSSV, nitrogen rich strip sensor (NDVI) value; FSV, field sensor (NDVI) value; INSEY, in-season estimated yield; DAE, number of days after emergence; YP0, yield potential (kg/ha); k, conversion factor for lb/ac to kg/ha = 1.121; NR, N rate (kg/ha); FNCS, fraction of N in cotton seeds after harvest (4/100), entered into this equation as 0.04; NUEF, N use efficiency fraction for cotton (50/100), entered into this equation as 0.50.
For this system, the NDVI sensor data needed to be collected on the go as the irrigation system moved along the field. The GreenSeeker handheld optical active sensor (Trimble Inc., Sunnyvale, CA, USA) was selected to measure NDVI. In order for the system to calculate the N rates using the handheld GreenSeekers, the GreenSeekers were integrated with an Arduino (Mega 2560 Rev 3, Arduino, Scarmangno, TO, Italy). Initially, data from each GreenSeeker were collected from the USB output of the sensor and converted into serial data format using a USB to serial interface module; however, these modules required external computer drivers for data translation. Therefore, a workaround solution was implemented in which the serial data output was extracted directly from the Greenseeker’s UART interface chip on the printed circuit board (PCB). The leads were soldered to the Transistor-Transistor Logic (TTL) signal lines, and the output data from the GreenSeeker were parsed on an Arduino microcontroller (Figure 1).

2.2.2. Flow Meter

The second piece of equipment that needed to be integrated with an Arduino for the system was the Raven flow meter (RFM-5, Raven Industries, Sioux Falls, SD, USA). This flow meter was used to monitor the flow of N coming from a supply tank for the fertigation system. A scaled simulator was developed in a laboratory setting to investigate the output signals from the flow meter (Figure 2). The Raven flow meter communicated flow rate through a variable frequency output signal which was recorded by an Arduino Mega microcontroller and converted into a corresponding flow rate value using Equation (6) which was derived from the flow meter’s documentation:
FR = f p
where FR, flow rate (L/s); f, output frequency of the flow meter (Hz); p, 136.8 (pulses/L).

2.2.3. Flow Control Valve

The simulator shown in Figure 2 was also used to develop the software control of the TeeJet flow control valve (344BPR-2, TeeJet Technologies, Wheaton, IL, USA). The TeeJet flow control valve regulated the flow rate of N being supplied to the fertigation system. A four-channel relay controlled by an Arduino microcontroller supplied 12 VDC power for the flow control valve.

2.2.4. Rotary Encoder

For the system to apply site-specific N rates throughout the field, the applicator system needed to locate its position accurately. Because of inherent inaccuracies associated with off the shelf GPS receivers, a rotary encoder with a non-slip roller wheel (GHW38, CALT, Beijing, China) was installed on one of the wheels of the irrigation system to measure its travel distance (Figure 3D) more accurately. The rotary encoder produced 2400 pulses per revolution. An Arduino microcontroller was used to count the pulses and calculate the distance traveled using a simple conversion factor shown in Equation (7). With the encoder mounted against the lateral’s wheel, one full rotation of the encoder’s nonstick wheel equals 0.41 m of travel.
d = e f
where d, distance traveled (m); e, output of the encoder (pulses); f, calibration factor to convert the output of the encoder to distance in meter = 1.7 × 10−4.

2.3. Development of the Microcontroller Software and Construction of the Fertigation System

The integrated components were used to construct the fertigation system on the lateral, and custom software was developed to control the automated N management decision making and application process.

2.3.1. N Fertigation System Components

The N Fertigation system was constructed first to aid in developing the electronic control components. The N supply system used in the previous iteration of the fertigation system was used for this project [14]. This system consisted of a roller pump driven by an electric motor that pumped N from a 1250 L tank mounted on the side of the lateral into a bladder tank shown in Figure 3F. The bladder tank was used to provide a steady pressurized supply of N to a flow control bypass system that consisted of a TeeJet flow control valve and Raven RFM5 flow meter (Figure 4). A 12V DC power supply was used to power the flow control valve. The flow control bypass system was regulated by a custom controller that read the N mainline’s flow rate and adjusted the flow rate using the flow control valve.
The N mainline supplied N to ten application zones across the length of the lateral using a 19.05 mm (3/4 in) CPVC mainline pipe. A pressure gauge was fitted to the mainline pipe after the bypass system to monitor the mainline pressure. In each zone, N was pumped through a 19.05 mm (3/4 in) check valve (18341, United States Plastic Corp., Lima, OH, USA) to prevent water from entering the application system if the pressure dropped below the irrigation system pressure. N application rates to each zone were controlled by pulsing 2-way, 24 VAC electric solenoid valves (USS-LSV00005, Joyfay International LLC, Cleveland, OH, USA). Each solenoid valve supplied N to four hoses fitted to the irrigation drop lines using a t-fitting (Figure 3A). A CP4916-15 orifice disk (TeeJet Technologies, Springfield, IL, USA) was placed at each t-fitting to regulate flow and pressure.

2.3.2. N Fertigation Electronic Control Components

The applicator consisted of three main electronic controlling components: a mainline controller, field nodes, and an encoder node (Figure 5). The mainline controller, field nodes, and encoder node all communicated wirelessly using XBee radios (XBS2C, Digi, Hopkins, MN, USA) mounted to an Arduino Mega using XBee shields (EK1185, Gikfun, Dongguan, GD, China). The mainline controller communicated with the field and encoder nodes through a custom call and response protocol using the Arduino serial ports to ensure no data were lost.
A 5V DC 4-channel relay module (4450182, Jbtek, Karol Bagh, ND, India) was used to open or close the flow control valve. For post-application analysis, a GPS module (MC035BK, GoGoRC, Kwun Tong, Hong Kong) was added to the controller to locate the system and denote what time the NDVI samples were taken during application. A real-time clock (RTC) module (DS3231, Maxim Integrated, San Jose, CA, USA) was added in case the GPS was not accessible. Lastly, a Micro SD card module (Shenzhen Haihongli Technology CO., LTD., Shenzhen, GD, China) was added to the controller to save the application data for post analysis (Figure 5C).
Field nodes were developed for each of the ten application zones to collect NDVI data and regulate the individual zones’ application rate through pulsing a solenoid valve. The field nodes and GreenSeeker sensors were fixed in the center of each application zone along the length of the lateral using custom-made mounting brackets (Figure 6). The sampling height of the GreenSeekers could be adjusted by sliding the mounting brackets in vertical direction, and a bolt was used to secure the sliding mount in place.
The field nodes were controlled using Arduino Mega microcontrollers. A 5V 1A linear voltage regulator (LM7805CT, onsemi, Phoenix, AZ, USA) was added to charge the GreenSeeker’s battery during application. A 5V DC 2-channel relay module (4332716883, Shenzhen Qunqi E-Commerce Co LTD, Shenzhen, GD, China) was used to control the solenoid valve and trigger the GreenSeeker to take a reading. A solderless RS232 serial adapter (Anmbest Control Technology Co., Ltd., Wuhan, HB, China) was fitted to the GreenSeeker USB power line, serial data wires, and trigger wire so the sensor could easily be secured to the Arduino while also securing a semi-weatherproof connection (Figure 3G). The field node components were placed inside a 3D printed housing and secured in place using silicon and super glue. A lid was constructed for the box that could be screwed shut and sealed with silicon (Figure 3H). A step-down transformer 24V AC power line was run along the length of the lateral to power the solenoid valves.

2.3.3. N Fertigation Control Program

The mainline controller was developed to calculate site-specific N rates and classify them into five flow rate categories (1.5, 4.5, 7.5, 10.5, and 13.5 mL/s) for each zone. The typical maximum N application rate prescribed using the Clemson University cotton N rate algorithm was approximately 15 mL/s, this was divided into five groups of 3 mL/s, and the target application rate was set for the center of each group. This kept minor changes in NDVI values from causing the system to constantly change application rates. Additionally, the controller would send the N rate data to the field nodes, automatically regulate the flow control bypass system, and save the application data to an SD card. These functions were accomplished with an Arduino microcontroller using a custom program developed for this study (Figure 7A).
User parameters were uploaded to the controller using an SD card that was also used to save the application data. The user parameter file contained a soil EC classification plot map generated using the Veris 3100 soil EC meter, the N rich strip values for the soil EC classifications, the lateral ground speed, and the number of days after emergence. This allowed the user to apply different irrigation rates and to apply the side-dress N on a day of their choosing.
The NDVI sensor data were collected by sending a general “read” signal to ensure readings were collected simultaneously. The NDVI data were collected from each node individually by sending a special character set (e.g., S1, S2, S3…S10) that corresponded to each specific node. Because each element of the system works independently, the controller was equipped to receive error codes and determine if a node was “offline”. Field nodes that were unresponsive would timeout after 700 ms, and an offline error code (e.g., 7) would be recorded instead of an NDVI value. If a field node was responsive but unable to collect NDVI data, an error code (e.g., 9) indicating such would be recorded. Lastly, if a field node was able to trigger an NDVI sample but the GreenSeeker displayed a proximity error, the field node could relay a proximity error code (e.g., 8) to the controller.
Once NDVI data were collected from each node, the controller would refer to the plot map of the field with the user parameters on the SD card to establish what soil EC classification each zone fell into and then calculate the N rate for each zone, respectively. The controller would then send a signal to all the nodes containing the mainline required flow rate and the prescribed flow rate for each zone. The controller would then apply N to the field by running a function to regulate the flow control bypass system to read the flow meter output frequency to calculate the mainline flow rate. The controller adjusts the mainline flow rate by relieving the system’s pressure through the flow control bypass system using the flow control valve. After the controller finished its application process, the program would start over; collecting NDVI data, prescribing and categorizing N rates, and applying N until the lateral moved to the end of the field.
The rotary encoder located the lateral as it travelled through the field and used an interrupt service routine to signal lateral movements that would frequently freeze the controller’s main program due to the interrupts coming during functions that could not be interrupted. For this reason, a rotary encoder node was developed using a separate Arduino Mega placed in the controller housing (Figure 3C). Upon request from the controller, the encoder node wirelessly transmitted the locational data to the controller, thus reducing the risk of the controller freezing during testing and application. Before the controller collected data from the field nodes, it would collect the distance traveled data from the encoder node, ensuring the N rate calculations for each node were associated with the soil EC data in each plot area. Once the encoder node received the controller’s signal, it converted the encoder value to distance traveled and reported the value. Otherwise, the encoder node would record the number of pulses sent from the rotary encoder using an interrupt function (Figure 7C).
A custom program was written for each field node to collect and send the field NDVI data upon request from the controller, calculate the percentage of flow needed to take from the mainline, and regulate the zone application rate using the solenoid valve (Figure 7B). The field node stayed in the pulsing solenoid function unless the controller commanded it to do something else. At the start of the program, the field node did not apply anything but was technically in the pulsing function. This was done by initializing the pulsing rate to zero.

2.4. N Fertigation System Testing

The ability of the system to collect GreenSeeker data, make application rate decisions, utilize the flow meter, and locate itself using the encoder was examined through a series of tests. The first evaluated element of the system was the ability to collect and use NDVI data to make management decisions; two methods were used. The first method for testing the system’s capacity to make accurate application rate decisions was by having the field nodes send simulated NDVI data to the controller in the laboratory. This allowed the controller to calculate N rates for a wide range of NDVI values (0.1–0.99), ensuring the system could properly use the Clemson sensor-based cotton N rate algorithm to prescribe N rates. The system calculated application rates for two soil EC classifications using simulated N- rich strip data uploaded using the SD card. The encoder value was set to be midway through the field, allowing the system to calculate N rates based on a soil EC classification plot map. This was done to verify that the system could determine the soil zone and prescribe different rates between zones. Application rates were calculated in MS Excel using the same NDVI data and compared to the controller calculated rates using a two-sample t-test. The fertigation system was also used to collect GreenSeeker data throughout the field to ensure that the system was capable of collecting NDVI data in a field environment. A single soil EC classification was uploaded to the controller for this test.
The Raven flow meter accuracy in the system was verified using the applicator system constructed on the experimental field’s lateral move irrigation system. Water was used to test the system in place of N. The five application rates (1.5, 4.5, 7.5, 10.5, and 13.5 mL/s) were uniformly applied in each zone for six minutes and collected using bottles attached to the irrigation drop hoses. The collected water volumes were totaled and converted into an application rate for the entire system. Application data saved by the system were compared to the measured data using linear regression and absolute mean error in MS Excel.
A series of tests were conducted to investigate its ability to locate itself, apply the correct N rate, and use the flow meter. The system’s ability to use the rotary encoder to estimate its distance travelled through the field was tested by moving the lateral through the field for 62 m while the controller saved the distance travelled values sent by the encoder node to the SD card. The accuracy of this method was determined by comparing the distance the lateral would have travelled based on its known speed compared to the distance travelled values recorded by the controller using linear regression analysis in MS Excel.

3. Results

The first element of the system that was evaluated was the system’s ability to collect and use the sensor data to automate the N management decision-making process and was assessed using two tests. Results for the two-sample t-test show that there was no significant difference (p > 0.05) between the controller and hand calculations for the kg/ha N application rate prescription (p = 0.99), the mL/s N application rate prescription (p = 0.99), or the mL/s N application rate classification (p = 1). The average application rates calculated by the mainline controller and by hand were 179.7494 kg/ha and 179.7491 kg/ha, respectively. The average flow rates calculated by the mainline controller and by hand for 25-S Urea were 13.82717 and 13.8273, respectively. The average classified application flow rate was 8.72 for both the mainline controller’s calculations and the hand calculations.
Figure 8 shows the full range of N application rate prescriptions (kg/ha) produced by the controller for two soil EC classifications. The typical range of NDVI values observed in the field is between 0.30 and 0.80. Figure 9 shows the controller’s N application rate prescription (mL/s) and the controller’s application rate classification (mL/s) for this range.
The second test used GreenSeeker data collected by the lateral over the entire length of the field to analyze the system estimation of N rates in a more practical environment while also testing if the system could collect GreenSeeker data in the field. Figure 10 shows the full range of NDVI values collected using the application system, as well as the prescribed and categorized N application rates for a single soil EC classification.
Another component of the system that was tested was the system’s ability to use the Raven flow meter to measure flow rates. Five application rates were applied uniformly across the lateral, and water volumes were collected for six minutes. The collected volumes for all zones were added together and converted into an actual flow rate to be compared to the data recorded from the controller. Absolute mean error and linear regression analysis were used to analyze the relationship between the measured flow rate of the system and the controller estimation of the flow rate using the flow meter. An absolute mean error of 7.12% was observed between the two sample populations. Figure 11 displays the linear regression analysis, which shows a slope of 1.0056 and an R2 of 0.95.
The third component of the system tested was the system’s ability to use the rotary encoder to estimate the distance travelled while the lateral moved through the field. This was tested by running the system through the field for 62 m and allowing the system to record the encoder values. Linear regression was used to compare the system’s estimation of the distance travelled to the distance the lateral should have moved over time. Figure 12 shows the linear regression analysis results, which show a slope of 0.9956 and an R2 of 0.99.

4. Discussion

The developed fertigation system was fully capable of using Arduino microcontrollers to collect NDVI sensor values from the GreenSeekers and automate the N management decision-making process. Additionally, the fertigation system was capable of using high precision flow metering and regulating equipment to monitor and regulate the N being supplied to the system using open-source components. These are useful tools for farmers and researchers working to develop nutrient management systems at a reduced budget and increase the availability of systems such as these to producers worldwide.
There are changes that could be made to improve the functionality of the system for a farmer’s use, and these could be done while continuing to use open-source, off-the-shelf components. The current code is limited to performing one task at a time. This means that the controller cannot monitor the system’s flow rate or pressure and collect NDVI data or calculate application rates at the same time. The same problem exists in the field nodes. For instance, when the field nodes are applying N, they cannot collect NDVI data simultaneously. A new version of multi-task Arduino code or Raspberry Pi can be used to enable multiple functions to run simultaneously, thus enabling the system to monitor the flow control bypass system while also collecting NDVI data and making decisions. This would also enable the addition of on-the-go error tracking for this system, which would be greatly beneficial for the operator as well. A dashboard or live visual interface could be added to a Raspberry Pi because of its preexisting visual interface capabilities. Displaying the errors live would also allow for an operator to pause an application and fix any problem on the go without having to find issues with an application in post-analysis.
The application system currently uses Application Transparent (AT) mode for the XBee radios, which means that the radios do not build a network. The radio module transmits the received serial data to a remote destination wirelessly. This also means that all the data sent throughout the system need to be managed very carefully to avoid losing data or locking up the system with incomplete messages. If the system was converted into Application Programming Interface (API) mode, communication could be conducted more quickly, using a data exchange protocol. It would also allow the system to build a sensor network, effectively eliminating any issues related to radio range and allowing the implementation of this system on a lateral of any length.
While the system performed well in these tests, the accuracy of the system needs to be further investigated. The system’s ability to vary rates across zones and the accuracy of the GreenSeekers data collected by the application system should be tested. These questions were not addressed in this study because the focus was on the development of the application system.

5. Conclusions

An automated variable-rate N fertigation based on real-time NDVI measurements for a lateral move irrigation system was developed in the Coastal Plain region of South Carolina. The system supplies N independently from the lateral’s flow of water so that a variable rate irrigation system could be used without affecting the flow of N. Although the developed system was integrated with a lateral move fertigation system, the methods, instrumentation, and software demonstrated in this paper can be adopted for automated N fertigation center pivot irrigation systems with some changes in distance calculation with the encoder. Automated, site-specific N fertigation using overhead irrigation systems would allow farmers to simultaneously reduce the environmental and financial impacts of N application without sacrificing yield quality or quantity. Most importantly, the ability to integrate high precision flow control and metering equipment, as well as handheld GreenSeekers with open-source platforms, such as Arduinos, gives farmers and researchers alike the ability to design and customize site-specific nutrient management equipment to meet their individual needs and increases the accessibility of these systems for producers everywhere. Future research must include testing the efficiency of the developed system for cotton in field conditions.

Author Contributions

Conceptualization, S.B., A.B.K., J.M.M., J.P., A.K. and M.M.; methodology, S.B., A.B.K., J.M.M. and A.K.; software, S.B.; validation, S.B., A.B.K. and J.M.M.; formal analysis, S.B.; investigation, S.B.; resources, A.B.K., J.M.M., J.P., A.K. and M.M.; data curation, S.B.; writing—original draft preparation, S.B.; writing—review and editing, S.B., A.B.K., J.M.M., J.P. and M.M.; visualization, S.B.; supervision, A.B.K.; project administration, A.B.K.; funding acquisition, A.B.K., J.M.M., J.P., A.K. and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This material is based upon work supported by NIFA/USDA, under project number SC-1700593. Technical Contribution No. 7046 of the Clemson University Experiment Station.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the general farm staff at the Edisto Research and Education Center for the support and resources provided throughout the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. USDA. USDA Coexistence Fact Sheets Cotton; USDA: Washington, DC, USA, 2015.
  2. USDA-ERS. Cotton Sector at a Glance; USDA-ERS: Washington, DC, USA, 2020.
  3. Almas, L.K.; Colette, W.A.; Warminski, P.L. Reducing irrigation water demand with cotton production in West Texas. In Proceedings of the 2007 Annual Meeting, Mobile, AL, USA, 4–7 February 2007. [Google Scholar]
  4. García-Vila, M.; Fereres, E.; Mateos, L.; Orgaz, F.; Steduto, P. Deficit irrigation optimization of cotton with AquaCrop. Agron. J. 2009, 101, 477–487. [Google Scholar] [CrossRef]
  5. Meisner, M.H.; Rosenheim, J.A.; Agkopoulos, I. A data-driven, machine learning framework for optimal pest management in cotton. Ecosphere 2016, 7, e01263. [Google Scholar] [CrossRef] [Green Version]
  6. Makhdum, A.H.; Khan, H.; Ahmad, S. Reducing cotton footprints through implementation of better management practices in cotton production; a step towards Better Cotton Initiative. In Proceedings of the Fifth Meeting of the Asian Cotton Research and Development Network, Lahore, Pakistan, 1 January 2011. [Google Scholar]
  7. Khan, A.; Tan, D.K.Y.; Munsif, F.; Afridi, M.Z.; Shah, F.; Wei, F.; Fahad, S.; Zhou, R. Nitrogen nutrition in cotton and control strategies for greenhouse gas emissions: A review. Environ. Sci. Pollut. Res. 2017, 24, 23471–23487. [Google Scholar] [CrossRef] [PubMed]
  8. Khan, A.; Najeeb, U.; Wang, L.; Tan, D.K.Y.; Yang, G.; Munsif, F.; Ali, S.; Hafeez, A. Planting density and sowing date strongly influence growth and lint yield of cotton crops. Field Crops Res. 2017, 209, 129–135. [Google Scholar] [CrossRef]
  9. Khan, A.; Munsif, F.; Akhtar, K.; Afridi, M.Z.; Ahmad, Z.; Fahad, S.; Ullah, R.; Khan, F.A.; Din, M. Response of fodder maize to various levels of nitrogen and phosphorus. Am. J. Plant Sci. 2014, 5, 47679. [Google Scholar] [CrossRef] [Green Version]
  10. Kablan, L.A.; Chabot, V.; Mailloux, A.; Bouchard, M.È.; Fontaine, D.; Bruulsema, T. Variability in corn yield response to nitrogen fertilizer in eastern Canada. Agron. J. 2017, 109, 2231–2242. [Google Scholar] [CrossRef] [Green Version]
  11. Hake, K.; Cassman, K.; Ebelhar, W. Cotton nutrition-N, P and K. Cotton Physiol. Today 1991, 2, 2. [Google Scholar]
  12. Porter, W.M. Sensor Based Nitrogen Management for Cotton Production in Coastal Plain Soils; Clemson University: Clemson, SC, USA, 2010. [Google Scholar]
  13. Hearn, A. Effect of preceding crop on the nitrogen requirements of irrigated cotton (Gossypium hirsutum L.) on a vertisol. Field Crops Res. 1986, 13, 159–174. [Google Scholar] [CrossRef]
  14. Williams, P.B.; Khalilian, A.; Marshall, M.W.; Maja, J.M.; Liu, H.; Park, D.; Nafchi, A.M. Development and Testing of a Variable Rate Nitrogen Application System through an Overhead Irrigation System. J. Water Resour. Prot. 2018, 10, 994. [Google Scholar] [CrossRef] [Green Version]
  15. Inman, D.; Khosla, R.; Westfall, D.; Reich, R. Nitrogen uptake across site specific management zones in irrigated corn production systems. Agron. J. 2005, 97, 169–176. [Google Scholar] [CrossRef]
  16. Socolow, R.H. Nitrogen management and the future of food: Lessons from the management of energy and carbon. Proc. Natl. Acad. Sci. USA 1999, 96, 6001–6008. [Google Scholar] [CrossRef] [Green Version]
  17. Wright, D.L.; Wiatrak, P.J.; Rhoads, F.M.; Marois, J.J. Nitrogen requirements for cotton in the Lower Coastal Plain. In Proceedings of the Beltwide Cotton Conference, Quincy, FL, USA, 4–6 January 2022. [Google Scholar]
  18. Williams, P.B.; Khalilian, A.; Marshall, M.W.; Maja, J.M.; Liu, H.; Park, D.; Nafchi, A.M. Cotton response to variable nitrogen rate fertigation through an overhead irrigation system. Agric. Sci. 2019, 10, 66–80. [Google Scholar] [CrossRef] [Green Version]
  19. Jones, M.; Farmaha, B.; Greene, J.; Marshall, M.; Mueller, J.; Smith, N. South Carolina cotton growers’ guide. Clemson Univ. Coop. Ext. Publ. EC 2019, 589, 4. [Google Scholar]
  20. Gerik, T.J.; Oosterhuis, D.M.; Torbert, H.A. Managing cotton nitrogen. Adv. Agron 1998, 64, 115–147. [Google Scholar]
  21. Wiatrak, P.; Khalilian, A.; Wallace, D.; Henderson, W.; Hallmen, R. Incorporating soil electric conductivity and optical sensing technology to develop a site-specific nitrogen application for corn in South Carolina. In Proceedings of the 2008 Southern Conservation Agricultural Systems Conference, Tifton, GA, USA, 29–31 July 2008; pp. 107–112. [Google Scholar]
  22. Zhang, L.-Z.; Spiertz, J.; Zhang, S.-P.; Li, B.; Van der Werf, W. Nitrogen economy in relay intercropping systems of wheat and cotton. Plant Soil 2008, 303, 55–68. [Google Scholar] [CrossRef] [Green Version]
  23. Vories, E.D.; Jones, A.S.; Sudduth, K.A.; Drummond, S.T.; Benson, N.R. Sensing Nitrogen Requirements forIrrigated and Rainfed Cotton. Appl. Eng. Agric. 2014, 30, 707–716. [Google Scholar]
  24. Sui, R.; Wilkerson, J.; Hart, W.; Wilhelm, L.; Howard, D. Multi-spectral sensor for detection of nitrogen status in cotton. Appl. Eng. Agric. 2005, 21, 167–172. [Google Scholar] [CrossRef]
  25. Muñoz-Huerta, R.F.; Guevara-Gonzalez, R.G.; Contreras-Medina, L.M.; Torres-Pacheco, I.; Prado-Olivarez, J.; Ocampo-Velazquez, R.V. A Review of Methods for Sensing the Nitrogen Status in Plants: Advantages, Disadvantages and Recent Advances. Sensors 2013, 13, 10823–10843. [Google Scholar] [CrossRef]
  26. Lukina, E.; Freeman, K.; Wynn, K.; Thomason, W.; Mullen, R.; Stone, M.; Solie, J.; Klatt, A.; Johnson, G.; Elliott, R. Nitrogen fertilization optimization algorithm based on in-season estimates of yield and plant nitrogen uptake. J. Plant Nutr. 2001, 24, 885–898. [Google Scholar] [CrossRef]
  27. Rogers, N.G. Sensor Based Nitrogen Management for Corn Production in Coastal Plain Soils; Clemson University: Clemson, SC, USA, 2016. [Google Scholar]
  28. Shanahan, J.; Kitchen, N.; Raun, W.; Schepers, J.S. Responsive in-season nitrogen management for cereals. Comput. Electron. Agric. 2008, 61, 51–62. [Google Scholar] [CrossRef] [Green Version]
  29. Khalilian, A.; Rogers, N.G.; Williams, P.B.; Han, Y.J.; Nafchi, A.M.; Maja, J.M.; Marshall, M.W.; Payero, J.O. Sensor-Based Algorithm for Mid-Season Nitrogen Application in Corn. Open J. Soil Sci. 2017, 7, 278–287. [Google Scholar] [CrossRef] [Green Version]
  30. USDA-NRCS. Web Soil Survey; USDA-NRCS: Washington, DC, USA, 2022.
  31. U.S.ClimateData. Climate Blackville-South Carolina. Available online: https://www.usclimatedata.com/climate/blackville/south-carolina/united-states/ussc0025 (accessed on 3 September 2021).
Figure 1. A modified GreenSeeker’s internal circuitry from two perspectives. (A) GreenSeeker PCB with leads soldered to the TTL signal lines and ground terminal, (B) The Greenseeker’s TTL TX signal lead soldered on the PCB.
Figure 1. A modified GreenSeeker’s internal circuitry from two perspectives. (A) GreenSeeker PCB with leads soldered to the TTL signal lines and ground terminal, (B) The Greenseeker’s TTL TX signal lead soldered on the PCB.
Agriengineering 04 00022 g001
Figure 2. The scaled applicator system used to integrate the Raven flow meter with an Arduino.
Figure 2. The scaled applicator system used to integrate the Raven flow meter with an Arduino.
Agriengineering 04 00022 g002
Figure 3. The systems’ major components: (A) the N injection point, (B) the mainline controller, (C) the mainline controller interior, (D) the rotary encoder mounted against the lateral’s wheel, (E) the flow control bypass system, (F) the N supply tank, (G) the modified handheld GreenSeeker, and (H) the field node.
Figure 3. The systems’ major components: (A) the N injection point, (B) the mainline controller, (C) the mainline controller interior, (D) the rotary encoder mounted against the lateral’s wheel, (E) the flow control bypass system, (F) the N supply tank, (G) the modified handheld GreenSeeker, and (H) the field node.
Agriengineering 04 00022 g003
Figure 4. The flow control bypass system mounted on the lateral’s mainline.
Figure 4. The flow control bypass system mounted on the lateral’s mainline.
Agriengineering 04 00022 g004
Figure 5. The major components of the electronic control system. (A) flow meter, (B) the flow control valve, (C) an SD card module, (D) a GPS receiver, (E) a four-channel relay module, (F) an XBee radio, (G) an Arduino Mega microcontroller, (H) a two-channel relay module, (I) a Trimble handheld GreenSeeker, and (J) a rotary encoder.
Figure 5. The major components of the electronic control system. (A) flow meter, (B) the flow control valve, (C) an SD card module, (D) a GPS receiver, (E) a four-channel relay module, (F) an XBee radio, (G) an Arduino Mega microcontroller, (H) a two-channel relay module, (I) a Trimble handheld GreenSeeker, and (J) a rotary encoder.
Agriengineering 04 00022 g005
Figure 6. GreenSeekers and nodes mounted on brackets along the length of the lateral move fertigation system.
Figure 6. GreenSeekers and nodes mounted on brackets along the length of the lateral move fertigation system.
Agriengineering 04 00022 g006
Figure 7. The program flow charts for the software written to wirelessly control the applicator system (A) the program flowchart for the mainline controller, (B) the program flow chart for the field nodes, and (C) the program flowchart for the encoder node.
Figure 7. The program flow charts for the software written to wirelessly control the applicator system (A) the program flowchart for the mainline controller, (B) the program flow chart for the field nodes, and (C) the program flowchart for the encoder node.
Agriengineering 04 00022 g007
Figure 8. Controller–calculated N prescriptions (kg/ha) vs. controller-recorded NDVI values. The vertical lines represent the typical range of NDVI values observed in cotton. The average application rate for the full range of NDVI values was 179.75 kg/ha for the soil EC classification 1 and 106.49 for the soil EC classification 2.
Figure 8. Controller–calculated N prescriptions (kg/ha) vs. controller-recorded NDVI values. The vertical lines represent the typical range of NDVI values observed in cotton. The average application rate for the full range of NDVI values was 179.75 kg/ha for the soil EC classification 1 and 106.49 for the soil EC classification 2.
Agriengineering 04 00022 g008
Figure 9. Controller–calculated N rates for soil EC classification 1 (black triangles) and soil EC classification 2 (grey triangles) and the controller-categorized N rates for soil EC classification 1 (black squares) and soil EC classification 2 (grey squares) vs. controller-recorded NDVI values. The average application rates for each rate classification of 0, 1.5, 4.5, 7.5, 10.5, and 13.5 mL/s were −10.107, 1.328, 4.433, 7.472, 10.622, and 29.329 mL/s for the soil EC classification 1. The average application rates for each rate classification of 0, 1.5, 4.5, 7.5, 10.5, and 13.5 mL/s were −14.024, 1.529, 4.644, 7.625, 11.477, and 33.279 mL/s for the soil EC classification 2.
Figure 9. Controller–calculated N rates for soil EC classification 1 (black triangles) and soil EC classification 2 (grey triangles) and the controller-categorized N rates for soil EC classification 1 (black squares) and soil EC classification 2 (grey squares) vs. controller-recorded NDVI values. The average application rates for each rate classification of 0, 1.5, 4.5, 7.5, 10.5, and 13.5 mL/s were −10.107, 1.328, 4.433, 7.472, 10.622, and 29.329 mL/s for the soil EC classification 1. The average application rates for each rate classification of 0, 1.5, 4.5, 7.5, 10.5, and 13.5 mL/s were −14.024, 1.529, 4.644, 7.625, 11.477, and 33.279 mL/s for the soil EC classification 2.
Agriengineering 04 00022 g009
Figure 10. Controller–prescribed N rates (mL/s) and the controller-categorized N rates (mL/s) vs. controller-recorded NDVI values.
Figure 10. Controller–prescribed N rates (mL/s) and the controller-categorized N rates (mL/s) vs. controller-recorded NDVI values.
Agriengineering 04 00022 g010
Figure 11. Controller-estimated flow rate using the flow meter vs. the actual flow rates measured by hand. The black triangles represent the recorded flow meter values while the grey line represents a 1:1 line.
Figure 11. Controller-estimated flow rate using the flow meter vs. the actual flow rates measured by hand. The black triangles represent the recorded flow meter values while the grey line represents a 1:1 line.
Agriengineering 04 00022 g011
Figure 12. Controller–estimated distance travelled using the rotary encoder vs. actual distance travelled. The black triangles represent the encoder values, while the grey line represents a 1:1 line.
Figure 12. Controller–estimated distance travelled using the rotary encoder vs. actual distance travelled. The black triangles represent the encoder values, while the grey line represents a 1:1 line.
Agriengineering 04 00022 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bell, S.; Koc, A.B.; Maja, J.M.; Payero, J.; Khalilian, A.; Marshall, M. Development of an Automated Linear Move Fertigation System for Cotton Using Active Remote Sensing. AgriEngineering 2022, 4, 320-334. https://doi.org/10.3390/agriengineering4010022

AMA Style

Bell S, Koc AB, Maja JM, Payero J, Khalilian A, Marshall M. Development of an Automated Linear Move Fertigation System for Cotton Using Active Remote Sensing. AgriEngineering. 2022; 4(1):320-334. https://doi.org/10.3390/agriengineering4010022

Chicago/Turabian Style

Bell, Stewart, A. Bulent Koc, Joe Mari Maja, Jose Payero, Ahmad Khalilian, and Michael Marshall. 2022. "Development of an Automated Linear Move Fertigation System for Cotton Using Active Remote Sensing" AgriEngineering 4, no. 1: 320-334. https://doi.org/10.3390/agriengineering4010022

APA Style

Bell, S., Koc, A. B., Maja, J. M., Payero, J., Khalilian, A., & Marshall, M. (2022). Development of an Automated Linear Move Fertigation System for Cotton Using Active Remote Sensing. AgriEngineering, 4(1), 320-334. https://doi.org/10.3390/agriengineering4010022

Article Metrics

Back to TopTop