# Mobile Coverage in Rural Sweden: Analysis of a Comparative Measurement Campaign

Authors: P. G. Sudheesh; Jaap van de Beek
Journal: Mobile Information Systems
 (2021)

---

## Abstract

Under the umbrella of 1G to 5G, different technologies have been used to provide mobile communication. Various technologies are being proposed to bring a person in remote area under coverage. However, a statistical analysis on what these users get from already existing technologies has not been carried out. We fill this gap by carrying out such a study using a measurement campaign, where we present a framework for analyzing mobile signal strength experienced at the user end. Measurements are taken throughout the Norrbotten county, the northernmost county in Sweden, using mobile phones recording various parameters at regular intervals. Based on measured signal strength, a coverage map has been made via inverse distance weighting (IDW) interpolation. Based on the coverage map, various analyses are carried out on signal strength over residential areas and roads of Norrbotten. Overall, we lay a framework to analyze and quantify the effect of signal strength on users.


---

## Body

## 1. Introduction

Earlier this year, 3rd Generation Partnership Project (3GPP) released the first versions of its 5G new radio (NR) standard, first as a non-stand-alone mode, and later the stand-alone mode. For a while now, it has become clear that this standard essentially has its focus on the urban areas. Use cases and scenarios invariably address dense network environments, small cells, and a large number of mobile users per unit area.Naturally, concerns have been raised by stakeholders representing rural values, needs, and interests. With ever-decreasing cell sizes, how willareal coverage be guaranteed is based not only on the old standards and systems but also on new essential 5G functionalities as network slicing, ultralow latencies, and massive bandwidth [1].The first 5G testbed in Norrbotten, Sweden, was tested in June 2019 [2] and, therefore, for majority of population, 5G network is still not accessible. Also, 5G handsets are also not so widely available in the market. Hence, a major portion of the population is still using 2G, 3G, and/or 4G network for mobile communication. Therefore, expanding coverage area and bringing more people under mobile coverage map require network planning and installing new base stations.However, installing terrestrial base stations (BS) and providing backhaul are subject to the revenue earned from the mobile users [3]. Therefore, installing a terrestrial base station and associated backhaul in “not so popular” areas is not attractive from an economical point of view. This in turn results in deployment of lesser number of base stations or even absence of base stations in such sparsely populated areas.As a result, the people in these areas are subjected to poor or absence of coverage. It is widely known that larger operators and stakeholders of telecom industry do not want to put up new base stations in nonprofitable areas. In other words, the telecom companies hold the right not to put up new base stations in places that are not profitable to them [4]. It is rather interesting that the notion of full coverage is limited only to the areas defined by the operators.To overcome these limitations and provide full coverage irrespective of population density, different solutions like positioning UAV BS [5] and new allocation spectrum [6] are considered. Swedish regulators have come up with new 700 MHz band that focuses on lesser coverage areas [6]. To this end, Swedish authorities and especially Norrbotten, the northernmost county in Sweden, wanted to identify and focus on the areas with no or poor coverage [4]. Neither the coverage map from operators nor the location of base stations was available to researchers. To this end, we generate a coverage map generated using crowd sourced measurements across Norrbotten.It is found in [7] that mobile coverage is directly linked to socioeconomic performance. Therefore, in this paper, we measure the coverage poverty experienced by the people. Several attempts have been done to calculate areas without coverage [8, 9]. Effect of mobile coverage on railway network in Norway is studied in [10]. These attempts, however, do not really give an idea of quality of service experienced by the users. Therefore, to quantify the quality of service, the coverage poverty experienced by (a) people living in their apartments and (b) cars plying through the roads of Norrbotten province of Sweden is mathematically calculated.The paper is organized as follows. A review of techniques used to analyze effect of signal strength in various conditions is provided in Section2. A detailed explanation of generating coverage map from crowd sourced measurements is provided in Section 3. Section 4 is devoted to result analysis, while conclusions are listed in Section 5.


## 2. Coverage Maps from Scattered Samples

In this section, we elaborate the methods by which the signal strength over the populated areas and roads are analyzed using measured signal strengths. There are different ways of generating a coverage map and analyzing the signal strength: (1) identifying BS and then calculating the signal strength using factors such as path loss [11], (2) finding base station location from manually measured signal strength at different locations [12–14], and (3) generating coverage map from manually measured signal strength values. Since the locations of BS are not disclosed to public, we resort to the third method where we use tens of millions of measurements to generate coverage map. With millions of measurements, it is more likely to get a predicted signal strength than resorting to earlier methods. The signal strength measurements are taken for different operators and radio access technology (RAT) at difference time instances, which are then interpolated to form rectangular zones depicting average signal strength over that area. However, as we are only relying on available signal strength to generate coverage map, we miss lot of area as unmapped region.


### 2.1. Generating Coverage Map

The availability of measurements is limited to GPS coordinates on roads in general. This is because measurements are usually carried out by mobiles that are kept in the vehicles. Therefore, practically it is impossible to measure signal strength throughout the area to be mapped and identify exact signal strength at those locations. Therefore, we use the signal strength from the available signal strength measurements to calculate signal strength in the rest of the areas.In order to generate coverage map of various RATs and operators, it is important to separate the measured signals. The measurement from mobile phones contains network type info, operator, and signal strength, based on which it is possible to filter the point to corresponding sets that represent different RATs and operators. Therefore, for each operator, sets of different RATs such asSLTE,SUMTS,SGSM can be made. Sets SLTE, SUMTS, and SGSM represent the sets of measurements under LTE, UMTS, and GSM, respectively. The subscripts to S represent different RATs.Interpolation techniques such as kriging and inverse distance weighting [15] are widely used for obtaining coverage area map using measured signal strengths. Although different versions of kriging such as ordinary kriging and universal kriging are popular for crowd sourced measurements, a computationally less complex algorithm is used here as we consider large geographical area. In [16, 17], the authors compare kriging with IDW algorithm and find that kriging provides slightly better results compared to IDW. However, it is mentioned in [17] that the computational complexity for kriging is much higher, which is not a great concern when the area under study is small. But, as the area under study is large, computationally simple IDW algorithm is preferred over kriging to generate coverage maps for all operators and RATs.IDW algorithm predicts the signal strength in locations where measured points are not available. IDW algorithm makes use of the fact that measurements in the near vicinity contribute better than farther points. Therefore, bigger number of measured signal strength makes the signal strength prediction more closer to a possible real measurement. It is interesting to note that the IDW algorithm produces an averaging effect on measurements that are isolated. However, in case of crowded measurements, the signal strength assigned to the polygon is affected not only by the measurements in it but also by the measurements near to it.The output of IDW algorithm is a raster map, made up of a grid of points. These points are located at the center of each pixel, which is the basic building block of the raster. To generate a coverage map, a signal strength value is assigned to each point in the grid, based on the available signal strength values. The points in grid and measured point are marked as blue and red dots, respectively, in Figure1. Although the IDW algorithm does not limit the number of points that are used in order to calculate the signal strength of point in the grid, since the points in the vicinity contribute more to the measurement, we restrict the search to a circle, where each point in the grid is the center of the search circle.Figure 1
Selecting points to run IDW algorithm.The points in the grid are spaced at distanced, which is also the length of the edge of each raster (coverage map) cell. Using Tobler’s rule [18], the grid spacing d is taken as(1)d=map scale2×1000,where map scale depends on where the points under study are spread. In effect, the length of edges of raster is dependent on the distance of right-most, left-most, top-most, and bottom-most points among the available points.To generate the raster, each pixel or spatial object must hold the value of an element from a setM. Therefore, an object O must have two properties:(a)
A spatial extent,Zi(b)
A functionfO that assigns a value from set M to each point in ZiConsidering the search radiusr, object O is defined as(2)O=Zi,MiΨr.Thus, objectO is a tuple Zi,MiΨr, where Mi∈M and Zi⊆R2. The function fO:Zi⟶M, where M=R,“NOSIGNAL”. In effect, the value fOxi,yi reflects the signal strength to be attributed to xi,yi∈Zi. The distance, d, between neighbouring points in a raster, xi,yi, is obtained from equation (1). The signal strength Mi for each point in raster at location xi,yi is given as(3)Mi=∑xj,yj∈SSjxi−xj2+yi−yj2−αΨrxj,yj∑xj,yj∈Sxi−xj2+yi−yj2−α,(4)Ψrxj,yjfOxj,yj=fOxj,yj,ifxi−xj2+yi−yj2<r,0,ifx0−xj2+y0−yj2>r,whereSj∈S, S=R,“NOSIGNAL” and dij is the Euclidean distance between measured point j and point i in the grid. α denotes the distance exponent or power index, which is taken as 2. The characteristic function Ψr is defined as in equation (4), where r is the search radius. Hence, after successful implementation of (3) on SLTE, SUMTS, and SGSM, we get MLTE, MUMTS, and MGSM representing the coverage map for LTE, UMTS, and GSM, respectively.


### 2.2. Coverage Specific Aspects

To generate coverage map of Telia, Tele2, and Telenor, for each LTE, UMTS, and GSM, IDW algorithm must be carried out to each ensemble of measurements dedicated to each of the above operators and RATs. The IDW applied to each of these groups of measurements is explained in Section2.1. A simple IDW with all the measurements, where ”NO SIGNAL” is also an entry, will generate wrong results and therefore needs to be treated separately.Here, we formulate a coverage specific approach to generate coverage map, by taking into account both real signal measurements and points with no signal measurements. The algorithm is represented in Figure2. The first step is to separate points with valid signal measurements and points without signal strengths and generate coverage map for each set. These rasters or grids of points are then converted to polygons. Now, the difference between coverage maps with valid measurements and no signal strengths is created. This output is then merged to the coverage map with valid measurements.Figure 2
Algorithm to generate coverage map.It is still possible to have polygons with more points without signal strength and less points with valid measurements. To address such situations, a final step is required. That is, the polygon is assigned ”NO SIGNAL,” when the number of points without valid measurements is more even if there is a valid signal strength assigned to it. In the absence of such situations, the signal strength value assigned to the final coverage polygon remains the same as the one from the coverage map generated from valid measurements.


## 2.1. Generating Coverage Map

The availability of measurements is limited to GPS coordinates on roads in general. This is because measurements are usually carried out by mobiles that are kept in the vehicles. Therefore, practically it is impossible to measure signal strength throughout the area to be mapped and identify exact signal strength at those locations. Therefore, we use the signal strength from the available signal strength measurements to calculate signal strength in the rest of the areas.In order to generate coverage map of various RATs and operators, it is important to separate the measured signals. The measurement from mobile phones contains network type info, operator, and signal strength, based on which it is possible to filter the point to corresponding sets that represent different RATs and operators. Therefore, for each operator, sets of different RATs such asSLTE,SUMTS,SGSM can be made. Sets SLTE, SUMTS, and SGSM represent the sets of measurements under LTE, UMTS, and GSM, respectively. The subscripts to S represent different RATs.Interpolation techniques such as kriging and inverse distance weighting [15] are widely used for obtaining coverage area map using measured signal strengths. Although different versions of kriging such as ordinary kriging and universal kriging are popular for crowd sourced measurements, a computationally less complex algorithm is used here as we consider large geographical area. In [16, 17], the authors compare kriging with IDW algorithm and find that kriging provides slightly better results compared to IDW. However, it is mentioned in [17] that the computational complexity for kriging is much higher, which is not a great concern when the area under study is small. But, as the area under study is large, computationally simple IDW algorithm is preferred over kriging to generate coverage maps for all operators and RATs.IDW algorithm predicts the signal strength in locations where measured points are not available. IDW algorithm makes use of the fact that measurements in the near vicinity contribute better than farther points. Therefore, bigger number of measured signal strength makes the signal strength prediction more closer to a possible real measurement. It is interesting to note that the IDW algorithm produces an averaging effect on measurements that are isolated. However, in case of crowded measurements, the signal strength assigned to the polygon is affected not only by the measurements in it but also by the measurements near to it.The output of IDW algorithm is a raster map, made up of a grid of points. These points are located at the center of each pixel, which is the basic building block of the raster. To generate a coverage map, a signal strength value is assigned to each point in the grid, based on the available signal strength values. The points in grid and measured point are marked as blue and red dots, respectively, in Figure1. Although the IDW algorithm does not limit the number of points that are used in order to calculate the signal strength of point in the grid, since the points in the vicinity contribute more to the measurement, we restrict the search to a circle, where each point in the grid is the center of the search circle.Figure 1
Selecting points to run IDW algorithm.The points in the grid are spaced at distanced, which is also the length of the edge of each raster (coverage map) cell. Using Tobler’s rule [18], the grid spacing d is taken as(1)d=map scale2×1000,where map scale depends on where the points under study are spread. In effect, the length of edges of raster is dependent on the distance of right-most, left-most, top-most, and bottom-most points among the available points.To generate the raster, each pixel or spatial object must hold the value of an element from a setM. Therefore, an object O must have two properties:(a)
A spatial extent,Zi(b)
A functionfO that assigns a value from set M to each point in ZiConsidering the search radiusr, object O is defined as(2)O=Zi,MiΨr.Thus, objectO is a tuple Zi,MiΨr, where Mi∈M and Zi⊆R2. The function fO:Zi⟶M, where M=R,“NOSIGNAL”. In effect, the value fOxi,yi reflects the signal strength to be attributed to xi,yi∈Zi. The distance, d, between neighbouring points in a raster, xi,yi, is obtained from equation (1). The signal strength Mi for each point in raster at location xi,yi is given as(3)Mi=∑xj,yj∈SSjxi−xj2+yi−yj2−αΨrxj,yj∑xj,yj∈Sxi−xj2+yi−yj2−α,(4)Ψrxj,yjfOxj,yj=fOxj,yj,ifxi−xj2+yi−yj2<r,0,ifx0−xj2+y0−yj2>r,whereSj∈S, S=R,“NOSIGNAL” and dij is the Euclidean distance between measured point j and point i in the grid. α denotes the distance exponent or power index, which is taken as 2. The characteristic function Ψr is defined as in equation (4), where r is the search radius. Hence, after successful implementation of (3) on SLTE, SUMTS, and SGSM, we get MLTE, MUMTS, and MGSM representing the coverage map for LTE, UMTS, and GSM, respectively.


## 2.2. Coverage Specific Aspects

To generate coverage map of Telia, Tele2, and Telenor, for each LTE, UMTS, and GSM, IDW algorithm must be carried out to each ensemble of measurements dedicated to each of the above operators and RATs. The IDW applied to each of these groups of measurements is explained in Section2.1. A simple IDW with all the measurements, where ”NO SIGNAL” is also an entry, will generate wrong results and therefore needs to be treated separately.Here, we formulate a coverage specific approach to generate coverage map, by taking into account both real signal measurements and points with no signal measurements. The algorithm is represented in Figure2. The first step is to separate points with valid signal measurements and points without signal strengths and generate coverage map for each set. These rasters or grids of points are then converted to polygons. Now, the difference between coverage maps with valid measurements and no signal strengths is created. This output is then merged to the coverage map with valid measurements.Figure 2
Algorithm to generate coverage map.It is still possible to have polygons with more points without signal strength and less points with valid measurements. To address such situations, a final step is required. That is, the polygon is assigned ”NO SIGNAL,” when the number of points without valid measurements is more even if there is a valid signal strength assigned to it. In the absence of such situations, the signal strength value assigned to the final coverage polygon remains the same as the one from the coverage map generated from valid measurements.


## 3. Practical Case Study

### 3.1. Measurements

In order to generate raster, we measure the signal strengths of various operators and RATs with GPS coordinates and time stamps. The measurement is made by using set of mobile phones, which are locked to specific operators and RATs. The set of mobile phones is kept in a bag as shown in Figure3, where each phone is locked to a single operator and RAT. To this end, three operators in Sweden are compared, where Tele2 and Telenor share the same infrastructure for 2G (GSM) and 4G (LTE). Therefore, each bag contains 2 phones to measure 2G signal and 3 phones to measure 3G signal and again 2 phones to measure 4G signal, which records measurements at regular intervals.Figure 3
Measurement setup by IQMTEL, Sweden.We consider three operators, Telia, Tele2, and Telenor, over the Norrbotten region, which is the northernmost province in Sweden. The mobile phones that measure signal strengths are enclosed in a bag and are placed in vehicles and boats plying through the area under study. The mobile phones collect data, namely, GPS coordinates, date, time, network type, signal strength, operator, cell identification code (CID), location area code (LAC), and eNode bid at regular intervals. It is worth noting that the measurements from various phones in the bag are not time-synchronized.The two GSM phones measure broadcast channel (BCCH) and the signal strength measured in GSM is the received signal strength indicator (RSSI). The 3G RAT used in the phones is universal mobile telecommunications system (UMTS), which is a wideband code division multiple access (WCDMA) system. Therefore, the mobile phones dedicated to measuring 3G signals measure the received power on one code over the primary common pilot channel (CPICH). This value is referred to as the received signal code power (RSCP). The 4G RAT is also known as the long-term evolution (LTE) standard, which uses time-frequency resource block with varying channel size (from 1.4 to 20 MHz). Therefore, the useful power is contributed by different resource blocks, which is measured as reference signal received power (RSRP).The 2G and 3G phones are locked to 2G and 3G RAT, respectively. This means that the mobile phone cannot perform handover operation from 3G to 2G when it has poor signal strength. However, the 4G mobile phones are not locked to RATs and as a result they switch RAT from LTE to UMTS and then from UMTS to GSM based on the measured signal strength. In other words, the phones that measure LTE signals provide RSRP values in fewer areas, even though there are RSRP values when phone switches to other RATs from LTE.The number of measurements received per municipality is given in Table1. In this paper, we study the distribution of signal strength over residential areas and roads. The population and road lengths associated with each of the municipalities are listed in the table. The number of measurements per municipality consists of three RATs and three operators. The population map considered for analysis is a 1 km × 1 km polygon shape file provided by [19]. The road data is given by [20], which provides information about number of vehicles moving through that particular road.Table 1
Measurements across municipalities of Norrbotten.
SI. no.MunicipalityNo. of measurementsArea (km2)PopulationRoad (km)1Älvsbyn7081911790.48256530.12Arjeplog88755814498.32903720.13Arvidsjaur18267546085.764901019.94Boden24092194277.128373875.95Haparanda37887016806.79851284.66Jokkmokk2089506933.25132796.37Kalix94006819314.416240652.78Luleå46032231850.3781607679Pajala182162220642.361571004.510Piteå38839382155.841932919.611Gällivare30706802913.918023976.612Kiruna23158002515.723233730.113Överkalix3251238103.73409521.714Övertorneå6712103246.14582556.1


### 3.2. Framework to Analyze the Quality of Service

With the measurements in hand, this section lays the framework to analyze the distribution of signal strength over residential areas and roads. For this purpose, we use two sets of polygon files, population and road, as a base layer that extracts signal strength from the superimposed coverage map. While the first helps us to study the effect in static users, the latter provides insight into the coverage for mobile users. Further, we define different levels from LTE, UMTS, and GSM signal strength, with which the quality of service enjoyed by the users can be easily identified. The algorithm used for classification of user levels is mentioned in Figure4. The fact that LTE, UMTS, and GSM offer better service in the respective order is used for forming the algorithm. While RSRP >−80 dBm is an excellent signal strength in LTE, −90 dBm < RSRP < −80 dBm offers poorer performance than the former [21]. Therefore, we assign the former as level 1 and the latter as level 2. Following the pattern, we assign higher levels to the LTE and lower levels to UMTS [21] as the data rate decreases in UMTS compared to LTE. This pattern is repeated in GSM as well [22].Figure 4
Assigning signal levels.Note that the points labelledNo Signal and NULL are not considered while assigning levels till level 12, where No Signal represents a situation where UEs did not receive any signal, whereas NULL considers a scenario where polygon does not have a signal strength to map to. Level 13 is used to represent the signal strength that is lesser than the minimum required signal strength for a reliable communication. A scenario where the UE cannot make call in any of the RAT is considered in level 14. The difference between level 13 and level 14 is that connection is not possible in level 14, while connection may be possible with possible chance of disconnection in level 13 [21, 23]. For example, if the device under test has LTE signal strength below −100 dBm, UMTS signal strength below −85 dBm, and GSM signal strength below −110 dBm, the signal strength values in mobile phones fall below the minimum recommended signal strength [21, 23]. Such a situation is possible in level 13. Meanwhile, if all devices are not getting any signal from the base station, we assign level 14, which is the worst situation among all. Level 15 does not represent the quality of the signal, as it houses the regions where the measurement campaign was not performed.


#### 3.2.1. Framework to Analyze Effects in Residential Areas

In this subsection, we study the effect of signal strength on users in residential areas. Figure5 demonstrates the extraction of signal strength from the coverage map. The upper layer is the coverage map, which holds the signal strength value. The lower layer, which is population map in this contest, extracts signal strength from the overlapping coverage map and adds the signal strength value Mk to each polygon of the population or road map.Figure 5
Extracting signal strength from coverage map.In order to extractMk from an overlapping coverage map, each polygon uses the settings in Table 2, which is executed using QGIS, a free geographic information system (GIS) software. That is, the polygon in population or road map considers those signal strengths from the overlapping coverage map, when the polygon in coverage map intersects with it. In case the population or road map polygon intersects with multiple polygons of the coverage map, the algorithm selects one random value among them, thereby creating a one-to-one mapping with each polygon in population or road map and signal strength value.Table 2
Feature extraction.
SI. no.OptionSetting in algorithm1Geometric predicateIntersect2Mapping (join type)One-to-one3Fields to addSignal strength


#### 3.2.2. Framework to Analyze Effects in Roads

Since the roads can stretch to few kilometers and assigning signal strength to each line is impossible, we split the lines representing roads to segments of 100 m or less, for absorbing information from overlapping coverage map. For example, a road whose length is 1.050 km is divided into 10 units, each of 100 m and one unit of length 50 m. The contribution of the unit of lengthx, where x<100m, is taken as x/100 while calculating total number of segments in a particular level. Signal strength value is extracted from overlapping coverage map of each RAT, SLTE, SUMTS, and SGSM, and is assigned to each line segment. As a result, in addition to signal strength value, the line segment object also carries the length of each segment. The line segment extracts signal strength value from overlapped coverage map M with the settings shown in Table 2. As a result, each line segment object q is a tuple Zq,Mq,lq, where lq is the length of the line segment q.Furthermore, to analyze the performance of various operators to users on road, we define two variables,λq and ϕq. The basic logic here is that the coverage poverty experienced by 100 cars in 1 km is the same as that of 1 car in 100 km. Therefore, to quantify the coverage poverty, we multiply the number of cars moving through segment q with lq. Also, to incorporate smaller road segments, we normalize the product by length. Therefore, we define λ for segment q:(5)λq=number of cars throughq×lq100.Since different roads have speed limits, we must add this factor into consideration, as the vehicle stays on road for more time in poor coverage areas. Therefore, we define another parameter:(6)ϕq=λq×36001000,where the numerator specifies the speed in km/hr and to convert it to m/s and we multiply the term by 3.6. These two variables quantify the coverage poverty of absence of coverage to the users.


## 3.1. Measurements

In order to generate raster, we measure the signal strengths of various operators and RATs with GPS coordinates and time stamps. The measurement is made by using set of mobile phones, which are locked to specific operators and RATs. The set of mobile phones is kept in a bag as shown in Figure3, where each phone is locked to a single operator and RAT. To this end, three operators in Sweden are compared, where Tele2 and Telenor share the same infrastructure for 2G (GSM) and 4G (LTE). Therefore, each bag contains 2 phones to measure 2G signal and 3 phones to measure 3G signal and again 2 phones to measure 4G signal, which records measurements at regular intervals.Figure 3
Measurement setup by IQMTEL, Sweden.We consider three operators, Telia, Tele2, and Telenor, over the Norrbotten region, which is the northernmost province in Sweden. The mobile phones that measure signal strengths are enclosed in a bag and are placed in vehicles and boats plying through the area under study. The mobile phones collect data, namely, GPS coordinates, date, time, network type, signal strength, operator, cell identification code (CID), location area code (LAC), and eNode bid at regular intervals. It is worth noting that the measurements from various phones in the bag are not time-synchronized.The two GSM phones measure broadcast channel (BCCH) and the signal strength measured in GSM is the received signal strength indicator (RSSI). The 3G RAT used in the phones is universal mobile telecommunications system (UMTS), which is a wideband code division multiple access (WCDMA) system. Therefore, the mobile phones dedicated to measuring 3G signals measure the received power on one code over the primary common pilot channel (CPICH). This value is referred to as the received signal code power (RSCP). The 4G RAT is also known as the long-term evolution (LTE) standard, which uses time-frequency resource block with varying channel size (from 1.4 to 20 MHz). Therefore, the useful power is contributed by different resource blocks, which is measured as reference signal received power (RSRP).The 2G and 3G phones are locked to 2G and 3G RAT, respectively. This means that the mobile phone cannot perform handover operation from 3G to 2G when it has poor signal strength. However, the 4G mobile phones are not locked to RATs and as a result they switch RAT from LTE to UMTS and then from UMTS to GSM based on the measured signal strength. In other words, the phones that measure LTE signals provide RSRP values in fewer areas, even though there are RSRP values when phone switches to other RATs from LTE.The number of measurements received per municipality is given in Table1. In this paper, we study the distribution of signal strength over residential areas and roads. The population and road lengths associated with each of the municipalities are listed in the table. The number of measurements per municipality consists of three RATs and three operators. The population map considered for analysis is a 1 km × 1 km polygon shape file provided by [19]. The road data is given by [20], which provides information about number of vehicles moving through that particular road.Table 1
Measurements across municipalities of Norrbotten.
SI. no.MunicipalityNo. of measurementsArea (km2)PopulationRoad (km)1Älvsbyn7081911790.48256530.12Arjeplog88755814498.32903720.13Arvidsjaur18267546085.764901019.94Boden24092194277.128373875.95Haparanda37887016806.79851284.66Jokkmokk2089506933.25132796.37Kalix94006819314.416240652.78Luleå46032231850.3781607679Pajala182162220642.361571004.510Piteå38839382155.841932919.611Gällivare30706802913.918023976.612Kiruna23158002515.723233730.113Överkalix3251238103.73409521.714Övertorneå6712103246.14582556.1


## 3.2. Framework to Analyze the Quality of Service

With the measurements in hand, this section lays the framework to analyze the distribution of signal strength over residential areas and roads. For this purpose, we use two sets of polygon files, population and road, as a base layer that extracts signal strength from the superimposed coverage map. While the first helps us to study the effect in static users, the latter provides insight into the coverage for mobile users. Further, we define different levels from LTE, UMTS, and GSM signal strength, with which the quality of service enjoyed by the users can be easily identified. The algorithm used for classification of user levels is mentioned in Figure4. The fact that LTE, UMTS, and GSM offer better service in the respective order is used for forming the algorithm. While RSRP >−80 dBm is an excellent signal strength in LTE, −90 dBm < RSRP < −80 dBm offers poorer performance than the former [21]. Therefore, we assign the former as level 1 and the latter as level 2. Following the pattern, we assign higher levels to the LTE and lower levels to UMTS [21] as the data rate decreases in UMTS compared to LTE. This pattern is repeated in GSM as well [22].Figure 4
Assigning signal levels.Note that the points labelledNo Signal and NULL are not considered while assigning levels till level 12, where No Signal represents a situation where UEs did not receive any signal, whereas NULL considers a scenario where polygon does not have a signal strength to map to. Level 13 is used to represent the signal strength that is lesser than the minimum required signal strength for a reliable communication. A scenario where the UE cannot make call in any of the RAT is considered in level 14. The difference between level 13 and level 14 is that connection is not possible in level 14, while connection may be possible with possible chance of disconnection in level 13 [21, 23]. For example, if the device under test has LTE signal strength below −100 dBm, UMTS signal strength below −85 dBm, and GSM signal strength below −110 dBm, the signal strength values in mobile phones fall below the minimum recommended signal strength [21, 23]. Such a situation is possible in level 13. Meanwhile, if all devices are not getting any signal from the base station, we assign level 14, which is the worst situation among all. Level 15 does not represent the quality of the signal, as it houses the regions where the measurement campaign was not performed.


### 3.2.1. Framework to Analyze Effects in Residential Areas

In this subsection, we study the effect of signal strength on users in residential areas. Figure5 demonstrates the extraction of signal strength from the coverage map. The upper layer is the coverage map, which holds the signal strength value. The lower layer, which is population map in this contest, extracts signal strength from the overlapping coverage map and adds the signal strength value Mk to each polygon of the population or road map.Figure 5
Extracting signal strength from coverage map.In order to extractMk from an overlapping coverage map, each polygon uses the settings in Table 2, which is executed using QGIS, a free geographic information system (GIS) software. That is, the polygon in population or road map considers those signal strengths from the overlapping coverage map, when the polygon in coverage map intersects with it. In case the population or road map polygon intersects with multiple polygons of the coverage map, the algorithm selects one random value among them, thereby creating a one-to-one mapping with each polygon in population or road map and signal strength value.Table 2
Feature extraction.
SI. no.OptionSetting in algorithm1Geometric predicateIntersect2Mapping (join type)One-to-one3Fields to addSignal strength


### 3.2.2. Framework to Analyze Effects in Roads

Since the roads can stretch to few kilometers and assigning signal strength to each line is impossible, we split the lines representing roads to segments of 100 m or less, for absorbing information from overlapping coverage map. For example, a road whose length is 1.050 km is divided into 10 units, each of 100 m and one unit of length 50 m. The contribution of the unit of lengthx, where x<100m, is taken as x/100 while calculating total number of segments in a particular level. Signal strength value is extracted from overlapping coverage map of each RAT, SLTE, SUMTS, and SGSM, and is assigned to each line segment. As a result, in addition to signal strength value, the line segment object also carries the length of each segment. The line segment extracts signal strength value from overlapped coverage map M with the settings shown in Table 2. As a result, each line segment object q is a tuple Zq,Mq,lq, where lq is the length of the line segment q.Furthermore, to analyze the performance of various operators to users on road, we define two variables,λq and ϕq. The basic logic here is that the coverage poverty experienced by 100 cars in 1 km is the same as that of 1 car in 100 km. Therefore, to quantify the coverage poverty, we multiply the number of cars moving through segment q with lq. Also, to incorporate smaller road segments, we normalize the product by length. Therefore, we define λ for segment q:(5)λq=number of cars throughq×lq100.Since different roads have speed limits, we must add this factor into consideration, as the vehicle stays on road for more time in poor coverage areas. Therefore, we define another parameter:(6)ϕq=λq×36001000,where the numerator specifies the speed in km/hr and to convert it to m/s and we multiply the term by 3.6. These two variables quantify the coverage poverty of absence of coverage to the users.


## 3.2.1. Framework to Analyze Effects in Residential Areas

In this subsection, we study the effect of signal strength on users in residential areas. Figure5 demonstrates the extraction of signal strength from the coverage map. The upper layer is the coverage map, which holds the signal strength value. The lower layer, which is population map in this contest, extracts signal strength from the overlapping coverage map and adds the signal strength value Mk to each polygon of the population or road map.Figure 5
Extracting signal strength from coverage map.In order to extractMk from an overlapping coverage map, each polygon uses the settings in Table 2, which is executed using QGIS, a free geographic information system (GIS) software. That is, the polygon in population or road map considers those signal strengths from the overlapping coverage map, when the polygon in coverage map intersects with it. In case the population or road map polygon intersects with multiple polygons of the coverage map, the algorithm selects one random value among them, thereby creating a one-to-one mapping with each polygon in population or road map and signal strength value.Table 2
Feature extraction.
SI. no.OptionSetting in algorithm1Geometric predicateIntersect2Mapping (join type)One-to-one3Fields to addSignal strength


## 3.2.2. Framework to Analyze Effects in Roads

Since the roads can stretch to few kilometers and assigning signal strength to each line is impossible, we split the lines representing roads to segments of 100 m or less, for absorbing information from overlapping coverage map. For example, a road whose length is 1.050 km is divided into 10 units, each of 100 m and one unit of length 50 m. The contribution of the unit of lengthx, where x<100m, is taken as x/100 while calculating total number of segments in a particular level. Signal strength value is extracted from overlapping coverage map of each RAT, SLTE, SUMTS, and SGSM, and is assigned to each line segment. As a result, in addition to signal strength value, the line segment object also carries the length of each segment. The line segment extracts signal strength value from overlapped coverage map M with the settings shown in Table 2. As a result, each line segment object q is a tuple Zq,Mq,lq, where lq is the length of the line segment q.Furthermore, to analyze the performance of various operators to users on road, we define two variables,λq and ϕq. The basic logic here is that the coverage poverty experienced by 100 cars in 1 km is the same as that of 1 car in 100 km. Therefore, to quantify the coverage poverty, we multiply the number of cars moving through segment q with lq. Also, to incorporate smaller road segments, we normalize the product by length. Therefore, we define λ for segment q:(5)λq=number of cars throughq×lq100.Since different roads have speed limits, we must add this factor into consideration, as the vehicle stays on road for more time in poor coverage areas. Therefore, we define another parameter:(6)ϕq=λq×36001000,where the numerator specifies the speed in km/hr and to convert it to m/s and we multiply the term by 3.6. These two variables quantify the coverage poverty of absence of coverage to the users.


## 4. Data Analysis

### 4.1. Signal Strength across Residential Areas in Municipalities

Histogram of signal strength is analyzed for various RATs and operators. For the ease of comparison, we consider LTE by Telia in Lulea and Pajala, a comparatively densely populated municipality and a sparsely populated municipality, respectively. Further, the histogram of GSM signal strength is plotted and compared with that of LTE.Figure6 shows histogram of GSM signal strength, where the height of each bar shows the percentage of people falling in the respective slot. Both municipalities have maximum number of people in −70 dBm to −60 dBm range for GSM. Also, the range is the same for most of the municipalities.Figure 6
GSM coverage in residential areas. (a) Coverage in residential areas of Lulea. (b) 2G Coverage in residential areas of Pajala.
(a)(b)From Figure7, we can see that the behaviour in LTE is similar to that of GSM, but majority of people fall in −100 dBm to −90 dBm range. Most of the municipalities show similar behaviour. It is visible that there is a shift in mostly used signal strength in GSM and LTE. That is, from Figure 4, it is known that −70 dBm to −60 dBm range in RSSI is an excellent GSM signal strength, and −100 dBm to −90 dBm LTE signal strength offers only fair service to the users. On the contrary, the quality of service offered by LTE for signal strength from −100 dBm to −90 dBm is still better than GSM signal strength of −70 dBm to −60 dBm.Figure 7
LTE coverage in residential areas. (a) 4G coverage in residential areas of Lulea. (b) 4G Coverage in residential areas of Pajala.
(a)(b)


### 4.2. Signal Strength across Roads in Municipalities

Figure8 shows the histogram of signal strength for Telia GSM. For fair comparison of histograms generated for residential areas, we consider the same municipalities, Lulea and Pajala. Unlike residential areas, where both municipalities have maximum number of people in −70 dBm to −60 dBm range for GSM, roads of these municipalities have maximum number of people in −90 dBm to −80 dBm and −80 dBm to −70 dBm ranges, respectively. It is observed that this shift is common for most of the municipalities. Another interesting observation is that the sparsely populated municipality, Pajala, has better signal strength in most of its roads, despite the low population density.Figure 8
GSM coverage on roads. (a) Coverage on roads of Lulea. (b) Coverage on roads of Pajala.
(a)(b)The LTE signal strength in roads of Pajala is reported in Figure9. Compared to Figure 8, a decrease in signal strength is observed. Although there is a dip in average signal strength in LTE, the performance at the user end is still better than that of GSM. Unlike the histograms of residential areas, where the absence of measurement (“No Data”) is almost zero, the percentage of roads without measurement is bigger. This is mainly due to the mode of measurement, where the measurements are taken by phones placed in mostly public utility vehicles.Figure 9
LTE coverage on roads. (a) Coverage on roads of Lulea. (b) Coverage on roads of Pajala.
(a)(b)


### 4.3. Effect of Traffic Pattern

In this subsection, we quantify the effect of 3G signal strength on users on roads. Figure10(a) plots histogram of λ, where λ=λq, for various 3G operators in Jokkmokk. With this approach, the best operator can be identified. Similarly, Figure 10(b) plots histogram of ϕ, where ϕ=ϕq, for various 3G operators in Jokkmokk. While the latter gives the idea about temporal signal quality experienced by an average user in car, the former shows the number of cars experiencing the signal strength. However, both figures have more or less the same pattern.Figure 10
Effect of traffic pattern. (a) Cars× 100 m. (b) Cars × sec.
(a)(b)


### 4.4. Quantifying the Coverage Poverty

It is quite clear that different municipalities get affected in different levels. A low signal strength in highly populated municipality will have huge impact compared to a less populated area. Figure11 analyzes average GSM signal strength of Telia on total population of municipalities. While the position of circles represents the mean signal strength and percentage of people without service, the radius of the circles in scatter plot shows the population or total road length in respective municipalities. Figure 11(a) reports the percentage of population without signal versus mean signal strength with GSM. It is clear that most of the users in residential areas enjoy fair-quality GSM signal. Most of the populated municipalities have decent signal quality. Also, a very small percentage (less than 1.5) of residential area is not covered. However, on analyzing signal strength on roads, in Figure 11(b), a decline in signal strength is observed. Although majority of municipalities get most of their roads under coverage area, there are municipalities like Övertorneå, where considerable fractions of roads are not covered.Figure 11
Impact on % of roads/population in Telia. (a) Residence area. (b) Roads.
(a)(b)


### 4.5. Signal Strength in Urban, Suburban, and Rest of the Areas

The effect of signal strength to users in residential areas in urban, suburban, and rest of the areas is reported in Figure12. The term rest of the areas is used to denote residential areas that do not fall under urban or suburban category. The CDF of population is plotted against GSM signal strength. While GSM offers best signal strength in urban areas, the medium and least signal strengths are observed in suburban areas and rest of the areas, respectively. A fraction of population is under no coverage area of Tele2 and Telenor in suburban areas and rest of the areas. However, despite the region classification, Telia offers 1–4 dB better signal to users. It is also important to note that the results in Figure 12 are based only on the available data. We have not considered areas without measurements for analysis.Figure 12
Signal strength in urban, suburban, and rest of the areas.


### 4.6. Comparison with Another County

GSM signal strengths of two different counties in Sweden, Norrbotten and Östergötland, are compared in Figure13. We compare the signal strengths observed in urban and suburban areas of both counties. From Figure 13(a), it is clear that, despite of the county, the signal strength in urban areas remains the same. However, from 13(b), it can be observed that the user gets 5 to 10 dB signal strength reduction, when the user moves from urban to suburban area.Figure 13
Coverage in Norrbotten and Östergötland. (a) Östergötland (urban). (b) Norrbotten (suburban).
(a)(b)


### 4.7. Effect of Roaming

It is observed that some areas near Finland are not getting service by Tele2 and Telenor. Although these areas are in Sweden, Tele2 and Telenor switch the call to another operator from Finland, forming large areas without measurements from Tele2 and Telenor. Telia, on the other hand, offers connectivity from Sweden. Figure14 shows the histogram of 3G signal strength in roads of Övertorneå. It is seen that 70% of roads and 60.2% of roads from Tele2 and Telenor are under ”No Data” bin. This is because there is no measurement from Tele2 and Telenor in majority of roads and Tele2 and Telenor transfer the call to Elisa, a Finnish operator.Figure 14
Histogram of 3G signal strength in roads of Övertorneå.


### 4.8. Correlation with 700 MHz Spectrum’s Focus Areas

Recently 700 MHz spectrum has been allocated to certain areas of focus [6]. Figure 15 presents the correlation of Telia’s coverage map developed using algorithm presented in Figure 2 and focus areas of new 700 MHz spectrum, represented by rectangles. Different colours are assigned to each of the locations based on the signal quality, from level 1, which is the best, to level 14, which is the worst. A major portion of the focus area does not intersect with our map due to lack of measurements and is shown in dashed lines. On comparing the focus areas provided with our coverage map, we found that 44.1% of the areas that are intersected with our map experience poor signal (below level 10). Around 54.7% of areas that are intersected with our map experience below-average signal (level 8–9).Figure 15
Correlation with 700 MHz spectrum’s focus area map.


## 4.1. Signal Strength across Residential Areas in Municipalities

Histogram of signal strength is analyzed for various RATs and operators. For the ease of comparison, we consider LTE by Telia in Lulea and Pajala, a comparatively densely populated municipality and a sparsely populated municipality, respectively. Further, the histogram of GSM signal strength is plotted and compared with that of LTE.Figure6 shows histogram of GSM signal strength, where the height of each bar shows the percentage of people falling in the respective slot. Both municipalities have maximum number of people in −70 dBm to −60 dBm range for GSM. Also, the range is the same for most of the municipalities.Figure 6
GSM coverage in residential areas. (a) Coverage in residential areas of Lulea. (b) 2G Coverage in residential areas of Pajala.
(a)(b)From Figure7, we can see that the behaviour in LTE is similar to that of GSM, but majority of people fall in −100 dBm to −90 dBm range. Most of the municipalities show similar behaviour. It is visible that there is a shift in mostly used signal strength in GSM and LTE. That is, from Figure 4, it is known that −70 dBm to −60 dBm range in RSSI is an excellent GSM signal strength, and −100 dBm to −90 dBm LTE signal strength offers only fair service to the users. On the contrary, the quality of service offered by LTE for signal strength from −100 dBm to −90 dBm is still better than GSM signal strength of −70 dBm to −60 dBm.Figure 7
LTE coverage in residential areas. (a) 4G coverage in residential areas of Lulea. (b) 4G Coverage in residential areas of Pajala.
(a)(b)


## 4.2. Signal Strength across Roads in Municipalities

Figure8 shows the histogram of signal strength for Telia GSM. For fair comparison of histograms generated for residential areas, we consider the same municipalities, Lulea and Pajala. Unlike residential areas, where both municipalities have maximum number of people in −70 dBm to −60 dBm range for GSM, roads of these municipalities have maximum number of people in −90 dBm to −80 dBm and −80 dBm to −70 dBm ranges, respectively. It is observed that this shift is common for most of the municipalities. Another interesting observation is that the sparsely populated municipality, Pajala, has better signal strength in most of its roads, despite the low population density.Figure 8
GSM coverage on roads. (a) Coverage on roads of Lulea. (b) Coverage on roads of Pajala.
(a)(b)The LTE signal strength in roads of Pajala is reported in Figure9. Compared to Figure 8, a decrease in signal strength is observed. Although there is a dip in average signal strength in LTE, the performance at the user end is still better than that of GSM. Unlike the histograms of residential areas, where the absence of measurement (“No Data”) is almost zero, the percentage of roads without measurement is bigger. This is mainly due to the mode of measurement, where the measurements are taken by phones placed in mostly public utility vehicles.Figure 9
LTE coverage on roads. (a) Coverage on roads of Lulea. (b) Coverage on roads of Pajala.
(a)(b)


## 4.3. Effect of Traffic Pattern

In this subsection, we quantify the effect of 3G signal strength on users on roads. Figure10(a) plots histogram of λ, where λ=λq, for various 3G operators in Jokkmokk. With this approach, the best operator can be identified. Similarly, Figure 10(b) plots histogram of ϕ, where ϕ=ϕq, for various 3G operators in Jokkmokk. While the latter gives the idea about temporal signal quality experienced by an average user in car, the former shows the number of cars experiencing the signal strength. However, both figures have more or less the same pattern.Figure 10
Effect of traffic pattern. (a) Cars× 100 m. (b) Cars × sec.
(a)(b)


## 4.4. Quantifying the Coverage Poverty

It is quite clear that different municipalities get affected in different levels. A low signal strength in highly populated municipality will have huge impact compared to a less populated area. Figure11 analyzes average GSM signal strength of Telia on total population of municipalities. While the position of circles represents the mean signal strength and percentage of people without service, the radius of the circles in scatter plot shows the population or total road length in respective municipalities. Figure 11(a) reports the percentage of population without signal versus mean signal strength with GSM. It is clear that most of the users in residential areas enjoy fair-quality GSM signal. Most of the populated municipalities have decent signal quality. Also, a very small percentage (less than 1.5) of residential area is not covered. However, on analyzing signal strength on roads, in Figure 11(b), a decline in signal strength is observed. Although majority of municipalities get most of their roads under coverage area, there are municipalities like Övertorneå, where considerable fractions of roads are not covered.Figure 11
Impact on % of roads/population in Telia. (a) Residence area. (b) Roads.
(a)(b)


## 4.5. Signal Strength in Urban, Suburban, and Rest of the Areas

The effect of signal strength to users in residential areas in urban, suburban, and rest of the areas is reported in Figure12. The term rest of the areas is used to denote residential areas that do not fall under urban or suburban category. The CDF of population is plotted against GSM signal strength. While GSM offers best signal strength in urban areas, the medium and least signal strengths are observed in suburban areas and rest of the areas, respectively. A fraction of population is under no coverage area of Tele2 and Telenor in suburban areas and rest of the areas. However, despite the region classification, Telia offers 1–4 dB better signal to users. It is also important to note that the results in Figure 12 are based only on the available data. We have not considered areas without measurements for analysis.Figure 12
Signal strength in urban, suburban, and rest of the areas.


## 4.6. Comparison with Another County

GSM signal strengths of two different counties in Sweden, Norrbotten and Östergötland, are compared in Figure13. We compare the signal strengths observed in urban and suburban areas of both counties. From Figure 13(a), it is clear that, despite of the county, the signal strength in urban areas remains the same. However, from 13(b), it can be observed that the user gets 5 to 10 dB signal strength reduction, when the user moves from urban to suburban area.Figure 13
Coverage in Norrbotten and Östergötland. (a) Östergötland (urban). (b) Norrbotten (suburban).
(a)(b)


## 4.7. Effect of Roaming

It is observed that some areas near Finland are not getting service by Tele2 and Telenor. Although these areas are in Sweden, Tele2 and Telenor switch the call to another operator from Finland, forming large areas without measurements from Tele2 and Telenor. Telia, on the other hand, offers connectivity from Sweden. Figure14 shows the histogram of 3G signal strength in roads of Övertorneå. It is seen that 70% of roads and 60.2% of roads from Tele2 and Telenor are under ”No Data” bin. This is because there is no measurement from Tele2 and Telenor in majority of roads and Tele2 and Telenor transfer the call to Elisa, a Finnish operator.Figure 14
Histogram of 3G signal strength in roads of Övertorneå.


## 4.8. Correlation with 700 MHz Spectrum’s Focus Areas

Recently 700 MHz spectrum has been allocated to certain areas of focus [6]. Figure 15 presents the correlation of Telia’s coverage map developed using algorithm presented in Figure 2 and focus areas of new 700 MHz spectrum, represented by rectangles. Different colours are assigned to each of the locations based on the signal quality, from level 1, which is the best, to level 14, which is the worst. A major portion of the focus area does not intersect with our map due to lack of measurements and is shown in dashed lines. On comparing the focus areas provided with our coverage map, we found that 44.1% of the areas that are intersected with our map experience poor signal (below level 10). Around 54.7% of areas that are intersected with our map experience below-average signal (level 8–9).Figure 15
Correlation with 700 MHz spectrum’s focus area map.


## 5. Conclusions

This paper has presented a framework for analyzing mobile signal strength experienced by users. Based on measured signal strength, a coverage map has been made via IDW interpolation. Various analyses are carried out on signal strength over residential areas and roads of Norrbotten. Further, measurements are compared to those of Östergötland and it was found that both municipalities have almost similar measurements.By analyzing coverage across all 14 municipalities of Norrbotten, in contrast to the suspicion that rural areas have poor signal strength, we found that 2G and 4G provide satisfactory results. However, 3G fails to provide coverage in some areas. This is worse in some areas, resulting in more than 50% of areas to be outside coverage area at some places. These areas are mostly near the Finland border, which results in the fact that roaming and additional chargesmay be applicable to the 3G user in these areas.



---

*Source: 8869534-2021-01-05.xml*