# A Model with Leaf Area Index and Trunk Diameter for LoRaWAN Radio Propagation in Eastern China Mixed Forest

Authors: Yin Wu; Genwei Guo; Guiyun Tian; Wenbo Liu
Journal: Journal of Sensors
 (2020)

---

## Abstract

Internet of Things (IoT) is a very promising technology in forest engineering, especially for the environment and plant growth monitoring. LoRa Wide Area Network (LoRaWAN) is a prevailing choice for the Forestry IoT owing to its low-power and long-range ability. Real-world deployment and network optimization require accurate path-loss modeling, so the LoRaWAN radio channel in the forest is needed to be intensively studied. However, most of the subsistent propagation models do not involve specific forestry environmental parameters. In this paper, two parameters related with the trees are considered: the leaf area index and the tree trunk diameter. Due to the time-changing characteristics of these two items (from spring to winter), an empirical model has been developed through extensive measurement campaigns: Firstly, the channel measurement platform is designed based on a real scene of mixed forest. Secondly, the fading characteristics of the channel transmission for LoRa nodes are tested, and the corresponding model is presented and evaluated. Lastly, an energy harvesting LoRaWAN is deployed and operated in a sampled forest region of Eastern China for environment monitoring based on our propagation model. The results show that 433 MHz LoRa path loss in the mingled forest could be precisely predicted by our proposed model. Moreover, network coverage and energy consumption optimization of the LoRa nodes could be performed, which enables the perpetual development of reliable forestry evolution monitoring system.


---

## Body

## 1. Introduction

The multispecies forest is an important silvicultural regime for sustainable forest management. Due to its vertical stratification structure, canopy plasticity, and niche separation, the mixed forest could utilize canopy space more efficiently to improve forest productivity and nutrient cycle more effectively than the natural forest. Along with its increasing area, mixed forest management faces different challenges [1]: competitions for natural resources of different tree species, different tending measures of growth stages, and improving pest control technologies. So it is urgent to develop and exploit reasonable resource monitoring systems for the prediction and management of mixed forest [2, 3]. Low-power wide-area networks (LPWANs) are a promising technology because of their high applicability, low energy consumption, and self-organizing mechanism. An overview and comparison of different emerging LPWAN technologies have been presented in [4–6] and listed in Table 1. In particular, LoRaWAN is one of the most prospective LPWAN technologies, gaining greater interest from the research and industrial communities. It has also been increasingly adopted in the Internet of Things (IoT) in the forest [7–10].Table 1
Comparison of LPWAN technologies.
LoRaNB-IoTSigfoxZigBeeWiFiFrequency868 MHz (EU); 915 MHz (USA); 433 MHz (Asia)Depends on the frequency licensed to LTE868 MHz (EU); 915 MHz (USA); 433 MHz (Asia)868 MHz (EU); 915 MHz (USA); 433 MHz (Asia); 2.4GHz2.4 GHz and 5 GHzTransmission distanceLong (1-2 km in urban, 20 km in rural)Long (10 km in rural)Long (40 km in rural)Short (10-100 m)Short (50 -100 m)Data transfer rate0.3-50 kbps160 kbps-250 kbps100 bps250 kbps2.4 GbpsTheoretical battery life10 years/AA battery10 years/AA battery10 years/AA battery1 year/AA battery1 year/AA batteryNetwork delay1 sec-3 sec6 sec-10 sec20 sec-25 secLess than 1 secLess than 1 secNetworking modeGateway-basedTelecommunication operator-basedGateway-basedGateway-basedGateway-basedLoRaWAN links every node via wireless communication, just as Figure1 shows. When they are deployed in the forest, the radio signal may experience attenuation as a result of diffraction, reflection, scattering, and absorption effects caused by various obstacles of trees and shrub. Mixed forests are challenging radio propagation environments, and there is a need to understand better the propagation characteristics in order to design efficient wireless systems that can operate inside them. However, the modeling of forest radio propagation is relatively complicated as different environmental parameters (such as type of vegetation, density, height of trees, and weather conditions) can have strong effects.Figure 1
The diagram of LoRaWAN in a mixed forest.Researchers have developed several well-known empirical models to study the radio channel characteristics in the forest: Meng et al. performed experimental path loss modeling for a near ground radio wave propagation in a plantation at 240 and 700 MHz over a large forest depth [11], the proposed LITU-R mode mainly considers the lateral wave effect and shows better performance for the prediction of foliage loss over a large foliage depth (up to 5 km). Azevedo et al. found that the signal attenuation is dependent on the vegetation density, defined by the product of the tree density and the average diameter of the trunks [12]. Anderson et al. presented the results of a UWB measurement campaign and comprehensive parametric analysis for four different forest densities across multiple diverse locations in Maryland of Virginia. They demonstrated that UWB propagation in a forest environment heavily depends on the forest density, antenna type, and forest configuration [13]. Palaios et al. studied the near-ground propagation characteristics in a typical central European forest in UHF-bands around 485 MHz [14]; the result showed that the forest type has a significant effect to the propagation and thus mixed-type of forest requires its own model.On the other hand, the characterization of the LoRa radio channel also got few achievements recently: Benaissa et al. studied LoRa path loss and temporal fading (868 MHz) for dairy cows in barns [15]. Hosseinzadeh et al. assessed the indoor propagation performance of LoRa technology and identified the model that best describes the process [16]. El Chall et al. investigated the LoRaWAN radio channel in both indoor and outdoor environments at urban and rural locations in Lebanon [17]. Sanchez-Iborra et al. presented a comprehensive performance evaluation of LoRaWAN under different environmental conditions [18]. A point-to-point LoRa path loss characterization has been researched in the urban, forest, and coastal environment [19]. Especially Sardar et al. analyzed some environmental factors in the forestry LoRa communication such as different antenna heights and packet lengths [20].Based on the above literature, we find that most of the existing empirical models concentrated on single species of trees such as cedars and pines. Small Acacia confusa or cryptomeria, glass, and shrub wood were rarely taken into account, much less the mixed forest. Secondly, most models are based on the radio frequency and transmission distance, while the environmental parameters were scarcely considered. The prior models mainly divided the trees into in-leaf and out-of-leaf situations. However, the leaf area of trees changes with different seasons, and the cross-sectional area of tree trunk would also influence the radio propagation in its growth stage. Thirdly, the randomly distributed trees in a mixed forest are spatially different from the artificially planted trees that are arranged in rows and columns, whereas sub-GHz LoRa for this kind of forest has not yet been investigated up to now. Finally, since on-field measuring devices are always energy-constrained, it is likely that novel energy-harvesting (EH) IoT technologies will find their use in the forest monitoring [21–23]. So to combine the energy-efficient radio transmission power control mechanism with the accurate propagation model would be essential and inevitable for the EH-wireless sensor system in the forest.In this work, practical LoRaWAN sensor nodes are used to characterize the path loss effect at 433 MHz band in the mingle forest environment at different heights from the ground. The measurement result which compared with the Okumura-Hata-Rural path loss model [24] and the LITU-R model [11] shows a significant difference. Hence, the Leaf Area Index (LAI) that can reflect the density of leaves at different growth stages is introduced as one characteristic parameter, and the trunk Diameter at Breast Height (DBH) that represents cross-section area of trees as the other. An empirical model comprising these two new parameters is then developed, and we compared it with the existing models to demonstrate its accuracy among the LoRaWAN nodes deployed in a mixed forest district in Nanjing city of Eastern China. The measured data, as well as the developed path loss models, can be used for efficient planning and deployment of EH-WSN in mixed forest environments. To the best of our knowledge, this is the first study on the analysis of propagation characteristics of LoRa in a mixed forest.The rest of this paper is organized as follows: Section2 describes the measurement campaign and presents the practical sensor nodes used and the investigated scenarios. Data analysis and model construction are discussed in Section 3. Section 4 tests the proposed model in a real field and compares it with existing models. Conclusions and future studies are given in Section 5.


## 2. Experimental Setup and Measurement Campaign

### 2.1. Mixed Forest Environment Description

We hold the modeling procedure at one typical mixed forest in Nanjing Forestry University, Jiangsu Province, Eastern China (118°48′55.98″E, 32°04′59.01″N). This demonstration forest base was built in 1999 and has a length of 300 m and a width of 200 m. Pinus massoniana, Ilex chinensis, Quercus variabilis, Lindera glauca, Photinia beauverdiana, and Rubus corchorifolius were the dominant variety with an average crown of approximately 2.75 m. The plants were approximately 1.3-11.9 m high with trunk DBH of approximately 0.11-0.56 m.The modeling measurement was mainly carried out at four growth stages including spring germination stage, summer leaf expanding stage, autumn leaf falling stage, and winter senescence stage in 2018. A variation of the same tree at the four growth stages is shown in Figure2: In early March, the tree came up and had sparse leaves (Figure 2(a)); next, the leaves grew rapidly in late June (Figure 2(b)); and the LAI and DBH also reached the peak value; then the tree was in leaf falling stage by September (Figure 2(c)); finally, the leaves faded away and only branches existed in Winter (Figure 2(d)). These different morphological characteristics could definitely affect the wireless signal propagation that passed through and nearby.Figure 2
The variation of the same tree at four different growth stages in the sample forest.
(a)
Spring(b)
Summer(c)
Autumn(d)
Winter


### 2.2. Measurement Devices

The experimental platform used in our measurement campaigns is depicted in Figure3. An Arduino UNO with Dragino LoRa shield expansion board was used as the sensor node, powered by a 3.7 V, 3000 mAH rechargeable lithium battery [25]. It has an integrated LoRa SX1278 transceiver and runs LoRaWAN protocol. One node without associated sensors that connected to a laptop was used as the sink, and it was able to receive LoRa frames from −20 to −141 dBm depending on the LoRa BW and SF [26]. Both sink and sensor nodes used an omnidirectional dipole antenna of 3 dBi gain. Upon the reception of each frame, the sink provided the received signal strength indicator (RSSI), the signal-to-noise ratio (SNR), and the payload message. These received parameters were recorded on the server side on the laptop for further analysis and processing. The quality of reception can be monitored in real-time with an MQTT Web client application.Figure 3
The DBH and PAR sensor node amounted in trees.
(a)
DBH LoRa node(b)
PAR LoRa nodeThe LoRa nodes were configured to send a packet every interval. The packets also included a sequential number in order to identify the packet loss. The transmit power was set to 14 dBm, SF to 12 (to achieve the best receiver sensitivity), and the BW to 125 kHz, using the three default channels (433.175, 433.375, and 433.575 MHz). The packet payload was 37 bytes including 13 bytes MAC header. The measurement setup parameters are listed in Table2. During all the measurements, the position of the sink was fixed, whereas the sensor nodes could be moved to different locations.Table 2
Measurement setup parameters.
ParameterValueSensor nodeSX1278SinkLG01-P, Dragino IoT kitTX power14 dBmFrequency433 MHzBandwidth (BW)125 KHzAntenna3 dBi omnidirectionalSpreading factor (SF)12Coding rate4/5Payload length37 bytesTime interval10 secSensor node antenna height0.5 m, 1.3 m, 2.0 m, 2.5 m, 3.0 mSink antenna height1.3 mAs can be seen from Figure3, there were two particular designed forestry sensors: DBH and Photosynthetically Active Radiation (PAR). The DBH sensor mainly consisted of a high accuracy sliding film potentiometer, a circular flexible elastic wheel rim, a fixed position deck, and a signal processing bridge circuit. After tightened installation onto the tree trunk, the potentiometer’s slider could point out the exact value of its resistance, with which we can deduce the DBH. The PAR sensor was designed for the measurement of LAI, which was based on the reference system [27, 28], it included two parts: the testing components and the reference component. The testing PAR sensors should be placed around the crown to measure the luminous flux which passed through, where the reference sensor could be set up in a clear field to provide contrast and evaluation. Both kinds of sensors make up an effective and indispensable supply to the proposed radio propagation analysis of IoT in the forest.


### 2.3. Practical Measurement Setup

Throughout the experiments, the intention was to place the sink node at the center of designated deployment field and collect RSS readings at eleven different distances (i.e., 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100 meters) and along twelve different 30 degrees separated radials. Therefore, an area of200mx200m is ideally needed to carry out the real field deployment experiment.The antenna of the sink node was placed 1.3 m above the ground, i.e., height of DBH. RSSIs were then measured by the transmitting node with 5 heights from bottom to top at heights of 0.5, 1.3, 2, 2.5, and 3 m. The transmitting and receiving antenna were maintained erect in all of the tests. The top and side views of the measurement site are shown in Figure4. Note that as a result of extensive randomly scattered leaves, branches and bushes, the signal strength may change considerably even when there is a little variation in the position of transmitting antenna. So for each measurement point, 100 samples were collected, providing a sufficiently large dataset for approximating important statistical properties of the signal, such as mean and variability. Hence, the RF measurement of every 660 measuring point is an average of each 100 RSSI samples.Figure 4
Measurement site: the dashed lines represent the testing points, while the center dot is the sink.
(a)
Vertical top view(b)
Side viewBesides, DBH and LAI results have also been collected and stored in the laptop which is connected with the sink. These two kinds of dedicated sensor nodes sent their data once a day to the sink for subsequent data analysis. They were all mounted on some specific trees selected by forest cultivation criterion [29]. We have held up the whole measurement campaign in four growth stages of the mixed forest once each; the average LAI and DBH values of the selected trees were used as the exact novel parameters in the next path loss model analysis.


## 2.1. Mixed Forest Environment Description

We hold the modeling procedure at one typical mixed forest in Nanjing Forestry University, Jiangsu Province, Eastern China (118°48′55.98″E, 32°04′59.01″N). This demonstration forest base was built in 1999 and has a length of 300 m and a width of 200 m. Pinus massoniana, Ilex chinensis, Quercus variabilis, Lindera glauca, Photinia beauverdiana, and Rubus corchorifolius were the dominant variety with an average crown of approximately 2.75 m. The plants were approximately 1.3-11.9 m high with trunk DBH of approximately 0.11-0.56 m.The modeling measurement was mainly carried out at four growth stages including spring germination stage, summer leaf expanding stage, autumn leaf falling stage, and winter senescence stage in 2018. A variation of the same tree at the four growth stages is shown in Figure2: In early March, the tree came up and had sparse leaves (Figure 2(a)); next, the leaves grew rapidly in late June (Figure 2(b)); and the LAI and DBH also reached the peak value; then the tree was in leaf falling stage by September (Figure 2(c)); finally, the leaves faded away and only branches existed in Winter (Figure 2(d)). These different morphological characteristics could definitely affect the wireless signal propagation that passed through and nearby.Figure 2
The variation of the same tree at four different growth stages in the sample forest.
(a)
Spring(b)
Summer(c)
Autumn(d)
Winter


## 2.2. Measurement Devices

The experimental platform used in our measurement campaigns is depicted in Figure3. An Arduino UNO with Dragino LoRa shield expansion board was used as the sensor node, powered by a 3.7 V, 3000 mAH rechargeable lithium battery [25]. It has an integrated LoRa SX1278 transceiver and runs LoRaWAN protocol. One node without associated sensors that connected to a laptop was used as the sink, and it was able to receive LoRa frames from −20 to −141 dBm depending on the LoRa BW and SF [26]. Both sink and sensor nodes used an omnidirectional dipole antenna of 3 dBi gain. Upon the reception of each frame, the sink provided the received signal strength indicator (RSSI), the signal-to-noise ratio (SNR), and the payload message. These received parameters were recorded on the server side on the laptop for further analysis and processing. The quality of reception can be monitored in real-time with an MQTT Web client application.Figure 3
The DBH and PAR sensor node amounted in trees.
(a)
DBH LoRa node(b)
PAR LoRa nodeThe LoRa nodes were configured to send a packet every interval. The packets also included a sequential number in order to identify the packet loss. The transmit power was set to 14 dBm, SF to 12 (to achieve the best receiver sensitivity), and the BW to 125 kHz, using the three default channels (433.175, 433.375, and 433.575 MHz). The packet payload was 37 bytes including 13 bytes MAC header. The measurement setup parameters are listed in Table2. During all the measurements, the position of the sink was fixed, whereas the sensor nodes could be moved to different locations.Table 2
Measurement setup parameters.
ParameterValueSensor nodeSX1278SinkLG01-P, Dragino IoT kitTX power14 dBmFrequency433 MHzBandwidth (BW)125 KHzAntenna3 dBi omnidirectionalSpreading factor (SF)12Coding rate4/5Payload length37 bytesTime interval10 secSensor node antenna height0.5 m, 1.3 m, 2.0 m, 2.5 m, 3.0 mSink antenna height1.3 mAs can be seen from Figure3, there were two particular designed forestry sensors: DBH and Photosynthetically Active Radiation (PAR). The DBH sensor mainly consisted of a high accuracy sliding film potentiometer, a circular flexible elastic wheel rim, a fixed position deck, and a signal processing bridge circuit. After tightened installation onto the tree trunk, the potentiometer’s slider could point out the exact value of its resistance, with which we can deduce the DBH. The PAR sensor was designed for the measurement of LAI, which was based on the reference system [27, 28], it included two parts: the testing components and the reference component. The testing PAR sensors should be placed around the crown to measure the luminous flux which passed through, where the reference sensor could be set up in a clear field to provide contrast and evaluation. Both kinds of sensors make up an effective and indispensable supply to the proposed radio propagation analysis of IoT in the forest.


## 2.3. Practical Measurement Setup

Throughout the experiments, the intention was to place the sink node at the center of designated deployment field and collect RSS readings at eleven different distances (i.e., 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100 meters) and along twelve different 30 degrees separated radials. Therefore, an area of200mx200m is ideally needed to carry out the real field deployment experiment.The antenna of the sink node was placed 1.3 m above the ground, i.e., height of DBH. RSSIs were then measured by the transmitting node with 5 heights from bottom to top at heights of 0.5, 1.3, 2, 2.5, and 3 m. The transmitting and receiving antenna were maintained erect in all of the tests. The top and side views of the measurement site are shown in Figure4. Note that as a result of extensive randomly scattered leaves, branches and bushes, the signal strength may change considerably even when there is a little variation in the position of transmitting antenna. So for each measurement point, 100 samples were collected, providing a sufficiently large dataset for approximating important statistical properties of the signal, such as mean and variability. Hence, the RF measurement of every 660 measuring point is an average of each 100 RSSI samples.Figure 4
Measurement site: the dashed lines represent the testing points, while the center dot is the sink.
(a)
Vertical top view(b)
Side viewBesides, DBH and LAI results have also been collected and stored in the laptop which is connected with the sink. These two kinds of dedicated sensor nodes sent their data once a day to the sink for subsequent data analysis. They were all mounted on some specific trees selected by forest cultivation criterion [29]. We have held up the whole measurement campaign in four growth stages of the mixed forest once each; the average LAI and DBH values of the selected trees were used as the exact novel parameters in the next path loss model analysis.


## 3. Measurement Results and Empirical Model

In this section, we present the results of measurements conducted in the sample Mixed Forest of Nanjing Forestry University as described in Section2. Then, we derive Pass Loss (PL) models from four different environments (spring, summer, august, and winter). At last, the proposed PL model is analyzed and discussed.The RSSI values collected in the mixed forest environment were all converted into PL values which represent the signal attenuation on the path between the transmitter and the receiver antennas. As for the PL, value is a unique property of the transmission environment. Generally, the relationship between PL, transmission power, and RSSI is given as:(1)PLdB=PTdBm−PRdBm+GTdB+GRdB,where PL refers to the path loss, PT refers to the transmitting power, PR is the receiving power (PR=RSSI+SNR), GT is the transmit antenna gain, and GR is the receive antenna gain. Table 3 just presents the PL results in the spring stage measurement when transmitting node’s height is 0.5 m.Table 3
Path loss measurement result (dB) in the spring stage with transmitting height 0.5 m.
Radial no.Degree5 m10 m20 m30 m40 m50 m60 m70 m80 m90 m100 m1054.560.178.681.984.290.589.8101.3105.4110.7115.023052.264.382.185.281.793.297.397.1110.5112.2118.336055.166.975.379.786.994.198.895.6112.4123.8109.949056.968.273.783.888.297.796.9103.4107.6121.1117.1512049.869.579.989.778.486.985.4106.8103.2126.7114.6615051.663.383.586.582.588.188.7109.799.3117.5107.3718056.959.285.483.489.689.499.4103.598.8111.2109.7821059.365.789.190.392.392.598.5101.9104.7107.5119.2924053.667.880.684.985.185.891.097.6108.1106.4110.81027048.862.677.281.186.894.394.793.4102.6109.3133.51130055.270.379.888.388.993.693.9108.2104.9104.6128.41233058.464.878.482.684.687.795.5110.1103.8109.9116.7Then, we calculate the average PL values related to each testing point including the very corresponding 5 heights. Figure5 illustrates the average PL values versus the log of distance. Moreover, it illustrates the corresponding linear regression model lines.Figure 5
Average path loss and linear regression lines for the four-stage environments.
(a)
Spring stage(b)
Summer stage(c)
August stage(d)
Winter stageIt can be clearly pointed out that the variations among PL values which at the same distances are relatively large. In the mixed forest environment, such variations mainly occur as a result of the absence of Light of Sight communication between the transmitter and receiver, resulting in great signal fluctuation. In some instances, multiple waves of the transmitted signal may arrive at the sink, causing multipath fading. In other instances, only a small scattered portion of the transmitted signal may be received. Meanwhile, the season changes of mixed forest have a great influence on the PL values, the growth, and decay of branches and leaves directly transform the transmission routes of LoRaWAN.Also, we can see that the average PL values closely follow the linear regression line, which indicates a dependency between the PL values and the log of distance. As a matter of fact, the linear regression provided in Figure5 produces a reliable coefficient of determination (i.e., 0.8913<R2<0.985), so the log-normal model is suitable for estimating the PL:
(2)PLLNdB=PLd0+10nlog10dd0+Xσ,where n is the path loss exponent that indicates the rate at which the signal attenuates with the distance, PLd0 is the path loss at a known reference distance d0 in the far-field, Xσ denotes a zero-mean Gaussian random variable with standard deviation σ, which reflects the variation of the received power around the average.Figure6 shows the variations of attenuation index n and PLd0 according to different heights. For the heights above 1.3 m, n increases along with the height, especially when in the summer time, it is very reasonable due to the lush branches and leaves. The communication path always passes through the crown when the height is over 2 m. However, the attenuation of 0.5 m is more or less larger than the 1.3 m except for the winter stage, we mainly attribute the phenomena to the bushes and thick undergrowth of the mixed forest. PLd0 shows the opposite trends as the testing nodes’ height increase, the values decline mostly because of the obstructions between transmitter, and receiver increase simultaneously. Furthermore, according to the different seasons, these two curves show a similar variation tendency, which we should thoroughly analyze. In other words, we may try to combine specific forest parameters together.Figure 6
Average attenuation indexn and PLd0 of the four-stage environments.
(a)(b)Here, we record the average DBH and LAI values with respect to the distance between transmitting nodes and sink. These two particular sensor nodes were located half-randomly based on the forestry cultivation rules [29]. We can see that the average DBH values almost keep the same from spring to summer, then it has a rapid growth until winter in Figure 7(a). However, Figure 7(b) shows that the LAI values mainly change along with the seasons, while not with the distance. The two items could perfectly represent the characteristic of the sample mixed forest.Figure 7
Average DBH and LAI of the four-stage environments.
(a)(b)The last parameter of the log-normal model, represented byXσ, can be defined by the standard deviation. Figure 5 reveals that a good approximation of the values can be obtained by using the mean. Here, we considered that Xσ=6.94 dB.So to precisely evaluate the parametersn and PLd0, we introduce joint polynomials to combine the DBH and LAI together. The formulae that are shown below elaborate the three constants of the linear fitting (xn1,xn2,xn3) with distance (dM, from sender to receiver) and height (hG, nodes’ vertical position above the ground). These parameters were determined by curve fitting of the above Figures 6 and 7:
(3)n=xn1⋅Φ+xn2,(4)Φ=xn3⋅10⋅log10VDBH−8.63+VLAI,(5)xn3=0.036⋅dM−240.61,(6)xn2=−0.5169⋅hG2+2.9773⋅hG+1.485,(7)xn1=−0.0931⋅hG2+0.8575⋅hG−0.04392.Similarly, we can have:(8)PLd0=xPL1⋅Ψ+xPL2,(9)Ψ=xPL3⋅log10VDBH−VLAI,(10)xPL3=0.045⋅dM+0.32,(11)xPL2=41.8,(12)xPL1=−5.1856⋅hG.VDBH and VLAI are just the values of DBH and LAI. Thus, a new empirical model could be created by using equations ((3)–(12)) to estimate the LoRaWAN path loss in a typical mixed forest.


## 4. Evaluation and Discussion

### 4.1. Model Validation

We carried out the experimental evaluation at a different mixed forest in Eastern China’s Nanjing city in 2019, it is located at the west side of Zi Jin Shan Mountain and has a length of 170 m and a width of 120 m. The tree species structure is similar with our modeling forest, except for the heavy brush approximate 0.6 m high. Also, the DBH and LAI values are higher due to the well-protected natural environment.The testing method was similar to that in the modeling experiments. Here, the validation measurement was executed at 5 different heights (0.75 m, 1.3 m, 1.75 m, 2.25 m, and 2.75 m) along one straight line between the transmitter and receiver. On the same time, two forestry parameters (LAI and DBH) at each validation point could be collected conveniently. So we implemented the whole validation from March to July, and the testing positions were chosen with an extensive representativeness for different growth stages as much as possible.The comparison of measurement results and the proposed new model is shown below in Figure8. 95% of the R2 values are larger than 0.85, the least value being 0.81. Most of the Root Mean Square Error (RMSE) values were between 2 and 4, as shown in Table 4.Figure 8
Comparisons between the measurements and the new model at 8 selected points.
(a)
LAI=0.67, DBH=7.8(b)
LAI=1.44, DBH=8.5(c)
LAI=2.02, DBH=9.1(d)
LAI=2.85, DBH=12.4(e)
LAI=3.56, DBH=15.6(f)
LAI=4.38, DBH=12.8(g)
LAI=5.16, DBH=11.5(h)
LAI=5.23, DBH=13.3Table 4
Validation results at different heights for 8 testing points.
Antenna height (m)(a)(b)(c)(d)(e)(f)(g)(h)0.750.930.910.880.950.930.900.940.962.871.893.843.152.692.412.953.371.30.960.960.840.930.870.860.910.953.912.162.962.372.193.583.401.771.750.940.930.950.960.900.920.880.991.723.742.062.893.064.073.762.552.250.920.900.970.940.950.970.850.972.493.813.413.863.342.612.792.722.750.880.890.850.870.910.930.820.843.142.762.822.632.993.674.112.09∗The first line denotes the R2 values, and the second line indicates the RMSE values at each height.We can see from the results that the proposed model matches the measurements quite well. In general, the presence of high density of trees can lead to a bad quality reception: the higher value of LAI and DBH always means that there are more obstructions of radio scattering, reflections, diffractions, and absorptions. Hence, this nonline of sight PL values increase sharply along with the transmission distance and antenna height.Then, the proposed new model has also been compared with some existing models at testing point (d), as shown in Figure9: The results for the empirical model Okumura-Hata-Rural are represented by a dash-dot line, it can be observed that Okumura–Hata model underestimates the PL values and presents a slower rising tendency; By contrast, the result of LITU-R model in dotted line shows smaller error since it was developed to follow data for higher dissipation. However, they all show relatively large deviations compared with our proposed models. So to evaluate the different path loss models thoroughly, we use three parameters: the RMSE, the Mean Absolute Error (MAE), and the Mean Absolute Percentage Error (MAPE). From Table 5, we can observe that the best prediction is achieved by using the proposed path loss model.Figure 9
Comparisons between the new model and two existing models (LAI=2.85, DBH=12.4).Table 5
The RMSE, MAE, and MAPE values between the proposed model and existing models.
LocationOkumura-Hata-ruralLITU-RProposed modelRMSEMAEMAPE (%)RMSEMAEMAPE (%)RMSEMAEMAPE (%)Height 0.75 m11.6610.1710.708.847.038.113.152.792.96Height 1.30 m9.788.569.446.935.495.872.371.962.17Height 1.75 m15.4113.7914.1213.3511.2412.632.892.022.79The primary cause is that none of the existing models could estimate the path loss at any height without any LAI or DBH. Furthermore, the value of RMSE, MAE, and MAPE in Okumura-Hata-Rural and LITU-R highly vary from one height to another, therefore, a slight variation of antenna or forestry parameters would negatively affect the prediction of the path loss in these models. Our proposed new model outperforms the existing models with the lowest RMSE, MAE, and MAPE regardless of the height of the antenna. Thus, the predicted values of our model are closer to real observations than the existing models.


### 4.2. Field Measurement on Energy-Harvesting LoRa Nodes

In this section, we evaluate the performance of the proposed radio propagation model based on the wind energy harvesting LoRaWAN which has been mounted in the sample mixed forest at Nanjing Forestry University [22]. We take the wind-powered sensor nodes serving as transmitters (Figure 10 shows the node’s mainboard, while the sensors could be attached on when in practice). An unlimited power supply sink is responsible for receiving messages and returning the packet delivery ratio (PDR) every day. The transmitting rate is 250 bps, and the payload length is 50 bytes. The selected nodes are located at 3 different positions: (i) LAIi=1.99, DBHi=13.82; (ii) LAIii=3.01, DBHii=10.78; (iii) LAIiii=2.63, DBHiii=15.84. The corresponding communication distances and antenna heights are (i) di=37m, hi=1.0m; (ii) dii=55m, hii=1.3m; (iii) diii=64m, hiii=2.0m.Figure 10
Wind energy harvesting LoRa node (mainboard).We validate the effectiveness of our propagation loss calculation module by using two comparisons. The existing fixed transmission power scheme is implemented as a reference, which is referred to as the Fixed-TX scheme. One of the wind-harvesting LoRaWAN nodes works under our proposed model at position (iii), while the other nodes located at positions (i) and (ii) be set with the fixed transmission powerTX=7dBm and TX=14dBm, respectively. The node (iii) could compute its transmitting power by equation (1) based on the RSSI value efficiently. All three nodes start with an initial energy of 1 Joule, and the average wind speed during the experiment is 4.4 m/s.Table6 presents the average PDR of total of 10,080 packets (ten packets per minute during 7 days) and the residual energy of three nodes in the final time. We can see that the proposed method achieves the lowest energy consumption compared with the other transmission power strategies no matter the distance. Compared with the Fixed-TX-14 dBm scheme, our proposed module achieves energy reduction around 30%. Compared with the Fixed-TX-7 dBm scheme, our proposed module improves the PDR around 5%. The wireless path loss and the background noise are always time-varying, the fixed transmission power definitely cannot guarantee the optimal transmission energy efficiency. The proposed PL model is effective in dynamically calculating the optimal transmission power along with the lowest energy cost.Table 6
The residual energy and PDR of three nodes.
Average testing resultProposed PL modelFixed-TX-14 dBmFixed-TX-7 dBmSupercapacitor voltage (volt, residual energy)3.1942.3722.468PDR (%)1009795


### 4.3. Rain Attenuation of LoRaWAN in the Forest

Rain and fog is quite common in the forestry district. Many studies indicate that the propagation of the radio signal is affected by rain and there are some remarkable models deduced to analyze its influence [30, 31]. However, most existing researches focus on the rain attenuation by wireless signals with frequencies above 10 GHz. So in this part, we investigate the influence of rain attenuation within the proposed forestry LoRaWAN for the first time. The experiment setup is just the same as part 4.2 and we record the variations of RSSI in position (ii) every 10 seconds when rainfall comes. The measurement data presented in Figure 11 are filtered by a moving average window with 30 samples. It is noticeable that during rainfall periods there is an obvious attenuation of PL.Figure 11
Variation of RSSI vs. rainfall in position (ii) from 2019.05.01 to 2019.05.08.
(a)(b)Figure11(a) shows a significant impact on path loss of LoRa signal when it is rainy during a 7-day consecutive observation, every increasement of PL values corresponds to a rainfall event accurately. Although the rain rate which has been recorded by a standard rain gauge [32] does not exceed the extremely heavy level, a nearly 8 dB variation of PL has occurred. Figure 11(b) shows a detailed changing process of rain attenuation, PL value increases sharply when the rain begins to fall, and it returns to the original state relatively slow after the rain stops. The main reason of fading may be the water molecular which is adhered on the antenna of LoRa nodes and pervaded in the transmission space, it would change the propagation characteristic and the antenna property drastically.According to reference [33], the rain attenuation can be represented by aRPb, where RP is the rain rate with a given period P while a and b are depending on the radio wave frequency and the rain temperature. Hence, we calculate this item with the collected rainy data and deduce the result as RP is set to 3.86 mm/hr based on the conventional rain information from local weather station, a and b just are 2.11 and 0.55, respectively. So the updated general PL model would be:
(13)PLLNdB=PLd0+10nlog10dd0+Xσ+aRpb.Furthermore, we also apply the updated model on the wind-harvesting LoRaWAN node to verify its effectiveness. It is located in position (iii) as well and the average wind speed is 4.9 m/s during the rainy period. The other two nodes located at positions (i) and (ii) are set the same as part 4.2, and the corresponding results are shown in Table7.Table 7
The residual energy and PDR of three nodes during rainy time.
Average testing resultProposed PL modelFixed-TX-14 dBmFixed-TX-7 dBmSupercapacitor voltage (volt, residual energy)3.0312.4472.584PDR (%)1009390It can be seen that the residual energy of the updated model decreases slightly otherwise the two Fixed-TX nodes get higher due to the larger wind speed. However, PDR of the updated scheme keeps 100%, while the comparative two nodes drop significantly. This indicates that the updated PL model matches the real rain attenuation efficiently and works stably in energy neutral mechanism.


## 4.1. Model Validation

We carried out the experimental evaluation at a different mixed forest in Eastern China’s Nanjing city in 2019, it is located at the west side of Zi Jin Shan Mountain and has a length of 170 m and a width of 120 m. The tree species structure is similar with our modeling forest, except for the heavy brush approximate 0.6 m high. Also, the DBH and LAI values are higher due to the well-protected natural environment.The testing method was similar to that in the modeling experiments. Here, the validation measurement was executed at 5 different heights (0.75 m, 1.3 m, 1.75 m, 2.25 m, and 2.75 m) along one straight line between the transmitter and receiver. On the same time, two forestry parameters (LAI and DBH) at each validation point could be collected conveniently. So we implemented the whole validation from March to July, and the testing positions were chosen with an extensive representativeness for different growth stages as much as possible.The comparison of measurement results and the proposed new model is shown below in Figure8. 95% of the R2 values are larger than 0.85, the least value being 0.81. Most of the Root Mean Square Error (RMSE) values were between 2 and 4, as shown in Table 4.Figure 8
Comparisons between the measurements and the new model at 8 selected points.
(a)
LAI=0.67, DBH=7.8(b)
LAI=1.44, DBH=8.5(c)
LAI=2.02, DBH=9.1(d)
LAI=2.85, DBH=12.4(e)
LAI=3.56, DBH=15.6(f)
LAI=4.38, DBH=12.8(g)
LAI=5.16, DBH=11.5(h)
LAI=5.23, DBH=13.3Table 4
Validation results at different heights for 8 testing points.
Antenna height (m)(a)(b)(c)(d)(e)(f)(g)(h)0.750.930.910.880.950.930.900.940.962.871.893.843.152.692.412.953.371.30.960.960.840.930.870.860.910.953.912.162.962.372.193.583.401.771.750.940.930.950.960.900.920.880.991.723.742.062.893.064.073.762.552.250.920.900.970.940.950.970.850.972.493.813.413.863.342.612.792.722.750.880.890.850.870.910.930.820.843.142.762.822.632.993.674.112.09∗The first line denotes the R2 values, and the second line indicates the RMSE values at each height.We can see from the results that the proposed model matches the measurements quite well. In general, the presence of high density of trees can lead to a bad quality reception: the higher value of LAI and DBH always means that there are more obstructions of radio scattering, reflections, diffractions, and absorptions. Hence, this nonline of sight PL values increase sharply along with the transmission distance and antenna height.Then, the proposed new model has also been compared with some existing models at testing point (d), as shown in Figure9: The results for the empirical model Okumura-Hata-Rural are represented by a dash-dot line, it can be observed that Okumura–Hata model underestimates the PL values and presents a slower rising tendency; By contrast, the result of LITU-R model in dotted line shows smaller error since it was developed to follow data for higher dissipation. However, they all show relatively large deviations compared with our proposed models. So to evaluate the different path loss models thoroughly, we use three parameters: the RMSE, the Mean Absolute Error (MAE), and the Mean Absolute Percentage Error (MAPE). From Table 5, we can observe that the best prediction is achieved by using the proposed path loss model.Figure 9
Comparisons between the new model and two existing models (LAI=2.85, DBH=12.4).Table 5
The RMSE, MAE, and MAPE values between the proposed model and existing models.
LocationOkumura-Hata-ruralLITU-RProposed modelRMSEMAEMAPE (%)RMSEMAEMAPE (%)RMSEMAEMAPE (%)Height 0.75 m11.6610.1710.708.847.038.113.152.792.96Height 1.30 m9.788.569.446.935.495.872.371.962.17Height 1.75 m15.4113.7914.1213.3511.2412.632.892.022.79The primary cause is that none of the existing models could estimate the path loss at any height without any LAI or DBH. Furthermore, the value of RMSE, MAE, and MAPE in Okumura-Hata-Rural and LITU-R highly vary from one height to another, therefore, a slight variation of antenna or forestry parameters would negatively affect the prediction of the path loss in these models. Our proposed new model outperforms the existing models with the lowest RMSE, MAE, and MAPE regardless of the height of the antenna. Thus, the predicted values of our model are closer to real observations than the existing models.


## 4.2. Field Measurement on Energy-Harvesting LoRa Nodes

In this section, we evaluate the performance of the proposed radio propagation model based on the wind energy harvesting LoRaWAN which has been mounted in the sample mixed forest at Nanjing Forestry University [22]. We take the wind-powered sensor nodes serving as transmitters (Figure 10 shows the node’s mainboard, while the sensors could be attached on when in practice). An unlimited power supply sink is responsible for receiving messages and returning the packet delivery ratio (PDR) every day. The transmitting rate is 250 bps, and the payload length is 50 bytes. The selected nodes are located at 3 different positions: (i) LAIi=1.99, DBHi=13.82; (ii) LAIii=3.01, DBHii=10.78; (iii) LAIiii=2.63, DBHiii=15.84. The corresponding communication distances and antenna heights are (i) di=37m, hi=1.0m; (ii) dii=55m, hii=1.3m; (iii) diii=64m, hiii=2.0m.Figure 10
Wind energy harvesting LoRa node (mainboard).We validate the effectiveness of our propagation loss calculation module by using two comparisons. The existing fixed transmission power scheme is implemented as a reference, which is referred to as the Fixed-TX scheme. One of the wind-harvesting LoRaWAN nodes works under our proposed model at position (iii), while the other nodes located at positions (i) and (ii) be set with the fixed transmission powerTX=7dBm and TX=14dBm, respectively. The node (iii) could compute its transmitting power by equation (1) based on the RSSI value efficiently. All three nodes start with an initial energy of 1 Joule, and the average wind speed during the experiment is 4.4 m/s.Table6 presents the average PDR of total of 10,080 packets (ten packets per minute during 7 days) and the residual energy of three nodes in the final time. We can see that the proposed method achieves the lowest energy consumption compared with the other transmission power strategies no matter the distance. Compared with the Fixed-TX-14 dBm scheme, our proposed module achieves energy reduction around 30%. Compared with the Fixed-TX-7 dBm scheme, our proposed module improves the PDR around 5%. The wireless path loss and the background noise are always time-varying, the fixed transmission power definitely cannot guarantee the optimal transmission energy efficiency. The proposed PL model is effective in dynamically calculating the optimal transmission power along with the lowest energy cost.Table 6
The residual energy and PDR of three nodes.
Average testing resultProposed PL modelFixed-TX-14 dBmFixed-TX-7 dBmSupercapacitor voltage (volt, residual energy)3.1942.3722.468PDR (%)1009795


## 4.3. Rain Attenuation of LoRaWAN in the Forest

Rain and fog is quite common in the forestry district. Many studies indicate that the propagation of the radio signal is affected by rain and there are some remarkable models deduced to analyze its influence [30, 31]. However, most existing researches focus on the rain attenuation by wireless signals with frequencies above 10 GHz. So in this part, we investigate the influence of rain attenuation within the proposed forestry LoRaWAN for the first time. The experiment setup is just the same as part 4.2 and we record the variations of RSSI in position (ii) every 10 seconds when rainfall comes. The measurement data presented in Figure 11 are filtered by a moving average window with 30 samples. It is noticeable that during rainfall periods there is an obvious attenuation of PL.Figure 11
Variation of RSSI vs. rainfall in position (ii) from 2019.05.01 to 2019.05.08.
(a)(b)Figure11(a) shows a significant impact on path loss of LoRa signal when it is rainy during a 7-day consecutive observation, every increasement of PL values corresponds to a rainfall event accurately. Although the rain rate which has been recorded by a standard rain gauge [32] does not exceed the extremely heavy level, a nearly 8 dB variation of PL has occurred. Figure 11(b) shows a detailed changing process of rain attenuation, PL value increases sharply when the rain begins to fall, and it returns to the original state relatively slow after the rain stops. The main reason of fading may be the water molecular which is adhered on the antenna of LoRa nodes and pervaded in the transmission space, it would change the propagation characteristic and the antenna property drastically.According to reference [33], the rain attenuation can be represented by aRPb, where RP is the rain rate with a given period P while a and b are depending on the radio wave frequency and the rain temperature. Hence, we calculate this item with the collected rainy data and deduce the result as RP is set to 3.86 mm/hr based on the conventional rain information from local weather station, a and b just are 2.11 and 0.55, respectively. So the updated general PL model would be:
(13)PLLNdB=PLd0+10nlog10dd0+Xσ+aRpb.Furthermore, we also apply the updated model on the wind-harvesting LoRaWAN node to verify its effectiveness. It is located in position (iii) as well and the average wind speed is 4.9 m/s during the rainy period. The other two nodes located at positions (i) and (ii) are set the same as part 4.2, and the corresponding results are shown in Table7.Table 7
The residual energy and PDR of three nodes during rainy time.
Average testing resultProposed PL modelFixed-TX-14 dBmFixed-TX-7 dBmSupercapacitor voltage (volt, residual energy)3.0312.4472.584PDR (%)1009390It can be seen that the residual energy of the updated model decreases slightly otherwise the two Fixed-TX nodes get higher due to the larger wind speed. However, PDR of the updated scheme keeps 100%, while the comparative two nodes drop significantly. This indicates that the updated PL model matches the real rain attenuation efficiently and works stably in energy neutral mechanism.


## 5. Conclusion

IoT in forest environments can support applications such as fire alarm, environment surveillance, pest diagnosis, and tree growth monitoring. Large-scale deployment of LoRaWAN in forest environments requires accurate characterization of the propagation channel utilizing practical sensor nodes. However, most of the existing studies in forest environments use signal generators instead of practical sensor nodes, which, in turn, may lead to inaccurate models and, as a result, poor decision making during large-scale deployment of LoRa nodes. Inaccurate models may also result in the poor energy efficiency of the network as well as inaccuracy in localization and coverage control applications.In this work, empirical path loss models for LoRaWAN in the mixed forest of Eastern China are researched using practical sensor node measurements at 433 MHz. New proposed model based on the DBH and LAI with different heights is developed and compared with the existing Okumura-Hata-Rural model and LITU-R path loss models to demonstrate its accuracy. Moreover, the performance of wind energy harvesting based LoRa nodes equipped with our model has also been investigated. Results show that the RMSE, MAE, and MAPE values of the new model were almost all smaller than the existing models, thus the proposed model best matches the real field application status of LoRaWAN in mixed forest.The measured data, as well as the proposed models, can be useful for efficient planning and deployment of LoRaWAN in forest environments, even during the rainy season. In addition, a future study will include using the same approach to investigate the behavior of other LPWAN protocols like LTE-M and NB-IoT in forest environments.



---

*Source: 2687148-2020-06-22.xml*