Abstract
Nowadays, ubiquitous network access has become a reality, thanks to Unmanned Aerial Vehicles (UAVs) that have gained extreme popularity due to their flexible deployment and higher chance of lineofsight links to ground users. Telecommunication service providers deploy UAVs to provide areal network access in remote rural areas, disasteraffected areas, or massiveattended events (sport venues, festivals, etc.), where full setup to provide temporary wireless coverage would be very expensive. Of course, a UAV is batterypowered with a limited energy budget for both mobility aspect and communication aspect. An efficient solution is to allow UAVs switching their radio modules to the sleep mode in order to extend the battery lifetime. This results in temporary unavailability of communication feature. Within such a situation, the ultimate deal for a UAV operator is to provide a costeffective service with acceptable availability. This would allow meeting some target quality of service while having a good market share granting satisfactory benefits. In this article, we exhibit a new framework with many interesting insights into how to jointly define the availability and the access cost in UAVempowered flying access networks to opportunistically cover a target geographical area. Yet, we construct a duopoly model to capture the adversarial behavior of service providers in terms of their pricing policies and their respective availability probabilities. Optimal periodic beaconing (advertising the presence of the UAV) is to be addressed, given the UAVs with limited battery capacity and their recharging constraints. A full analysis of the game, both in terms of equilibrium pricing and equilibrium availability, is derived. We show that the availabilitypricing game exhibits some nice features as it is submodular with respect to the availability policy; whereas, it is supermodular with respect to the service fee. Furthermore, we implement a learning scheme using best response dynamics that allows operators to learn their joint pricingavailability strategies in a fast, accurate, and distributed fashion. Extensive simulations show convergence of the proposed scheme to the joint pricingavailability equilibrium and offer promising insights into how the game parameters should be chosen to efficiently control the duopoly game.
1. Introduction
In the recent few years, Unmanned Aerial Vehicles (UAVs), or drones, have attracted a lot of attention, since they present advantages of ubiquitous usability, high maneuverability, and lowcost deployment. UAVs, equipped with navigation systems and smart sensors, are currently being deployed for surveillance operations, rescue missions, and rapid ondemand communication. Along with the maturity of the UAVs technology and relevant regulations, they have become an important market, where a UAVs worldwide deployment is expected. For instance, the registered number of UAVs in use in the U.S. exceeds 200,000 just in the first 20 days of January 2016, after the Federal Aviation Administration (FAA) started requiring owners to sign up [1]. Indeed, UAVs can be deployed fruitfully to support cellular communication systems when terrestrial infrastructure networks are damaged. Moreover, it supports wireless communication in exceptional scenarios such as hard to reach rural areas, festivals or sporting events, and emergency situations where the terrestrial base stations installation may be too expensive. In this context, UAVs as flying base stations have been reported as a promising approach that can boost the capacity and coverage performance of existing wireless networks. Using UAVs as aerial base stations provides several benefits such as high maneuverability and robustness, flexible deployment, and efficiency on demand telecommunication and also mobility [2].
Thanks to their high mobility and air location, UAVs are likely to have good communication channel since the UAVground link is more likely to have LineofSight (LoS) links [3]. Furthermore, UAVs provide unlimited telecommunication applications for their advantage to efficiently sense information about the surrounding environment. They are now a key technology for intelligence, recognition, search, inspection tasks, surveillance, public services, and so on. Another key feature of UAVs is the Internet of Things (IoT) which is a technological revolution that has made the leap from conceptual to actual [4, 5]. IoT enables devices to exchange data and interoperate within the internet infrastructure, which affords ubiquitous connectivity while reducing the transmission cost and providing extended range for lowpower communication. Technically, UAVs play an important role in the Internet of Things’ vision. It offers flexible deployability and reprogramming in mission possibilities to deliver numerous operating solutions and services, for IoT partners.
While deploying drones as flying base stations offers numerous advantages, still a number of economic and technical challenges arise. These challenges comprise not only technical and physical issues but also effective pricing and availability of management features. Yet, for UAVs service providers to overcome these relevant challenging issues, availability or beaconing scheduling becomes a very important but largely unexplored topic. For this, comprehensive modeling and performance analysis of UAV’s setups has become extremely attractive. In this article, we provide a UAV’s joint pricing and availability problem, considering their limited battery capacity. We focus on the scheduling of availability periods as a key enabler that can potentially lead to better energy efficiency and then, to a better longterm availability (lifetime) as well as a satisfactory quality of service. More precisely, we are interested in the joint pricingavailability problem for the UAVbased network market. Then, we design a new scheme that jointly considers and solves the UAV’s pricing and energyefficient issue.
In order to address these issues, we examine the availability and pricing interactions between service providers as a noncooperative game. We propose a duopoly setup, where a finite network of mobile unmanned aerial vehicles is deployed according to a homogeneous Poisson Point Processes (PPP) and serving a number of ground IoT devices. The UAVs are moving according to a Random Way Point (RWP) model. As a first proposal, we derive the coverage probability expression as well as the service probability for each UAV in the current scenario. In order to achieve the maximum system performance in terms of pricing policy and energy efficiency, our proposal introduces the Nash equilibrium analysis. Furthermore, we propose a learning automaton to derive the joint priceavailability equilibrium. Finally, we provide extensive numerical simulations to highlight the importance of taking priceavailability as a joint decision parameter and provide thereby important insights/heuristics into how to set them.
The remainder of this article is organized as follows. In Section 2, we give an overview of the proposed UAVs duopoly system model and strategic pricingavailability. In Section 3, we present the availability game analysis, existence, and uniqueness of NE solution and also define the game submodularity characteristics. Furthermore, the pricing game is provided in Section 4; the property of supermodularity was discussed mathematically as well analytically. Section 5 considers a joint availabilitypricing approach with numerical learning implementations, and we further gave the impact on some important parameters on the learning process. Finally, Section 6 concludes the paper.
2. Related Work
A significant existing literature investigates interesting features on UAVs technological performance and reducing cost. For instance, the authors in [6] study energyefficient UAV communication via optimizing the UAV’s trajectory; for this, they consider a UAV trajectory with general constraints, under which the UAV’s flight radius and speed are jointly optimized. An important study of the UAV’s obstacle avoidance plan is conducted in [7], where the writers propose algorithm based on Iterative Regional Inflation by Semidefinite in order to improve the fighting efficiency. A useful derivativefree one, Sliding Mode Control (SMC) theorybased learning algorithm, for the control and guidance of a UAV is proposed in [8]. Robustness, in realtime applications and finitetime convergence, was proven through intensive simulations. Another learning approach is treated in [9], where the authors provide a thorough review on deep learning methods for UAVs, including the most relevant developments as well as their performances and limitation. An important issue which is IoT applications is discussed [10]. This paper investigates the efficient deployment and mobility of multiple unmanned aerial vehicles, in order to enable reliable uplink communications for the IoT devices with a minimum total transmit power. Another IoTrelated work is in [11], and it presents the features that a UAV system should implement for drone traffic management, namely, the air traffic control network, the cellular network, and the Internet. In [12], the authors introduce different classifications of Internet of Things with examples of utilizing IoT technologies. They refer to the concept of globalscaled IoT, where using unmanned aerial vehicle offers a large coverage and more flexibility, thanks to the mobility approach. The UAVbased IoT concept has been further studied in the paper [13], which provides a comprehensive survey on the UAVs and the related issues. It highlights potential for the delivery of UAVbased valueadded IoT service from height.
3. AvailabilityPricing Game among UAV Operators
3.1. The Flying Access Network Framework
3.1.1. Coverage Probability
Consider a circular geographical area with a radius R within which a number N of wireless users are deployed according to a homogeneous Poisson point processes with density of u (number of users per m^{2}). In this area, a set of UAVs moving randomly according to a random waypoint mobility model are used as aerial mobile base stations to provide wireless service for the ground IoT devices (Figure 1). The drones are belonging to different operators and are engaged to provide an effective coverage for mobile IoT users. In this setting, the drones have the same characteristics such as altitude h, total available bandwidth, and maximum transmit power.
The coverage probability is formally defined as follows:where is the expression for a device user that connects to the UAV, is the received signal power from the UAV including fading and path loss, N is the noise power, I is the total interference power from existing transmitters, and β refers to the threshold.
According to [14], for an uplink user located at , the coverage probability is given by the following equation:where r and ϕ are the radius and angle in a polar coordinate system, and the UAV is located at the center of the area of interest. , respectively, are the lineofsight, respectively, nonlineofsight probabilities.
The average coverage probability is computed by taking the average of over the area with the radius R:where is the path loss exponent over the userUAV link, η is an additional attenuation factor due to the NLoS connection, and is the UAV transmit power.
The airtoground signal propagation is almost affected by obstacles in the surrounding environment. Thus, a common approach for channel modeling between the UAV and downlink users is based on probabilistic direct lineofsight (LoS) and nonlineofsight (NLoS) links. Whereby each link may occur with a specific probability which depends on the propagation environment, the UAV and users locations, and the elevation angle (here, θ).
In this scenario, we suppose designing a UAVbased communication system in a remote rural area where full setup to provide temporary wireless coverage would be very expensive. Due to its higher altitude and the low average height of building in hard to reach rural areas, UAV has higher probability of direct lineofsight links to ground users. The probability of the NLoS component occurrence is therefore significantly lower than that of the LoS component. For simplicity, we ignore the small scale NLoS, and as evident, the coverage probability formula would be significantly simplified as follows:
3.1.2. Service Probability
Given the limited capacity and recharging difficulties of a UAV battery, periodic beaconing is one of the most challenging problems that need to be addressed. Each UAV i is periodically sending beacons of duration advertising this presence to mobile IoT devices. The beacon/idle cycle is periodically repeated every time slot T during a time window: (Figure 2). UAVs should define their beaconing periods strategically in order to maximize their encounter rate with the IoT ground users. However, they should avoid battery depletion resulting from maintaining useless beaconing in the absence of contact opportunities.
The first encounter follows an exponential distribution with a random parameter λ. In order for a UAV i to encounter first the ground IoT destination at time , the following conditions must hold:(i)The UAV i has to be beaconing at .(ii)All potential encounters with ground destination, happening before the time instant , by the other UAVs need to be unsuccessful. In other words, the encounters need to happen while UAV i’s competitors are inactive.
Consequently, the successful contact probability is given by the following equation:
For mathematical calculations, we focus on the case of two operators ( and ), and this minimizes mathematical complication but still allows us to analyze the important features of operators strategies. We define the probability of beaconing while encountering for the first time the destination within [0; m]:
For a , the probability of being idle is given by the following equation:
The probability that encounters first the ground destination without accounting for its state (probing/idle) is expressed as follows:
And finally, we define the probability that encounters first the ground destination without accounting for its state:
3.2. The Game Formulation
3.2.1. Customers Behavior
Most telecommunication markets are oligopolistic, this means, dominated by a certain number of service providers called oligopolists. In order to maximize benefits, the competition between service providers is becoming increasingly tough. From an oligopolist’s point of view, increasing market share is the most important objective. Thus, service providers are supposed to define the optimal pricing policy, and the best QoS (represented by availability) requested in order to attract more and more customers. Because customers are rational by nature, they are likely to register to an operator rather than the others or to stay at no subscription state according to their own satisfaction (service price, quality of service, and quality of experience).
In this paper, we assume a duopoly telecommunication market where just two service providers compete against one another to provide service for IoT ground users. Each UAV i belonging to SP i picks its availability duration represented by the periodic beaconing time chosen within the interval and a service fee per data unit . Formally, SPs are intelligent individuals in conflict, and this situation is involving rational decisions.
The noncooperative game theory provides the appropriate tools for determining optimal behavior in such competitive environments. It provides mathematical methodology for modeling and analyzing interactions between intelligent rational players in conflicts of interest. This branch of game theory deals largely with how intelligent individuals interact with one another in an effort to achieve their own goals. Players are selfish, and each individual player makes decisions independently without forming alliances.
We capture the pricingavailability interaction among two UAV operators as a normal form game. Each SP seeks to choose the level of service allowing a reasonable battery lifetime while maximizing its monetary profits, given the action of the competitor. The noncooperative game is formally defined as follows:where(i) is the set of players (UAV operators)(ii) is the strategy space of both UAV operators(iii) is the vector of payoffs defined as the difference between their respective profits and costs.
Increasing the experienced market demand is the most important objective of each operator. Then, from an operator perspective, the question is to define the best beaconing duration and the best pricing policy that attracts a high portion of subscribers. Ideally, one UAV would like to simultaneously maximize the proposed fee and minimize the total consumed energy. However, these are two conflicting objectives, since ground users are likely to be rational. Inspired by the logit model [15], we model the probability that a given IoT device is registering with the service provider i which can also be seen as the market share of a UAV_{i}. The demand function follows an exponential distribution with the parameter μ. This last parameter is called the constant hazard or temperature, which defines the users’ rationality, so that, the higher is μ, the less users are rational:
The average number of customers that are registering to UAV i is . Because UAVs compete in quality of service and price, it seems to be obvious that the experienced demand depends on these parameters. However, it depends also on the environment variables.
3.2.2. Service Providers’ Strategic Decision
We are interested in the problem of optimal pricingbeaconing in the fixedwing UAVs that should avoid battery depletion resulting from maintaining useless activity duration in the absence of contact with the ground. We propose a new design paradigm that considers both the UAV’s pricing and energyefficient problems.
The UAV is periodically sending beacons advertising the presence to mobile users on the ground. From the UAV perspective, there are tradeoffs between the strategic parameters τ and energy consumption. On one hand, as the beacon duration increases, the coverage chance as well as the service probability increase. On the other hand, total energy consumption is mainly proportional to the activity period duration. The same applies for the service price, where increasing the pricing policy value may increase revenues, but in the same time, it negatively affects the experienced market demand. Here, the strategy space and the reward are common knowledge; however, the chosen beaconing period and taxed fee are not since decisions are taken simultaneously. The utility of each UAV SP is a function of its own strategy as well as the decisions of the other:where is the dissipated energy for a :where is the normalized throughput of the user u served by a , is the energy cost per slot for sending beacons, is the reception energy, is the transmission energy cost, is the unit energy for remaining switching the transceiver status, and is the set of ground IoT devices served by the given drone i.
4. Availability Game with Fixed Fee
In literature, the pricing scheduling is almost the key parameter that affects straightly operators’ revenues. However, the realism of this assumption is sometimes questioned [16]. Indeed, similar market price may affect the customers’ perceptions about the UAVs’ availability and then their loyalty. Due to a rough competition, availability has started to become an important strategic tool for operators to expand their market share. We examine energyefficientpricing tradeoffs within a noncooperative game framework. This incites operators to make efforts in providing better QoS while fixing optimal periodic advertising beacons. The proposed approach is realistic enough to analyze the interplay of availability on service providers’ revenues. The problem is then, to determine for each UAV, the optimal availability (beaconing period) while rival behaves optimally, and this represents exactly the Nash equilibrium status.
Recall that the Nash equilibrium strategy is the best response of each UAV to the expected best response behavior of the other.
Proposition 1. The availability game admits at least one Nash equilibrium.
Proof 1. To show the existence of an equilibrium, the sufficient condition is the quasiconcavity of the utility function. We notice that the strategy set is a convex, close, and compact interval. Furthermore, the utility function is continued with respect to . In addition, we show that the second derivative with is negative. After calculation, the second derivative is expressed as follows:whereIt has been easily proven that the second derivative of the service probability with respect to is of negative sign, and therefore, for the second derivative to be negative, it could be sufficient to make sure that the following condition fulfills:This condition seems logical since the demand is positive, and each operator requires a certain demand threshold as an incentive to deploy a drone resource; this means that for a UAV operator to be motivated to coverage a certain area, the minimum required demand equals at least the cost of energy that is going to waste. From this, the utility function of a UAV i belonging to an operator i is quasiconcave, and then admits a local optimum point at which is it its best response against strategies of other operators.
Proposition 2 (uniqueness of availabilitytermed Nash equilibrium). The availability game is dominance solvable [17] and consequently satisfies Rosen’s diagonal strict concavity condition. Henceforth, it has a unique Nash equilibrium.
Proof 2. The dominance solvability condition is given by the following equation:whereThe dominance solvability conditions are verified and consequently the concave game satisfies Rosen’s diagonal strict concavity condition for the uniqueness of the Nash equilibrium.
Interactions among UAVs in this noncooperative game have a very attractive property, if one UAV reduces its availability (beaconing period) and the other UAV also has an interest in decreasing its own. Informally, this type of games led us to conclude that it is about submodular games.
Lemma 1. The game with N mobile devices and two UAVs service providers is submodular.
Proof 3. To prove the submodularity of the game, we compute the secondorder mixed derivative of with respect to two beaconing durations and . With the mixed derivative expression given in (19), it is straightforward to conclude that it is effectively about a submodular game.
With submodularity, the best response of a UAV is a decreasing function of the other UAV availability strategy (beaconing period). Consequently, this game possesses unique equilibria that the best response dynamic learning scheme can reach with probability 1.
5. Pricing Game with Fixed Availability Period
Despite the several advantages of availability competition, still the price is a more common strategical decision in telecommunication business. This motivates us to investigate the pricing game with the fixed availability parameter. We argue in this section that operators set prices and let the consumers decide on their experienced demand. That is, each operator is maximizing its profit believing that its competitor is doing the same. Conceptually, this situation is a typical Bertrand–Nash equilibrium (see papers in [18, 19]). We turn first to mathematically analyze the equilibrium existence.
Proposition 3. The pricing game admits at least one Nash equilibrium.
Proof 4. The strategy set represented in is a convex, closed, and compact interval. In addition, the utility function is continuous with . Furthermore, we show that the second derivative with respect to has a negative sign. After calculation,The secondorder derivative is expressed as follows:which is clearly negative. Hence, the profit function of each UAV is quasiconcave with respect to its own fee and then the NE existence.
Proposition 4 (uniqueness of pricebased Nash equilibrium). The pricing game has a unique Nash equilibrium.
Proof 5. The utility function satisfies the dominance solvability conditions and consequently satisfies Rosen’s conditions, which is defined as follows:A remarkable property of this game is that when one player picks a higher action, the others have incentive to follow and do the same. Roughly, this introduces the notion of a supermodular game characterized by strategic complementarities.
Lemma 2. The game with N mobile devices and two UAVs service providers is supermodular.
Proof 6. To prove the supermodularity of the game, we compute the secondorder mixed derivative of with respect to its proposed fee and that of its competitors ( and ). With the following mixed derivative expression, which is clearly of a positive sign, we argue that it is a supermodular game:Supermodularity tends to be analytically appealing. An interesting property is that it behaves well under various learning rules. Furthermore, convergence towards the Nash equilibrium point is guaranteed under best response dynamics.
6. Joint PricingAvailability Game
6.1. Insights into RealWorld Implementation: Fully Distributed Learning
Learning is a fundamental component of intelligence and cognition; it is defined as the ability of synthesizing the acquired knowledge through practice, training, or experience in order to improve the future behavior of the learning agent. In this section, we introduce a learning scheme [20, 21] aiming to understand the behavior of users during the interactions and to reach the eventual convergence towards Nash Equilibrium. The goal is that UAV’s service providers learn their own payoffs and determine the NE joint availabilitypricing strategy. To accomplish this, the best response algorithm is a suitable class; it leads to a Nash equilibrium status by exploiting fast monotony of best response functions. In fact, the NE concept is implicitly based on the assumption that players follow best response dynamics until they reach a state from which no player can improve his utility by changing unilaterally his strategy [22].
Algorithm 1 summarizes the best response dynamics that each service provider performs in order to converge to the game Nash equilibrium. The best response strategy of a player is defined as its optimal one, given the strategies of the other players. At a time iteration t, each service provider chooses the best joint availabilitypricing strategy against the opponents strategies chosen at the previous round . It is worth to notice that convergence of the best response dynamics is granted for smodular games [23].

The best response algorithm exhibits enormous advantages as it offers an accurate and fast convergence to the unique Nash equilibrium point. We say UAVs learn to play an equilibrium, if after a given number of iterations, the strategy profile converges to an equilibrium status.
6.2. Numerical Investigations
The pricingavailability best response learning results are provided using MATLAB software. Parameters used in the simulation are set based on typical values such as , and . Here, we will analyze the impact of the various parameters especially market temperature, coverage probability, and encounter rates on the Nash equilibrium price as well as the provided QoS (measured by beaconing periods) for both competitors.
Under this setup, Figure 3 illustrates the behavior of service providers under the best response algorithm and the convergence towards the Nash equilibrium price with time iterations. Here, we consider symmetric operators with randomly generated initial points. A smaller temperature μ value means that consumers are more and more rational, so that the operators’ NE price tends towards lower values. Indeed with , the NE price is absolutely zero . However, when μ becomes greater (consumers are less rational), the NE price arises and reaches its maximal value with . The same applies to Figure 4, where we depict the convergence of Algorithm 1 towards Nash equilibrium beaconing periods with different values of μ. It is worth noting that the best response algorithm converges accurately and rapidly to the Nash equilibrium point. This feature is required in efficiently designing such a mechanism in realistic systems. With , the is zero; however, is increasing with μ. By linking that with the NE price analysis, we can say that, with an increasing availability periods, the operators propose greater prices. This shows that higher quality of service (linked to network availability) incures additional cost.
We depict in Figures 5 and 6 the individual best response prices, respectively, beaconing periods for two asymmetric UAVs. Each UAV has a different coverage probability value (may be due to high parameter, for example). Counterintuitively, Figure 5 shows that the NE price does not depend on the coverage probability. Even with different coverage probabilities, the two competitors converge towards the same equilibrium price; this last value is increasing with UAV’s coverage which is quite intuitive.
Figure 6 depicts convergence to the equilibrium beaconing time under different values of coverage probabilities and different rational levels. It is worth noticing that, the UAV with greater coverage converges to a higher beaconing NE value, thus, offers a better QoS (availability service). We exclude from this analysis the case with a small value of the parameter μ, where the two operators converge to null values. Yet, the greater is the coverage (likely UAV high), the longer is the NE beaconing availability advertising periods.
Figure 7 illustrates the impact of both μ and the encounter rate on the NE price on each UAV i in an asymmetric setup (one UAV has an encounter rate advantage over the other). We remark that, even with different encounter rates, the NE prices are the same; however, the NE with a minimum temperature is greater than that of the symmetric case.
We plot in Figure 8 the beaconing best response learning evolution for two inhomogeneous UAVs (one UAV has an encounter rate advantage over the other), and the NE beaconing periods are straightforward affected by μ as well as λ. Thus, with a greater μ and λ, the best response converge to a high value of ; however, the NE beaconing period for a UAV with smaller λ is always lower.
In order to study the evolution impact of the parameter μ on the learning results, Figures 9 and 10 illustrate the NE price and NE beaconing period versus temperature parameter for two symmetric UAVs and different IoT population sizes, respectively. The results in Figure 9 show that with temperature parameter smaller than a certain threshold (for example, for ), the NE price is in its minimum value, and from this point, the price is increasing until the maximum price value () (for example, μ = 6 for N = 50). In this figure, we argue that, with a greater population density, the price tends to converge faster to the maximum value. These findings impact obviously the beaconing period evolution illustrated in Figure 10. With a greater population, the NE beaconing increases and reaches the equilibrium status in a smaller time window.
Pricing and beaconing learning evolution under different coverage probability values are illustrated in Figures 11 and 12, respectively. In Figure 11, the NE price is minimal with low coverage, and from a certain coverage probability value (nearly 0.6), the best strategy tends to . However, Figure 12 shows clearly the linear impact of on NE beaconing strategies. This last are increasing with , that confirms results in Figure 5.
Figures 13 and 14 illustrate the NE price with coverage probability and NE beaconing period for two symmetric UAVs and different IoT population sizes, respectively. These two figures are provided to study the impact of coverage probability on the learning results while considering . Figure 13 shows clearly the “natural” impact of population density on the service price, with a greater population that the operators can set greater prices. Furthermore, an interesting finding in Figure 14 is that, with a greater population, the NE beaconing duration is maximized, confirming thereby the results in Figure 13. It is worth mentioning how fast and accurate the proposed algorithm is for the convergence to the unique Nash equilibrium, which provides numerous insights into how to implement it and measure its effectiveness.
7. Conclusion
In this article, we deal with the pricingavailability interaction among adversarial unmanned aerial vehicles acting as flying base stations. We construct a theoretic framework based on the noncooperative game theory and characterize the equilibrium strategies for each UAV, both in terms of equilibrium fees and equilibrium availability probability. Furthermore, a special feature is obtained as the availability game with fixed price is submodular while the pricing game is supermodular under fixed availability. Next, we check that a simple iterative best responsebased algorithm allows exploring the unique Nash equilibrium of the game. Performance evaluation at equilibrium results allows UAV’s service providers to optimize their energy consumption while maximizing their monetary revenues.
As a future work, we are generalizing our proposal while considering heterogeneous mobility motions. Furthermore, some field experiments are also envisioned.
Data Availability
We hereby declare availability of all data related to our submission.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was conducted within the framework of Project funded by the Moroccan Ministry of Higher Education and Scientific Research and the National Centre for Scientific and Technical Research.