Understanding of Spatial Urban Growth Pattern and Associated Landscape Dynamics in Congolese Mining Cities Bordering Zambia from 1990 to 2023

1. Introduction

Global urbanization is a multifaceted process influenced by demographic, economic, and social factors [1]. By 2018, over half of the world's population resided in urban areas, a figure projected to reach 68% by 2050 [2]. Urbanization rates vary globally, with rapid acceleration observed in developing regions such as Africa, fueled by migration and population growth [3]. While Africa remains predominantly rural, it is expected to undergo significant urbanization by 2050 [4,5], with countries like Nigeria, the Democratic Republic of the Congo (DR Congo), Ethiopia, and Tanzania projected to contribute significantly to this global urban population increase by 2050 [6,7].

Rapid urbanization in Africa induces profound changes in urban landscapes [8] and presents challenges such as housing shortages and infrastructure strain, underscoring the need for effective urban planning and governance [9]. This is particularly crucial in small and medium-sized border cities, necessitating tailored policies for sustainable development [10]. These cities, often rife with uncertainties and complexities, concentrate extreme poverty, particularly among those with limited options [11,12]. Located strategically near national borders, they experience rapid population growth due to cross-border activities and economic potential, requiring careful consideration and strategic interventions [13]. The unplanned expansion of border cities presents significant urban planning challenges, surpassing local authorities' planning capacity, leading to unsustainable urban sprawl affecting land use, infrastructure, and the environment [14].

In Central-Southern Africa, the Katangese Copperbelt Area (KCA), situated in southeastern DR Congo bordering Zambia, exemplifies this trend. Cities in this area reflect shared histories and cultural exchanges driven by mining activities and trade. Proximity to national borders transforms these cities into dynamic economic hubs fostering cross-border trade and investment [15,16]. Migration, primarily motivated by employment opportunities, significantly contributes to rapid urbanization and demographic growth, intensifying the demand for housing and infrastructure [17,18]. Efficient transport infrastructure plays a crucial role in supporting commercial activities and regional connectivity, while industrial policies influence growth through the establishment of industries and economic zones, driving urban expansion and economic development [19,20]. However, these cities face urban planning and governance challenges, leading to rapid, unregulated urbanization resulting in a lack of green spaces, and inadequate infrastructure and services [21]. This issue is particularly acute in the DR Congo due to rapid population growth, poor governance, and widespread poverty.

Congolese border cities like Kipushi, Kasumbalesa, Mokambo, and Sakania confront notable challenges concerning urban development and governance, notably housing provision. Mining activities historically shaped Sakania and Kipushi, attracting workers from the city of Lubumbashi (the provincial capital) and beyond post-2002 mining sector liberalization [22,23]. Kipushi developed residential areas for Congolese army personnel, while Kasumbalesa, pivotal for cross-border exchanges, expanded without planning due to recent demographic explosions. Mokambo's growth stems from the road network and a border post, fostering economic opportunities. These agglomerations expand spatially but face largely unplanned, informal urbanization, marked by land use conflicts [24,25]. Social ghettos exacerbated social disparities, and visual landscape degradation lead to biodiversity loss, as evidenced in the city of Lubumbashi, where 3.6 km² of green space is lost annually [26]. Yet, urban trees, green belts, and peri-urban forests are essential for sustainable development.

Systematic mapping and quantitative assessment of spatial dynamics in Congolese border cities are essential for understanding urban development implications. Remote sensing and landscape ecology techniques offer a robust framework for monitoring spatio-temporal landscape changes [27,28]. Additionally, landscape ecology explores landscape structure and ecosystem interactions [29]. The synergy between these approaches fosters a comprehensive understanding of border urbanization and associated landscape dynamics, facilitating evidence-based policies for sustainable growth and resource management [30].

The objective of this study is to conduct an in-depth spatial analysis of border urban dynamics and associated landscape changes, between the DR Congo and Zambia. The hypothesis being tested is that due to rapid population growth and economic opportunities, border towns are experiencing accelerated urban expansion, particularly in medium-sized cities with dynamic economies. This disrupts landscapes by increasing heterogeneity and promotes the expansion of savannas due to forest fragmentation.

2. Materials and Methods

2.1. Study Area: Congolese Cities Bordering Zambia

The cities of Kipushi, Kasumbalesa, Mokambo, and Sakania are strategically located in the southeastern region of the DR Congo (Figure 1), within the Katangese Copperbelt Area (KCA) [31].

Situated between 10° and 12° S and 26° and 29° E, this region falls under the Cw climate classification, characterized by a distinct rainy season from November to March and a dry season from May to September, with transitional months in April and October [32]. The region receives an annual precipitation of approximately 1200 mm and has an average annual temperature of 20 °C [32]. The topography is predominantly flat with scattered hills typical of the Katanga landscape, and the soils are primarily ferralsol-type [33]. Rapid urbanization and human activities have fragmented the original miombo woodland, transitioning it into savannahs [34-36]. These cities are experiencing continuous population growth, with estimates of 396,000 in Kipushi, 700,000 in Kasumbalesa, 47,000 in Mokambo, and 55,000 in Sakania [37]. Economic activities in this region are diverse, focusing on mining, agriculture, and cross-border trade, particularly significant due to proximity to the Zambian border. Trade between the DR Congo and Zambia is dominated by the exchange of minerals, especially copper and cobalt. The DR Congo mainly exports these minerals to Zambia, while Zambia provides the DR Congo with manufactured and agricultural products, which transit through these border cities. Kasumbalesa serves as a major border city with established customs infrastructure; Mokambo features a dry port and critical road connections; Sakania is influenced by its mining history and industrial potential, while Kipushi is recognized for its mining legacy and ongoing urban revitalization efforts.

2.2. Methodology

2.2.1. Data

Landsat images with a 30-meter spatial resolution and less than 5% cloud cover, spanning from 1990 to 2023, were utilized for this study. The selected intervals (1990-2000, 2000-2010, and 2010-2023) enable a decadal analysis, capturing long-term trends and spatial dynamics changes. The 1990-2000 period was characterized by political turbulence and conflicts in the Eastern region. The 2000-2010 era witnessed mining sector liberalization (2002), the first electoral cycle (2006), subsequent infrastructure modernization, and the global financial crisis (2008). The 2010-2023 period saw additional electoral cycles (2011 and 2018), provincial restructuring (2015), and a change in political regime (2019). The selection of Landsat images specifically during the dry season minimizes cloud cover, ensuring clear visibility for accurate interpretation of landscape features [23]. This approach is crucial for precise mapping and quantification of spatio-temporal landscape changes [38].

2.2.2. Classification

Using the WGS-84 reference ellipsoid, the Landsat images were georeferenced in the UTM Zone 35S coordinate system, which corresponds to the study region [23]. Additionally, radiometric normalization was applied to minimize reflectance variations between Landsat image dates. Global corrections for atmospheric effects were also performed using the dark object subtraction method [39]. Consecutively, a false-color composite of selected Landsat images was meticulously constructed by combining the mid-infrared, near-infrared, and red bands, to enhance discrimination among diverse vegetation types [40,41]. For precision and clarity, distinctive land cover units were methodically identified and assigned unique codes across different scenes. To establish a solid foundation for subsequent analyses, Regions of Interest (ROIs), representing training areas, were meticulously delineated for each study year, during the dry season. The selection of ROIs was guided by sampling polygons, strategically positioned to avoid transition zones, thereby minimizing the impact of the mixel phenomenon, and enhancing subsequent analysis accuracy [42,43]. To refine analytical capabilities, these carefully crafted ROIs were utilized to construct a comprehensive model for training the Random Forest classifier under Google Earth Engine. This ensemble approach, incorporating multiple decision trees, provided a robust and adaptable foundation for subsequent land cover classification [44].

The methodology for Landsat image classification, using Google Earth Engine (GEE) and the Random Forest algorithm, represented a rigorous approach implemented to characterize five land cover types: forest, savanna, agriculture, built-up & bare soil, and other land cover (Table 1).

To evaluate the accuracy of the obtained classifications, we followed the best practices recommendations of Olofsson et al. [45]. An unbiased surface estimators and estimated uncertainty was constructed by collecting a sample of reference observations from change maps between 1990 and 2023. This process relied on truth points collected in each land cover class. Samples were randomly stratified according to a 9-strata map for each period, including 5 stable strata (forest, savannah, agriculture, and built-up & bare soil) and 4 relevant change strata for each period (forest lost, savannah gain, agriculture gain, and built-up & bare soil gain). The sample size was determined using Cochran's method [46], with 800 points sampled for each period (1990-2000, 2000-2010, and 2010-2023). Based on the proportion of each stratum, 250 points were assigned to strata occupying more than 40%, 150 points to strata occupying between 10 and 40%, and 100 points to strata occupying less than or equal to 10%. Subsequently, QGIS 3.26.1 software was used to calculate the error matrix, expressed in terms of estimated surface proportions [43]. Measurement accuracies, including overall accuracies, user's and producer's accuracies, were also automatically generated using the same software.

2.2.3. Quantifying Urban Landscape Pattern Changes

To quantify human impact on landscape morphology, several metrics have been calculated [47,48]. First, the class area which refer to the relative extent of specific land cover types within a defined landscape, helping in understanding the composition of the landscape by identifying the predominant land cover matrix. The number of patches plays a crucial role in assessing landscape fragmentation since a high number of patches indicates fragmentation and scattered distribution, whereas a low number suggests infilling or aggregation of patches. Furthermore, Shannon's diversity was calculated to measure the landscape heterogeneity [49], while a disturbance index, defined as the ratio between the cumulative area of anthropogenic land cover in the landscape and the forest area [50], was calculated to assess the level of landscape anthropization. However, urban sprawl was quantified using the Urban Expansion Intensity Index (UEII), which quantifies the rate and magnitude of urban growth within a specified area over time. This index provides a comprehensive view of urbanization intensity and its impacts [51]. Additionally, the Largest Patch Index (LPI), defined as the ratio between the largest patch area and the class area while providing information on the fragmentation of a land cover following its reduction [49], was also calculated. Finally, the average Euclidian distance to the nearest neighbor was calculated. This metric indicates the average distance between each patch in the landscape and its closest neighboring point, offering an understanding of spatial pattern dispersion [43].

Between two dates, changes in patch number and class area indicate distinct spatial transformation processes in landscape dynamics, between two dates [40,47,52]: Decreases in both patch number and class area indicate attrition, while increased class area with decreased patch number suggests aggregation. Unchanged patch numbers with increased class area signify enlargement, whereas growth in both class area and patch number reflects the creation of new patches. Dissection is marked by reduced class area and increased patch number, often with linear disruptions causing minimal area loss. Fragmentation, on the other hand, involves patch increase accompanied by significant class area loss. To differentiate between fragmentation and dissection, the ratio of total areas at different time points was examined, with a ratio above 0.75 indicating dissection dominance and a ratio at or below 0.75 suggesting prevalent fragmentation [53].

3. Results

3.1. Classification Accuracy and Mapping

Table 2 presents the accuracy performance of supervised classifications of Landsat images, using the Random Forest classifier, for the years 1990, 2000, 2010, and 2023. The overall accuracy exceeds 90 % for each analyzed period, underscoring the reliability in distinguishing between different land cover types. The user's and producer's accuracy ranging from 93 % to 100 %, further confirm the high quality of the classifications, indicating minimal errors in classifying the various land cover types. Additionally, applying a 95 % confidence interval to estimate the stratified area of each land cover class across the different periods reveals a margin of error below 5 %. This low uncertainty enhances the credibility of the results, suggesting that the area estimates for each class are reliable and precise.

The visual analysis reveals the growth of built-up and bare soil is observed in different directions: eastward in the city of Kipushi, northward and southwestward in the city of Kasumbalesa, eastward in the city of Mokambo (C), and westward in the city of Sakania (D). Additionally, an important decrease in forest cover, particularly around the city of Kipushi (A) and the city of Kasumbalesa (B), were observed. Concurrently, savannas are gradually encroaching upon former forested areas around all cities studied, indicating a shift in ecosystem dynamics. There is also a noticeable increase in agricultural lands surrounding all four cities, underscoring the significant impacts of human activities on the landscape. The reduction in forest cover, coupled with the expansion of built-up and bare soil areas, agriculture, and savannas, highlights the extensive footprint of deforestation and urbanization on local landscapes.

Figure 2. Mapping spatial land cover dynamics in Kipushi (A), Kasumbalesa (B), Mokambo (C), and Sakania (D) landscapes from 1990 to 2023 using supervised classification of Landsat images with the Random Forest classifier.

Figure 2. Mapping spatial land cover dynamics in Kipushi (A), Kasumbalesa (B), Mokambo (C), and Sakania (D) landscapes from 1990 to 2023 using supervised classification of Landsat images with the Random Forest classifier.

3.2. Landscape Composition Dynamics

Between 1990 and 2023, there was a marked decrease in forest cover surrounding Kipushi, Kasumbalesa, Mokambo, and Sakania. Initially, forests dominated these regions, but by 2023, their extent had drastically diminished, with less than 15% remaining in Kipushi, 10% in Kasumbalesa, 30% in Mokambo, and 15% in Sakania. Conversely, savannas multiplied fourfold on average, built-up and bare soil increased eightfold, and agriculture surged by a factor of 46 across the four border cities. These changes signify a significant reshaping of the land cover in the region (Figure 3). Furthermore, an analysis using the Shannon Diversity Index indicated shifts in landscape composition, with Kasumbalesa, Mokambo, and Sakania displaying increased diversity between 1990 and 2023, suggesting a transition toward more heterogeneous landscapes. This transformation underscores the diversification of land cover characteristics and types within these border cities (Figure 4).

3.3. Intensity of Urban Expansion, Sprawl and Associated Landscape Anthropization

During the 1990-2000 decade, the urban expansion intensity index for Kipushi, Kasumbalesa, Mokambo, and Sakania was 0.08, 0.14, 0.02, and 0.03, respectively. These results indicate that Kasumbalesa experienced the most significant urban growth intensity during this period. From 2000 to 2010, there was a substantial increase in the urban expansion intensity index for all border cities studied, with values of 0.32, 0.29, 0.04, and 0.21 for Kipushi, Kasumbalesa, Mokambo, and Sakania, respectively. Compared to the previous decade, Kipushi saw a more pronounced urban expansion intensity during this period. Lastly, from 2010 to 2023, the urbanization intensity index continued to rise, reaching values of 0.54, 0.12, and 0.21 for Kasumbalesa, Mokambo, and Sakania, respectively. For the city of Kipushi, the value decreased fourfold during this period. These results indicate a more important urban expansion intensity in Kasumbalesa, placing it at the forefront in terms of spatial urban growth compared to the other studied cities (Figure 5).

However, the border cities of Kipushi, Kasumbalesa, Mokambo, and Sakania, along with their surrounding areas, experienced a progressive anthropization of their landscapes (Table 3). In 1990, the largest forested patch covered approximately 80% of the class area, with an average distance between patches of less than 70 meters. However, this largest patch area decreased by approximately 10%, while the average distance between forest patches multiplied nearly threefold by 2023. Concurrently, an inverse trend was observed for agriculture, savannas, as well as built-up and bare soil. Over the 33 years from 1990 and 2023, the average distance between their patches decreased, while the size of their respective largest patches increased (Table 3), suggesting a trend towards fragmentation and degradation of landscapes.

Additionally, during the period 1990-2023, a significant increase in disturbance index was observed within the landscapes of the studied border cities (Figure 6). In Kipushi, this index multiplied by 14, increasing from 0.5 in 1990 to 6.5 in 2023. In Kasumbalesa, an even more pronounced increase was noted, with a multiplication factor of 53, raising the index from 0.2 in 1990 to 11.8 in 2023. For Mokambo, the index multiplied by 10, increasing from 0.3 in 1990 to 2.8 in 2023. Lastly, in Sakania, the index multiplied by 37, reaching 7 in 2023 compared to an initial value of 0.2 in 1990 (Figure 6). These increases reflect a significant intensification of human activity in these regions over this 33-year period.

3.4. Analysis of Landscape Spatial Pattern Dynamics

Applying the decision tree algorithm of Bogaert et al. (2014) on the data in Table 4, our findings underscore that between 1990 and 2000, the forest experienced an increase in PN and a decrease in CA for the cities of Kipushi, Kasumbalesa, Mokambo, and Sakania, suggesting a spatial transformation process of dissection (tobs0.83>0,75). Conversely, the savannas, agriculture, and built-up and bare soil recorded an increase in both PN and CA, indicating a creation process.

During the period from 2000 to 2010, fragmentation was the dominant spatial transformation process for the forest for all studied border cities (tobs0.58<0,75), since the decrease of CA was accompanied by an increase in PN. Savannas exhibited an increase in CA and a decrease in PN, suggesting an aggregation-type spatial transformation process. Agriculture, and built-up and bare soil were characterized by an increase in PN and CA, indicating creation as the spatial transformation process. From 2010 to 2023, the forest showed an increase in PN and a decrease in CA for the border cities of Kipushi, Mokambo, and Sakania, indicating a process of fragmentation (tobs 0.38<0,75). However, around the border city of Kasumbalesa, there was a simultaneous decrease in both PN and CA of forest, indicating attrition as the spatial transformation process. Agriculture exhibited a process of patch creation, as the increase in CA resulted from an increase in PN. Finally, for savannas, and built-up and bare soil, there was noted a decrease in PN alongside an increase in CA, suggesting an aggregation-type spatial transformation process.

4. Discussion

4.1. Methodology

The use of Landsat images, despite its spatial resolution of 30 meters, remains pertinent due to its long-term availability, global coverage, and ability to capture multispectral data [54]. This characteristic allows for the analysis of regional-scale land cover changes, which is crucial for monitoring urban expansion evolution and associated landscape changes over time [54,55]. Among the available land cover classification algorithms, the Random Forest was selected for its capability to handle large volumes of multispectral data and produce detailed land cover maps. By employing a set of decision trees, this algorithm categorizes image pixels based on their spectral, spatial, and contextual characteristics, facilitating effective differentiation between various land cover types [56]. The combination of Landsat imagery and a Random Forest classifier proved effective in assessing urban expansion and associated landscape changes in the New Juaben Municipality, Ghana) [57]. Moreover, incorporating landscape metrics offers quantitative assessments of changes in urban landscape patterns and associated ecological impacts. Parameters like fragmentation, connectivity, urban patch shape, and landscape diversity serve as crucial tools for comprehending the dynamics of urban expansion, growth patterns, and their environmental repercussions) [58]. However, it's important to note that Landsat data primarily provides information on the physical environment and may not capture directly the social, economic, and political factors influencing landscape dynamics, which could be a limitation. Nevertheless, aligning Landsat image dates with significant events such as economic developments or political changes allowed to link observed spatial changes with these events, enhancing the understanding of how socio-political and economic factors impact urban spatial expansion, as demonstrated by Cabala et al. [31] in the KCA study. However, expanding the assessment of land cover changes through predictive modelling for future scenarios would be crucial for anticipating potential environmental impacts, including biodiversity loss, particularly in the context of climate change and unforeseen urban growth. Nevertheless, analyzing the period from 1990 to 2023 enables the identification of long-term landscape changes, such as urban area expansion, and shifts of ecosystems, offering insights into potential future developments. Similarly, studying land cover pattern dynamics is important as it enables deductions about ecological processes. For instance, while Li et al. [59] deduced a threat for biodiversity conservation in Shangai (China) through an increase in the cover of green space planted with exotic species.

4.2. Urban Expansion Intensity and Associated Landscape Dynamics

The period of 1990s-2001 was characterized by a slow urban expansion, along the RD Congo-Zambia border, due to political conflicts, leading to economic and social instability. Moreover, political tensions can disrupt commercial activities and hinder economic growth, thereby limiting employment opportunities and urban development [60]. Additionally, during conflict periods, priorities often shift towards addressing immediate security and political stability issues, relegating urbanization projects to the background [61]. Conversely, we found an acceleration in the urbanization extent during 2000-2010, which correspond to the period of progressive mining liberalization, favouring the increased of foreign and domestic investments, and stimulating economic growth and creating jobs. This economic growth also resulted in an increased demand for labor, attracting a growing population to urban areas near mining zones [62,63]. Khoji et al. [23] and Cabala et al. [31] findings corroborate the acceleration of the urbanization of main agglomerations within the KCA during the same period. Furthermore, mining companies could finance the construction of housing for mine workers, promoting the expansion of border cities [64].

The decade 2010-2023 is characterized by a new phase of rapid urbanization of the Congolese border cities due to political stabilization and economic recovery after the global financial crisis. Indeed, the increasing demand for mineral resources attracted new investments and labor, thereby stimulating urban growth. Yet, the increasing demand for housing among new citizens is driving uncontrolled urban expansion, extending beyond any form of government control. Indeed, Congolese cities are generally surpassing their own limits and encroaching on adjacent rural areas. Consequently, given the expensive urban lifestyle prevalent in many urban centers, a significant portion of the population in border cities opts for areas that maintain their rural essence, where land resources remain relatively affordable [65]. Concurrently, during this decade, the roads connecting the city of Lubumbashito Sakania via Mokambo and Kasumbalesa, have all been asphalted. This process has had a significant impact on the spatial expansion of Mokambo, notably by promoting the purchase of plots with traditional houses and their transformation into modern houses for customs officers and their families. This pattern aligns with Arimah's [66] findings regarding infrastructure's role in enhancing the prosperity of African cities. However, this trend of urban modernization can induce the loss of female local knowledge, particularly regarding the painting of traditional houses.

Unfortunately, the spatial urban expansion observed in border cities between the DRC and Zambia is largely driven by self-construction, leading to urban sprawl as revealed by the increase in CA and PN simultaneously, as well as UIIE. This self-construction is due to inadequate urban planning, gaps in land management, corruption, and the influence of political and economic interests on urban decisions [67]. This process often leads to the development of informal settlements with excessive low built-up density and limited access to basic services, as illustrated by Groupe Huit [68] in the city of Lubumbashi. Additionally, the self-construction favour land speculation, amplify socio-economic disparities, as well as food insecurity due to persistent agricultural activities displacement [69], despite the general trend of increase in the CA of agriculture between 1990 and 2023. Indeed, most farmers lack ownership titles, making their land vulnerable. On the other hand, land speculation makes the conversion of agricultural land into buildable land economically profitable in the short term and with lower risk [70]. Consequently, the spatial expansion of border cities leads to the exurbanization of some producers, transforming urban agriculture into peri-urban, and then rural agriculture. This shift is exacerbated by limited land availability and the prioritization of infrastructure over agricultural activities in (peri-)urban areas [71].

The expansion of built-up areas is accompanied by a decline in forest cover and an increase in savannahs around the studied border cities. Urbanization often leads to direct deforestation for building, road infrastructures, and other urban installations. Trees are cut down to make space, reducing forest cover, as demonstrated by Cabala et al. [35] in the Lubumbashi plain and Bamba et al. [72] in the Kisangani region in the DR Congo. Additionally, urban growth can favour the exploitation of surrounding natural resources for fuelwood or other forest products. Yet, excessive exploitation leads to a significant decrease in forest cover [41,73]. Moreover, urban expansion can stimulate the conversion of forested lands into agricultural lands to compensate for the loss of agricultural lands encroached upon by urbanization to meet the growing food needs of the urban population, as observed in Freetown [70]. However, the felling of trees promotes sunlight penetration, a crucial element for the development of herbaceous vegetation dominating savannahs [74], justifying the progression of savannahs in the surrounding areas of the studied border cities. Furthermore, the fragile fertility of agricultural lands in these regions often leads to their rapid abandonment after 2 to 3 years of cultivation, favouring their colonisation by herbaceous vegetation.

Our research findings affirm the impact of city size on deforestation, a notion supported by previous research [55,72]. The rapid urban expansion observed in Kasumbalesa can be attributed to the benefit from significant commercial flows due to its strategic position at national borders. This dynamic commercial activity creates employment opportunities, attracts investments, and stimulates economic growth, leading to rapid urban expansion. In contrast, territorial capital border cities or mining towns (i.e. Kipushi and Sakania) exhibited different characteristics since they may face specific challenges. These cities are more focused on a specific industry, which can limit their ability to diversify their economy and attract a more varied population. Similarly, territorial capital border cities are more focused on administrative and governmental functions, which can impact their economic dynamism and attractiveness to investors and migrants. Moreover, when examining the environmental impact of Kisangani, a city of considerable economic and demographic influence, in contrast to Ubundu, a city of moderate importance, Bamba et al. [72] discovered a more significant deforestation trend around Kisangani.

However, our results underline urban densification in Kipushi that can be attributed to the stable supply of electricity and water in the area. When infrastructures such as electricity and water are consistently available, it creates an environment conducive to the extension of former buildings. But, the recent trend of populations moving from Lubumbashi to settle in Kipushi could justify the urban sprawl observed in this city. This population movement can be attributed to the increasing insecurity in Lubumbashi. Additionally, job opportunities in mining industry may attract workers and their families. In contrast, the situation in Sakania is different due to the nature of its mining workforce. While the town hosts a significant number of mine workers, it is important to note that many of these workers come from other regions of the country and do not have their families and residence in Sakania.

4.2. Implications for Regional Urban Planning

Our findings highlight a rapid urban expansion and a trend towards urban sprawl in the studied border cities. Yet, urban sprawl contributes to community fragmentation, increased traffic congestion, and diminished quality of life through reduced green spaces and environmental degradation. It results also in escalated costs for infrastructure and public services, heightened transport inefficiencies, and reduced productivity associated with extended commuting times. To mitigate these effects, sustainable urban planning strategies are imperative. These encompass the implementation of urban densification policies aimed at optimizing the utilization of existing spatial resources, promotion of mixed-use urban development to minimize travel distances, and targeted investments in efficient public transportation infrastructure. Additionally, preserving green spaces and adopting urban growth management policies are pivotal for enhancing urban resilience and fostering balanced and sustainable urban development.

Our findings reveal that urban expansion leads to the disappearance of forests, posing threats to biodiversity conservation. Countless adverse effects of deforestation can be observed in the region, as the gradual disappearance of numerous non-timber forest resources [75], the noticeable reduction in wildlife [32], and the reduction in albedo, resulting in fewer rainy days [76]. Also, preserving agricultural land is crucial for ensuring food security in the area given its susceptibility to urban expansion. It is crucial to adopt environmental protection policies such as creating protected areas, strengthening land use regulations, and promoting sustainable agricultural practices. For illustration, the establishment of the Bururi Forest Reserve (near the city of Bururi, Burundi), with active ecoguards, has facilitated the regeneration of forest resources [77]. Concurrently, in the plateau des Batéke in Kinshasa (DR Congo), agroforestry practices that involve combining Acacia trees with maize crop enable households to enhance their yield and obtain wood for charcoal production in close proximity [78]. Certainly, reforestation solutions exist near the city of Kipushi [75]. However, despite ongoing progress, it has become evident that this positive human influence on the landscape through the planting of exotic species will not be sufficient to offset the loss of forest [79] while posing the risk of further ecosystem degradation [80]. It is crucial to utilize native species in reforestation projects and to enhance collaboration between traditional authorities and city managers to make more land available on the outskirts of cities for reforestation projects [81]. Additionally, there is a need to enhance the capacity of technical public service agents in terms of methodology and technical expertise to better supervise these types of reforestation projects and preserve reforested land in the long term [82]. Finally, the promotion of urban forestry in cities where green spaces are neglected and surrounding forests are disappearing is also an essential strategy for improving quality of life and the environment. By reducing air pollution, providing shade, lowering temperatures, and encouraging community participation through the planting of fruit trees, this approach creates healthier, more resilient, and more livable urban environments [83,84]. This integrated method is crucial for addressing the environmental and social challenges of contemporary urban areas, thereby promoting sustainability and collective well-being.

5. Conclusions

The study delves into the intricate spatial urban growth and associated landscape dynamics of Kipushi, Kasumbalesa, Mokambo and Sakania along the border with Zambia, an area of variable cross-border trade. Employing a robust methodology that integrates remote sensing, GIS, and landscape ecology analysis tools has facilitated a comprehensive mapping and quantification of landscape dynamics over time. The findings affirm that rapid population growth is primary driver of substantial landscape changes within these cities, leading to the transformation of natural landscapes into sprawl urban zones. However, there is the variance in the extent of urban expansion among cities, with those experiencing lower customs activity levels demonstrating a more restrained spatial urban growth. This discrepancy underscores the influential role of economic factors in shaping border urban dynamics, highlighting the intricate interplay of factors driving urban development in these regions. Additionally, the study has identified a marked acceleration in urban spatial growth since the early 2000s, indicating a swift transformation of these border agglomerations. This leads to an ecosystem shift, marked by a decline in forest cover, previously dominant in the landscapes of the cities studied in 1990, except for Kipushi, replaced by emerging and merging savannah patches in the landscape in 2023. The anthropization of the landscape of these cities and their heterogeneity are further intensified over time by the continuous creation of agricultural patches. These observations lead to the conclusion that border urban areas in the DRC-Zambia region are undergoing significant spatial changes driven by a complex interplay of economic, demographic, and infrastructural factors. Understanding these dynamics is crucial for informed urban planning and policy formulation. Sustainable development strategies must be devised to address challenges such as managing urban sprawl, promoting economic resilience, and fostering cross-border cooperation in border regions.