Propagation Analysis for Wireless Sensor Networks Applied to Viticulture

Wireless sensor networks have been proposed as a solution to obtain soil and environment information in large distributed areas. Themain economic activity of the São Francisco Valley region in theNortheast of Brazil is the irrigated fruit production.The region is one of themajor agricultural regions of the country.Grape plantations receive large investments andprovide goodfinancial return. However, the region still lacks electronic sensing systems to extract adequate information from plantations. Considering these facts, this paper presents a study of path loss in grape plantations for a 2.4GHz operating frequency. In order to determine the position of the sensor nodes, the research dealt with various environmental factors that influence the intensity of the received signal. It has been noticed that main plantation aisles favor the guided propagation, and the vegetation along the secondary plantation aisles compromises the propagation. Diffraction over the grape trees is the main propagation mechanism in the diagonal propagation path. Transmission carried out above the vineyard showed that reflection on the top of the trees is the main mechanism.


Introduction
The term Precision Agriculture (PA) refers to a set of technologies used to understand the changes that occur during the plantation cycles [1].These tools should allow better monitoring of the stages of agricultural production.In this context, Wireless Sensors Networks (WSNs) have been proposed as a solution to perform this monitoring [2,3].The emergence of protocols and low-cost devices enable the deployment of the technology in those environments, mainly by allowing the supervision of large areas.
This research has been developed in the São Francisco Valley region, Northeast of Brazil.The main economic activity is irrigated fruit growing for export, and the production of the region reaches the foreign market, especially Europe, where sanitary and quality standards are strict [4].As mentioned, WSNs are candidates to improve the production process, as they provide new ways to manage the harvesting areas, along with the fact that they do not need transmission wires.Therefore, the cost of deployment of WSNs may be cheaper, and more adequate, than wired networks.
One of the possible applications of this research is to collect soil moisture data.With the humidity data it is possible to develop an automatic irrigation system taking into account areas of soil that are most in need of water, avoiding the uniform irrigation of the plantation.The proposed idea can be seen in Figure 1, which represents nodes arranged in a moisture map.
However, before the implementation of a communication system, it is important to know the path loss due to different antennas settings, related to the presence of obstacles, for example.Factors such as the transmission power, vegetation density, and the occurrence of line-of-sight also influence the project [5,6].Due to the lack of specific models, or parameters for the existing models in this type of application, this article presents an experimental study of the path loss found in grape plantations based on a in loco measurement campaign.Classical literature models were chosen and the parameters values were obtained empirically.From these values, discussions regarding the results have been performed.The objective of the research was to analyze the communication channel, predict the path loss, and evaluate the effects of the propagation phenomena on the signal intensity.
Based on the experiments performed in the same farm during several weeks, it has been noticed that the main plantation aisles favor the guided propagation.On the other hand, there was a considerable amount of vegetation along the secondary plantation aisles which compromised the propagation.Measurements taken following a diagonal path, in relation to the main aisle, have shown that diffraction by the grape trees is the main propagation mechanism.Transmission carried out above the vineyard, with both transceivers above it, has shown that reflection on the top of the trees is the main mechanism.
This paper is organized as follows.Section 2 presents related works in the area of WSNs in agriculture.Section 3 presents theoretical concepts regarding the propagation mechanisms.Section 4 shows the environment under study.The measurements are discussed in Section 5. Conclusions and future works are presented in Section 6.

Wireless Sensor Networks and Agriculture
The literature presents some works that deal with WSN applied to the agriculture, which motivated this research.The system proposed in [3] was developed for large-scale agricultural environments, using static sensor nodes with two goals: to collect information from newborn animals and to collect measurements of soil moisture.The nodes automatically read moisture data from the soil, typically in intervals of one minute, for each node.Then, the data is aggregated at the base to generate an updated profile for the entire pasture.
In [7], software and hardware for deployment of a WSN in Washington State have been developed.There are two networks, one to measure regional temperature and one to monitor plantation freezing.In the regional network, a master node is configured with multiple repeaters to build a backbone.The proposed network monitors the plantation and uses star topology.The biggest problems were encountered when the energy management and the switches were damaged during storms.
The Lofar Agro Project [8] has been developed to monitor microclimates in plantations.The purpose of monitoring is to combat the Phytophthora infestans bacteria that attack potato.The authors comment that humidity is an important factor in the development of diseases.A sensor network monitors, in addition to the humidity, the temperature and whether or not the leaves are wet, which are the main factors related to the Phytophthora infestation.The communications network has a weather station that measures light, air pressure, precipitation, wind direction, and strength.
Studies were performed in a cornfield [9] to evaluate the propagation characteristics with the IRIS Crossbow Inc. nodes.The influence of antenna height, the different directions within the planting area, the weather conditions (sunny day and rainy day), and transmitted power were considered.All factors were evaluated for the worst case: plantation with maximum height, maximum width of the stems, and long, thick leaves.It was possible to identify the conditions for increasing the coverage area by the sensor nodes.
A WSN was implemented to estimate the amount of water in the crop from the observed signal attenuation [10].The nodes are deployed with the objective of collecting the signal.From the mean value of the power measured by the sensor nodes, the amount of water is estimated from the propagation model proposed by the authors.The advantage of this method is that the sensors are not required to obtain data.Some studies using ZigBee are presented in [11][12][13].In [11] an irrigation control system to optimize water usage is presented.In [12], a framework is proposed to monitor humidity, temperature, and lighting to control water consumption in agricultural areas.A WSN was developed in [13] using two modifications to the grid topology to improve network performance.

Theoretical Background
The WSNs suffer from various types of transmitted signal degradation.In wireless systems, the signal is attenuated in the path between the transmitter and receiver, according to the characteristics of the environment in which the electromagnetic wave propagates.These features include terrain deformations, obstacles between the antennas, and vegetation density.In addition, several propagation mechanisms influence the signal propagation, such as reflection, scattering, refraction, and diffraction [14].

Classical Propagation Models.
The path loss is caused by the dissipation of power radiated by the transmitter, as well as by the effects of the propagation channel [15].It can be defined as the power loss between the transmitter and receiver in a radio communication link.
There are two ways to model the path loss [16]: a priori modeling and measurement-based modeling.A priori models are appropriate for predicting path loss in places where it is difficult to perform measurements.The path loss is predicted using only a priori knowledge.The measurementbased models are related to methods to collect data and predict values of signal strength at locations where data cannot be collected.

Free Space Propagation.
In free space, it is possible to devise a model with a flexible path loss exponent [16] for ideal isotropic antennas.The power (in dBm) at the receptor is given by in which   is the transmitted power,   and   are the transmitter and receiver gains, respectively,  is the wavelength of the transmitted signal,  is the distance between transmitter and receiver, and  is the system loss.The parameter  is adjusted using the least squares or least absolute residual regression models, for example, [17].

Ground Reflection Model.
The ground reflection model, or the plane earth model [18,19], is an a priori model.It considers a direct, and a reflected, ray with no obstacles between the transmitter and the receiver.The path loss is given by [19] in which   and   are the transmitter and receiver gains, respectively, ℎ  and ℎ  are the heights of the antennas of the transmitter and the receiver, respectively, and  is the distance between transmitter and receiver.

Log-Distance Model.
The log-distance model attenuation parameter  is obtained empirically.It is used for environments with shadowing, and the received power is given by [20] in which ( 0 ) is the average power measured at a reference distance, usually taken at 1 meter,  is the distance loss exponent, and  is a Gaussian random variable with mean  and variance  2 .

Propagation between Two Parallel Planes.
The parallel plane guide is formed by two metal planes.There is a dielectric material in the middle with permittivity .The electromagnetic wave propagates inside the structure in the  direction (see Figure 2).Depending on the permeability and permittivity it is possible to have the value of the received power, inside parallel planes, very close to the received power in free space for the same distance; that is, the power along the parallel planes is a function of  2 .

Empirical Foliage Models.
There are some models based on measurements to predict losses in different types of vegetation and at different frequencies.The exponential decay model proposed by Weissberger [21] is used when the path between the transmitter and the receiver is full of dense vegetation, composed of trees with low humidity.The loss due to vegetation is given by [21] in which  is the frequency, in GHz, and  is the height of the trees, in meters.An ITU-R recommendation was developed from measurements performed primarily in the UHF band.The distance between the transmitter and the receiver must be shorter than 400 m, to allow the signal to propagate without the formation of a lateral wave.The excess attenuation is given by in which  is the frequency and  is the distance.
According to these models, it is possible to verify that the excess loss due to vegetation can be represented by the expression [22]  (dB) =     . ( From measurements taken in the environment, the parameters , , and  can be obtained using optimization techniques to reduce the least square error.

Link Budget.
The link budget [14] is used to obtain the signal strength at the receiver, considering all path losses between the transmitter and receiver.The received power, accounting for all gains and losses, is given by in which   is the received power in dBm,   is the power of the transmitter,   is the transmitter gain,  path is the total path loss, and   is the gain of the receptor.

Environments
The plantation studied covers an area of approximately 10 hectares and its average height is 2 meters.For this structure, tree trunks and wires are used to support the leaves and grapes in order to uniformize the height of the crop.This structure is called a vineyard, as can be seen in Figure 3. Four height and position terminal settings have been identified within the plantation, which result in different dominant propagation mechanisms.The first case is when the terminals are below the vegetation top and positioned in the type A aisle.The second case occurs when the terminals are below the vegetation top, positioned in the type B aisle.When the terminals are below the top, and there are obstacles between the transmitter and the receiver, the third case is characterized.Finally, the fourth case occurs when both terminals are above the crop.Figure 6 presents the superior view in which type A aisle is also used as a reference axis.
The type A aisle is characterized, mainly, by the following features: Type A aisle is presented in Figure 3.The distance between the vines is around 3 m, and there is a larger amount of vegetation, since the vegetation along type B aisle grows to a height of 90 cm, approximately (Figure 4).There are also ground elevations along the rows.Furthermore, there is a large quantity of grapes throughout the path, as shown in Figure 5.There are areas where the vegetation is about 5 cm height.
The diagonal environment is similar to the environment described for the type B aisle, except for the tree trunks positioned between the two devices.There is a lot of vegetation and there are elevations along the rows.There are grapes between nodes and there exist areas, where the vegetation is about 5 cm high.The vegetation that covers most of the area is supported by wire, by trunks of vines and wooden stakes.

Methodology.
The experiments described in this section were performed in the Minuano farm, located 70 km from Petrolina city, in the state of Pernambuco, Brazil.The climate is classified as tropical semiarid, with low rainfall rates throughout the year, favoring the production of grapes.Besides the scarcity and irregularity of rainfall, the rains occur mainly in the summer, and there is strong evaporation, as a result of the high temperatures.The average temperature during the experiments was 36 ∘ C.   A pair of XBee ZB Pro S2B nodes was used.The devices have +18 dBm of output power, −102 dBm sensitivity specified by the manufacturer, and the operating frequency is 2.4 GHz.The antenna is one dipole type.Both devices presented the same specifications, except that one was connected to a notebook, so that the data could be stored in a database.This terminal was configured as a coordinator.
One of the terminals (endpoint) was kept at the same position, and the other one (coordinator node) was strategically positioned along the plantation.At each point, the coordinator node requires one Received Signal Strength Indicator (RSSI) package of the endpoint.The measurement procedure was as follows: (1) The coordinator (transmitter) is connected to a notebook and positioned at a distance of 1 meter from endpoint.
(2) The coordinator sends one package to the endpoint, requiring it to calculate the RSSI.
(3) The endpoint transmits the package with the computed value.
(5) This procedure is repeated at different distances, until the power is low enough, or the plantation fence is reached.
(6) Data are analyzed using Gnuplot and MATLAB.
(7) One propagation model is chosen and adjusted for each case.

Measurements.
The measurements taken along type A aisle were performed with the nodes positioned in the middle of the aisle.The receiver was installed at a height of 1 m.First, the transmitter was installed at a height of 20 cm above the ground.One hundred signal intensity measurements were made at intervals of 50 meters.The experiment was repeated with the transmitter positioned at heights of 40 cm, 60 cm, and 1 m.These measurements were taken during the "in leave" plantation stage.
Other measurements were performed for different stages of the plantation, at a height of 1 m.The transmitter was placed along type A aisle when the plantation did not have any leaves, only the wires that supported them covered the top of the plantation.In a second round of measurements, data were collected during the "in fruit" plantation stage.
At a distance of 400 m, the end of the plantation was reached.However, comparing the power value at that distance with the sensitivity specified by the manufacturer, it is noticed that the range is larger than 400 m.Any aisles long enough to perform the measurements beyond 400 m have not been found in this plantation.
Along type B aisle, measurements were performed with the nodes positioned in the middle of the aisle.The enddevice was kept at a fixed position at a height of 1 m.The transmitter was installed at different heights, such as in type A aisle, 20 cm, 40 cm, 60 cm, and 1 m.It has been found in preliminary experiments that the range, in this case, is about 390 m.One hundred signal strength measurements were made at intervals of 25 m instead of 50 m to better observe the signal decay behavior.
As in type A aisle, measurements were limited, due to the size of the plantation, to a distance of 200 m.Using a direct comparison between the power at the distance of 200 m and the device sensitivity (−102 dBm), one realizes that the range goes beyond 200 m.However, no type B aisle long enough was found to perform measurements beyond 200 m.
For the diagonal case, the measurements were performed following some angles, as shown in Figure 6.The enddevice was kept in a fixed position with a height of 1 m.The transmitter is installed, first, at a height of 20 cm above the ground.One hundred signal strength measurements were made following the dashed lines forming 15, 30, 45, and 75 degrees (lines 1, 2, 3, and 4 of Figure 6, resp.) with type A aisle.The experiment was repeated with the transmitter at heights of 40 cm, 60 cm, and 1 m.
In some positions, when the signal strength was around −97 dBm, the communication was not established.It was necessary to make several attempts until the data could be collected.Furthermore, in this case, the range is smaller than type A and type B aisles.This is an expected result because there are obstacles between the transmitter and the receiver.
Measurements performed with the coordinator above the plantation were performed, as illustrated in Figure 7.The end-device was installed at a height of 3 m, positioned in the middle of the aisle.Then, the coordinator was moved along type A aisle at a height of 3 m.

Least Absolute Residual Method.
The least absolute residual (LAR) method was used to adjust the curves.It minimizes the sum of the deviation modules.The advantage of this method is that it is more robust against outliers [17].
Considering   the power and   the distance, for the unidimensional case, the method consists of determining  0 parameter in so that the absolute sum of the errors is minimized, The value of  0 that minimizes the sum is the sample median, which explains the decrease of the outliers contribution.For linear regression, 0 and  1 parameters are calculated using linear programming to minimize [17] min

Interpretation of Measurements
In order to measure the goodness of the model adjustments, the Root Mean Square Error (RMSE) has been used [17]: in which  is the number of points in the sample, and   and x are the measured and predicted values, respectively.The determination coefficient  2 is also used as a metric in the study of regression and is given by [23] The total variance is the sum of the squares of the differences between the measured value and the average of the measured values.The explained variance is the sum of squares of the differences between each predicted value and the average of the measured values.This measure indicates how good is the regression, that is, how much the regression line approaches the collected data.The coefficient  2 is in the range [0, 1].The closer to one, the better the regression line describes the behavior of the data.

Type A Aisle.
Figure 8 shows the measurements of type A aisle, for Friis free space and plane earth models without correction factors.The vertical axis represents the average received power in dBm and the horizontal axis represents the distance between antennas.
The models were chosen for the first analysis because there are no obstacles between the terminals.It is observed, however, that the signal strength decay does not follow the plane earth model.One can verify this statement by averaging the differences between the value predicted by the plane earth model and the measured value.Thus, the average gain is equal to 7.87 dB.Table 1 presents the RMSE values.
From the curves and RMSE values, it seems that the geometry of the plantation favors the guided propagation of the wave.This fact is evidenced mainly by the RMSE value for free space that is lower than the plane earth model.Graphically, it is also possible to notice this behavior.However, a dependency on the height is identified, caused by the reflected ray at  = 100 m (limit of the first Fresnel zone).In addition, two adjustments of the theoretical models were  made using the LAR method in MATLAB.The exponents of distance for both free space and the plane earth model were adjusted to minimize the RMSE.Changing the plane earth model, since the decay is not proportional to  −4 , in which  is the adjusted distance exponent.The values obtained for the exponents are presented in Table 2, with the RMSE values for each case.The values of   ,   , and   are provided by the XBee manufacturer and are equal to 18 dBm, 2.1 dB, and 2.1 dB, respectively.It can be seen from Table 2 that the Friis model, with the adjusted exponent, has an exponent close of two, with little improvement of the RMSE.Regarding the plane earth model, the RMSE decreases, indicating a better fit, but the decay is not proportional to the fourth power of the distance.Figure 9 presents the curves with the adjusted exponents.
Figure 2 shows the results for measurements done for different stages of the plantation at the height of 1 m.Both "in fruit" and "in leave" cases are characterized by the waveguide effect, although, mainly because of spreading caused by the  water in the fruits, the attenuation is greater for the "in fruit" case.When the wave guide is unset, that is, when there are no leaves on top of the plantation, the signal power decays much faster and the plane earth model does not fit well.As the measurements do not have a logarithmic behavior, the linear regression was used instead; the following expression has been obtained: This result can be explained by spreading and absorption of the electromagnetic field caused by the wires.Table 3 presents the exponents of the models.It is important to highlight the difference of 0.36 between "in fruit" and "in leave" adjustments.The corresponding plot of the received power (dBm), as a function of the distance, is presented in Figure 10.

Type B Aisle.
As type B aisle is characterized by a higher density of vegetation, scattering and absorption cause losses that cancel the canalization effect.In Figure 11, it can be  seen that the decay is faster for short distances and, as the transmitter moves away from the receiver, the power decays slowly, indicating a logarithmic behavior.Therefore, the log-distance model was used with some modifications.The attenuation is given by in which  is the distance,  0 is the reference distance taken at 1 meter, and  is obtained from the LAR method;  is a Gaussian random variable with zero mean.One modification is the substitution of ( 0 ) by   +   +   ,  =   +   +   −   =  −   (17) in which  is given by As this work does not address the small-scale fading, the  term does not appear in the adjusted equations.Therefore, the received power is given by The value of  is obtained from the LAR method implementation provided by MATLAB and is equal to 4.537.The resulting RMSE is equal to 3.567, and  2 equals 0.9399.The curve obtained is shown in Figure 11, with the distance in a linear scale.
Another option to adjust the curves is to use an equation that includes the height, so that the model provides more parameters and lower RMSE.Start from the following expression: in which ( 0 ) is the power at the reference distance (1 meter), (ℎ  ) is expressed as a function of the receiver distance,  is a constant, and  is the distance.The values of ( 0 ) + (ℎ  ) and  are obtained using the LAR method.
Considering ℎ  = 1 m the reference height, it is possible to obtain the following formula, for ℎ  = 1: The adjusting equation chosen for (ℎ  ) was in which  and  are obtained using the LAR method and are equal to −12.07 and −2.41, respectively.It was possible to find a complete equation, given by  in which ( 0 ) is the average received power at the reference distance.The graph in linear scale for the distance is shown in Figure 12.Table 4 presents the RMSE and  2 values.

Diagonal.
First, a comparison between the Weissberger model and the recommendation of the ITU-R for the UHF band regarding all performed measurements along the diagonal has been performed.The equations for this case, including the excess loss by vegetation, are in which   is the attenuation in free space,  is the frequency, and  is the distance.The plot of the measurements and the model are presented in Figure 13.Although the measurements made by Weissberger have been taken in different environments, the proposed equation presents a reasonable fit in this case.However, in order to reduce the RMSE value, the parameters , , and  of the excess attenuation were adjusted.In addition, the exponent of log-distance model and a linear fit have been determined.The curves for the adjusted models are presented in Figure 14.Table 5 presents the exponents, the RMSE, and  2 values obtained.

Terminals above the Vineyard.
When the two terminals are above the plantation, the electromagnetic field is reflected on the top of the plantation structure formed by the metallic wires and leaves.Thus, the plane earth model is used, considering the heights of the nodes equal to 1 m, instead of 3 m (3 m for the total height from the ground, minus 2 m height of the piles).In this situation, until a distance of 100 meters is reached, the propagation occurs in free space, since the clearance of the route is larger than 0.6, in which  is the radius of the first Fresnel zone.For measurements taken up to a 100 meters, the gap relative to the top structure is approximately 0.6; that is, the plane earth and free space model leads to the same attenuation value.For this distance, the main propagation mechanism is the plane earth with the reflected ray on the top of the plantation.The RMSE for the plane earth model, with the adjusted exponent, equals 5.64.The value of the exponent is equal to 4.017.The graph of Figure 15 shows the measurements collected with the terminals above the plantation following type A aisle and Friis and plane earth models.

Conclusions and Future Works
This paper presents a study of path loss in grape plantations for a 2.4 GHz operating frequency.The main contributions of this research are the identification and interpretation of the propagation phenomena in the studied environment, as well as the methodology adopted to obtain the data parameters.The researched environment has characteristics that favor certain propagation mechanisms depending on the situation.
The region was divided into four areas characterized by different dominant propagation mechanisms.Type A aisle

2 InternationalFigure 1 :
Figure 1: Representation of sensors arranged in a humidity map that presents data by color.

Figure 2 :
Figure 2: Waveguide diagram composed of two parallel planes.
(i) The aisle is formed by the vines and stems which are evenly spaced on both sides.(ii) The aisle is almost entirely covered by the leaves of the vines.(iii) Vegetation along the aisle has a maximum height of 40 cm.(iv) There are areas where the ground vegetation has 5 cm height.

Figure 3 :
Figure 3: Photo of the vineyard used for the experiments.

Figure 4 :
Figure 4: Photos of type B aisle vegetation, in the Vineyard used for the tests.

Figure 5 :Figure 6 :
Figure 5: Photo of grapes along type B aisle.

Figure 7 :
Figure 7: Diagram which represents the position of sensor nodes above the plantation.

Figure 8 :
Figure 8: Received power (dBm) as a function of distance, and Friis free space and plane earth models for type A aisle (linear scale for distance).

Figure 9 :
Figure 9: Received power (dBm) as a function of distance for Friis free space and plane earth models with adjusted exponents for type A aisle (linear scale for distance).

= 4 )Figure 10 :
Figure 10: Received power (dBm) as a function of distance for different stages (in leave, in fruit, and no fruit) of the plantation (linear scale for distance).

Figure 11 :
Figure 11: Received power (dBm), as a function of distance and adjusted logarithmic equation (linear scale for distance).

Figure 12 :
Figure 12: Received power (dBm), as a function of distance and the adjusted linear equation (linear scale for distance).

Figure 13 :
Figure 13: Received power (dBm), as a function of distance, Weissberger and ITU-R equations (linear scale for distance).

Figure 14 :
Figure 14: Received power (dBm), as a function of distance, modified Weissberger model, and linear and log-distance adjustments (linear scale for distance).

Figure 15 :
Figure 15: Received power (dBm) as a function of distance, Friis free space, and plane earth models for the case in which the terminals are above the plantation (linear scale for distance).

Table 1 :
RMSE for Friis free space and plane earth models for type A aisle.

Table 2 :
Friis free space and plane earth models RMSE with adjusted exponents for type A aisle.

Table 3 :
RMSE of the models used for different plantation states along type A aisle.

Table 4 :
RMSE and  2 of the adjusted linear equation.

Table 5 :
Parameters, RMSE values and  2 for Weissberger, and logarithmic and linear equations.