Evaluation of LoRa Network Performance for Water Quality Monitoring Systems

1. Introduction

Industry 4.0, uniting both the physical and digital realms, has revolutionized the industry value chain. The transformation-constructed smart machine Internet-based notorious Internet of Things (IoT) enables data exchange between machines. IoT reduces human interference by introducing a remote monitoring system that can optimize real-time connection values. Remote monitoring systems improve the reliability of data collection, analysis, and operations [1]. The implementation of IoT in autonomous monitoring has transformed the potential of various industries, including the set of environmental applications of Smart Earth technology [2]. The increase in IoT demand in Smart Earth applications for its efficiency and precision services while operating at low cost and energy is parallel to its strategy for ecosystem management.

The Anthropocene epoch contributed to the World Health Organization’s (WHO) instigation and formulated the Sustainable Development Goals (SDG) blueprint to protect the planet and promote a healthy lifestyle. Objective 3 of the SDG progresses towards healthy lives, including the provision of target 6 of the SDG, protection of water ecosystems, and management of sanitation of water resources from pollution and scarcity [3]. Consequently, Malaysia has adopted Integrated Water Resources Management (IWRM) and Integrated River Basin Management (IRBM) programs to regulate water resources, predominantly rivers, because they are the main water resources [4]. However, only 47% of the water resources in Malaysia are classified as clean, while the remaining 43% are slightly polluted, and the remaining 10% are polluted due to poor management of water resources. Even during the pandemic Covid-19 between the years 2020 and 2021, the rivers of Sungai Gong and Sungai Semenyih in Selangor were polluted by hazardous chemicals and agricultural waste [5]. The unscheduled water supply due to water shortage affects the economy and productivity, as it drives the water plant to cease operations for cleaning purposes.

Therefore, monitoring the water quality conditions based on physical, chemical, and biological characteristics is the primary operation to support IWRM and IRBM [6,7]. The analysis of qualitative data on water quality characteristics was supervised to classify water bodies based on the Water Quality Index (WQI) and National Water Quality Standards (NWQS) in Malaysia [8]. The conventional method for monitoring the water quality requires a period of laboratory analysis. Thus, the assistance of autonomous water quality monitoring can provide an analysis of water quality properties that can vary over time [9]. The modification of autonomous water quality monitoring with IoT implementation acts as a quick decision-maker for water environmentalists based on real-time data. This reduces human labor and increases the reliability of data [10].

The employment type of the wireless network, which is the main element in IoT operation, necessitates the purpose of the system, type of experimental data, and location of the study. The invention of long-range (LoRa) networks in the wireless network community has become a game changer in IoT applications [11]. LoRa is the derivation of Chirp Spread Spectrum (CSS) technology to operate at long-range using an air interface [12]. This technology allows a low-type data rate to be transmitted within geographical and urban areas by manipulating its vast network coverage at a low power consumption and low cost. Relating to the ability of LoRa to operate with a large volume of sensors with low-type data rate output, this research is intended to study the design of a continuous water quality monitoring system with the LoRa approach performance evaluation.

The remainder of this paper is organized as follows. Section 2 discusses related work on the implementation of wireless networks as IoT solutions in water quality monitoring systems and reviews LoRa technology and its applications. The parameters, development, and design of autonomous water quality monitoring are elaborated in Section 3, including descriptions of the LoRa network parameters to be evaluated. The results of the experiment are explained and analyzed in Section 4, and lastly concluded and are referred to for future reference in Section 5.

2. Literature Review

2.1. Water Quality Monitoring Based With IoT Application with Different Wireless Network

The classification of water through water quality monitoring serves a broad purpose in water resource conservation based on the assessment of pollution control, assuring the sanctuary of biological ecology and health status for living organisms, and evaluating the current water quality trends, efficiency of waste management, and water treatment. The development of a continuous water quality monitoring system that delivers real-time data overcomes the drawbacks of manual methods in improving the accuracy of water quality data analysis.

Continuous water quality monitoring using IoT-based architecture deployed Wireless Sensor Networks (WSN). It dedicates sensors to collecting water quality data that require a wireless network as a bridge to communicate with end devices that act as a portal to the IoT platform. The significance of determining the types of wireless networks to be implemented in a water quality monitoring system reflects the capabilities and limitations of the application objectives. Effective IoT execution includes balancing the power, cost, efficiency, and maintenance [13].

Adopting various types of wireless networks into the water quality monitoring system, IoT has been proposed and evaluated by previous researchers. A Personal Area Network (PAN)-type network is a solution to IoT technology, offering its competencies to operate in closed proximity applications while being cost-and energy-effective. Moldovaeba et al. and Wang et al. used Zigbee as their PAN radio frequency module for continuous water quality monitoring [14,15]. Data were transmitted through the GPRS from the GSM version, extending the communication range. Signified by water quality data rates, the merger of two networks for extension unnecessarily consumes additional energy and cost, as the SIM card is required for communication activation.

Wi-Fi is another notorious network technology in IoT applications owing to its high data rates and low-cost antennas. Dilshad et al. and Naj et al. implemented Wi-Fi as a communication network for water quality monitoring system [16,17]. The vast amount of data transmitted improved the accuracy of water quality analysis with low latency. However, the limited range of Wi-Fi requires additional Wi-Fi extenders for river environments, which increases the utilization of operational costs and power.

Joseph et al. subscribed to the 3G network as a communication medium for water quality monitoring applications that transmit data to a remote database through Short Message Service (SMS) [18]. The use of a cellular network extended from the GSM version is secure and has low latency. However, it is more suitable for IoT applications with high data throughput and overwhelmingly low-type data sensors. In Stockholm, Ericsson collaborated with municipal councils and universities to deploy IoT systems with massive sensors and narrowband NB-IoT to monitor water quality throughout Stockholm City [19]. This system can increase the efficiency of sensors in longer and larger areas while delivering direct feedback for water quality classification. Nevertheless, both networks require expensive reserved channels to exploit implementation costs.

2.2. LoRa Overview

LoRa is a network solution categorized as low-power, wide-area networking (LPWAN) that allows long-range communication between devices designated for low-bit-rate data applications operated on batteries. The development of LoRa aims to exclude repeaters, improve network capacity, and support a large number of devices [20]. The LoRa solution was composed of two layers of LoRa and LoRaWAN that were significant to each other.

2.2.1. LoRa and Radio Modulation

The radio modulator-based CSS modulation technique in the physical layer (PHY) process gain in the radio frequency channel increases the data signal in the total signal spectrum. It maintains the signal data rate while reducing the output power for the transmission. CSS employs the orthogonal spreading factor (SF) method to synchronize transmission within a similar channel frequency. However, this affects the frequency offsets within the transmission and receiver, which are diminished in the decoder. Furthermore, it increases resistance to the Doppler effect. In addition, the SF method optimizes the power levels and data rates of each node to preserve the battery life and improve the efficiency of the receiver to decode signals [21]. The signal can be demodulated to an even lower signal-to-noise ratio (SNR) below the noise floor between -7.5 dB and -20 dB.

The integration of SF with the bandwidth (BW) and coding rate can deduce the expected data rate (DR). The relationship between SF and BW, deducing DR, is given by (1).

$$DR=SF*\frac{BW}{2^{SF}}$$

(1)

The chirp $2^{SF}$ used in LoRa doubles the bandwidth correlated to the symbol, and the bit rate increases the transmission rate. The relationship between BW and double SF and the duration of package transmission ($T_S$) is given by (2).

$$T_S=\frac{2^{SF}}{BW}$$

(2)

The generated diverse data transmission differs from the difference in the spread factor, as the LoRa data transmission can be determined every second. The symbol rate is denoted by $R_s$ expressed as in equation (3).

$$R_s=\frac{BW}{2^{SF} }$$

(3)

The sync word encoder in the LoRa frame differentiates LoRa networks that use the same frequency band. The LoRa frame, formatted with a short header, consists of information on bytes, Cyclic Redundancy Check (CRC), and Coding Rate (CR), as illustrated in Figure 1.

The influence of bandwidth on decoder sensitivity causes a change in CR that impacts the Packet Error Data (PER) when quick erupted interference occurs [21]. The effective bit rate ($R_b$) which is significant for the correlation of CR, BW, and SF in Equation (1), is given by Equation (4).

$$R_b=SF\times \frac{BW}{2^{SF}}\times CR=\left(SF\right)R_s$$

(4)

The LoRa in the PHY layer manages the radio frequency signals within the transmission line, including the frequencies, modulation waveforms, and power levels. Reflecting the condition of non-line-of-sight (NLOS) and multipath propagation factors in urban environments requires authentication in wireless medium characterization, such as path loss, shadowing, and multipath fading. Furthermore, the signal values obtained by the IoT platform must adapt to the sensitivity of the transceivers.

Equation (5) defines the corresponding losses and gains of the received power, $P_{rx}$at the receiver:

$$P_{rx}=P_{tx}+G_{tx}-L_{tx}-L_{rx}-L_{fs}-M+G_{rx}$$

(5)

Where $P_{tx}$ is the transmission power, and the gains at the transmitter and receiver are denoted as $G_{tx}$ and $G_{rx }$ respectively. The $L_{tx}$ and $L_{rx}$are the expressions of the loss antenna and cable connectors at the transmitter and receiver, respectively, including the loss from environmental influence, $L_{fs}$. M is the fading margin between the sensitivity of the receiver and the signal strength level.

In free space propagation, attenuation occurs as the power spreads over the areas. The power flux, which is the signal strength, is calculated as shown in (6).

$$P_d=\frac{P_{tx}}{4\pi d^2}$$

(6)

The effective area of an isotropic antenna, $G_e$ is used in LoRa, as shown in equation (7).

$$G_e=\frac{{\lambda }^2}{4\pi }$$

(7)

Equation (8) defines received power in free space.

$$P_{rx}=P_d{G}_e$$

(8)

$$P_{rx}=P_{tx}\times \frac{{\lambda }^2}{4\pi }$$

(9)

Path loss is generally defined as in equation (10):

$$L_{fs}=P_{tx}-P_{rx}$$

(10)

Substitute equation (9) in (10):

$$L_{fs}=20{\log }_{10}\left(4\pi \right)+20{\log }_{10}d-20{\log }_{10}\lambda$$

(11)

In free space, as λ (in km) = 0.3/f turn equation (11) into equation (12):

$$L_{fs}=32.44+20{\log }_{10}d+20{\log }_{10}f$$

(12)

The propagation path loss model of the earth plane entailed attenuation from the earth plane influence in Equation (13).

$$L_p=40{\log }_{10}d-20{\log }_{10}{h}_t-20{\log }_{10}{h}_r$$

(13)

Although there is still a non-definite path loss model for LoRa evaluation, the empirical formulation of the Okumura-Hata model, which can derive the propagation loss in an urban environment, is given by (14).

$$L_H=69.55+26.16{\log }_{10}f-13.82{\log }_{10}{h}_T-a\left(h_R\right)+(44.9-6.55{\log }_{10}{h}_T){\log }_{10}d$$

(14)

where $h_R$and $h_T$are the heights of the receiver and transmitter, respectively; f is the operating frequency (MHz); d is the distance between the transmitter and receiver; and $a\left(h_R\right)$is the correction parameter that is affected by the type of environment. $a\left(h_R\right)$for a large city environment (refer to (15)).

$$a\left(h_R\right)=3.2{[{\log }_{10}(11.75\bullet h_R)]}^2-4.97$$

(15)

Following the limitation of the Okumura-Hata model, which is applicable to land mobile services with a frequency range of 100-1500 MHz, distances of 1-20 km, and base station height of 30-200 m, Petajajarvi et al. obtained the maximum range for different transmitted powers with a constant variable height from the Okumura-Hata model [22]. Furthermore, the received power obtained was as low as -120 dBm for the two transceivers, resulting in a lower signal translation. The analysis by Oliveira et al. on the LoRa communication ranges impacted from the surroundings resulted in maximum communication ranges of 5.6 km and 2 km in rural and urban environments, respectively [23]. Furthermore, the received power obtained was as low as -120 dBm for the two transceivers resulting in a lower signal translation.

2.2.2. LoRaWAN

LoRaWAN is an open network protocol standardized by the LoRa Alliance located at the MAC layer to secure bi-directional communication that known as ALOHA. The ALOHA protocol allows uplink transmission at low power, followed by two short downlinks that receive the windows. LoRaWAN utilizes a star topology to connect end-devices, gateways, and network servers for information exchange. LoRaWAN applies the ADR method to the network server to initiate the increment of the spread factor, which allows a higher delivery ratio between the end-nodes and gateway when a lower spread factor fails.

LoRaWAN offers two types of activation for security verification and service quality authorization: over-the-air activation (OTAA) and Activation by Personalization (ABP). LoRaWAN reflects the requirements of the application address, downloading, and energy limitations. It scheduled three classes for downlink communication: Class A, Class B, and Class C. The default Class A is supported by all LoRaWAN end-devices, whereas Classes B and C are optional supports. Class B mode scheduling opens downlinks from periodic beacons and receives the network windows with regulated latency. The Class C mode keeps the receivers of end devices open except during transmission, resulting in fading of latency [24]. A summary of the class mode applied by the LoRaWAN is shown in Figure 2.

Ertürk et al. found that the performance of LoRaWAN on successful packet transfer is affected by high volumes of end-nodes after simulating the Class A mode on the MAC layer using MATLAB [25].

Delobel et al. analyzed the delay of confirmed downlink frames in class and concluded that the gateway was influenced by a sudden data outburst, causing interference in radio transmission [26].

2.2.3. Implementation of LoRa in Application from the Previous Study

Since the LoRaWAN protocol was officially publicized in 2015, it has been studied, analyzed, and implemented in various IoT applications. The performance of LoRa in an urban environment was evaluated using propagation measurement by Parades et al., who were able to predetermine the interference that occurs during transmission [27]. Rizzi et al. used the Okumura-Hata model to calculate the path loss in an urban environment, requiring a reduction in the signal-to-noise ratio with a higher antenna gain to increase the transmission range [28].

A flood warning system based on a Wireless Sensor Network (WSN) developed by Leon et al. used an ultrasonic sensor alerting a user on Twitter with LoRa as a data exchanger and found no loss of packet data, but it had a 1 cm margin error of the water level [30]. Codeluppi et al. built an ad-hoc-level module with a LoRaWAN architecture, and transmitted sensor data to provide advanced information to farmers [31]. The results differed based on its location in the greenhouse and vineyard, and there were also few samples owing to the lack of solar radiation during rainy days.

Abbasi et al. utilized LoRa in smart grid applications and found that LoRa is sensitive to an increase in the number of nodes, thereby increasing the energy consumption [32]. Furthermore, the employment of a directional antenna approximately fixed the Data Extraction Rate (DER). Prakosa et al. monitored the soil state based on IoT with LoRa for smart agriculture applications in rural areas and concluded that the large coverage influenced by high SF increased delays in data transfer [33]. The recent epidemic of the SARS-Cov 2 virus initiated Lousado and Antunes to monitor the health of elderly people using LoRa perceived that high repetition was influenced by constant connection and LoRa capability was limited owing to its low coverage in that region [34].

Different applications and scenarios result in differences in the LoRa performance. The LoRa network can be modified and enhanced for effective utilization based on its application.

3.1. Water Quality Monitoring System Built

The type of water quality sensor was selected based on water characteristics regulated by the WHO. The robustness and reliability of the water quality sensors embedded in the water quality monitoring station are based on the specifications and datasheets supplied. The water quality sensors that were used are listed in Table 1.

The water sensors were connected to a microcontroller board (TTGO ESP32-LoRa32) at an operating frequency of 868/915 MHz, which is suitable for the Malaysian region. The end-nodes deployed the Frequency Shift Keying (FSK) modulation mode, equipped with a data rate of 1.2 kbps to 300 kbps. The antenna installed at the water station was of the whip 915 MHz omnidirectional type. Table 2 and Table 3 describe the LoRa end-node parameter setup with different SF for the antenna gain of 2.1 dBi and 3 dBi.

The prototype water quality monitoring station was powered by a polycrystalline solar panel of 6 W/6 V, as shown in Figure 4

Water stations P1, P2, and P3 were located at three different locations in both rivers at the Universiti Malaya, Sungai Anak Air Batu and Sungai Pantai, at distances of 117 m, 1560 m, and 566 m, respectively, from the gateway, as shown in Figure 5.

3.2. LoRa Gateway Configuration

The Raspberry Pi HAT was assembled using an LPWAN concentrator module deployed as the LoRa gateway. The gateway module includes a GPS module and a Heat Sink. The package is based on the Semtech transceiver concentrator, and the S1257/58 front-end chirps allow the management of packets from the scattered endpoints.

The gateway operated with the support of the global license-free frequency AS920-923, which corresponds to the location of Malaysia within the Asian region. The antenna gain used for the gateway was 5.8 dBi for a larger spread and better perception in an urban environment.

The gateway was integrated with Balena.io and connected to Thing Stack V3.16.2, using the latest Semtech Packet Forward protocol. The water quality data and performance of LoRa on the Received Signal Strength Indicator (RSSI), signal-to-noise ratio (SNR), downlink, uplink, and packet data can be directly monitored from the server, The Thing Stack.

The system does not require a full-time alert for the end nodes and allows them to fall into the sleep mode. Its purpose is to minimize the power consumption of the end nodes and prolong battery life.

The gateway was positioned on the rooftop of Block C Faculty of Engineering at an estimated height of 25 m from the ground, as shown in Figure 6.

3.3. Water Quality Monitoring Dashboard

The Thing Stack server was programmed to interconnect with the IoT analytic platform, ThingSpeak. ThingSpeak was trained as a cloud to observe, record, and store daily water quality monitoring activities. The recorded water quality data were extracted from the cloud and analyzed for water quality status classification.

An Android APK application was developed using MIT Inventor Developer integrated with ThingSpeak as the GUI. Users can acquire and monitor the water quality status of Sungai Anak Air Batu and Sungai Pantai daily. The display includes pH, Total Dissolved Solids (TDS), Dissolved Oxygen (DO), temperature, and NTU values from the water stations at Sungai Pantai and Sungai Anak Air Batu. The GUI of the Android version of the water quality monitoring dashboard is shown in Figure 7.

3.4. Evaluation of LoRa Performance

Several network performance indicators were evaluated to verify the implementation of LoRa in a continuous water quality monitoring design. This section discusses the network performance that was evaluated in this study.

First, the commonly inspected network characteristic RSSI estimates the strength of the radio signal received during transmission to determine the capability of the device to listen to a signal. The value of the RSSI is affected by the location, Line of Sight (LOS) variable, n, and distance between the node and receiver, d, as shown in equation (16). In addition, the background noise energy influences signal strength, causing communication errors. The Expected Signal Power (ESP), referring to equation (17), eliminates the impact of the noise effect on the RSSI in equation (16) and acts as an indicator to measure the exact received power.

$$RSSI=-10n{\log }_{10}d+A$$

(16)

$$ESP=RSSI+SNR-{10\log }_{10}\left(1+{10}^{0.1SNR}\right)$$

(17)

The SNR in (17) differentiates the received signals within the spectrum from unfeasible background signals. The SNR is simplified in (18).

$$SNR=\frac{P_{Sn}}{P_n}$$

(18)

$$SNR \left(dB\right)=10{\log }_{10}\frac{P_{RX}}{P_n}$$

(19)

Where $P_{Sn }$is the signal power, and $P_{Sn}$is the noise power. A higher signal power can demodulate the noise signal, thereby eliminating the potentially corrupted signals. The corrupted signals below the noise floor initiate retransmission, which requires airtime in wireless networks and hence degrades throughput and latency. Nevertheless, the LoRa technology can demodulate the minimum signal below the noise floor based on its spreading factor.

Next, the observed network characteristic, packet loss, is usually affected by errors in the lack of signal strength at the receiver, severe system noise, overload of network nodes, network congestion, hardware issues, and distances. Degrading communication from packet loss reduces throughput, security communication, and inadequate data transmission, thereby causing loss in data encryption.

In this study, LoRa is a wireless network that applies at UDP protocol. Although it tolerates packet loss, the sender cannot detect whether packets have been received. The server, The Thing Stack, applies the Semtech UDP Packet Forwarder. The monitoring of packet loss was observed and recorded through the gateway log. Packet loss can be defined by Equation (20).

$$PL \left(\%\right)=\frac{P_{st}-P_{rc}}{P_{st}\times 100}$$

(20)

Where PL is the packet loss defined over time with $P_{st}$, sent packet and $P_{rc}$, the received packet.

In this study, the conclusion from the K-means clustering result of water quality properties will be integrated with packet loss to assist in detecting patterns in network traffic data. This pattern allows the deduction of various possible factors during transmission, such as hardware failure, network congestion, and signal interference.

One of the influential factors for all network performance is path loss because it determines the energy loss within the transmitter and receiver. The general path loss is given by Equation (21).

$$L_p= P^t+G^r+G^t+P^r$$

(21)

Where$L_p$, the path loss, which is a combination of $P^t$, $G^r,G^t,$ and $P^r$, which are the antenna transmission power, receiver power gain, transmitter power gain, and RSSI, respectively. The path loss with wavelength or frequency is denoted by (22).

$$L_{fs}={\left(\frac{4\pi d}{\lambda }\right)}^2$$

(22)

$$L_{fs}=10\log \left(\frac{{\lambda }^2}{{4\pi }^2{d}^2}\right)$$

(23)

The carrier frequency in this study, 915 MHz, changes the value of λ as in equations (24) and (25).

$$\lambda =\frac{v}{f}$$

(24)

$$\lambda =\frac{3\times {10}^8}{915\times {10}^6}=0.3279m$$

(25)

In a free-space damping propagation application, the usual logarithm extended from (23) and (12) amalgamates with the energy spread and antenna fault, as given by equation (26).

$$L_{fs}=20{\log }_{10}d+91.67-G_{tx}-G_{rx }$$

(26)

where $G_{tx}$ and $G_{rx}$are the transmitter antenna gains, and the antenna receiver gain includes feeder losses. The increment in distance increased the loss in free space by 6 dB. Reflection, refraction, and penetration of radio waves affected by the attenuation structure are also significant for the impact of energy losses on the budget link. Referrings to equation (14) in the previous section, the Okumura-Hata module is used for propagation loss in an urban environment.

The last network performance analyzed is the received power because its quantification links the performance from the link budget. Temporary attenuation during transmission requires the fade margin to specify its power level performance. The minimum power level is determined using the link budget related to the received sensitivity. The determination of the received sensitivity corresponds to the minimum power received by the receiver node to decode the transmitted bit promptly. The receiver sensitivity, $R_x$applied Equation (27) tolerance to thermal noise.

$$R_{xs}=-174+10{\log }_{10}BW+NF+SNR$$

(27)

Where BW is the bandwidth and NF is the noise factor, referring to equation (27). The received power declines as it passes through the channel over a distance and in the environment. The received sensitivity in LoRa can be below -130 dBm to allow interpretation of the lower signal. The received power accumulates all losses and gains, as mentioned in Section 2, referring to equation (5). Furthermore, setting the value of SF also affects the sensitivity, as it is used to set the data transfer against the range. The setup bandwidth of the receiver operated at 922 MHz and the noise level generated in this study also affected the receiver power level.

4. Result and Discussion

The monitored water quality data were collected and analyzed to validate the compatibility of LoRa as a linkage for autonomous water quality monitoring systems. This segment also observed and analyzed the trend of the LoRa performance based on the network characteristics selected in the previous section.

4.1. Water Quality Data

Users can access the water quality status from the application dashboard, as shown in Fig. 6 in the previous section. The trend of water quality data for each water quality property for the station, as mentioned in the previous section, can be previewed on the Thing Speak monitor dashboard, as shown in Figure 8.

For further analysis, water quality data were extracted from Thing Speak, which acts as a cloud for data storage. Water quality data were analyzed based on WQI and NWQS by Malaysia for pH, Dissolved Oxygen (DO), turbidity, temperature, and Total Dissolved Oxygen (TDS). Table 4 shows the results of the K-means pre-processing of raw water quality data for 1 month’s observation. Moreover, K-means allows for the detection of abnormalities in unsupervised water quality data during data collection. The observation relating the K-means results of water quality data and network performance is further discussed in Section 4.2.2.

Based on the tabulated data in Table 4, the status of water quality at each station can be classified based on the NWQS standard guide. The pH values for all water stations had similar ranges, with near-neutral pH values of 7. Referring to the WQI and NWQS, the pH of all stations qualified to be included in Class Ⅰ because its range remained within 6.5–8.5. For turbidity, which measures the salinity of water, the value obtained from station P1 was slightly different from P2 and P3, with an estimated difference of 30 NTU. However, the water status from all stations classified was included in Class Ⅲ because the turbidity standard tabulated at 0 NTU-50 NTU is above Class ⅡB.

The standard river water temperature in Malaysia ranges from 22.0°C to 31.7°C, depending on the seasonal effect. The tabulated data in Table IV show that the temperature of river water from all locations falls within the standard range; hence, the class of water temperature based on the WQI and NWQS remains Class Ⅰ. However, the water temperature at station 2 was 2°C higher than that at the other stations. Therefore, according to NWQS Malaysia, an additional 2°C from the current average temperature at that time classifies the water as class Ⅱ.

Although the total dissolved solids contained the water from all water stations fluctuated marginally from each other, their quality remained within Class Ⅰ because the range still fell below 500 mg/l. The last water quality property monitored in this study, dissolved oxygen, was constant at all locations, with a total volume of 8 mg/l. Based on the WQI and NWQS guidelines, the total volume of dissolved oxygen in the river water in this study is suitable for aquatic organisms and is labeled as Class Ⅰ.

The conclusions from the analyzed water quality properties can be used to classify water into water quality classes based on the WQI and NWQS benchmarks. The class of river at each water station at Sungai Anak Air Batu and Sungai Pantai, Universiti Malaya , was summarized in a pentagonal shape as five parameters of water quality were monitored in this study, as shown in Figure 9. The shape and blue color of the pentagon effectively visualizes the summarization of real-time monitoring of water quality to support the study of sampling site data [7].

Based on the summary pentagonal shape shown in Figure 9, the water bodies at all locations exhibited a similar trend. The full pentagonal shape classifies the water body as Class Ⅰ, but the trend of water in the figures above for all stations had an imperfect shape, as it was affected by the water quality class of turbidity. Water bodies at station 2 were affected by the temperature and turbidity class. In conclusion, the water bodies at Sungai Pantai and Sungai Anak Air Batu at Universiti Malaya are classified as class ⅡB, which are suitable for recreational purposes but are undrinkable. Furthermore, conventional treatments are required for the water supply.

The absolute conclusion based on the WQI calculation of its water quality properties subindex (SI), as shown in (25), cannot be defined as the presence of water quality elements such as Ammoniacal Nitrogen (AN), Biochemical Oxygen Demand (BOD), and Chemical Oxygen Demand (COD), which cannot be obtained with real-time methods and require laboratory analysis.

While 0 ≤ WQI ≤ 100,

$$WQI=\left(0.22\mathrm{*}SIDO\right)+\left(0.19\mathrm{*}SIBOD\right)+\left(0.16\mathrm{*}SICOD\right)+\left(0.15\mathrm{*}SIAN\right)+\left(0.16\mathrm{*}SISS\right)+\left(0.12\mathrm{*}SipH\right)$$

(28)

Where SIDO is the subindex of dissolved oxygen, SIBOD is the subindex of biochemical oxygen demand, SICOD is the subindex of chemical oxygen demand, and SIAN and SIpH are the subindices of ammoniacal nitrogen and pH, respectively.

The result of river classification in Universiti Malaya from the early prediction in this study is significant compared to a previous study by Gafri et al., who concluded that the study area contained a moderate level of pollution due to urban activities [7].

Nonetheless, early prediction from real-time systems allows water environmentalists to periodically monitor the health of water bodies as an early step toward conserving clean river water.

4.2. LoRa Network Performance

The importance of monitoring LoRa network performance during water quality data transmission corresponds to the efficiency of LoRa as a communication bridge for water quality monitoring in real-time. The analysis of the network performance for the radio coverage of LoRa is discussed in this section based on the setup of different transmitter antennas, as discussed in the previous section.

4.2.1. RSSI and SNR

The recorded RSSI and SNR from The Thing Stack were measured through standard deviation to observe its spread variance, which was affected by differences in the antenna gain, spread factor (SF), distances, and NLOS. Table 5 list the maximum minimum and mean RSSI for each LoRa node at all water stations.

The measurement of RSSI determines the quality of the received signal. The minimum acceptable range value of the RSSI for LoRa was -120 dBm, and a low RSSI indicated that the signal was weak. Table 5 shows that the lowest RSSI recorded was -118 dBm at the second node, P2, which was located farthest from the receiver for SF of 7 and 9. The lowest recorded RSSI was still above the minimum KPI for the RSSI of LoRa and wireless networks. The differences in the antenna gains and spread factors influenced the trend of the RSSI for each distance. Figure 10 and Figure 11 show the mean, maximum, and minimum trends of RSSI based on each antenna gain.

As shown in Figure 10 and Figure 11, the trends of the minimum and maximum values from the mean RSSI value are spread over certain distances. The increment of antenna gains by 0.9 dBi from 2.1 dBi based on Fig. 11, affects the increment of the maximum and minimum RSSI for all SF from all locations.

The use of a 3 dBi antenna increases the minimum value from the farthest location of node P2 at 1560 m from the receiver with the lowest SF of 7. The minimum and maximum values obtained were -116 dBm, and -110 dBm respectively. The RSSI of the nearest end-node at 117 m from the receiver also shows its best performance antenna gain of 3 dBi, as it reached above -90 dBm even when implied with SF 7.

Furthermore, the spread value of the RSSI was affected by distance, urban environmental influences, weather, interruption of other radio frequencies, and other variables. It can be concluded that the RSSI value of LoRa decreases as the distance increases with the effect of NLOS [37]. In addition, the higher gain of the antenna strengthens the received signal while reducing interference during transmission.

Although the SNR of LoRa can operate below the noise floor, the higher the SNR during data transmission, the higher is the quality of the signal transferred. The SNR collected simultaneously with the RSSI shown in Table 6 falls below the noise floor for all the minimum, maximum, and mean values from all locations for each implied spread factor.

The accumulated SNR sample value iss slightly different from the implied SF value for each distance. The capability of LoRa to allow the signal to travel below the noise floor constructed the KPI of the LoRa demodulator to translate the received signal for SF of 7, 9, and 12 as -7.5 dB, -12.5 dB, and -20 dB, respectively. The recorded SNR value stays above the KPI of the LoRa demodulator except for station P2, which was located 1560 m away from the receiver when SF 7 was implied for both antennas gain of 2.1 dBi and 3 dBi with its mean values of -13.41 and -12, respectively. Although the signal can travel below the noise floor, it cannot be interpreted by the demodulator; therefore it can affect the packet data. Further aspects of the SNR are also monitored in Figure 12 and Figure 13, which illustrate the minimum, maximum, and mean over distances.

Based on Figure 12 and Figure 13, antenna gains of 3 dBi increased all SNR values from station P1 above the noise floor, even with the implementation of SF 7. Nonetheless, SF 7 is inapplicable for station P2, as mentioned previously, because its SNR value did not surpass its KPI of -7.5 dB. Meanwhile, at water station P3, the use of SF 7 and an antenna gain of 3 dBi exceeded its KPI, with a minimum SNR value of -6 dB. Furthermore, the dispersion trend of the SNR for an antenna gain of 3 dBi was smaller and more stable compared with the implementation of an antenna gain 2.1 dBi.

The stability of the SNR trend supported by an antenna gain of 3 dBi can improve the optimization of the transmission and reduce undesirable signals. The increase in antenna gain also increases the signal power to overcome the NLOS and temperature factors [38]. High and stable SNR values can reduce packet loss, thus increasing throughput.

Relating the RSSI and SNR is significant to the expected signal power, as shown in equation (17). ESP information estimates the average desired signal in a noisy environment.

4.2.2. Packet Loss

The execution of the RSSI and SNR corresponds to the dropped packet data. Table VII lists the daily average packet losses for each period in which the end node transmits data to the receiver.

Referring to Table 7, packet loss from both station P1 and P3 remained below 10% for all conditions, whereas the highest average packet data dropped by 22% from location P2, which was embedded with an antenna gain of 2.1 dBi and SF of 7.

The occurrence of high packet loss was heavily affected by a low SNR, which could not be demodulated by the receiver, and the NLOS effect in the urban environment.

The effects of high packet loss degrade the network throughput and increase latency, thereby influencing the accuracy of the water quality data. In addition, the increment in SF can increase the latency, thereby affecting the time delay and packet loss, depending on the density of the environment [33]. However, the packet loss from station P2 was reduced with changes in SF and antenna gain, reaching a maximum of 5.4%.

Nonetheless, absolute prevention of packet loss in any network condition is impossible, but it can be minimized through preventive measures. Sudden changes in data packets support the prediction of network abnormalities.

Anomalies in water quality data caused by high packet loss corrupt and reduce the accuracy of the generated result. A previous method by Jáquez et al. detected an anomaly with a scalar-type program using a Boolean value [35]. They detect errors in the data owing to malfunctioning hardware from the sensor and the water pump.

In this study, packet loss was integrated with the results of K-means from water quality data to predict anomalies in the data. The observation of K-means from the water quality data integrated with packet loss from all water stations is shown in Figure 14, Figure 15 and Figure 16.

While the trend of K-means of water quality properties at water station 3 was stable, turbidity and TDS trends at station 1 increased and decreased respectively, in the middle of October 2022. However, its packet data operates with a stable trend and a low packet loss rate. Hence, it can be concluded that the water bodies at station 1 varied in salinity and dissolved content. Water bodies at Station 2 showed stable fluctuations in turbidity and TDS throughout the experimental timeline. There was a high packet loss at Station 2 recorded between November and December 2022 because of the affected transmitted packet data caused by the seasonal change in the monsoon season in Malaysia. Observations of packet loss and water quality data allow environmentalists to inspect and maintain the condition of water sensors and network performance. By using packet loss as an additional input feature, the proposed solution can better distinguish between true and false positives, resulting in an improved overall performance. The integration analysis resulted in enhancements in the water quality monitoring and network performance analysis.

Although, the trend of packet loss from station 2 had a high loss compared to the other stations, the packet loss rates were still within acceptable ranges for low packet data. A previous study by Wei et al. concluded that an increase in communication distance increases the packet loss rate [37].

Thus, monitoring packet loss in this study is important for troubleshooting the network performance and ensuring the efficiency of operational tasks.

4.2.3. Path Loss

Identifying path loss for wireless network applications in urban environments with possible power attenuation is important for defining the budget link. The path losses in this study were analyzed using the propagation of free space (FSDP), referring to (26), and the Okumura-Hata (Oku-Hata) model, referring to (14). The path losses are presented in Table 8.

Because the FSDP method includes antenna gain and feeder losses, the path loss with an antenna gain of 3 dBi showed a slight decline from each end node compared to the implementation of an antenna gain of 2.1 dBi. However, there were no differences in path loss when the antenna gain changed with the Oku-Hata method because its significant elements were the height of the transmitter and receiver, frequency, and distances. However, the Oku-Hata method results in a higher path loss than the FSDP method because it accounts for the effects of diffraction, penetration, and reflection of multiple obstacles in the exterior realm. The obtained FSDP for distance with an estimation of 1 km in this study was significant, as calculated by Ali et al. for a similar omnidirectional antenna used for both the receiver and transmitter [39].

Nevertheless, the use of the FSDP model can provide information for modifications to reduce the average free-space path loss.

4.2.4. Received Power

The computation of the transmission power, gains, and losses influencing the amount of received power has a significant effect on the total budget link. correlating to the link budget, the findings of the received power and received sensitivity are significant in concluding the value of the fading margin. The result of the fading margin allows the design specification to ensure the system performance level.

This parameter was not widely discussed in a previous study of the LoRa monitoring system, as it is known that the receiver sensitivity of LoRa can operate below -130 dBm. Hence, it allows the interpretation of lower signals; however it is crucial to determine the minimum strength to be detected and processed for transmission quality assurance. Table 9 shows a comparison of the total received sensitivity and received power.

The received power gain from all locations can be perceived by the gateway, because the received sensitivity is distributed at a minimum of -129.1 dBm. The received power analyzed from all locations remained below the received sensitivity values. The low value of the received sensitivity in this study allowed the receiver to interpret the weak signals transmitted from the LoRa nodes at each water station. The use of a higher antenna gain in this design system allows it to maintain a minimum fading margin with an estimation of 10 dB from all distances.

5. Conclusions

The objective of this study is to monitor water quality data in an urban environment using the LoRa network as the primary communication method. Real-time analysis of water quality data allows for a quick decision-making system for water environmentalists to prevent pollution outbreaks. Therefore, evaluating the LoRa network performance for developing an autonomous system is signified by its heterogeneous operating conditions in an urban environment. A detailed analysis of the RSSI, SNR, packet loss, path loss, and received power demonstrated the LoRa compatibility for independent operation. Changes in the network properties of LoRa, such as SF and antenna gain, heavily impacts the performance of the LoRa network. The increase in SF and antenna gain increases the distance, receiver sensitivity, and ability to overcome obstacles within the transmission line. The results can assist different configurations for environmental monitoring with an IoT-based LoRa network in urban settings. Furthermore, the large distance between the end-nodes can reduce the redundancy of packet data. The integration of observations of water quality data from the K-means method and packet loss can be used to analyze the anomaly deduction on instruments, changes in water quality properties and network performance. Although LoRa is known for its capability to operate below -130 dBm, enhancement of its network properties, such as bandwidth, can improve its efficiency in terms of its and data type. Finally, this study demonstrates the real-world performance of LoRa as a water quality monitoring communication system, highlighting its various network properties, within an urban area on campus.