Mapping Rural Mobility in the Niakhar Area (Fatick Region, Senegal)

Despite the development of information and communication technologies, mobility still remains an important factor in the development of the social capital, by promoting social interactions. However, rural African populations face more difficulties getting around than city dwellers. In addition, studies on mobility in Africa are relatively weak. In rural African areas, this ignorance of mobility is even more important. It is necessary to identify the characteristics and the spatial structure of mobility flows before starting any reflection on solutions, particularly in terms of spatial allocation of transport resources. Our main objective is to map the current movement flows in the Niakhar area. We chose to use the cartographic approach to analyze mobility flows. The data come from a survey carried out in 2018, as part of the French Institute of Research for Development (IRD) Niakhar Atlas project. For the data analysis, we used 4 methods: dominant flux method, net flux calculation method, directional mean statistics and areal interpolation. The results revealed differences in mobility behavior depending on the season and in space, and the existence of three dominant centers: Ngayokhème, Diohine and Toucar. They also made it possible to highlight the main directions of the flows and the scalar fields of mobility.


Introduction
The development of means of transport has increased the mobility capacity of populations and exchanges at the global level [1]. At the same time, the unequal distribution of resources (infrastructure, equipment and basic services) in the areas has increased the mobility needs of populations [2]. Despite the development of information and communication technologies, mobility still remains an important factor in the development of social capital, by promoting social contacts [3]. It allows people to access the equipment, services and jobs available in their area. As a structuring factor of territories, mobility makes it possible to account for discriminating processes, the relationship of populations to space and to reveal the level of socio-spatial integration.
However, studies on mobility in Africa are relatively weak. In rural African areas, that weakness is even greater. However, rural populations face more difficulties getting around than city dwellers. These mobility constraints are linked to various factors, including isolation, the state of communication channels, the lack of organized transport system, the often high cost of transport and the very low motorization of households. Those obstacles, combined with the prevailing poverty in these rural areas [4], considerably hinder the mobility capacity of rural populations. It is important to know the characteristics and the spatial structure of mobility flows before starting any reflection on solutions, particularly in terms of spatial allocation of transport resources. Indeed, the Niakhar area is a heterogeneous space from all points of view [5] and this heterogeneity can generate differences in mobility behavior.
Several studies have been carried out [6] [7] [8] [9] [10] to analyze the mobility. Admittedly, these studies are important, since they have made it possible to understand the socio-economic behavior of mobility and to reveal the spatial distribution of populations' movements, but they are insufficient. Indeed, they do not sufficiently take into account the spatial dimension, which is nevertheless very important for understanding mobility, the spatial allocation of transport equipment and services, and the fight against spatial disparities. The scale of data analysis is often arbitrary. The zoning used does not capture the process underlying the data, since the zones are often different in size and shape. So the results may simply be zoning artefacts. This makes their use for demand modeling and transport planning problematic. Moreover, concerning the specific case of the Niakhar area, so far, there is no study concerning the mapping of mobility flows.
In view of all these observations, we have chosen to use the cartographic approach to analyze rural mobility flows in the Niakhar area.
Our main objective is to map the current mobility flows in the Niakhar area.
This involves in particular revealing the dominant flows and centers, the attractiveness and flow fields. To achieve this goal, we have followed four steps. The first consisted in determining the dominant flows and centers, the second in measuring attractiveness. In the third, we measured the polarization. The fourth step allowed us to determine the flow fields.

A. Ndonky et al. Journal of Geographic Information System
In this study, we are interested in rural mobility, considered as the set of current trips made in rural areas to carry out activities. Indeed, everyday travel is necessary to accomplish social (wedding/baptism, visits from friends/relatives…) and economic (shopping/sales, work…) activities. This type of travel is becoming more and more important with the development of production tracks and weekly markets in rural areas, particularly that of Niakhar. In addition, it should be noted that long-term displacements in the Niakhar zone are outside the Niakhar zone, that is to say in other regions of Senegal, in particular that of Dakar. However, we are interested in movement within the Niakhar area. Mobility in our discussion associates all the travel flows of individuals during the last month preceding the day of the survey, the reasons for these trips and the means of transport used.
The data come from a survey carried out in 2018, as part of the Niakhar Atlas project of the IRD. For data analysis, we used 4 methods. These are firstly the dominant flux method proposed by [11] and the net flux calculation method [12], used respectively to measure dominant fluxes and attractiveness. Next, the directional mean statistic [13] [14] was used to measure flux polarization (flux vector field). Finally, the use of area interpolation [15] made it possible to determine the scalar flux field.   Figure 1). Sahelo-Sudanese type, its climate is characterized by the alternation of two seasons (rainy season and dry season). In 2013, its population was 44,994 inhabitants, with a density of 221.6 inhabitants/km 2 [16]. This population is very young, because in 2012, young people under the age of 20 represented 53% of the total.

Data
The data used come from the survey on rural mobility carried out by IRD in 2018. Our survey plan was done taking into account the differences in social and cultural characteristics of the population of the study area and integration of the specific geographic (spatial) dimension of this area. In fact, we can note a social and spatial heterogeneity of the Niakhar zone and an interlocking of administrative meshes. In total, two criteria (constraints) guided the methodological choices: respect for the socio-demographic diversity of the Niakhar area, taking into account the heterogeneity of the rural space. Only the use of the finest mesh proposed by the administrative grid of the Niakhar zone (hamlet allowing an Journal of Geographic Information System exhaustive statistical treatment of the Niakhar zone) makes it possible to better understand this heterogeneity. The hamlet is the smallest administrative division unit in the Niakhar zone located below the village mesh. It is often homogeneous from a socio-demographic and spatial point of view, due to its small size.
The number of people to be surveyed is based mainly on the proportion of people aged 6 and over moving outside their village of residence. As this proportion is not known, to calculate the size of the sample, we will use the following formula: with: n = sample size of working people; D = sounding effect fixed at 2.5; e = Margin of error fixed at 5%; t2: 95% confidence level corresponding to 1.96.
The calculation gives a minimum total sample size of 884 individuals. To increase the power of the test, we increased the total sample size to 900.
The target population includes individuals aged 6 and over. Given the absence of an updated sampling frame, this could not be identified by random selection before the survey. Faced with this situation, we carried out a multi-stage cluster survey. Each hamlet is defined as a cluster. Thus, a random selection at the following degrees: hamlet, household and individuals, was carried out according to rigorous procedures and guaranteeing the random nature of the selection. In the first degree, we randomly drew the hamlets, 90 in total. In each hamlet, 10 households were selected and in each household, only one individual was chosen for the survey. Figure 2 illustrates the spatial distribution of the selected hamlets. Journal of Geographic Information System

Data Processing and Analysis Strategy
The statistical processing of the data collected is done using CS Pro and Stata software. Geoprocessing, mapping and spatial analysis were done by ArcGIS and Quantum GIS software.
Our analysis includes 4 major steps: determination of the dominant flows and centers, measurement of attractiveness, mapping of the vector fields of observed flows and finally mapping of the scalar fields of estimated flows.
To determine the dominant fluxes and centers, the dominant flux method proposed by [11] was adopted. This method consists in representing at most one flow at the start of each spatial entity and offers the possibility of locating the main direction of the flows towards the dominant (center) zone. A zone A is said to be "dominant" of zone B when A receives the greatest flux from B and A receives more flux than B [17]. A zone A is said to be "dominated" by zone B when A transmits the greatest flow towards B and B receives more flow than A.
Thus the use of flow classes (thresholds) becomes more relevant than proportional variations to make the flow map more readable. The great difficulty with which this method is confronted is the choice of the threshold beyond which the flow must be represented. The solution proposed by [18] is to take the mean of the flows as the threshold value. Another solution, proposed by [19], is to choose Journal of Geographic Information System a certain percentage of dominant flows. Due to the strong asymmetry of the data in favor of low flows, we adopted the solution proposed by [19], by retaining a threshold of 5%.
However, this method of dominant flows does not allow the measurement of attractiveness, which is nevertheless very essential for the description of mobility flows in the space.
To measure attractiveness, an attractiveness index was calculated for each village by making, for each of them, the balance/total volume of flows ratio. In order to summarize the information in a simple and efficient way, we have combined on a single map the total volume of flows exchanged and the attractiveness index. We can then locate the most important poles, attractive centers and repulsive places in space. Admittedly, the two methods make it possible to describe the mobility flows, but they do not offer the possibility of apprehending the direction of the mobility flows, which is however very important for characterizing the movements in the space. The use of directional mean statistics made it possible to measure the polarization of mobility flows. This method is based on Tobler's idea according to which each flow can be represented by a vector linking an origin point and a destination point and comprising a norm proportional to the number of flows exchanged [20]. The vector sum of the flows originating in a place produces a vector with an angle that corresponds to the average direction of these flows. The norm of this vector gives additional information on the intensity of the spatial orientation and the polarization of the fluxes emitted by the mean direction. The equation for the Linear Mean Directional (LMD) statistic is as follows: with θ i = the directions of a series of polylines from the same origin [21].
That statistics is used to calculate the circular variance which indicates the degree of variation of the individual directions around the mean direction. Its value is between 0 and 1. The smaller it is (close to 0), the more the input vectors have very similar directions and the higher the level of polarization. When it is larger (close to 1), then the input vectors have very different directions and the polarization level low. The circular variance equation (CirVar) is as follows: with n = sample size [21].
We used the inverse of this variance to calculate the polarization index. The higher the index value, the stronger the polarization. So reading it will be easier. The methods used so far have produced interesting results. However, they do not allow the representation of scalar flux fields, which are very important in the analysis of movements in the space.
To map the scalar flux fields, we often use point interpolation methods. We can note, for example, the weighted average (or inverse distance) method, trend surface analysis, kriging [22] [23]. Certainly these methods are interesting for representing the specific phenomenon on a continuous surface. However, they have limitations with regard to the zonal layout data, as is the case in this study. For all these reasons, we chose to use the area interpolation method. The latter makes it possible, from the zonal implantation data, to produce a prediction surface (smoothed data) and errors, as well as to re-aggregate the data at the scale of the target polygons [15].
The overall approach for choosing the methods used is summarized in Figure  3. This figure highlights the complementarity of the methods used. So, the combination of all these methods allows to produce rich and complementary results.

Types of Centers and Dominant Flows
The map of the dominant wintering centers and flows ( Figure 4) is produced        Figure 9 and Figure 10 show the vector fields of movements observed respectively in wintering and in the dry season. In general, we can observe that Toucar, Ngayokhème and Diohine, the largest villages in the area, polarize the most displacement flows within the Niakhar area. The results also show that the flows are more polarized in wintering than in the dry season ( Figure 10 and Figure 11). The most polarized flows come from the villages located on the outskirts of the area. We can note, for example, the flows coming from the villages of Godel, Khassous, Sob, Kothiokh. Those flows are all the more polarized as we are in wintering. On the other hand, flows from villages located in the center of the zone are less polarized.

Scalar Fields of Estimated Emitted Flows
To check the quality of the prediction, we use variographic analysis and cross-validation statistics. Variographic analysis involves changing the parameters of the model until 90% of the empirical covariances (blue crosses) fall within the confidence interval (in red) ( Figure 12). Thus in the dry season, as in wintering, 91.12% of these covariances fall within the confidence interval ( Figure 12 and Figure 13), which shows the good quality of the prediction (well-adjusted interpolation model).
The cross-validation statistics of the spatial interpolation model give a weighted root mean square close to 1.2 in the dry season and 1.1 in wintering, values close to 1 (ideal value for a well-fitted model). This confirms the results of the variographic analysis.
Since the predicted (estimated) data is valid, we will now present the results. Figure 14 and Figure 15 show the scalar fields of the emitted flows. They

Discussion
Our study revealed the main poles, secondary poles of attraction and repulsive zones, the main sources of emission of intra-zone travel flows in Niakhar. Our results made it possible to highlight the principal directions and the scalar fields of flow (estimated flows in a continuous space). This type of study is a first at the scale of this area, making it difficult to compare its results with those of similar studies in the same area. Nevertheless, we can compare its results with those of studies done elsewhere.
Like other works carried out in various contexts, such as that of [17] on Cameroon, Belgium and Czechoslovakia, by [3] on the agglomeration of Dakar (Senegal), our study highlighted the poles of attraction, the sources of emission flows and the main directions of mobility flows. Unlike the results of the study by [3] which showed a very strong polarization of mobility flows, this study revealed a weaker polarization. This difference is explained by the difference in geographic contexts. Indeed, the agglomeration of Dakar is a peninsula where travel possibilities are limited in certain directions. In addition, this agglomera-tion is strongly dominated by a city center which polarizes all trips. On the other hand, the Niakhar area is an open rural environment where the mobility possibilities exist in all directions. In addition, there is not a hyper dominant center that polarizes most of the movement flows.
The main poles of attraction are Toucar, Diohine and Ngayokhème. Those localities also constitute the largest villages which concentrate more population, infrastructure, equipment and services, as well as certain administrative functions such as the Ngayokhème village, capital of the commune of the same name [24]. We can therefore say that the geography of mobility flows is determined by the spatial distribution of infrastructure, basic equipment, economic and social activities, as well as administrative functions.
The results revealed differences in travel practices in the space depending on the season: the volume of travel is greater in the dry season than in wintering and the flows more polarized in wintering than in the dry season. The greater volume of travel in the dry season can be explained by the fact that during this period the populations are done with their field work and have more time to move. Moreover, mobility conditions are easier in the dry season than in the rainy season. In fact, during the rainy season, some localities are cut off from the rest of the area due to the presence of water on the traffic axes. In addition, some roads are in very poor condition. The strongest polarization in wintering may be due to the fact that the possibilities of movement are more limited, because of the difficulties, or even the impossibility of moving on certain traffic axes during this period.
Certainly, we can observe that the level of polarization of travel flows is greater in the outlying villages and weaker in the villages located in the center. This difference should be put into perspective. It is largely linked to the geographical location of these villages. In fact, the villages located on the outskirts emit part of the flows to the neighborhood located outside the study area, but these were not taken into account because of methodological choices (only intra-zones are analyzed). On the other hand, in the villages located in the center, all the flows of movements in the neighborhood were taken into account, whatever the direction.
The advantage of the dominant flow method is that it allows the synthesis of flows. But it has limits, because the dominant flows depend both on the spatial structure of the flows and on the division of the zones, and this can cause difficulties for the comparison of the results. The area interpolation method has the advantage of being able to reveal the spatial variations of the flows within the villages, of being free from the centroids of the irregular meshes to produce the continuous surface. On the other hand, this method eliminates abrupt variations (local pockets of high value of the phenomenon inside the mesh), since it is based on the continuous progression of the density of the workforce [25].
Our study presents some originality compared to those carried out so far in the study area. Indeed, it made it allowed us to map the main mobility flows and to highlight the poles of attraction from a large volume of flows. The use of di-rectional statistics helped us to measure the levels of polarization of flows and thus identify the main directions of intra-zone movements. The representation of flows in the form of a continuous surface, from data of areas of irregular shapes and different sizes (villages), also constitutes another originality, because it produced results revealing spatial disparities in mobility, at a spatial resolution very fine.
With the demographic dynamics characterized by the strong growth of the rural population and the neoliberal policies, the spatial disparities of resources become increasingly large even in rural areas. In the specific case of the Niakhar area, large villages tend to be well equipped with basic infrastructure, services and equipment, while smaller localities are under-equipped [24]. Consequently, the inhabitants of the smallest villages are often far from areas of concentration of infrastructure, equipment and services, and forced to travel long distances in order to benefit from the resources located in the dominant villages. This is what the results of our study reveal. Rural mobility is therefore indicative of the functional division of space and resource allocation strategies in rural space.
The spatial characteristics of mobility serve as a basis for defining the spatial allocation of transport infrastructure and equipment [26]. Hence the importance of knowing these characteristics before initiating any reflection on the spatial allocation of transport resources. From this point of view, the results of our study are of definite use, since they can help us highlight the poles of attraction, the sources of emission of the flows, the main directions of movement, as well as the spatial variations in the distribution of emitted fluxes at very fine spatial resolution. Thus we can contribute to the fight against socio-spatial inequalities in mobility practices and access to resources in the rural area.

Conclusions
The cross-study of dominant centers and flows, vector fields and scalar fields of mobility flows, combining the dominant flow method, directional statistics and the area interpolation method confirms the relevance of our methodological choices, because it has made it possible to better reveal the spatial structure of intra-rural mobility flows, in particular the poles of attraction, the main directions and the scalar fields of flow. This is a significant advance compared to studies carried out so far in the same area.
Our results raise the debate on the relationships between spatial resource allocation policies and spatial segregation in rural areas. With strong demographic growth in rural areas and the lack of regulation of the spatial allocation of basic equipment and infrastructure, rural areas are increasingly segregated, the need for mobility is growing. In this context, spatial and social disparities will tend to increase, to the detriment in particular of populations located in small, isolated and poor villages. We are therefore in a situation which is all the more worrying because the means of collective transport are limited and the public transport systems are poorly organized in rural areas.

A. Ndonky et al. Journal of Geographic Information System
The production of these results is an encouraging first step for the analysis of the spatial structure of intra-rural mobility flows. But improvements need to be made to the methods used to refine the results. Thus, for future studies, socio-demographic characteristics, the location of basic infrastructure and equipment should be taken into account and put in relation to mobility flows. It is also important, to improve the measurement of the scalar field of mobility flows, to introduce a control variable, such as the level of attractiveness measured from the density of infrastructure, services and basic equipment.
Since our approach has been validated in the case of Niakhar, it could, for comparison purposes, be applied to other rural areas in Senegal or Africa.