Impact of the environment and the topology on the performance of hierarchical body area networks

Personal area networks and, more specifically, body area networks (BANs) are key building blocks of future generation networks and of the Internet of Things as well. In this article, we present a novel analytical framework for network performance analysis of body sensor networks with hierarchical (tree) topologies. This framework takes into account the specificities of the on-body channel modeling and the impact of the surrounding environment. The obtained results clearly highlight the differences between indoor and outdoor scenarios, and provide several insights on BAN design and analysis. In particular, it will be shown that the BAN topology should be selected according to the foreseen medical application and the deployment environment.


Introduction
Recent advances in ultra-low power sensors have fostered the research in the field of body-centric networks, also referred to as body area networks (BANs) [1][2][3][4].In these networks, a set of nodes (called sensors) is deployed on the human body.They aim at monitoring and reporting several physiological values, such as blood pressure, breath rate, skin temperature, or heart beating rate.Most of the time, sensing is performed at low rates but, in emergency situations, the network load may increase in seconds.Therefore, an in-depth analysis of the network outage, throughput, and achievable transmission rate can give insights on the maximum supported reporting rate and the corresponding performance.
In [5,6], we have considered a preliminary link-level performance analysis of BANs with centralized topologies.In the current study, we extend this approach, integrating the propagation channel characteristics and the impact of the hierarchy in a general network-level performance analysis framework.All considered networks will have hierarchical topologies, i.e., the sensor nodes will not be directly connected to a central controller.The modeling of the BAN channel has recently been thoroughly investigated [7][8][9][10][11].The main findings on the body radio propagation channel can be summarized as follows.First, the average value of the power decreases as an exponential function of the distance.However, unlike classical propagation models, where the received power is a decreasing function of the distance of the form d -a , the authors of [12,13] show that a law of the form 10 gd (g <0) characterizes more accurately body radio propagation.Second, the propagation channel is subject to distinct propagations mechanisms with respect to the location of the sensors on the body.More precisely, on-body propagation and reflections from the environment act jointly to create a particular propagation mechanism that is specific to BANs.
This article addresses the development of a specific framework for the accurate evaluation of the impact of the body-specific propagation and network topology.Our results are derived by means of the link throughput analysis, this metric being a traditional measure of how much traffic can be delivered, per time unit, by the network [14,15].Therefore, our analysis is expedient to understand the level of information which could be collected and processed in body-related applications (e.g., health or fitness monitoring).Furthermore, since energy is critical in the design of autonomous BANs in the context of medical applications [16][17][18], an accurate evaluation of the impacts of the BAN topology and transmission rate on the energy consumption is of fundamental interest.
The slotted ALOHA multiple access scheme [19] was recently proposed by the IEEE 802.15.6 working group as one of the reference medium access control (MAC) schemes for the wireless body networks in the context of the narrowband communications [20].In particular, in each time slot, the nodes are assumed to transmit independently with a certain fixed probability [21].This approach is supported by the observations in [ [22], p. 278] and [21,23], where it is shown that the traffic generated by nodes using a slotted random access MAC protocol can be modeled by means of a Bernoulli distribution.In fact, in more sophisticated MAC schemes, the probability of transmission at a node can be modeled as a function of general parameters, such as queuing statistics, the queue-dropping rate, the channel outage probability incurred by fading [24], the adaptation of the sampling to rate to patient's condition [25], the MAC strategy used [26], etc.Since the impact of these parameters is not the focus of the this study, the interested reader is referred to the existing literature [27][28][29] for further details.
The principal contributions of this article can be summarized as follows.First, a comprehensive and detailed analytic framework for BAN performance evaluation is developed, obtaining closed-form expressions for the link probabilities of outage in the context of multi-user communications.This framework encompasses the effect of the environment, the topology, and the traffic intensity.Next, different topologies, corresponding to various medical applications, are characterized in terms of achievable throughput.Finally, the performance of each topology is discussed, and practical insights are given on how to instantiate a real-life BAN with respect to the application demands and propagation context.Furthermore, throughout this entire article the indoor and the outdoor environments are treated separately and properly compared.
The remainder of this article is organized as follows.In Section 2, the propagation mechanisms are introduced and characterized.In Section 3, the conditional success probability of a link transmission for a node, given the transmitter-receiver and interferer-receiver distances, is derived.In the same section, the minimum required transmit power, over a given link, in the absence of any interfering node is computed in both indoor and outdoor environments.Then, in Section 4, the tree topologies analyzed in this article are presented, and the traffic model is discussed.Finally, in Section 5, an extensive performance analysis, in terms of network throughput and energy consumption, is performed.Section 6 concludes this article.

Propagation mechanisms
In order to build an accurate model for the on-body propagation, a Rohde & Schwartz ZVA-24 vector network analyzer was employed to capture the complex-valued frequency-domain transfer functions between 3 and 7 GHz, with a frequency step of 50 MHz.Omnidirectional Skycross SMT-3T010M ultra-wide band antennas were used during the entire measurement campaign.Their small-size (13.6 mm × 16 mm × 3 mm) and low profile characteristics precisely match the body sensor requirements.These antennas were separated from the body skin by about 5 mm to ensure a return loss value S 11 ≤ -9 dB.Finally, lowloss and phase-stable cables interconnect all components, and the IF-bandwidth was set to 100 Hz to enlarge the dynamic range to about 120 dB.
The experimental scenario is presented in Figure 1 and can be described as follows.The measurements were carried out around the 94 cm of the waist of a man (1m87, 83 kg) whose body is in a standing position, arms hanging along the side.The transmit antenna is placed around the body at a distance d from the receive antenna, which is located at the middle axis of the torso.

A. Measurements
First, the diffraction mechanism is analyzed by gradually shifting the transmitter around the body.The spatial values of the power are extracted from seven different sites separated by 4 cm each.For each level, the transmitter is also shifted one level below and one level above, and the observed measures are averaged.Second, the impact of the reflections off the surrounding environment was investigated for five positions of the transmitter around the body.Repeated measures are taken by positioning the human body on a rectangular grid of 7 × 7 position, each separated by 4 cm.This procedure is performed for a set of 20 locations in a standard office room with a surface of about 20 m 2 .
The baseband frequency response at the receiver was then converted into the delay domain using an inverse discrete Fourier transform [30].Next, a Hamming window was applied to reduce the side lobes up to -43 dB for the second lobe.The resulting complex impulse response allows a description of the BAN channel with a delay resolution of up to 0.25 ns.As shown in Figures 2 and 3, the different multipath and scattering mechanisms are well distinguished as a function of time.More precisely, the diffraction around the body is followed by the reflections off the environment.Both propagation mechanisms can be efficiently separated by applying a rectangular time gating at 7 ns.Finally, the narrowband power of the two distinct contributions mechanism is estimated by integrating the complex values of the temporal taps over each sub-channel.
The conclusions of this extensive measurement campaign, also highlighted in [13], can be summarized in three points.Firstly, there is propagation through the body.However, when high transmission frequencies are considered, the attenuation undergone by these waves is relevant and the corresponding contribution can be neglected.
A second mechanism corresponds to guided diffraction around the body.This mechanism is consistent with a surface wave propagation, and its properties depend on the body specific characteristics.
Finally, the last propagation contribution comes from the surrounding environment.More precisely, the third propagation mechanism originates from reflections off the body limbs (arms and legs) and the surrounding objects (walls, floor, and ceiling).Obviously, this mechanism is observed only in an indoor environment.
Based on an extensive measurement campaign, we now present accurate statistical models corresponding to the propagation mechanisms described above.

B. On-body propagation (guided diffraction)
As previously emphasized in [12,31], the average received power (in dB scale) is the following linearly decreasing function of the distance: where P(d) is the instantaneous received power (dimension: The average received power, in linear scale, can then be expressed as follows: where where L 0 is a function of L ref , d ref , and g. a In Figure 4a, the loss L is shown as a function of the distance, considering narrowband transmissions at 5 GHz.More precisely, in Figure 4a, experimental measurements (circles) and their linear interpolation (solid line) are shown.Finally, using (3) in (2) one obtains While expression (4) characterizes the average value, it does not provide insights on the instantaneous distribution of the received power.In [31], it has been experimentally observed that the on-body propagation channel is characterized by slow large-scale fading (i.e., shadowing).More precisely, the instantaneous received power at distance d can be expressed as follows: where X is a random variable (RV) which depends on the channel characteristics.As shown in [32] and confirmed by our measurements, X has a log-normal distribution b with parameters μ and s, where s dB typically ranges from 4 to 10 dB, μ dB is the average path loss on the link (dimension: [dB]).Since the loss is accounted for by the term L(d), it follows that μ dB = 0 dB, and the cumulative distribution function (cdf) of X reduces to the following: with the following corresponding probability density function:

C. Reflections off the environment
The second significant propagation mechanism originates from the multiple reflections off the environment.A substantial measurement campaign has shown that the contribution of the environment can be considered, on average, as an additive, constant power when the transmission distance is significant (i.e., when d >25 cm).The obtained results are shown in Figure 4b, and the power received by means of reflections from the surrounding environment is shown as a function of the distance.It can be observed that when d >25 cm, the value of the loss is, on average, around -78 dB.More precisely, for d >25 cm, the average value of the received power can be expressed in logarithmic scale as follows: where P is the transmit power and L (env) dB −78dB.Alternatively, the average received power can be expressed in linear scale as where L (env) = 10 L (w) dB /10.
Our measurement campaign has shown that the propagation channel can be accurately characterized as narrowband Rayleigh block fading.Therefore, the instantaneous received power P env has the following exponential distribution [33]:

D. A unified BAN propagation model
The combination of the two propagation mechanisms presented in Sections 2-B and 2-C allows to derive a unified propagation model for a generic BAN.It can be observed that the degree of importance of each mechanism depends on the distance between transmitter and receiver.More precisely, in close proximity, the dominant propagation mechanism is the on-body propagation described in Section 2-B.Above the cross-over distance d cross ≈ 25 cm, the contribution of the environment becomes dominant, and the second propagation mechanism, presented in Section 2-C, is the only relevant one.Therefore, a unified propagation model can be characterized as follows: • in an outdoor environment, the average received power can be computed using (4) (i.e., E[P(d)] ∝ P10 γ d ) and the instantaneous received   power is determined by the log-normal fading channel model given by ( 5); • in an indoor environment, -if d ≤ d cross , the average received power can be computed using (4) (i.e., E[P(d)] ∝ P10 γ d ) and the log-normal fading in ( 5) is used; -if d > d cross , the average received power is approximately constant (i.e., E[P(d)] = PL (env) ) and the instantaneous received power, owing to a Rayleigh faded channel model, has the distribution given by ( 8).
In Figure 5, the average path loss is shown as a function of the distance.In particular, the overall (unified) path loss can be expressed as follows:

Link-level performance analysis
In this article, we consider multi-user communicationsin a BAN, all sensors need to transmit to a central controller and, in this sense, the scenario at hand can be interpreted as a multi-user scenario.The transmission over a link of interest is denoted with the subscript "0."Besides the intended transmitter, other nodes may be interfering.Depending on their distance to the receiver, the interfering nodes will be denoted differently.More precisely, • in an indoor scenario, the interferers located at distances shorter than d cross are referred to as "closerange interferers," their number is indicated as N close , and the generic node will be denoted with a subscript i ∈ N close {1, 2, . . ., N close }; • in an indoor scenario, the interferers located at distances longer than d cross are referred to as "far-range interferers," their number is indicated as N far , and the generic node will be denoted with a subscript j ∈ N far {1, 2, . . ., N far }; • in an outdoor scenario, the number of interferers is indicated as N out , and the generic node will be denoted with a subscript k ∈ N out {1, 2, . . ., N out }.
The transmission state of the a node at time t is characterized by the following indicator variable: 1 if the node is transmitting at time t, 0 if the node is silent at time t.
Assuming slotted transmissions (i.e., t can assume multiples of the slot time), a simple random access scheme is such that, at each time slot, a node transmits with probability q A transmission in a given link is successful if and only if the signal-to-noise and interference ratio (SINR) at the receiver is above a certain threshold θ.This threshold value depends on the receiver characteristics, the modulation format, and the coding scheme, among other aspects.The SINR at the receiving node of the link is given by where P 0 (d 0 ) is the received power from the link source located at distance d 0 , N 0 is the power noise spectral density, B the channel bandwidth, and P int is the total interference power at the link receiver, i.e., the sum of the instantaneous received powers from all the undesired transmitters.More precisely, in an indoor environment, one has and, in an outdoor environment, one has Finally, as typical in the context of BANs, we assume that all nodes use the same transmit power, i.e., P i (0) = P j (0) = P k (0) = P 0 (0), ∀i, j, k.

A. Link probability of success with short-range transmission in indoor scenarios
The link probability of success for a required threshold SINR value θ in the context of a short, indoor, log-normal faded link is equal to In the Appendix, it is shown that where {c m } n m=1 and {a m } n m=1 , n being an integer determined by the expansion accuracy, are suitable coefficients.By using the function ζ(•; •) and recalling expression (10) for the interference power, the link probability of success ( 12) can be written as follows: where, in the last passage, we have used the fact that the RVs {Λ i , Λ j , P env , X i } are independent.The term in the third line of the expression at the right-hand side of ( 13) can be further expressed as The final integral expression in ( 15) can be numerically computed.The term in the third line of expression (13) can be expressed as follows: Finally, by inserting ( 15) and ( 16) into ( 13), the link probability of success can be given by the expression in (17).

B. Link probability of success with long-range transmission in indoor scenarios
The Rayleigh-faded channel model applies to indoor links with length d >d cross .In this scenario, E[P(d)] ≈ P env (for both the intended transmitter and interferers) and the link probability of success can be expressed as follows: It can be observed that the terms in the second and third lines at the right-hand side of (20) are similar to (15) and (16).Therefore, by using the same derivation of Section 3-A, with P 0 L(d 0 ) replaced by P 0 L (env) , one has and By inserting ( 21) and ( 22) into (20), one obtains the final expression (18) for the probability of successful transmission on the link.

C. Link probability of success in outdoor scenarios
In these scenarios, the links are subject to log-normal fading, and exponential power decreases.The link probability of success can simply be derived using the derivation in Section 3-A, setting N out = N close and N far = 0 (this does not mean that there are not far interferers, but that their propagation model is simply the same of close interferers).Therefore, the computation of the link probability of success P (outdoor) is straightforward from (17) and the final expression is given in (19).

D. Minimum transmit powers
The first terms in the sum at the right-hand side of (17) and the first multiplicative term at the right-hand side of (18) correspond to the link probabilities of success in a noise-limited regime, i.e., when no interferers are present.In fact, setting N close = N far = 0 (i.e., P int = 0) in ( 17) and ( 18), the probabilities of successful link transmission reduce to Therefore, if a threshold link probability of success equal to P th ∈ (0, 1) is required, the minimum required transmit power in an indoor scenario can be written as follows c : where N 0 has been expressed as Tk b , with T being the room temperature (dimension: [K]) and k b = 1.38 × 10 - 23 J/K being the Boltzman's constant, and B being the transmission bandwidth.
In an outdoor scenario, by setting N out = 0 in (19), the probability of successful link transmission reduces to Considering a threshold link probability of success equal to P th , the minimum required transmit power becomes In Figure 6, the minimum required transmit power P 0 for a successful link transmission in an indoor scenario with a ZigBee equipment (B = 5 MHz, θ = 5 dB), operating at T = 300 K and with log-normal fading characterized by s dB = 8 dB, is shown as a function of the distance, considering various values of the required link probability of success of P th .As expected, once the link distance overcomes the critical value around 25 cm, the required transmit power becomes constant.The dashed region corresponds to the typical operational region.In Figure 7, the minimum required transmit power for an outdoor scenario is shown as a function of the distance.The system parameters are set as in Figure 6.It can be observed that, unlinke an indoor scenario, in an outdoor scenario, the minimum required transmit power is an increasing function of the distance (in fact, there are no reflections from surrounding objects).
On the basis of the results presented in Figures 6 and  7, the following observations can be made.The value of P th plays a limited role on the minimum transmit power.If the transmit power is constrained by energy concerns, then only short-range communications (some tenths of centimeters) will be possible: therefore, a multi-hop network architecture is to be preferred.Finally, in an indoor environment, as seen from Figure 6, the reflections from the surrounding environment make the minimum transmit power become constant when d ≥ 25 cm.
In the remainder of this study, we will consider only interference-limited BANs, i.e., scenarios where condition ( 23) is satisfied.Formally, this is equivalent to assuming that N 0 B ≪ P int .

A. BAN tree topologies
In [35], a preliminary performance analysis of BANs with star topologies was carried out.Indeed, these topologies are well suited for medical applications since they exhibit low-power consumption [36] and can perform application-specific data aggregation [37][38][39].However, in order to limit the transmit power, the use of tree (hierarchical) BAN topologies is appealing.
In Figure 8, an illustrative tree topology is presented.It can be observed that, in a generic situation, multiple hierarchical levels have to be considered because of the existence of multiple measurement clusters.Each cluster has a cluster-head, which collects the data from its sensors (and its own data) and transmits them to the final sink.We assume that the links in each cluster are short (i.e., each cluster is in a regime of close-range interferers) and the links from the cluster-head to the coordinator are long (i.e., there is a regime of far-range interferers).However, the proposed framework is applicable to any type of tree architecture.
In this article, we will focus on the impact of the tree clustering on the throughput and energy consumption.More precisely, in Figure 9, three two-level (i.e., 3-tier) hierarchical topologies with 16 nodes are presented.

B. Medical applications of the tree topologies
The three topologies shown in Figure 9 are generic and suitable for a range of medical applications [40,41].More precisely, "Configuration A" refers to a multi-sensor site where highly dense clusters of nodes are deployed.This is representative of medical scenarios where intense monitoring, in a few areas of interest, is needed.Relevant medical applications are mobile EEG (ElectroEncephaloGraphy) or post-operative monitoring of localized critical health conditions.
The second configuration ("Configuration B") is more balanced and corresponds to multiple monitoring sites distributed over the body.Two typical BAN scenarios are encompassed: (i) redundant acquisitions of local physiological signals (for safety reasons) and (ii) multiple independent sensing devices, each having its own relay node (i.e., ECG (ElectroCardioGraphy) combined with limbs monitoring and motion sensors).Relevant medical applications include stroke or Parkinson's disease monitoring (through a combination of EEG, accelerometers, and a gyroscope), and cardiac arrest or ischaemic heart disease monitoring (through a combination of an ECG and a mechanoreceptor).The third configuration ("Configuration C") is representative of a generic sensing scheme where multiple sensors are networked and distributed all over the body without local clustering.In this sense, it is representative of a star topology, as each intermediate relay is connected to a single sensing unit.A relevant medical application is given by a wearable vest with multiple sensors across it (each node may measure local blood pressure, collect electrical signals for ECG, and measure local accelerations).

C. Multi-hop traffic model
In this study, we consider a slotted communication model, where T slot (dimension: [s]) denotes the duration of each slot.It is important to distinguish between data generation and data transmissions at the sensors.Data generation, in real applications, depends on the quantity to be measured; data transmission depends on the communication system design.We now show clearly that generation and transmission cross-influence each other.
Let us first model data generation.For the sake of simplicity, we assume that, in each slot, a sensor can generate at most one packet, and we denote by l [0, 1] the corresponding probability of packet generation.In other words, the number of packets generated by a sensor in a slot is a Bernoulli RV with parameter l, i.e., l can also be interpreted as the average number of packets generated in a slot.Therefore, l/T slot represents the average number of generated packets per unit time (dimension: [s -1 ]).Finally, denoting the packet length as L (dimension: [b/pck]) and the (fixed) transmission data-rate as R b (dimension: [b/s]), the packet duration is T pck , ≜ L/Rb (dimension [s]).For stability reasons, it has to hold that Given specific transmission technology (which determines R b ) and communication protocol, e.g., Zigbee (which determines the percentage of overhead in a transmitted packet), it is possible to determine the maximum payload per slot by imposing T pck = T slot .Denoting L payload < L the length of the payload and by r samp- med the sampling rate of the medical sensor (dimension: [b/s]), the following condition has to hold: In other words, the above inequality shows clearly that the communication/networking technology has an impact on the features of the (medical) sensors.We remark that a careful analysis of the transmission probabilities of the (medical) sensors will more likely lead to different values of l (and q) for each node, depending on the type of physiological constant and the congestion at the relays, among others.This analysis goes beyond the scope of this article and is the subject of future research.However, whatever the used sensors, it is possible to derive the equivalent value of l and, therefore, rely on the proposed framework.
At this point, we model data transmission.Under the considered assumption of slotted ALOHA MAC protocol, a simplified model for the MAC protocol, a sensor has probability q of transmitting a packet in a slot.Obviously, this makes sense only if the node has a packet to transmit.Moreover, for stability reasons, it has to hold that λ ≤ q.
In fact, the condition l > q would be equivalent to assuming that the sensor generates, per time unit, more packets than those it can actually transmit.In this case, there would be an overflow at the sensor, and packets would be lost.On the other hand, assuming l < q is meaningless as well, it is impossible that the transmission probability of a sensor node is higher than its generation probability (what would it transmit?).Therefore, in the considered simplified model, it follows that l = q, i.e., the generation and transmission processes coincide.Note also that, according to this model, q is equal to the per-node load (defined as the average number of packets generated during an interval equal to the duration of a packet transmission).Therefore, the network load G (adimensional) is simply equal to q • N tot , where N tot denotes the total number of sensor nodes in the BAN.
Let N be the set that consists of N leaf sensor nodes connected to a given relay node (i.e., the set of sensor nodes per cluster, excluding the relay).In half-duplex communications, a node transmits if and only if (i) it has data to send or (ii) it has no data to send but acts as a relay for other nodes.We denote by q leaf the probability that a leaf node has data to send.Obviously, q leaf = q, i.e., the probability that data are present and ready to be sent.A relay node will transmit if it gets data from a leaf (event denoted as "relay") or has to send sensed information (event denoted "data"), i.e., in the last passage, we have exploited the fact that the events "data" and "relay" are independent.By definition, P{data} = q.Obviously, the probability of relaying depends on (i) the probability of having data present at any node and (ii) their successful reception at the relaying node.Therefore, 1eaf→relay .(26) According to the assumption at the end of Section 3, a transmission is successful if the channel is not in an outage, i.e., if the (instantaneous) SINR exceeds a certain threshold θ.Therefore, P 1eaf→relay = P{SINR > θ} on the considered link.Since (i) all links in a cluster are, on average, equal and (ii) q leaf = q, one has Finally, P{data} = q, from (25), one has where P 1eaf→relay can be either (17) or (19), in indoor or outdoor scenarios, respectively.
In Figure 10, the probability of transmission of a relay node is shown as a function of the probability of transmission of a single node, considering various values for the number N 1 of nodes in a first-level cluster (i.e., leaves of the collection tree).It can be observed that when q ≤ 0.5, the value q relay depends on the relaying (i.e., q relay ≥ q since it accounts for the traffic of the leaves plus the traffic generated by the relay) and, when q >0.5, it is dominated by the relay probability of sending the data itself (i.e., q relay ≈ q since the relay transmits its data and prohibits reception of the ones from the leaves).
Finally, in multiple-tier topologies (more complex than the 3-tier considered in this article), the same approach can be applied to compute the probability of transmission of any node acting as a relay at a given hierarchical level of the network.In the considered 3-tier topologies, this approach can be straightforwardly applied to evaluate q sink , i.e., the probability of transmission from the sink (e.g., through a 3G connection).

Network-level performance analysis
The main simulation parameters are set as follows.With reference to the topologies in Figure 9, the distances between a leaf and its relay and between a relay and the sink are 10 and 30 cm, respectively.The SINR threshold is set to θ = 5 dB.The fading power of the lognormal propagation model is s dB = 8 dB.These values correspond to typical, multi-kpbs sensor nodes.

A. Performance metrics
In the following, we will consider two key performance metrics: (i) the link-level throughput, and (ii) the energy consumption rate.
Regarding the link-level throughput, a transmission will be successful if and only if a transmission link is not in an outage.This corresponds to requiring that the (instantaneous) SINR of the link is above the threshold θ, i.e.,

P = P{SINR > θ}
where P is either equal to P (indoor) far , P (indoor) close , or P (outdoor) depending on the scenario at hand.The prob- abilistic link throughput [42] (adimensional) of a node is defined as follows: • in the full-duplex communication case, it corresponds to the product of (i) P and (ii) the probabil- ity that the transmitter actually transmits (i.e., q); • in the half-duplex communication case, it corresponds to the product of (i) P, (ii) the probability that the transmitter actually transmits (i.e., q), and (iii) the probability that the receiver actually receives 1q).
The probabilistic link throughput can be interpreted as the unconditioned reception probability which can be achieved with a simple automatic-repeat-request scheme with error-free feedback [43].For the slotted ALOHA transmission scheme under consideration in the context of BAN (where transmissions can typically be organized in a full-duplex way), the probabilistic throughput is τ q × P.
More specifically, in our analyses, three different nodes are considered: leaves, relays, and the sink.The corresponding throughput metrics are τ relay = q relay × P relay→ sink .
since it is supposed that the sink acts as a special device and can communicate with the external equipments with probability equal to one (e.g., it is used to store data on a memory card or to send these data by means of a reliable transmission technology, such as, for example, 3G).
The second performance metric of interest is the energy consumption rate.The average energy consumed by the network in a slot, denoted as E (dimension: [J]), can be expressed as E = q N 1eaves E TX + q relay N relays (E TX + E RX ) + q sink E RX where E TX and E RX are the energies (dimension: [J]) consumed by single-packet transmission and reception acts, respectively.
In most wireless systems, E TX ≈ E RX , and one has where E E/E TX is denoted as energy depletion rate (adimensional) and corresponds to the ratio of the energy consumed by the network in a slot and the energy consumed to transmit a single packet.
We now provide the reader with a performance analysis of all the three hierarchical topologies deployed in outdoor and indoor scenarios.

B. Outdoor Scenarios
In Figure 11, the per-node throughputs at the three hierarchical levels (i.e., leaf nodes, relaying nodes, and central sink) are shown as functions of the sensors' probability of transmission q.The three subfigures refer to the three topologies in Figure 9.It can be observed that the three considered topologies lead to very different performances (in terms of throughput) for the leaves, the relays, and the sink.
In the Configuration A (Figure 9a) the throughputs of the leaves, relays, and sink are increasing functions of q for small values of q and reach the maximum values at q ≃ 0.15.For q >0.15, τ leaf starts decreasing.In contrast, τ relay remains approximately constant for q (0.15, 0.45): in fact, the packets' losses at the leaves' are compensated by data transmitted by the relay node itself, so that the overall value of τ relay tends to remain stable.It can also be observed that the throughput at the sink remains approximately constant for q (0.15, 0.45), and its value is close to the maximum achievable throughput of any slotted ALOHA system without lossy links, which is e -1 ≈ 0.37 [44].Since the throughput at the sink can be interpreted as the overall network throughput, it can be concluded that the network Configuration A yields an excellent channel utilization at the sink node.Regarding the operational region of this configuration, it can be seen from Figure 11a that for q ≥ 0.6 one has τ leaf ≃ 0, i.e., the load is too high, and the leaves tend to be disconnected from the network, i.e., packets from the leaves are no longer relayed and successfully transmitted to the sink.Consequently, for q ≥ 0.6, the throughputs at the relays and at the sink tend to decrease.
The performance with the second topology-referred to as Configuration B (Figure 9b)-is analyzed in Figure 11b.It can be observed that the throughput at the relays is, for small values of q, an increasing function of q and reaches a maximum value at q ≈ 0.1.Beyond this value, the throughput at the relays tend to rapidly decrease.On the other hand, τ leaf is an increasing function till q ≈ 0.3: this is because the number of sensors per cluster (3) is smaller than the number of relays (4) and, therefore, the throughput at the leaves continues to increase even if the throughput at the relays starts decreasing.Unlike Configuration A, in Configuration B, the maximum throughput at the leaves is higher than the maximum throughput at the relays.As the number of relays is relatively large, they tend to interfere with each other and, therefore, the throughput of the sink reaches a maximum at q ≈ 0.1 and, then, decreases.It can be observed that the maximum throughput at the sink with Configuration B is approximately equal to that of Configuration A. However, unlike Configuration A, in Configuration B, there is no interval of q where the throughput at the sink tends to remain constant.In other words, this configuration does not support, at network level, a larger interval of values of q.
The last scheme-denoted as Configuration C (Figure 9c)-is highly centralized.Its performance is investigated in Figure 11c.As each cluster contains only one leaf node, τ leaf is the highest.On the other hand, the relays interfere with each other while communicating to the sink and, therefore, τ relay remains very low (its maximum value is around 0.05).As a consequence, τ sink , after reaching its maximum value for q ≈ 0.08 (similarly the previous configuration), tends to decrease to zero much faster than in Configuration B. Note that the maximum value of the throughput at the leaves is close to the maximum value of the throughput at the sink.Finally, note that for q ≥ 0.5, even if τ leaf is high, τ sink is basically zero: in other words, no data transmitted by the leaves can be successfully transmitted by the sink to an external controller (e.g., through 3G communications).

C. Throughput in indoor scenarios
In Figure 12, the per-node throughputs at the various hierarchical levels are presented for the three topologies of interest.As a first, general, observation, it is seen that the per-node throughputs are much lower in indoor scenarios than the corresponding ones, as shown in Figure 11, in outdoor scenarios.This can be explained by the presence of a reflections off the limbs and the surrounding objects.Indeed, the initial antenna gain (at d = d ref ) is about L ref = -57.42dB, and this value is not very different from the gain of the environment, i.e., L (env) dB = −78dB.Therefore, short links are less affected (since the received signal power is much stronger than the reflected power), while longer links are more likely to suffer significant interference from the reflected waves.This was not the case in outdoor scenarios where distant nodes did not contribute to the interference thanks to the high path loss coefficient (i.e., g = 124 dB/m).
Regarding Configuration A, it is seen from Figure 12a that the leaves can support a wide range of values of q (i.e., the throughput is non-zero for any value q (0, 0.6)).As anticipated in the description of the results in outdoor scenarios, the relays effectively accumulate the leaves' and their own data, guaranteeing the highest  throughput almost for all values of q-for very low values of q, τ relay < τ sink .However, the last links (i.e., the relayto-sink links) are subject to strong interference because of the reflections off the surrounding environment, and the sink throughput is much lower than in the outdoor scenario.More precisely, the throughput reaches a maximum at q ≈ 0.05 and becomes insignificant for q ≥ 0.3.The performance of Configuration B is presented in Figure 12b.Since the tree is more balanced than in Configuration A (i.e., there are less leaves and more relays), the performance observed at the leaves is better in terms of throughput.However, the increase of the amount of relay nodes and the fact that these are more subject to environment interference (since these are considered as long links) make the throughput decrease significantly.Finally the throughput at the sink remains limited, compared to the outdoor scenario, for the reasons described previously in analysis of the Configuration A.
The third configuration-namely, Configuration C-is shown in Figure 12c.In this configuration, the throughput at the leaves is significant.This could have been expected by taking into account the facts that (i) the links are shorts and, therefore, nearly not subject to interference; and (ii) the amount of concurrent transmissions remains Since there are numerous relay nodes, the throughput at the relays is very low, because of the presence of multiple access collisions.Furthermore, the reflections off the environment reduce drastically the throughput at the sink when the relay probability of transmission increases.It can be observed that the value of τ sink rapidly reaches a maximum for q ≈ 0.05, before decreasing rapidly for increasing values of the parameter q.

D. Energy depletion rate
First, regarding a BAN deployed in an outdoor environment, the energy consumption rate E is presented in Figure 13 as function of q and for the three configurations of interest.It can be observed that the energy consumption rates of the three configurations present clearly different profiles.More precisely, Configuration A outperforms Configuration B which, in turn, is more energy efficient than Configuration C. Also, it can be observed that this remains true for any value of the node probability of transmission q.
The energy consumption rate in indoor scenarios is shown in Figure 14.It is noticeable that the values of the energy consumption in Figures 13 and 14 are approximately the same.Also, the relative performances   of the three configurations in indoor scenarios are the same as in the outdoor case: Configuration A presents the lowest energy consumption rate, whereas Configuration C is the most energy consuming.

E. An illustrative comparison with TDMA-based architectures
In this article, we refer to the IEEE 802.15.6 standard and the slotted ALOHA MAC.Owing to its random access strategy, this MAC protocol could be of less interest if the traffic is high.In that case, a centralized, time division multiple access (TDMA) access might be more interesting and appealing.In this subsection, we provide the reader with an illustrative comparison between the slotted ALOHA and TDMA schemes, focusing on the delay exhibited by both strategies-in fact, the throughput of a TDMA scheme is equal to 1.
In a slotted ALOHA system, the delay is directly related to the amount of (re-)transmissions needed to send (or forward) a packet.More precisely, the average number of transmissions can be written as where Q r is the probability of a successful communication realized after exactly r transmissions (i.e., r -1 unsuccessful transmissions and a successful transmission).It can be expressed as In other words, r is a geometric RV with parameter P s , and one obtains where the link probability of success P s depends on the scenario (i.e., outdoor or indoor, position of the nodes, etc.) and is given by relations ( 17), (18), and (19).Finally, the delay, expressed in number of time slots, is In a TDMA system, the delay D TDMA depends on the amount of time to wait before a dedicated slot takes place.In a generic approach, it can be supposed that a relay node will allocate exactly one slot per sensor to receive its data and a slot per sensor for the forwarding uplink.In Figure 15, the slots' allocation (i.e., chronogram) is presented for each relay (noted R 1 , R 2 , . . . ,R N 2 ).The time slots have a fixed duration of T s (dimension: [s]).As can be seen from Figure 15, each relay needs a frame of N 1 time slots to collect the (possibly generated) packets from its N 1 leaves.It then needs another frame (N 1 time slots) to forward them to the sink.At this point, it needs to remain idle for (N 2 -2) frames, as the sink is busy collecting the packets from the other relays.This corresponds to assuming the same transmission rates at leaves and relay, and the same TDMA-based approach as at first and second layers.Therefore, the distance between two consecutive slots assigned to a given leaf is equal to N 1 • N 2 slots: when a leaf generates a packet, it needs to wait a number of slots between 0 (its slot is the current one) and N 1 • N 2 -1 (its slot just passed).d As each number of slot has the same probability, the average delay (expressed in time slots) experienced by a given leaf node is The above derivation for D TDMA represents an "average" scenario where a node generates at most a packet in an interval equal to N 1 • N 2 time slots: this corresponds to requiring that l ≤ 1/(N 1 • N 2 • T s ).In this case, expression (32) for the delay is correct.If, on the other hand, l >1/(N 1 • N 2 • T s ), then it means that as soon as a leaf has transmitted a packet to its relay, it is likely to generate shortly a new packet, which will wait a period longer than (32).In general, it can be stated that In Figure 16, the average delay (expressed in contention time slots) incurred by a leaf to reach its relay is shown as a function of the transmission probability q.All curves refer to an indoor scenario.In the TDMA scheme, we consider the average expression e (32).As expected, it can observed that, when the node probability of transmission is low, the slotted ALOHA significantly outperforms the TDMA scheme.However, for increasing probability of transmission, i.e., for increasing traffic load, there exists a critical threshold above which the TDMA scheme is to be preferred.
To summarize, as a TDMA-based scheme has a throughput τ = 1, it becomes very attractive for values of q beyond the maximum of the considered slotted ALOHA system, as the latter becomes unstable, i.e., the value of the delay D ALOHA ∞ since P s 0. In scenarios with low reporting rate, the slotted ALOHA scheme is to be preferred.

Conclusions
In this article, we have presented an analytic framework for the evaluation of the link probability of success in interference-limited BANs subject to fading.This analytic derivation is based on novel experimental measurements which highlight two characteristic propagation mechanisms in BANs deployed in indoor scenarios: onbody propagation, and propagation through reflections from the environment.
Regarding the impact of the topology, three configurations have been analyzed.These analyses showed significantly different performances, in terms of per-node throughput and energy consumption rate.It can be concluded that a decentralized topology in an outdoor environment presents the best tradeoff between the throughputs at the leaves, the relays, and the sink.Moreover, the shape of the throughput curves shows that it is very stable, i.e., any increase or decrease of the node generation rate does not significantly impact the per-node throughput, and it is the most efficient in terms of energy consumption rate.In an indoor environment, the balanced tree seems to be more suitable.Indeed, it presents a higher throughput for the sink, the relays, and the leaves at the same time.
In conclusion, multi-user BANs deployment and operation should take into account the specificities of environment and adapt the routing algorithms and clustering strategies accordingly.In the particular context of an indoor environment or when the traffic load is high, narrowband-shared spectrum techniques are not performant, and other MAC schemes, such as TDMA should be considered instead, even if they are more complex to implement.

Figure 1
Figure 1 Possible positions of a transmitter-receiver pair in a BAN.
[W]) at distance d (dimension: [m]), P is the transmit power (dimension: [W]), d ref is a reference distance (dimension: [m]), L ref is the gain at the reference distance (adimensional, in dB), and g is a suitable constant (dimension: [m -1 ]).For instance, typical experimental values for these parameters are d ref = 8 cm, L ref = -57.42dB, and g = -124 dB/m [31].

Figure 2 Figure 3
Figure 2 Power delay profile as a function of the time in an indoor environment and for d ≤ 25 cm (body front).

Figure 4
Figure 4 Propagation loss as a function of the distance: (a) onbody propagation, and (b) propagation through reflections off the environment.In both cases, experimental results (circles) and their linear (or piece-wise linear, in case (b)) interpolations (solid line) are shown.

Figure 5
Figure 5 Generic propagation model (on-body and environment reflection superimposed).

Figure 6
Figure 6  Minimum transmit power as a function of the distance in an indoor scenario.The dashed region is the operational region of a BAN.

Figure 7 Figure 8
Figure 7 Minimum transmit power as a function of the distance in an outdoor scenario.The dashed region is the operational region of a BAN.

2 N 1 = 4 N 1 = 8 qFigure 10
Figure 10 Terminal probability of transmission of a relay node as a function of a single terminal probability of transmission and for N Î {2, 4, 8} neighbor nodes and P leaf→relay = 1.

Figure 11
Figure 11 Per link throughputs for the three considered strategies in an outdoor environment.(a) Configuration A: seven leaves at each of the two relays (N 1 = 7, N 2 = 2); (b) Configuration B: three leaves at each of the four relays (N 1 = 3, N 2 = 4); and (c) Configuration C: one leaf at each of the eight relays (N 1 = 1, N 2 = 8).

Figure 12
Figure 12 Per-node throughputs (at each level) for the three considered topologies in an indoor environment.(a) Configuration A: seven leaves at each of the two relays (N 1 = 7, N 2 = 2); (b) Configuration B: three leaves at each of the four relays (N 1 = 3, N 2 = 4); and (c) Configuration C: one leaf at each of the eight relays (N 1 = 1, N 2 = 8).

Figure 13 Figure 14
Figure13 Energy consumption rate E of an outdoor hierarchical BAN with 16 nodes as a function of the leaves' probability of transmission and for different hierarchical configurations.

Figure 17
Figure 17 The function ζ(z; s) as a function of z, considering various values of s (in dB).

Table 1
Coefficients for the approximation of the ζ function