Feasibility of Backscatter Communication Using LoRAWAN Signals for Deep Implanted Devices and Wearable Applications

This paper presents a method for low data rate transmission for devices implanted in the body using backscattered Long Range (LoRa) signals. The method uses an antenna loaded with a switch that changes between two load impedances at the rate of a modulating oscillator. Consequently, the LoRa signal transmitted by a LoRa node is reflected in the adjacent channels and can be detected with a LoRa gateway tuned to the shifted channels. A prototype developed to operate at Medical Implant Communication Service (MICS) and the Industrial Scientific and Medical (ISM) 433 MHz band is presented. The prototype uses a commercial ceramic antenna with a matched network tuned to the frequency band with high radiation efficiency. The effect of the coating material covering the antenna was studied. Simulated and experimental results using a phantom show that it is feasible to read data from deep implanted devices placed a few meters from the body because of the high sensitivity of commercial LoRa receivers.


Introduction
Long Range Wide Area Network (LoRaWAN) is one of the best wireless technologies for the Internet of Things (IoT) [1,2]. Its popularity is due to the robust modulation scheme based on chirp spread spectrum (CSS) modulation [3][4][5]. LoRa receivers exhibit an impressive sensitivity which, depending on the spread factor used, can be as low as −148 dBm (Semtech SX1276 [6]). As a result of this sensitivity, radio links of several kilometers long can be designed [7]. LoRa also provides advantages in both blocking and selectivity and can achieve high interference immunity with low energy consumption and low-cost oscillators. LoRa is based on a low power modulation designed to send small data packages over long distances, operating on a battery. It can demodulate signals which are 7.5 dB to 20 dB below the noise floor (see Table 13 in [6]. The performance of wearable and implantable devices is often limited by the lifetime of batteries. A considerable part of the power budget of these devices is spent by the radiofrequency circuits (oscillators, receivers, and transmitters) on wireless communication [8]. Recently, backscatter communications have been raising a great interest in a variety of applications including implanted devices [9][10][11][12]. Backscatter communication requires zero transmission power and works by backscattering the RF signal sent by an external transmitter [13,14]. This type of communication is used in Radio Frequency Identification (RFID) systems. Various backscattering systems compatible with commercial wireless transmitter signals have been proposed in the literature [15][16][17]. The simplest backscatters (see [14] for a survey of ambient communications) generate On-Off Keying (OOK) or Frequency Shift Keying (FSK) signals from an unmodulated carrier and are encoded by a custom receiver packets (with a spread factor of 12 and a bandwidth of 125 kHz) with the backscatter enabled (and the modulating oscillator enabled) and disabled. The measurement was made with an RTL-SDR software defined radio (SDR) connected to a monopole antenna using SDR Sharp software 1 m from the backscatter. The central channel is the transmitted channel, but when the backscatter is enabled, it reflects the signal shifted from the central carrier to the oscillation frequency (fosc = 300 kHz in this experiment).    The backscattered field of an antenna → E S is a function of the antenna load Z L . It can be split as a sum of two terms [21,22]: The first term is the structural mode → E est , which is independent of the load, and the second ( is the antenna mode (also known as tag mode).
→ E est is the scattered field when the tag is connected to a reference load Z L = Z a *, where Z a is the antenna's impedance. Γ is the power reflection coefficient given by [23]: When the load is the conjugate impedance of the antenna impedance, all the incident power is transferred to the load and the antenna only reflects the structural mode. The structural mode arises from the current induced on the antenna's conducting surface by the incident wave. It does not depend on the load but depends on characteristics such as antenna type, geometry, and material. Thus, the structural mode is independent of the load reflection coefficient, whereas the antenna mode is proportional to it.
The reflection coefficient can be modulated by switching the antenna load between two impedance states (Z ON , Z OFF ). As a first approximation, the reflection coefficient can be approximated by a square waveform with amplitude ∆Γ = Γ ON − Γ OFF . Therefore, this periodic signal with period 1/f osc can be developed in a Fourier series and the spectrum is given as an infinite sum of delta functions located at the harmonics: where c n are the Fourier coefficients, and f c is the input frequency of the incident signal that illuminates the tag. For a square waveform with duty cycle δ c , the coefficients c n are given by: Sensors 2020, 20, x FOR PEER REVIEW 3 of 26 to a reference load ZL = Za*, where Za is the antenna's impedance. Γ is the power reflection coefficient given by [23]: When the load is the conjugate impedance of the antenna impedance, all the incident power is transferred to the load and the antenna only reflects the structural mode. The structural mode arises from the current induced on the antenna's conducting surface by the incident wave. It does not depend on the load but depends on characteristics such as antenna type, geometry, and material. Thus, the structural mode is independent of the load reflection coefficient, whereas the antenna mode is proportional to it.
The reflection coefficient can be modulated by switching the antenna load between two impedance states (ZON, ZOFF). As a first approximation, the reflection coefficient can be approximated by a square waveform with amplitude ΔΓ = ΓON−ΓOFF. Therefore, this periodic signal with period 1/fosc can be developed in a Fourier series and the spectrum is given as an infinite sum of delta functions located at the harmonics: where cn are the Fourier coefficients, and fc is the input frequency of the incident signal that illuminates the tag. For a square waveform with duty cycle δc, the coefficients cn are given by: where Γavg is the average power reflection coefficient between ON and OFF states. Using (3), the backscattered field can be written as: Therefore, the first term in (5) results in the non-modulated term and depends on the structural mode of the tag. The second term represents the modulated side-bands that are a function of the antenna mode. Coefficients of higher amplitude in the Fourier expansion (n = ±1), are which result in the components at the frequencies fc ± fosc.
The result is that the LoRa channel at transmission frequency fTX is also reflected in the adjacent (4) where Γ avg is the average power reflection coefficient between ON and OFF states. Using (3), the backscattered field can be written as: Therefore, the first term in (5) results in the non-modulated term and depends on the structural mode of the tag. The second term represents the modulated side-bands that are a function of the antenna mode. Coefficients of higher amplitude in the Fourier expansion (n = ±1), are which result in the components at the frequencies f c ± f osc .
The result is that the LoRa channel at transmission frequency f TX is also reflected in the adjacent channels with central frequencies given by f TX ± f osc . Therefore, a LoRa receiver tuned to the frequency of one of these shifted channels at f TX ± f osc can demodulate the transmitted packed if the backscatter is enabled. Enabling and disabling the modulating oscillator can transmit low-rate data for telemetry purposes. LoRa gateways are intermediaries that allow sensing devices to transmit data to the cloud. Modern LoRa gateways are based on chipsets that can simultaneously receive LoRa messages on multiple frequency channels. Therefore, it is possible to demodulate the signal shifted by the backscatter to another frequency channel. The proof of concept uses the ADG902 (Analog Devices Inc., Norwood, MA, USA) CMOS switch [24] and the LTC6907 (Analog Devices Inc., Norwood, MA, USA) [25] oscillator, and the 0900AT43A0070 antenna (Johanson Technology, Camarillo, CA, USA) [26], which can be tuned to the desired frequency with some adjustments. LTC6907 is a low-power clock generator for which the oscillator frequency f osc can be adjusted to between 40 kHz and 4 MHz with a resistor that has a current consumption between 36-100 µA at 400 kHz to 1 MHz. The LTC6907 has 1% frequency accuracy from 0 • C to 70 • C, which is enough for this application. Figure 2 shows the system integrated into a small-size PCB. In addition, the LoRa receiver can measure the frequency offset in the received packet and therefore can be adjusted to minimize the effects of oscillator drifts. The experiments use the Semtech SX1278 LoRa transceiver (Semtech Corp In order to show the modulation of LoRa packets, Figure 3 shows the measured spectrogram (spectrum as a function of the time) of the modulated LoRa channel at 406 MHz continuously sending packets (with a spread factor of 12 and a bandwidth of 125 kHz) with the backscatter enabled (and the modulating oscillator enabled) and disabled. The measurement was made with an RTL-SDR software defined radio (SDR) connected to a monopole antenna using SDR Sharp software 1 m from the backscatter. The central channel is the transmitted channel, but when the backscatter is enabled, it reflects the signal shifted from the central carrier to the oscillation frequency (fosc = 300 kHz in this experiment).    In order to show the modulation of LoRa packets, Figure 3 shows the measured spectrogram (spectrum as a function of the time) of the modulated LoRa channel at 406 MHz continuously sending packets (with a spread factor of 12 and a bandwidth of 125 kHz) with the backscatter enabled (and the modulating oscillator enabled) and disabled. The measurement was made with an RTL-SDR software defined radio (SDR) connected to a monopole antenna using SDR Sharp software 1 m from the backscatter. The central channel is the transmitted channel, but when the backscatter is enabled, it reflects the signal shifted from the central carrier to the oscillation frequency (f osc = 300 kHz in this experiment).   The modulated backscattered power at the sideband is proportional to the differential radar cross-section RCS dif given by [23]: where λ is the wavelength, G a is the antenna gain, and m is a modulating factor that can be obtained from (4) as m = 4|c 1 | 2 /|∆Γ| 2 [27]. Assuming an ideal 50% duty cycle and a modulating square wave, m is 1/π 2 . The backscattered power can be maximized if the antenna resonates in the band (ImZ a = 0) and the two reflection coefficients are out of phase (e.g., close to short circuit and open circuit) with a phase difference of 180 • . In order to obtain the difference in the reflection coefficients (assuming that the antenna is well-matched to the reference impedance 50 Ω) between the two states ∆Γ, a prototype consisting of an SMA connector connected to the ADG902 switch is measured with the Vector Network Analyzer (VNA). Figure 4 shows the measured reflection coefficient when the switch is on the OFF state, which is loaded with a short circuit, and for the ON state, which is loaded with an open circuit.
A marker at 406 MHz shows that the two states have a phase difference close to 180 • . A loss of 1.5 dB with respect to the ideal short circuit and open circuit case is obtained due to the losses introduced by the switch. The modulated backscattered power at the sideband is proportional to the differential radar cross-section RCSdif given by [23]: where λ is the wavelength, Ga is the antenna gain, and m is a modulating factor that can be obtained from (4) as m = 4|c1| 2 /|ΔΓ| 2 [27]. Assuming an ideal 50% duty cycle and a modulating square wave, The backscattered power can be maximized if the antenna resonates in the band (ImZa = 0) and the two reflection coefficients are out of phase (e.g., close to short circuit and open circuit) with a phase difference of 180°. In order to obtain the difference in the reflection coefficients (assuming that the antenna is well-matched to the reference impedance 50 Ω) between the two states ΔΓ, a prototype consisting of an SMA connector connected to the ADG902 switch is measured with the Vector Network Analyzer (VNA). Figure 4 shows the measured reflection coefficient when the switch is on the OFF state, which is loaded with a short circuit, and for the ON state, which is loaded with an open circuit. A marker at 406 MHz shows that the two states have a phase difference close to 180°. A loss of 1.5 dB with respect to the ideal short circuit and open circuit case is obtained due to the losses introduced by the switch.

Data Transmission
The listen-before-talk approach is used to select a free channel. The LoRa gateway is configured to listen on the LoRa shifted channel. The procedure for data transmission is schematically shown in Figure 5. The protocol is inspired in the well-known asynchronous serial port protocol. The LoRa transmitter node is continuously sending packets, the duration of which (packet air time) depends on the spread factor selected and the number of bits per packet (see Section 2.5). For example, for a

Data Transmission
The listen-before-talk approach is used to select a free channel. The LoRa gateway is configured to listen on the LoRa shifted channel. The procedure for data transmission is schematically shown in Figure 5. The protocol is inspired in the well-known asynchronous serial port protocol. The LoRa transmitter node is continuously sending packets, the duration of which (packet air time) depends on the spread factor selected and the number of bits per packet (see Section 2.5). For example, for a packet with SF = 12 and 10 bits containing only the preamble data, its duration is 1 s (considering a small tolerance). The data are encoded in frames. Each frame consists of a start bit and N data bits with the sensor information and ends with a Stop bit. The backscatter reflects the packets enabling the internal oscillator. When the receiver in the gateway receives a packet corresponding to the start bit, it starts to decode the message after an interval between frames. Note that it does not need to receive a complete packet, it is only necessary to detect the presence of a packet. Therefore, it is enough to receive part of the preamble section of the packet. This feature simplifies the synchronization. The LoRa receiver can measure the frequency shift between the expected central frequency and the central frequency measured by the receiver. Therefore, small frequency offsets due to drift in the backscatter oscillator can be dynamically corrected between frames. The frame time is determined from the number of bits per frame (start, payload, and stop bits) multiplied by the interval between packets. Therefore, the frame time depends on the number of bits that will be sent in the application. It can be fixed or variable. In this last case, it is recommended to send some field indicating the total length of the frame. Additional protection such as parity bits, checksums, or CRC bits can be encoded within the data field. The implanted backscatter is assumed to update data cyclically. After the frame transmission, the backscatter changes to sleep mode to reduce power consumption. The backscatter initiates a new frame transmission after a period that is longer than the minimum interval between frames that must be a minimum of a bit duration. However, longer intervals between frames or transmission bursts help to increase the battery life. to drift in the backscatter oscillator can be dynamically corrected between frames. The frame time is determined from the number of bits per frame (start, payload, and stop bits) multiplied by the interval between packets. Therefore, the frame time depends on the number of bits that will be sent in the application. It can be fixed or variable. In this last case, it is recommended to send some field indicating the total length of the frame. Additional protection such as parity bits, checksums, or CRC bits can be encoded within the data field. The implanted backscatter is assumed to update data cyclically. After the frame transmission, the backscatter changes to sleep mode to reduce power consumption. The backscatter initiates a new frame transmission after a period that is longer than the minimum interval between frames that must be a minimum of a bit duration. However, longer intervals between frames or transmission bursts help to increase the battery life.

Phantom Design
In order to develop an antenna for a deep implant, a phantom material needs to be designed to simulate the muscle tissue for the desired frequency range. The dielectric properties of different tissues have been reported in the literature [28]. To this end, several phantom materials have been proposed for different frequency ranges and applications [29][30][31][32]. In our case, the dielectric properties must be designed to take into account the shift in the antenna frequency response introduced by the high dielectric constant of the tissues and to simulate the attenuation of the body and the reduction in the radiation efficiency of the antenna due to the high tissue conductivity. In our case, we are interested in finding a simple phantom for the UHF band (400-433 MHz). The dielectric constant of the muscle is about 57 and the conductivity is 0.8 at 406 MHz according to data provided by Gabriel et al. [28]. For frequencies higher than 100 MHz, the dielectric constant of the muscle is lower than that of water [30] so a phantom based on a saline water solution is considered because antennas can be used in it more easily than in a solid phantom and it is easy to manufacture. In [31], a mixture of

Phantom Design
In order to develop an antenna for a deep implant, a phantom material needs to be designed to simulate the muscle tissue for the desired frequency range. The dielectric properties of different tissues have been reported in the literature [28]. To this end, several phantom materials have been proposed for different frequency ranges and applications [29][30][31][32]. In our case, the dielectric properties must be designed to take into account the shift in the antenna frequency response introduced by the high dielectric constant of the tissues and to simulate the attenuation of the body and the reduction in the radiation efficiency of the antenna due to the high tissue conductivity. In our case, we are Sensors 2020, 20, 6342 7 of 26 interested in finding a simple phantom for the UHF band (400-433 MHz). The dielectric constant of the muscle is about 57 and the conductivity is 0.8 at 406 MHz according to data provided by Gabriel et al. [28]. For frequencies higher than 100 MHz, the dielectric constant of the muscle is lower than that of water [30] so a phantom based on a saline water solution is considered because antennas can be used in it more easily than in a solid phantom and it is easy to manufacture. In [31], a mixture of water with NaCl and sucrose was proposed to design a phantom at 2.45 GHz. The present study follows a similar procedure but adjusts the parameters to the 400-500 MHz band. The permittivity was measured using an open-ended coaxial probe [33] from Keysight (slim form probe 85070E/N1501A model, Keysight Technologies, Santa Rosa, CA, USA) with a Vector Network Analyzer (VNA) model Keysight PNA E8364C. Figure 6a shows the real part and imaginary part of the permittivity as a function of the mass fraction (percentage by weight) of sugar at 406 MHz. The imaginary part remains nearly constant whereas the real part decreases as a function of sugar. It can be shown that for a percentage of about 45% of sugar, the real part is close to the muscle value. Figure 6b shows the complex permittivity for a saline water solution with 45% of sugar and changing the quantity of NaCl at 406 MHz. Therefore, adding NaCl can increase the conductivity and adding sugar can reduce the real part. A solution consisting of 45% sucrose, 0.8% NaCl, and water (percentages are in weight) was chosen. Figure 7 compares the data reported by Gabriel et al. [28] and our measurements at ambient temperature (25 • C). Agreement was good in the 400-600 MHz band. The power transported by the wave decreases exponentially in accordance with Lambert's law. The attenuation constant introduced by the tissue can be computed using [34]: where ε c = ε − jε" is the complex relative permittivity, ε 0 and µ 0 are the vacuum permittivity and permeability, respectively, and f is the frequency. Considering the muscle complex permittivity, a typical attenuation of 165 dB/m or a power penetration depth of d p = 1/(2α(Np/m)) = 2.6 cm is introduced by the muscle at 406 MHz. Equation (7) was derived for incident plane waves in large bodies without considering reflections within the materials at the opposite interface with the air. Ayappa et al. [34] have shown that the above relationship only applies if the thickness of the tissues is larger than the critical length, which is about 5.4d p . Otherwise, computer simulation models are required to predict more realistic power distributions.
Sensors 2020, 20, x FOR PEER REVIEW 7 of 26 measurements at ambient temperature (25 °C). Agreement was good in the 400-600 MHz band. The power transported by the wave decreases exponentially in accordance with Lambert's law. The attenuation constant introduced by the tissue can be computed using [34]: where εc = ε' − jε'' is the complex relative permittivity, ε0 and µ0 are the vacuum permittivity and permeability, respectively, and f is the frequency. Considering the muscle complex permittivity, a typical attenuation of 165 dB/m or a power penetration depth of dp = 1/(2α(Np/m)) = 2.6 cm is introduced by the muscle at 406 MHz. Equation (7) was derived for incident plane waves in large bodies without considering reflections within the materials at the opposite interface with the air. Ayappa et al. [34] have shown that the above relationship only applies if the thickness of the tissues is larger than the critical length, which is about 5.4dp. Otherwise, computer simulation models are required to predict more realistic power distributions.

Antenna for Deeply Implanted Backscatter
Several implanted antennas have been proposed in the literature for ISM bands and MICS [35][36][37]. A ceramic (alumina) antenna with a superstrate was proposed in [35]. Chip antennas made with low temperature co-fired ceramics (LTCC) are available on the market and are extensively used in wireless modules (ISM band, Bluetooth, GPS, WLAN, etc.). The high permittivity of these materials enables the antennas to be smaller so that they can be integrated into small circuits. In the present study, a miniature commercial LTCC ceramic antenna (0900AT43A0070, Johanson Technology, Camarillo, CA, USA) designed for operation in the air at 900 MHz is used to operate as an implant in the 402-433 MHz band. The antenna is 7.7 mm long, 2 mm wide, and 0.8 mm high. Simulations and experiments were conducted to study the feasibility of commercial ceramic antennas for devices implanted in MICS and the 433 MHz band. The small size of the antenna (0.0715λ in the body) means that it is electrically small compared with the wavelength, so neither the pattern diagram nor the directivity will be greatly affected, and it has an omni-directional radiation pattern. Due to the high permittivity surrounding the antenna, the effective electrical length of the monopole decreases and changes the antenna impedance. This drawback can be overcome with a matching network that has discrete components (inductors and capacitors) as the manufacturer proposes for air but tuned to the impedance of the implanted antenna. However, the main drawback is that the antenna radiation efficiency is reduced by the lossy tissue material close to the antenna. This reduction in radiation efficiency, in turn, reduces the antenna gain and from Equation (6), we determine that the RCSdif and the backscattering power are also lessened. Implanted medical devices must be covered by a biocompatible material to prevent contamination [38] and metallic oxidation, and the possibility of a short-circuit due to the high conductivity of the tissues [39]. The coating material acts as a superstrate and the thickness and the material have a considerable influence on the radiation efficiency [35]. Increasing the thickness of this coating material or the air gap between the cover and the antenna improves the radiation efficiency [39] but also increases the size of the implanted device. Therefore, there is a constraint between radiation efficiency and size. In the present study, we have used rubber silicone (εr = 2.9, tan δ = 0.02) [40]. Figure 8 shows the reflection coefficient of the antenna connected to a short transmission line (width 1.85 mm) printed on 32 mil Roger RO4003C substrate (ε r = 3.55, tan δ = 0.0027) within the body with a 1 mm thick coating of silicone. The size of the ground plane was 10 mm × 10 mm. An image is shown in Figure 8a. This figure shows that the antenna is completely mismatched in the 400-433 MHz band. Therefore, a matching network is required. Assuming a wellmatched design, the realized antenna gain must be close to the antenna gain that does not take into

Antenna for Deeply Implanted Backscatter
Several implanted antennas have been proposed in the literature for ISM bands and MICS [35][36][37]. A ceramic (alumina) antenna with a superstrate was proposed in [35]. Chip antennas made with low temperature co-fired ceramics (LTCC) are available on the market and are extensively used in wireless modules (ISM band, Bluetooth, GPS, WLAN, etc.). The high permittivity of these materials enables the antennas to be smaller so that they can be integrated into small circuits. In the present study, a miniature commercial LTCC ceramic antenna (0900AT43A0070, Johanson Technology, Camarillo, CA, USA) designed for operation in the air at 900 MHz is used to operate as an implant in the 402-433 MHz band. The antenna is 7.7 mm long, 2 mm wide, and 0.8 mm high. Simulations and experiments were conducted to study the feasibility of commercial ceramic antennas for devices implanted in MICS and the 433 MHz band. The small size of the antenna (0.0715λ in the body) means that it is electrically small compared with the wavelength, so neither the pattern diagram nor the directivity will be greatly affected, and it has an omni-directional radiation pattern. Due to the high permittivity surrounding the antenna, the effective electrical length of the monopole decreases and changes the antenna impedance. This drawback can be overcome with a matching network that has discrete components (inductors and capacitors) as the manufacturer proposes for air but tuned to the impedance of the implanted antenna. However, the main drawback is that the antenna radiation efficiency is reduced by the lossy tissue material close to the antenna. This reduction in radiation efficiency, in turn, reduces the antenna gain and from Equation (6), we determine that the RCS dif and the backscattering power are also lessened. Implanted medical devices must be covered by a biocompatible material to prevent contamination [38] and metallic oxidation, and the possibility of a short-circuit due to the high conductivity of the tissues [39]. The coating material acts as a superstrate and the thickness and the material have a considerable influence on the radiation efficiency [35]. Increasing the thickness of this coating material or the air gap between the cover and the antenna improves the radiation efficiency [39] but also increases the size of the implanted device. Therefore, there is a constraint between radiation efficiency and size. In the present study, we have used rubber silicone (ε r = 2.9, tan δ = 0.02) [40]. Figure 8 shows the reflection coefficient of the antenna connected to a short transmission line (width 1.85 mm) printed on 32 mil Roger RO4003C substrate (ε r = 3.55, tan δ = 0.0027) within the body with a 1 mm thick coating of silicone. The size of the ground plane was 10 mm × 10 mm. An image is shown in Figure 8a. This figure shows that the antenna is completely mismatched in the 400-433 MHz band. Therefore, a matching network is required. Assuming a well-matched design, the realized antenna Sensors 2020, 20, 6342 9 of 26 gain must be close to the antenna gain that does not take into account the antenna mismatch. Figure 9 shows the gain, directivity, and radiation efficiency at 406 MHz as a function of the thickness of the silicone coating. The simulations were performed with the full-waveAnsys high-frequency structure simulator (HFSS) simulator including the model provided by the antenna manufacturer. We used a box of muscle tissue (ε r = 57, σ = 0.8 S/m) surrounding the antenna to simulate the effect of the body. The antenna directivity is not hardly affected by the thickness of the coating, but the radiation efficiency, and therefore the gain, increases with thickness. The efficiency fluctuates between 9.5% and 16%, which is considerably lower than in air (about 27%). Figure 10 shows the effect of the coating permittivity. This figure assumes the coating material to be 1 mm thick. The radiation efficiency increases slightly as the permittivity of the coating increases. However, it should be pointed out that coatings based on ceramic substrates with high permittivity are difficult to drill and cut.
Sensors 2020, 20, x FOR PEER REVIEW 9 of 26 account the antenna mismatch. Figure 9 shows the gain, directivity, and radiation efficiency at 406 MHz as a function of the thickness of the silicone coating. The simulations were performed with the full-waveAnsys high-frequency structure simulator (HFSS) simulator including the model provided by the antenna manufacturer. We used a box of muscle tissue (εr = 57, σ = 0.8 S/m) surrounding the antenna to simulate the effect of the body. The antenna directivity is not hardly affected by the thickness of the coating, but the radiation efficiency, and therefore the gain, increases with thickness. The efficiency fluctuates between 9.5% and 16%, which is considerably lower than in air (about 27%). Figure 10 shows the effect of the coating permittivity. This figure assumes the coating material to be 1 mm thick. The radiation efficiency increases slightly as the permittivity of the coating increases. However, it should be pointed out that coatings based on ceramic substrates with high permittivity are difficult to drill and cut.    In order to design the matching network with real materials and dimensions, a prototype with an SMA connector was manufactured (Figure 11a). The size of the board (7.5 mm × 21 mm) and ground plane (7.5 mm × 12 mm) together with the position of the discrete matching components is the same as in the implanted prototype ( Figure 2). The impedance of the antenna was measured with a VNA connecting a thru line (zero Ohm SMD resistance), after which the reference plane was shifted to obtain the antenna impedance. Then, an LC matching network was designed using the Smith chart. Figure 11b shows the measured reflection coefficient of the antenna in the phantom material before and after the matching network. The frequency band can be tuned by modifying the components of the matching network. A 56 nH inductance in series and a 0.5 pF in parallel were used to match the antenna at 406 MHz and there was less than −10 dB between 387.28 MHz and 469.3 MHz (82.1 MHz bandwidth or 19.2%). Table 1 compares the designed antenna with other printed implanted antennas proposed in the literature in the MICS band ( [41][42][43][44]) using different miniaturization techniques. The miniature antenna proposed in the present study achieves a high gain and radiation efficiency, and a bandwidth that is moderate but enough to cover the MICS band. It is also available commercially, which simplifies the commercialization of future implanted devices.

Link Budget
In this section, a link budget will be presented to calculate the maximum depth and coverage area at which the implanted backscatter can be read as a function of the system parameters. The received power can be calculated using the Friis transmission equation with the following parameters: the transmission power (PTx), the transmission antenna gain GTx and receiver antenna gain GRx, the wavelength λ, the differential RCS given by (6), the distance between the transmitter and

Link Budget
In this section, a link budget will be presented to calculate the maximum depth and coverage area at which the implanted backscatter can be read as a function of the system parameters. The received power can be calculated using the Friis transmission equation with the following parameters: the transmission power (P Tx ), the transmission antenna gain G Tx and receiver antenna gain G Rx , the wavelength λ, Sensors 2020, 20, 6342 12 of 26 the differential RCS given by (6), the distance between the transmitter and the body (d Tx ) and between the body and receiver (d Rx ), the body attenuation α(dB/m), and the implant depth (d): The parameters used in the simulations are listed in Table 2. The LoRa receiver sensitivity is a function of the noise figure of the receiver and the bandwidth BW. It is computed from the noise floor plus the required signal to noise ratio (SNR), which is a function of the Spreading Factor as shown in Table 3 [6,45]: S(dBm) = −174 + 10 log(BW) + NF(dB) + SNR(dB)   Table 3 assumes a typical BW of 125 kHz. Thus, sensitivity decreases if lower bandwidth is used in accordance with (9) (12 dB for the lower bandwidth compared to 125 kHz). Table 3 also shows the time on air and the equivalent bit rate computed from the equations given in [6,45] considering 10 bytes of payload with a preamble of 6 bytes and a coding rate of CR = 4/5.
However, due to multipath propagation and reflections on the ground and objects, in indoor environments, the power decays with exponential factors that are not those given in the Friis model considered in (8). The received power or RSSI can be expressed using the following empirical model, which is an extension to the backscatter channel of the empirical propagation model for indoor environments [41,46]: Sensors 2020, 20, 6342

of 26
In this model, the received power is modeled as an average term and a random term X, where d 0 is a reference distance (e.g., the midpoint between transmitter and receiver d 0 = (d Tx-Rx /2)) and P R (d 0 , d 0 ) is the average received power at this reference distance. The exponential decay factors n 1 and n 2 are a function of the environment and the antenna height but for indoor environments, they have values of about 2.5-3 [47,48]. From the Friis equation, P R (d 0 , d 0 ) is a function of transmitted power and the antenna gains. However, it is also a function of the antenna height and the diffraction. P R (d 0 , d 0 ) can be found experimentally by regressing the measured data. Nevertheless, in this section, we perform simulations by considering the following expression derived from the Friis transmission Equation (8): L obs is a term in dB that accounts for the attenuation for diffraction that depends on the antenna height [46]. The losses for ground diffraction can be neglected if the antenna height is higher than 0.6R 1 , where R 1 is the radius of the first Fresnel zone [46]: For a distance between the transmitter and receiver (backscatter) of 2 m, and antennas of the same height, the minimum antenna height is 0.36 m at 406 MHz. Therefore, for a typical antenna height of around 1-2 m, we can expect that the effect of ground diffraction will be small.
The random variable X takes into account the attenuation due to multipath propagation. It can be regarded as a log-normally distributed random variable X(dB)~N(0,σ), where σ in dB is the standard deviation and experimental results have given values between 2 and 4 dB for a single path [48] (in our case we expect it to be double this because of the backscatter channel).
The coverage probability (Prob) is the probability that the average received power or RSSI, P R (dBm), will be higher than the receiver sensitivity S. For a log-normal distributed channel, it is given by [46]: where erf is the Gauss error function. For σ = 8 dB and received power of 10 dB over the sensitivity (fading margin), the probability values are about 90%. Two cases were analyzed: monostatic and bistatic. In the monostatic case, the LoRa transmitter and receiver gateway are close together (d Tx = d Rx ). Figure 12 compares the maximum depth or distance inside the body that the backscatter can communicate for SF = 7 and SF = 12 with BW = 125 kHz. Figure 13 shows the maximum depth as a function of the transmitter to body distance and receiver to body distance in the bistatic case for SF = 7 and SF = 12 with BW = 125 kHz, respectively. In these simulations, the exponential decay factors (n 1 and n 2 ) were considered to be equal to 2.5. Depths can be greater than 10 cm at 4 m. In the bistatic case, depths can be greater if the transmitter is closer to the body than the receiver. These simulations have taken into account a fading margin due to multipath propagation of 10 dB. Results show the potential to read low rate data from implants using frequency-shifted LoRa backscatter. The number of packets per unit time increases if the spreading factor is low (e.g., SF = 7) as can be seen in Table 3 because the time on air is shorter than at high spreading factors. However, the sensitivity and the robustness of the system were better when spreading factors were high. In the bistatic case, to achieve better coverage, it is assumed that the transmitter will be located near the patient (e.g., close to the bed) and the receiver somewhere else in the room. Sensors 2020, 20, x FOR PEER REVIEW 14 of 26 Figure 12. Simulated maximum depth for an implanted backscatter as a function of the distance to the body using the parameters of Table 1 for the monostatic case with a spreading factor SF = 7 and 12. Propagation model parameters: n1 = n2 = 2.5, Fade Margin = 10 dB.
(a) (b) Figure 13. Simulated maximum depth for an implanted backscatter as a function of the distance of the transmitter and receiver from the body using the parameters of Table 2 for the bistatic case with a spreading factor SF = 7 (a) and SF = 12 (b). Propagation model parameters: n1 = n2 = 2.5, Fade Margin = 10 dB.
A simplified scenario shown in Figure 14 consists of the transmitter and the receiver and the implanted backscatter on the same plane. All these elements are used to estimate the coverage from the propagation model proposed. In this scenario, it is assumed that the transceivers are fixed and the backscatter can move around different positions. Figures 15 and 16 show the simulated coverage results for backscatter at a depth (d) of 5 cm and 10 cm, respectively. They depict the average received power computed with (10)-(11) considering the ideal free space exponential factors n1 = n2 = 2 and n1 = n2 = 2.5 as a function of the position of the backscatter. The monostatic and biostatic case (with a distance between the transmitter and receiver dTx-Rx = 1m) are shown. For devices implanted on the surface or at small depths (<5 cm), the read range area can be about 4 m × 4 m, whereas for deep implanted devices (10 cm), a read range area of about 1.8 m × 1.8 m can be covered.  Table 1 for the monostatic case with a spreading factor SF = 7 and 12. Propagation model parameters: n1 = n2 = 2.5, Fade Margin = 10 dB.
(a) (b) Figure 13. Simulated maximum depth for an implanted backscatter as a function of the distance of the transmitter and receiver from the body using the parameters of Table 2 for the bistatic case with a spreading factor SF = 7 (a) and SF = 12 (b). Propagation model parameters: n1 = n2 = 2.5, Fade Margin = 10 dB.
A simplified scenario shown in Figure 14 consists of the transmitter and the receiver and the implanted backscatter on the same plane. All these elements are used to estimate the coverage from the propagation model proposed. In this scenario, it is assumed that the transceivers are fixed and the backscatter can move around different positions. Figures 15 and 16 show the simulated coverage results for backscatter at a depth (d) of 5 cm and 10 cm, respectively. They depict the average received power computed with (10)-(11) considering the ideal free space exponential factors n1 = n2 = 2 and n1 = n2 = 2.5 as a function of the position of the backscatter. The monostatic and biostatic case (with a distance between the transmitter and receiver dTx-Rx = 1m) are shown. For devices implanted on the surface or at small depths (<5 cm), the read range area can be about 4 m × 4 m, whereas for deep implanted devices (10 cm), a read range area of about 1.8 m × 1.8 m can be covered.  Table 2 for the bistatic case with a spreading factor SF = 7 (a) and SF = 12 (b). Propagation model parameters: n 1 = n 2 = 2.5, Fade Margin = 10 dB.
A simplified scenario shown in Figure 14 consists of the transmitter and the receiver and the implanted backscatter on the same plane. All these elements are used to estimate the coverage from the propagation model proposed. In this scenario, it is assumed that the transceivers are fixed and the backscatter can move around different positions. Figures 15 and 16 show the simulated coverage results for backscatter at a depth (d) of 5 cm and 10 cm, respectively. They depict the average received power computed with (10)-(11) considering the ideal free space exponential factors n 1 = n 2 = 2 and n 1 = n 2 = 2.5 as a function of the position of the backscatter. The monostatic and biostatic case (with a distance between the transmitter and receiver

Results
A backscatter prototype has been designed as shown in the scheme in Figure 1. The prototype includes an antenna with a matching network, an RF switch (ADG902) used as a modulator, and a low-power oscillator (LTC6907). The system is controlled by a microcontroller (ATTiny85) [49] that can generate test messages. The backscatter is covered with 1 mm of heat shrink silicone (approximately 1 mm thick). This proof of concept prototype is externally biased with a battery pack connected by two wires that are also used to support the device. A photograph is shown in Figure  17. In order to check the correct modulation of the backscatter inside the body, the backscatter was inserted into a 20 cm diameter cylindrical vessel (see Figure 17). The vessel was filled with a mixture of saline solution with sugar described in the previous section. A signal generator was connected to a transmitting antenna and a spectrum analyzer was used as a receiver connected to a receiving antenna. Two whip monopoles tuned to 406 MHz were used in the generator and the spectrum analyzer. The distance from the Tx and RX to the vessel was 1 m. The frequency was swept and the backscattered power was received. Figure 18 shows the spectrum in which the sidebands at the modulating backscatter frequency at 406 MHz can be seen. The central peak at this carrier frequency was due to the coupling between the transmitter and receiver. The sideband peaks are spaced from the carrier at the same frequency at which the ADG902 is switching (300 kHz in this experiment). Figure 19 shows the measurement of the received power at the sideband as a function of the generator frequency. The maximum received power was around the 408 MHz band. It is observed that the backscatter covers the desired frequency band, which shows that the antenna was correctly tuned within the phantom.

Results
A backscatter prototype has been designed as shown in the scheme in Figure 1. The prototype includes an antenna with a matching network, an RF switch (ADG902) used as a modulator, and a low-power oscillator (LTC6907). The system is controlled by a microcontroller (ATTiny85) [49] that can generate test messages. The backscatter is covered with 1 mm of heat shrink silicone (approximately 1 mm thick). This proof of concept prototype is externally biased with a battery pack connected by two wires that are also used to support the device. A photograph is shown in Figure 17. In order to check the correct modulation of the backscatter inside the body, the backscatter was inserted into a 20 cm diameter cylindrical vessel (see Figure 17). The vessel was filled with a mixture of saline solution with sugar described in the previous section. A signal generator was connected to a transmitting antenna and a spectrum analyzer was used as a receiver connected to a receiving antenna. Two whip monopoles tuned to 406 MHz were used in the generator and the spectrum analyzer. The distance from the Tx and RX to the vessel was 1 m. The frequency was swept and the backscattered power was received. Figure 18 shows the spectrum in which the sidebands at the modulating backscatter frequency at 406 MHz can be seen. The central peak at this carrier frequency was due to the coupling between the transmitter and receiver. The sideband peaks are spaced from the carrier at the same frequency at which the ADG902 is switching (300 kHz in this experiment). Figure 19 shows the measurement of the received power at the sideband as a function of the generator frequency. The maximum received power was around the 408 MHz band. It is observed that the backscatter covers the desired frequency band, which shows that the antenna was correctly tuned within the phantom.   To test the reliability of the system, some experimental results were achieved using a LoRa transmitter (ESP32 TGNO with Semtech SX1278 transceiver) and another identical LoRa module used as a receiver. The central frequency of the transmitter was 406 MHz and the oscillation frequency was 300 kHz. The receiver was tuned to 406.3 MHz. The LoRa was configured to transmit and receive packets with SF = 12 and BW = 125 kHz for better sensitivity. Figure 20 shows an image of the setup used in the experiments.  Figure 21 shows the complementary distribution function (CDF) plot of the received signal strength indicator (RSSI) and the SNR in the shifted channel measured by the LoRa transceiver. In order to obtain the data, the receiver was tuned to the backscattered channel and was located at random points (uniformly distributed) in the laboratory within a coverage area of 4 × 4m. The distance from the LoRa transmitter to the vessel was 80 cm at a height of 1 m from the ground. The Measured Power(dBm) Figure 19.
First harmonic (f c + f m ) measured power as a function of the frequency with the backscatter immersed in the vessel described in Figure 17.
To test the reliability of the system, some experimental results were achieved using a LoRa transmitter (ESP32 TGNO with Semtech SX1278 transceiver) and another identical LoRa module used as a receiver. The central frequency of the transmitter was 406 MHz and the oscillation frequency was 300 kHz. The receiver was tuned to 406.3 MHz. The LoRa was configured to transmit and receive packets with SF = 12 and BW = 125 kHz for better sensitivity. Figure 20 shows an image of the setup used in the experiments.
Sensors 2020, 20, x FOR PEER REVIEW 18 of 26 Figure 19. First harmonic (fc + fm) measured power as a function of the frequency with the backscatter immersed in the vessel described in Figure 17.
To test the reliability of the system, some experimental results were achieved using a LoRa transmitter (ESP32 TGNO with Semtech SX1278 transceiver) and another identical LoRa module used as a receiver. The central frequency of the transmitter was 406 MHz and the oscillation frequency was 300 kHz. The receiver was tuned to 406.3 MHz. The LoRa was configured to transmit and receive packets with SF = 12 and BW = 125 kHz for better sensitivity. Figure 20 shows an image of the setup used in the experiments.  Figure 21 shows the complementary distribution function (CDF) plot of the received signal strength indicator (RSSI) and the SNR in the shifted channel measured by the LoRa transceiver. In order to obtain the data, the receiver was tuned to the backscattered channel and was located at random points (uniformly distributed) in the laboratory within a coverage area of 4 × 4m. The distance from the LoRa transmitter to the vessel was 80 cm at a height of 1 m from the ground. The  Figure 21 shows the complementary distribution function (CDF) plot of the received signal strength indicator (RSSI) and the SNR in the shifted channel measured by the LoRa transceiver. In order to obtain the data, the receiver was tuned to the backscattered channel and was located at random points (uniformly distributed) in the laboratory within a coverage area of 4 × 4m. The distance from the LoRa transmitter to the vessel was 80 cm at a height of 1 m from the ground. The backscatter was immersed in the center of a 20 cm diameter vessel with the phantom solution. We can observe that the RSSI is typically higher than −123 dBm (less than 10% of the locations receive levels below −123 dBm and SNR below −15 dB). Therefore, most locations were above the sensitivity of the receiver, and the packets were successfully backscattered.
Sensors 2020, 20, x FOR PEER REVIEW 19 of 26 backscatter was immersed in the center of a 20 cm diameter vessel with the phantom solution. We can observe that the RSSI is typically higher than −123 dBm (less than 10% of the locations receive levels below −123 dBm and SNR below −15 dB). Therefore, most locations were above the sensitivity of the receiver, and the packets were successfully backscattered. In order to show the coverage area, Figure 22 shows the heat maps of the measured RSSI for the transmitter located 80 cm from the vessel with the immersed backscatter prototype centered and immersed at a depth of 10 cm The backscatter prototype and the transmitter are fixed whereas the receiver is moved along different locations in the room (laboratory). Under these circumstances, data can be received within the laboratory (area 5 m × 10 m) and through walls provided the receiver is located in neighboring rooms on the same floor. This result is in accordance with Figure 13b which shows that devices implanted at a 10 cm depth can be detected up to 10 m from the receiver when the transmitter is closed (80 cm in this case). In order to show the coverage area, Figure 22 shows the heat maps of the measured RSSI for the transmitter located 80 cm from the vessel with the immersed backscatter prototype centered and immersed at a depth of 10 cm The backscatter prototype and the transmitter are fixed whereas the receiver is moved along different locations in the room (laboratory). Under these circumstances, data can be received within the laboratory (area 5 m × 10 m) and through walls provided the receiver is located in neighboring rooms on the same floor. This result is in accordance with Figure 13b which shows that devices implanted at a 10 cm depth can be detected up to 10 m from the receiver when the transmitter is closed (80 cm in this case). To emulate the functionality of an implanted sensor, the measurement of the temperature with an external sensor is transmitted using the encoding procedure described above. Figure 23 shows the received data obtained with frame transmission over 60 min. The transmitter was located 80 cm from the vessel surface with the backscatter immersed in the phantom liquid at the center. The receiver was located at a distance of 1.5 m. In this test, the better range SF configuration [3] was set up (SF = 12), keeping the bandwidth (BW) to 125 kHz to not increase drastically the time on air of the packet. To test the backscatter, a straightforward protocol (described in Section 2.2) to send a byte in a 10 bits frame (start bit plus 8 data bits and end of frame bit) was implemented. The average number of frames lost was 4%. These results pave the way for long-range wireless implanted devices based on backscattering communication in indoor environments.  To emulate the functionality of an implanted sensor, the measurement of the temperature with an external sensor is transmitted using the encoding procedure described above. Figure 23 shows the received data obtained with frame transmission over 60 min. The transmitter was located 80 cm from the vessel surface with the backscatter immersed in the phantom liquid at the center. The receiver was located at a distance of 1.5 m. In this test, the better range SF configuration [3] was set up (SF = 12), keeping the bandwidth (BW) to 125 kHz to not increase drastically the time on air of the packet. To test the backscatter, a straightforward protocol (described in Section 2.2) to send a byte in a 10 bits frame (start bit plus 8 data bits and end of frame bit) was implemented. The average number of frames lost was 4%. These results pave the way for long-range wireless implanted devices based on backscattering communication in indoor environments. To emulate the functionality of an implanted sensor, the measurement of the temperature with an external sensor is transmitted using the encoding procedure described above. Figure 23 shows the received data obtained with frame transmission over 60 min. The transmitter was located 80 cm from the vessel surface with the backscatter immersed in the phantom liquid at the center. The receiver was located at a distance of 1.5 m. In this test, the better range SF configuration [3] was set up (SF = 12), keeping the bandwidth (BW) to 125 kHz to not increase drastically the time on air of the packet. To test the backscatter, a straightforward protocol (described in Section 2.2) to send a byte in a 10 bits frame (start bit plus 8 data bits and end of frame bit) was implemented. The average number of frames lost was 4%. These results pave the way for long-range wireless implanted devices based on backscattering communication in indoor environments.

Discussion
The use of a frequency band (406-433 MHz) allows to increase the coverage because the propagation losses are lower than at 868/915 MHz or 2.45 GHz ISM band (free space attenuation decays with the square of the frequency) and besides, the walls [50] and indoor propagation [48] introduce lower attenuation at this frequency band. This helps to counter the fading effect caused by multipath propagation in indoor environments which, in backscattering systems, affects the up and downlinks [51][52][53]. The proposed system can easily use receiver diversity with multiple receivers to sniff shifted channels and prevent selective multipath fading. Channel hopping is not expected to provide significant improvements in the 406 MICS and 433 MHz ISM bands due to the small spacing between channels (less than the channel coherence bandwidth) so the channels are correlated [47,48]. However, to prevent interference or collision with another system, a listen-before-talk procedure (LBT) can easily be implemented in the LoRa transmitter so that the available channels can be listened to before the packets are sent. The parameters of the channel (frequency, SF, BW) can be transmitted to the receivers using a conventional LoRa link. This is different from the 2.4 GHz band, where channel hopping can be used to counter the effects of multipath propagation (e.g., in Bluetooth) [48].
This section compares the technology presented with other existing wireless technologies for implanted applications. Table 4 summarizes some examples of wireless technologies and compares such parameters as the carrier frequency, data rate, implant depth, or read range. These technologies can be classified according to whether are battery-less or not. Each technology has advantages and drawbacks depending on the requirements of the final application. Due to the higher losses of the body, Wireless Power Transfer (WPT) based on inductive links at low-frequency (LF) [54] and high-frequency (HF) [55,56] is the method preferred to power battery-less devices. Data are transferred by load modulation (Load Shift Keying-LSK) or by another wireless link at higher frequencies [57]. Depth and power depend on the size of the coils used [55], and the communication (read range) can be performed when the coils are close to the body. Another drawback is that it requires specialized readers and ASICs. In [58], the authors propose a method based on three coils to improve the depth and reading range using smartphones with NFC and commercial NFC IC with energy-harvesting capabilities.
When the carrier frequency increases, the communication uses far-field and the attenuation increases very quickly with the frequency. Passive UHF RFID systems for identification and sensor communication at small implanted depths (5-15 mm subcutaneous) have been proposed [59,60]. The communication is by backscattering and a commercial UHF reader can be used. The WPT link is limited by the sensitivity of the UHF RFID IC [60] (about −20 dBm for modern ICs such as Impinj Monza R6 and −10 dBm for AMS-SL900A with sensing capabilities). The reader transmitter power is limited by exposure regulations (SAR) [60] so when body losses are higher the read range is limited to 0.6-1.6 m depending of the depth of the implant (5 to 15 mm) [59].
Battery implanted devices are used when high performance and high data rate transfer are required. Commercial active implantable transceivers are available on the market. For example, the Microsemi ZL70103 transceiver [61] operates in the MICS band and the 433 MHz ISM band at a raw data rate of up to 400 kbps and with wake-up at 3.45 GHz. Its Rx sensitivity is −102 dBm at 200 kbps (with a packet error rate of 10%) and the Tx/Rx average current is 5 mA. Several implanted wireless devices have been proposed for endoscopy [62][63][64][65][66]. For this application to download video data, a high data rate is needed; therefore, active transceivers are used. However, the battery lifetime is limited to a few hours, which is enough for the exploration period required. Data rates from 1 Mbps [66] to 80 Mbps [65] have been achieved using different modulations such as BPSK [65], FSK [62], QPSK [63,64], or OFDM [66], at frequency bands ranging from 20 MHz to 925 MHz transmitting power between 1 mW and 6 mW. The price to pay is high power consumption and, therefore, the limited lifetime. It is not the best option for long term monitoring even if the data rate required is small. Zero power backscatter gives better power consumption because it does not transmit RF power, which extends the life of the battery although the link attenuation is double that of active transceivers due to the forward and backscattered channel. The proposed backscatter can be used with other modulations such as ASK (OOK) or FSK and the bitrate achieved (a few MHz) is high enough to transfer implanted video camera data. In this case, a CW tone illuminates the backscatter from the transmitter and the backscatter responds by modulating the carrier. A backscatter similar to that described in this study was used in [9] to achieve 30 Mbps working at 915 MHz at a depth of 6 cm using a loop antenna. In another study [67], an FSK with a backscatter achieved 1-5 Mbps at 600 MHz and a depth of up to 10 cm. In both cases, the reader antennas were located close to the phantom and the reader system was implemented using software-defined radio (SDR). The LoRa backscatter technique proposed in this study can achieve greater depths and read ranges outside the body because of the higher sensitivity of the LoRa modulation, but the main limitation is the low data rate. Another important point is that the proposed system does not require specialized readers because commercial low-cost LoRa transceivers can be used.

Conclusions
This paper shows the feasibility of communication based on receiving backscatter LoRa packets at a shifted channel for implanted devices. A proof of concept prototype has been developed. The prototype consists of an RF switch that presents two high reflective coefficients to an antenna controlled by a low-frequency modulating oscillator. A commercial ceramic antenna is tuned using a matching network at the MICS and ISM 433 MHz bands, which provides high miniaturization and high radiation efficiency. The backscatter can be interfaced with a low power microprocessor to send low rate data from sensors for purposes of telemetry. Although the backscatter can be used with other modulations such as ASK and FSK, the proposed method uses the reflection of LoRa packets to exploit the high sensitivity and robustness of LoRa modulation; therefore, there is no need to integrate high power-consuming FPGAs or other ICs. A complete link budget has been presented which shows the potentiality of the communication technique. The method proposed achieves great depths at distances of up to several meters from the body so that only one transmitter base station is enough to cover a typical room.