Edinburgh Explorer Urban Land-Use Dynamics in the Niger Delta: The Case of Greater Port Harcourt Watershed

: Cities in developing countries are urbanising at a rapid rate, resulting in substantial pressures on environmental systems. Among the main factors that lead to ﬂooding, controlling land-use change offers the greatest scope for the management of risk. However, traditional analysis of a “from–to” change matrix is not adequate to provide information of all the land-use changes that occur in a watershed. In this study, an in-depth analysis of land-use change enabled us to quantify the bulk of the changes accumulating from swap changes in a tropical watershed. This study assessed the historical and future land-use/land-cover (LULC) dynamics in the River State region of the Niger Delta. Land-use classiﬁcation and change detection analysis was conducted using multi-source (Landsat TM, ETM, polygon map, and hard copy) data of the study area for 1986, 1995, and 2003, and projected conditions in 2060. The key ﬁndings indicate that historical urbanisation was rapid; urban expansion could increase by 80% in 2060 due to planned urban development; and 95% of the conversions to urban land occurred chieﬂy at the expense of agricultural land. Urban land was dominated by net changes rather than swap changes, which in the future could amplify ﬂood risk and have other severe implications for the watershed. net changes (i.e., change in quantity). This study argues that accounting for net change only limits understanding of the sum of changes that occur within a landscape. The analysis of landscape persistence, gains, and losses is informative and vital to understanding the nature of changes that occur over time.


Introduction
Globally, cities are urbanising at unprecedented rates, resulting in profound and unintended impacts on landscapes and urban hydrology [1][2][3][4]. Urbanisation is defined as the concentration of people in urban settlements and the process of change in land-use occupancy, resulting from the conversion of rural lands into urban, suburban, and industrial communities [5,6]. In 2014, more than half (54%) of the world's population lived in cities [2,7] and by 2050, the urban population is projected to rise by 2.5 billion people, with 90% of this population anticipated to be concentrated in Asia and Africa, where 60% of this population is anticipated to live in slums [2]. There is now a wide consensus that cities in Asia and Africa are more vulnerable to the effects of flooding due to increased urbanisation, climate change, poor planning, weak regulations, and poor adaptive capacity [2,[8][9][10]. As a consequence, there is a critical need to understand the historical patterns of land-use change, the extent of future urbanisation and their resulting impacts on urban landscapes in these regions to inform flood risk management strategies.
Greater Port Harcourt (GPH) city in Nigeria is a compelling example of rapidly urbanising cities in Africa, with a delicate and dynamic wetland environment that is very sensitive to flooding [11]. Since 1965, Port Harcourt municipality has been the site of Nigeria's largest harbour and the centre of its petroleum activities [12,13]. From its initial development, the city's planning authority has struggled to cope with the physical expansion and overcrowding issues due to a continual rise in population [14,15]. Currently, their vision (developed in 2009) to expand, modernise, and transform the city into an internationally recognised destination is being implemented with a 50-year integrated urban master plan [14,16]. However, research has shown that both planned and unplanned developments can have significant negative effects on urban landscapes and could have severe flood risk implications [2,[17][18][19]. Therefore, understanding the extent and nature of the land-use/land-cover (LULC) changes using multi-temporal satellite data is vital for interpreting landscape dynamics and will aid in strengthening planning policy in the region.
Empirically, LULC change results have extensive use in flood modelling [20][21][22][23][24]. For example, Du et al. [20] used LULC data inputs to generate scenarios for assessing the effects of urbanisation in the Qinhui River Basin in China. Elsewhere, McColl and Aggett [25] utilised LULC data for forecasting flood impacts in Kittitas County in Washington, United States (USA). In terms of hydraulic modelling, Duan et al. [26] used remotely sensed data to produce flood maps of Chiang Mai province in Northern Thailand. For hydrological purposes, land-use data are applied as spatial data for generating curve number (CN, an index for determining run-off potential), percentage of impervious surfaces (PctImp), and Manning's N (roughness coefficient), which are all dependent on the land cover properties.
Satellite remote sensing is a method that is frequently used to detect land-use changes on landscapes [27], involving the application of multi-sourced and multi-date satellite imageries to estimate the difference in the extent of land cover. Remotely sensed data can be effectively used to quantify changes in a consistent and accurate way over a range of temporal and geographical scales [27][28][29][30]. Moreover, satellite remote sensing provides the most common source of data for mapping the change patterns that can, in most cases, be readily and repeatedly obtained in a digital and georeferenced format [31]. However, a limitation with this method is that different features can be difficult to identify if they appear the same to the sensor and may lead to misclassification; for example, artificial and natural grass in green light may display the same characteristics.
In the past, researchers have applied a range of change detection methods that can broadly be grouped into pre-classification and post-classification [17,27,[32][33][34]. Pre-classification involves the analysis of transformed images from two different dates. The methods used to perform this analysis include principal component analysis (PCA), band combination (BC), vegetation index differencing (VID), image regression, image rationing, and background subtraction techniques [31,34,35]. In contrast, the post-classification method employed in this study classifies multi-temporal images into thematic maps, and then implements a pixel-by-pixel comparison of the classified images.
The post-classification detection method offers a rational and more efficient means of estimating changes, because in addition to its traditional use for determining the extent (percentage) of change (for utilisation as an input for flood modelling), it can also be used to examine the nature of change i.e., to analyse the 'from-to' changes. However, the final accuracy of change determination depends on the quality of the classified image for each date [33]. As a result, several existing studies have provided limited insight into the processes that influence the pattern of land-use change and the component of change [36]. Beyond the extent of change and land-use conversions, this study was able to examine the other nature of change components, as well as the process of change, which could have a bearing on flooding in the watershed.
Presently, a number of studies have investigated the extent of land use change in the GPH watershed [20,22,37]. However, land-use detection studies that account for the processes that determine change through an analysis of the nature of this change are rare. In addition to accounting for the process of change, these studies neglected land persistence as a significant factor. According to Pontius et al. [36], total change in landscapes could result from two processes: swap changes (change in location) and net changes (i.e., change in quantity). This study argues that accounting for net change only limits understanding of the sum of changes that occur within a landscape. The analysis of landscape persistence, gains, and losses is informative and vital to understanding the nature of changes that occur over time. The concept of swap change is very important, as it gives information about the processes that determine the change. Accounting for the total change of landscapes is often overlooked in change detection studies [36]. This study addresses a gap in the existing body of knowledge relating to the study area, as it provides a more comprehensive interpretation of the sum of changes. Importantly, it explains the nature of change in the watershed that could have implications for flood risk.
In this study, supervised classification was applied in the image classification process and later used to derive the cross-tabulation matrix. The cross-tabulation matrices for 1986, 1996, and 2003 were extended to derive the gross gains, gross losses, and persistence by principal categories. Swap and net changes were derived from the extended cross-tabulation matrix.
Since land-use change resulting from the implementation of the GPH 50-year development plan can have serious implications for the watershed, the purpose of this paper is to understand the historical and future LULC dynamics in the Greater Port Harcourt watershed. The key questions are: What was the extent and nature of historical LULC changes? What will be the extent of future urban LULC changes due to implementation of the plan by 2060? What are the dominant forces and the trend of urban land-use change in the watershed? Lastly, does the application of multi-source data affect result accuracy?

Study Area
The study area under investigation is a lowland watershed made up of five major basins and 39 sub-basins spanning about 4820 km 2 and is in River State, which is the southernmost part of the Niger Delta in Nigeria ( Figure 1). This study area is a hydrologically sensitive coastal plain that is barely 20 m above sea level. It is relatively a flat terrain with a slope not greater than 3% [38], and is located between 4 • 42 north and 4 • 47 north latitude as well as 6 • 55 east and 7 • 08 east longitude in the equatorial region of the world.
Administratively, Greater Port Harcourt (GPH) city is an agglomeration of eight local communities surrounding the Port Harcourt city. GPH was comprised of Port Harcourt city, Obio-Akpo, and Okrika areas prior to the establishment of the Greater Port Harcourt Law in 1999. Currently, based on the newly implemented urban master plan, GPH city now also includes the Ikwerre, Oyigbo, Ogu/Bolo, Etche, and Eleme administrative areas [16]. Port Harcourt is one of two cities with the highest population in the Niger Delta; it has witnessed an exponential population growth with an annual growth rate of 3.0%, resulting in population eruption from about 180,000 people in 1963 to about two million in 2016 [39]. By 2020, the population of Rivers State is projected to rise from about three million to seven million people [40].
Development controls in Nigeria date back to the colonial era with the enactment of the Township Ordinance No. 29 of 1917 as the first attempt of introducing spatial orderliness into the land-use pattern in Nigerian cities, and later the Nigeria Town and Country Planning Ordinance (No. 4 of 1946) to provide for the development and improvement of different parts of the country through planning schemes initiated by the planning authorities [41][42][43]. However, those regulations only reflected the social and economic needs of the colonialists at the time.
Since independence in 1960, planning authorities have increasingly used development controls with the aim of achieving the efficient and effective use of land, as well as physical improvement across the country. The Nigerian Urban and Regional Planning Law, Decree No. 88 of 1992, and a subsequent amendment, the Urban and Regional Planning (Amendment) Decree No. 18 of 1992, were later enacted [41,43]. These latter land-use regulations generally set out limits, principles, and procedures for dealing with land-use issues such as overcrowding and congestion. Despite these statutory controls, issues of overcrowding and the construction of squatter settlements persist, and levels of compliance remain low [42]. The federal, state, and local governments share responsibility for land-use planning and development control, with the Federal Planning Commission taking responsibility for national physical development plans and sub-regional plans, the Rivers State Urban and Regional Planning Board controlling urban and regional plans in the state, and the local planning authorities controlling the town and rural plans [42,44].

Materials and Methods
This study adopted the post-classification change detection approach for analysing historical LULC changes. Four major steps were followed in analysing historical LULC changes, which were data acquisition, image enhancement, supervised classification, and change detection analysis. These steps can be grouped under the pre-processing, classification, and post-classification stages ( Figure 2). The analysis of the extent of future changes was estimated by overlaying the GPH master plan (digitised vector map) on the baseline map derived from the satellite data.

Data Acquisition
For this study, two multi-temporal and multi-spectral imagery data sets and a vector map were obtained to generate digital maps for analysing historical changes between 1986-2003. The three historical data sets obtained include one Landsat Thematic Mapper (TM) imagery for the year 1986; one polygon map for the year 1995; and one Landsat Enhanced Thematic Mapper plus (ETM+) map for the year 2003. These dates were selected to benchmark approximate decadal changes of the landscape. One limitation of this approach is the use of multi-sourced data that could increase uncertainty due to the varying spatial and spectral resolutions of the data set. This problem is common and often unavoidable in studies spanning over decades. Despite the limitation, a similar approach has been successfully followed in detecting land-use changes in prior studies [17,20,31].
The TM and ETM+ are Landsat satellite data (Table 1) were obtained free of charge from the United States Geological Survey Global Land Facility website (http://glcf.umd.edu/data/landsat/). TM sensor data were used because it is the only freely available data covering the spatial and temporal extent of the study area. The 1995 polygon map was obtained from the Rivers State's Ministry of Land and Housing (Riv-MoLH), with this data generated using the aerial survey method by the Riv-MoLH, and then converted to a polygon map. The polygon map was also applied because it was the only available cloud-free data for the 1995 time period. High cloud cover and satellite scanner issues hindered the use of more recent maps (2007 onwards) of the area. Originally, 2009 was selected as the baseline year as it benchmarks the commencement of the GPH development, but due to the unavailability of appropriate data for that year, 2003 was adopted as the baseline year. Although this could increase the uncertainty of the result, the effect may not be significant, given that the temporal difference is less than 10 years, and because the change in urban area is the main interest of this study. Ultimately, the criterion used for data selection was that maps must have 20% or less cloud cover.

Data Acquisition
For this study, two multi-temporal and multi-spectral imagery data sets and a vector map were obtained to generate digital maps for analysing historical changes between 1986-2003. The three historical data sets obtained include one Landsat Thematic Mapper (TM) imagery for the year 1986; one polygon map for the year 1995; and one Landsat Enhanced Thematic Mapper plus (ETM+) map for the year 2003. These dates were selected to benchmark approximate decadal changes of the landscape. One limitation of this approach is the use of multi-sourced data that could increase uncertainty due to the varying spatial and spectral resolutions of the data set. This problem is common and often unavoidable in studies spanning over decades. Despite the limitation, a similar approach has been successfully followed in detecting land-use changes in prior studies [17,20,31].
The TM and ETM+ are Landsat satellite data (Table 1) were obtained free of charge from the United States Geological Survey Global Land Facility website (http://glcf.umd.edu/data/landsat/). TM sensor data were used because it is the only freely available data covering the spatial and temporal extent of the study area. The 1995 polygon map was obtained from the Rivers State's Ministry of Land and Housing (Riv-MoLH), with this data generated using the aerial survey method by the Riv-MoLH, and then converted to a polygon map. The polygon map was also applied because it was the only available cloud-free data for the 1995 time period. High cloud cover and satellite scanner issues hindered the use of more recent maps (2007 onwards) of the area. Originally, 2009 was selected as the baseline year as it benchmarks the commencement of the GPH development, but due to the unavailability of appropriate data for that year, 2003 was adopted as the baseline year. Although this could increase the uncertainty of the result, the effect may not be significant, given that the temporal difference is less than 10 years, and because the change in urban area is the main interest of this study. Ultimately, the criterion used for data selection was that maps must have 20% or less cloud cover.

Image Processing
Image processing was carried out prior to change detection analysis. Images were manipulated in a geographical information system (GIS) environment in order to relate data to the true biophysical environment. A number of sub-processes were undertaken, including mosaicking, geometrical correction, map registration, atmospheric correction, and image enhancement to manipulate the selected scenes. Due to the position and scale of the maps, four map tiles were obtained and merged for each year. The maps were then spatially referenced to WGS_84_UTM_Zone_32N in the Universal Transverse Mercator (UTM) system.
Next, image enhancement was performed by adjusting the image values to highlight the distinct features for improving the visual interpretability of the maps. While the range of image enhancement techniques is broad, the composite colour generation technique was applied for enhancing images, because it makes the fullest use of the capabilities of the human eye. Natural colour band combination (3, 2, 1), characterised by short wavelengths, was applied for image classification of the years 1986 and 2003.

Supervised Classification
The image classification process involved the conversion of multi-band raster imagery into a single-band raster with categorical classes that relate to different LULC classes [45]. In this study, the supervised classification (a deterministic) method was applied. It was applied to the 1986 (TM) and 2003 (ETM+) satellite images, and involved clustering pixels into classes based on training data defined by the user [46]. Supervised classification was performed for the 1986 (TM) and 2003 (ETM+) satellite data. This method was preferred over the unsupervised classification method as it provides the advantage of understanding not just the extent of change or change magnitude, but also because it could be used to understand land-use conversions and change processes that have taken place through a pixel-by-pixel comparison using the matrix [47].
The Maximum Likelihood Classifier (MLC) method was utilised for image classification, which was achieved by generating signatures and training samples based on field knowledge and the existing map. The signature file store spectral signatures of LULC classes across the series of bands [45] and pure pixels were used as training samples. The procedure for classifying maps involved: clipping the classified and composite maps; reclassifying the existing classified (1995) map; identifying LULC classes; collecting training samples; generating signature files; executing the maximum likelihood classification tool; creating thematic maps; and reclassifying the map from six to five classes for the 1986 and 2003 maps, and from 10 to five classes for the 1995 map. Built-up areas were reclassified to the urban category; forest, thicket, and swamp forest were reclassified to the forest category; and plantation and cultivation were reclassified to the agricultural land category. High and low mangrove were reclassified to the mangrove category, while river was reclassified to the water body category.

Post-Processing
After classification, overall accuracy and kappa statistics were performed to validate the classification using confusion matrices. In total, 468 random points were used to compare classes in the truth and classified maps. The main steps in accuracy assessment involved: (1) random sampling, Urban Sci. 2018, 2, 108 7 of 24 extraction of reference (truth) and classified map data using the 'extract value-to-point' tool in ArcMap; and (2) comparison of the reference map to the classified map using an error matrix. Data were organised using frequency and pivot tools in ArcMap, and the overall accuracy was calculated as the total predicted pixels divided by the sum of the predicted pixels (of the major diagonal). According to Enaruvbe and Ige-Olumide [48], an overall accuracy of <50% is poor, 50-70% is moderate, while greater that 70% is considered good. Subsequently, the Kappa statistic was used to measure the magnitude of agreement between the reference and classified map [49]. Kappa is expressed as:

LULC Change Detection and Analysis
After classification, land-use change detection was performed to estimate the differences between two scenes. This was achieved in two ways, first by determining the extent of change (Tables 2-4), and then by determining the nature of the change (Tables 5-11). The extent or magnitude of change is used to describe the degree of expansion or reduction in the LULC size resulting from the classification [47], with negative values indicating a reduction in LULC size, and positive values indicating an increase in the size of the LULC class.
It is expressed as: where A = percentage of change (i.e., percentage increased or decreased); F = first date (date of the first imagery analysed); and I = reference date (date of the second imagery analysed) The nature of change was calculated from a transition matrix table, which consisted of the total and inter-category conversions of LULCs from 1986 to 2003. This was used to describe the type of "from-to" changes, transitions, and conversions that have occurred between the two dates [47,50]. The rationale for this analysis was to interpret the prominent transitions in a particular category from one year to the next.
Cross-tabulation was performed in the GIS environment. In the table, rows display the categories in Time 1, whereas the columns represent the categories in Time 2. Row totals at the right signify the proportion of the landscape by category in Time 1, and the column totals at the bottom represent the proportion of landscape by category in Time 2. As a first step, the off-diagonal entries were analysed, followed by the column and row total, and subsequently followed by the persistence, gross gains, and gross losses. This matrix was used to derive the total and inter-category gross gains, gross losses, and persistence, along with the absolute net change for the entire landscape and classes. Furthermore, as Pontius et al. [36] emphasised, other measures such as swap changes in the landscape as well as loss-to-persistence and gain-to-persistence ratios can be determined. These measures help define the processes that determine the transitions with the advantage of allowing determination of both the landscape and inter-categorical transitions.
Gross loss was derived by deducting persistence (diagonal entries) from the row total, and gross gain was calculated by deducting persistence from the column total. Total change was derived by summing the gross loss and gross gain, while absolute net change was calculated by subtracting the gross loss from the gross gain. By deducting the absolute net change from the total change, the swap change was derived. The loss-to-persistence ratio and the gain-to-persistence ratio can be computed using the approach in Pontius et al. [36].

Procedure for Estimating Changes to Future Urban Area
With the availability of the GPH Master Plan, the estimation of future changes in the urban area was performed using ArcMap 10.1, assuming that urban expansion would result mainly from the implementation of the master plan. Again, the study assumed that the conditions within the map area Urban Sci. 2018, 2, 108 8 of 24 would change due to master plan implementation, while the conditions of other LULC classes outside the GPH LULC map would largely stay the same.
First, the master plan was digitized and overlaid on the 2003 baseline map. The digitisation process was done, first by georeferencing the analogue maps (master plan) to an appropriate projected coordinate system, i.e., WGS_84_UTM_zone_32N, and later by digitising the LULC classes. The main urban classes that were reclassified to urban area from the master plan include built-up area, university premises, commercial light industrial, industrial, cemetery, a dumpsite, an Airforce base, and the international airport. Open spaces, active open spaces, and golf courses were also reclassified to agricultural land. Next, LULC classes from the GPH Master Plan were reclassified into urban, forest, agriculture, mangrove, and water based on percentage sealing and tree coverage. Riverine areas were reclassified to mangrove, while water body was reclassified to waterbody.  Tables 5-7). Kappa statistics calculated for the 1986 and 2003 classified maps were 0.54 and 0.50, respectively, which clearly indicate that the overall accuracy of both maps was reasonable, because they were above 70%, and the kappa statistics for both maps represent a moderate rating. area would change due to master plan implementation, while the conditions of other LULC classes outside the GPH LULC map would largely stay the same. First, the master plan was digitized and overlaid on the 2003 baseline map. The digitisation process was done, first by georeferencing the analogue maps (master plan) to an appropriate projected coordinate system, i.e., WGS_84_UTM_zone_32N, and later by digitising the LULC classes. The main urban classes that were reclassified to urban area from the master plan include built-up area, university premises, commercial light industrial, industrial, cemetery, a dumpsite, an Airforce base, and the international airport. Open spaces, active open spaces, and golf courses were also reclassified to agricultural land. Next, LULC classes from the GPH Master Plan were reclassified into urban, forest, agriculture, mangrove, and water based on percentage sealing and tree coverage. Riverine areas were reclassified to mangrove, while water body was reclassified to waterbody.  Tables 5-7). Kappa statistics calculated for the 1986 and 2003 classified maps were 0.54 and 0.50, respectively, which clearly indicate that the overall accuracy of both maps was reasonable, because they were above 70%, and the kappa statistics for both maps represent a moderate rating.

Extent of Historical Urban LULC Change
Row 1 in Table 2 shows the proportion and area of historical urban land use/land cover in the watershed between 1986-2003. The results indicate that significant urban expansion occurred between 1986-2003. The trend line in Figure 6 and historical maps in figures 3-5 further demonstrate that the historical change was progressive, while the bar chart indicates that the LULC changes between approximate successive decades were significant. For example, the urban LULC category increased in proportion from about 2.7% in 1986 to about 5.8% in 1995, and finally to around 11% in 2003 (see Tables 2-4). Table 2 and Figure 6 indicate that the landscape experienced a drastic change of about 309% between 1986-2003, because urban area expanded from about 134.5 km 2 in 1986 to about 550 km 2 in 2003. This means that the urban area expanded at a rate of about 24 km 2 /annum within the entire historical period studied. Although the percentage change in the later period (1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003) was slightly lower than that in the earlier period (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995), the results demonstrate that the expansion in the first and second time period were similar, showing that significant urban growth occurred in the past between 1986-2003, and that the growth was progressive.

Extent of Historical Urban LULC Change
Row 1 in Table 2 shows the proportion and area of historical urban land use/land cover in the watershed between 1986-2003. The results indicate that significant urban expansion occurred between 1986-2003. The trend line in Figure 6 and historical maps in Figures 3-5 further demonstrate that the historical change was progressive, while the bar chart indicates that the LULC changes between approximate successive decades were significant. For example, the urban LULC category increased in proportion from about 2.7% in 1986 to about 5.8% in 1995, and finally to around 11% in 2003 (see Tables 2-4). Table 2 and Figure 6 indicate that the landscape experienced a drastic change of about 309% between 1986-2003, because urban area expanded from about 134.5 km 2 in 1986 to about 550 km 2 in 2003. This means that the urban area expanded at a rate of about 24 km 2 /annum within the entire historical period studied. Although the percentage change in the later period (1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003) was slightly lower than that in the earlier period (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995), the results demonstrate that the expansion in the first and second time period were similar, showing that significant urban growth occurred in the past between 1986-2003, and that the growth was progressive.  Proportion of Urban Area (%)  Table 5 presents the results of the projected percentage change in urban area by 2060 resulting from implementation of the master plan. Based on the map in Figure 7 and chart in Figure 8, this study Urban Sci. 2018, 2, 108 12 of 24 projects that the urban LULC category would expand significantly by about 80% by 2060. Meanwhile, the stacked bar chart in Figure 9 indicates a progressive growth in the proportion of urban area, which is expected to rise from about 11% in 2003 to about 20% by 2060, meaning a double increase of the urban area. The bar chart displays a significant rise in the urban land category from about 550 km 2 in 2003 to about 988 km 2 in 2060, which represents respectively an 80% and 600% growth from the 2003 and 1986 calculated urban land cover. the proportion of urban area for the respective dates, while the line graph shows the trend and change in urban land-use area. Table 5 presents the results of the projected percentage change in urban area by 2060 resulting from implementation of the master plan. Based on the map in Figure 7 and chart in Figure 8, this study projects that the urban LULC category would expand significantly by about 80% by 2060. Meanwhile, the stacked bar chart in Figure 9 indicates a progressive growth in the proportion of urban area, which is expected to rise from about 11% in 2003 to about 20% by 2060, meaning a double increase of the urban area. The bar chart displays a significant rise in the urban land category from about 550 km 2 in 2003 to about 988 km 2 in 2060, which represents respectively an 80% and 600% growth from the 2003 and 1986 calculated urban land cover.

Urban Transitions
Based on the land-use conversion results (Tables 5-7), it was found that urban area increased significantly mainly at the expense of agricultural land. Results show that the conversion of all the other LULC categories to urban area caused a total of 549 km 2 in urban area between 1986-2003. The most prominent transition involves the conversion of agricultural land (about 422 km 2 ) to urban land. That means that about 93.3% of all conversions to urban land resulted from development on agricultural land. In contrast, the conversion of other LULC classes to urban areas were very low, such that the transition from mangrove, forest, and water to urban area were only about 15.3 km 2 , 13.3 km 2 , and 1.9 km 2 respectively.  Table 8 reveals that the prominent non-urban shifts between 1986-2003 were both from forest to agricultural land (about 321 km 2 ) and from agricultural land to forest (about 229.5 km 2 ), resulting in a net change from forest to agricultural land of 92.5 km 2 with the conversion of forest to agricultural land accelerating between 1995-2003. In addition, the conversion of 135 km 2 of forest to mangrove was another notable shift between 1995-2003. Generally, non-urban transitions were more significant in the latter period.

Analysis of Landscape Persistence and Components of Change
As shown in Table 8, all of the LULCs (i.e., the landscape) made a total gross gain of 1178.4 km 2 , which is about 25% of the watershed. A gross gain of one category is accompanied by a gross loss of another category, which means the total gross gain is equal to the total gross loss in a landscape. About 75% of the landscape persisted to change. In other words, about only 25% of the watershed transitioned from one LULC category in the entire historical period of 1986-2003. Importantly, the landscape experienced a total change of 48.8% and a total absolute net of 17%, while 31% of the entire watershed experienced a swap change. The urban area (with the highest intensity) had the greatest gain of about 452.8 km 2 between 1986-2003, followed by a gain in agricultural land of about 365.2 km 2 ( Table 8). At the same time, agricultural land experienced the highest gross loss of about 654.4 km 2 by a large margin, followed by mangrove (about 257.8 km 2 ). In summary, the entire landscape experienced more persistence than transition and more swap change than net change.

Inter-Category Net and Swap Change and Gain-to-Persistence Ratio
The results in Table 8 show that although agricultural land experienced the bulk of the total change, which was characterised by swap change rather than net change, in contrast, the bulk of the urban land changes was a net change. Agricultural land and forest lands experienced more swap changes between 1986-2003 (Table 8). This swap change was significant between 1986-1995 (Table 9), however it accelerated between 1995-2003 (Table 10), while urban land and water experienced only minor swap changes between 1986-2003 (Table 8). Moreover, Table 11 illustrates that the gain-to-persistence ratio is greater than one (>1) for the urban and forest categories, suggesting that urban and forest land had a higher tendency to expand than persist between 1986-2003. Urban land experienced the most prominent gain-to-persistence ratio (of about 4.7), demonstrating that the urban area was more prone to expanding than persisting when compared to other categories.

What Was the Extent of Historical LULC Changes
Based on the change detection analysis, this study reveals that the urban land cover experienced significant growth in extent of over 300% between 1986-2003, increasing from 135 km 2 to around 550 km 2 , signifying that urban area quadrupled in extent. In a recent study of the local area by Mmom and Fred-Nwagwu [37], the urban area increased by 86% between 1986-2007, while in another [48], it was reported that only 10% change occurred. One explanation for this disparity is differences in the scope (area) of the landscape studied in each study with the area covered by the previous studies situated at a municipal (smaller) scale ranging between 400-570 km 2 , while the spatial extent covered by this study was at a regional scale, covering the entire watershed of about 4800 km 2 . In summary, this study found that the urban area quadrupled in extent, but an important lesson is that the scale of the urban landscape that is studied can have an impact on the observed extent of change.
Compared to urban expansion documented elsewhere, changes in the urban extent measured as percentage change was found to be very significant (>100%) in this study, and other similar studies have also reported a dramatic change in urban land cover within the same time period. For example, Al Kuwari and Kaiser [31] discovered that urban expansion in the desert region of the western Nile Delta was very significant (4935%) i.e., from 23 km 2 in 1984 to 1158 km 2 in 2005. Another study by Shalaby and Tateishi [51] established that urban area expanded significantly by 666% in the northwestern coastal zone of Egypt. Likewise, the urban expansion of 309% that was experienced in the GPH watershed is considered significant. In contrast, Yuan et al. [52] found that only a 38% (<100%) change occurred in the Minnesota metropolitan area between 1985-2002. Hegazy and Kaloop [19] concluded that a less significant change (67%) took place between 1984-2000 in Mansoura and Talkha cities in Egypt. In short, previous studies have demonstrated great variability, with a range of drivers influencing the extent of land-use change. This study places the estimated land-use change to urban area between the upper and lower estimates of previous studies.
Urban land is usually a relatively small proportion of total land cover in many catchments [47,[53][54][55], which is consistent with the findings in this study. According to Allan [54], urban land exceeded only 5% of the catchment area in 29 river basins and 10% in only 10 of the 150 large basins studied. However, in this study area, urban area was greater than 10% by 2003, and is likely to reach 20% in 2060. However, even at low percentages, permeable surfaces in urban areas can have a disproportionate impact on flooding [54,56], which means that a significant change in the extent of urban land in the GPH watershed is expected to cause dramatic changes in peak discharge. The findings in this study are also consistent with other studies indicating that agricultural land occupies the largest proportion of land area in various developed catchments, whereas urban land occupies a much smaller portion [54].
Studies of LULC changes in other watersheds have shown variable results. Built-up areas increased by 100% (13.64 km 2 year −1 ) in the Upper Citarum Watershed, West Java Province in Indonesia between 1997-2005 [57], while the urban area in the Simly Watershed, Islamabad (Pakistan) increased by 80% (0.42 km 2 year −1 ) between 1992-2012 [58]. However, a study in the Weihe River Basin of northwest China [59] detected that the urban area increased by only 16% (15.5 km 2 /year −1 . A study analysing the multi-annual changes in urban area in Shijiazhuang City found that urban area increased by 81.5% between 1987-2001 [60], where urban growth was referred to as 'fast expansion' because the urban area grew at a rate of about 5.6 km 2 /year. In this study, it was found that the urban area increased at a rate of about 24 km 2 /year between 1986-2003, suggesting that the GPH watershed experienced a higher rate of urbanisation than in these other cities, and can be described based on Xiao et al. [60] study as "fast expansion". Identifying the driving forces of land-use change is valuable for enhancing the understanding of the nature, extent, and spatial pattern of urban expansion, and for delivering effective future land-use plans. Six (6) main categories of driving forces are frequently considered in studies: physical, socioeconomic, proximity, accessibility, neighbourhood, and policy factors [61][62][63]. Prior research on the study area has shown that socioeconomic factors (mainly population eruption) are the main driving force of land-use change [13,37,[64][65][66]. The population growth that was experienced in the area-from about 700,000 in 1991 to about one million inhabitants in 2006 [37]-coincides with the rapid land consumption for urban expansion that was reported in this study. A local social survey in 1989 revealed that migrants accounted for 72% of the city's population, where 65% of the migrant population originated from rural areas [13]. This suggests that rural-urban migration contributed considerably to the expansion of the city, which could be due to the petroleum-related economic activities in the area.
Like Nigeria, population growth due to rural-urban migration is also the main cause of urban expansion in India and China [4]. While the drivers of urban expansion in Africa and Europe are both socioeconomic in nature, the rate of urban expansion is much slower in Europe, and is more related to gross domestic product (GDP) growth than population growth [4]. According to the United Nations (UN), the rural population in this area will continue to decline until 2050 [67] due to rural-urban migration. Since the rise in urban population is predicted to continue in the future, urban sprawl is likely to continue, and the accompanying increase in impermeable surfaces could have serious flood risk and planning implications.
Other socioeconomic drivers of land-use change in the area include the high cost of available land and the desire for a traditional water-based living [68]. In the past, migrants from riverine communities have resorted to squatting in informal shelters on waterfronts in the south of the city, creating areas of high population density, whereas migrants from upland areas tend to build informal settlements in the north of the city [68], which helps to explain the spatial distribution of urban land in the fringes of the city. Physical factors also affected the distribution of land-use change in the study area, since features such as rivers as well as excessively ponded areas and marshlands condition where people live in the city [14]. Nonetheless, low-income groups still live in such areas of high flood risk.

What Is the Extent of Future Urban LULC Changes Due to the Implementation of the Plan by 2060?
Regarding future changes, this study revealed that a significant change in the extent of the urban area by 80% is expected by 2060, meaning that the new city's extent is projected to be twice the size of the city in 2003. This estimated change in size is comparable to trends in recent studies. For instance, one study projected that urban land is expected to increase by 27.5% between 2009-2030 [59]. Note that the built-up area is generally considered a parameter for quantifying urban sprawl, and can be quantified by measuring the changes in impervious surface [3,23,69]. Du et al. [20] predicted a rapid rise in impervious surface, from 23% to 31% between 2012-2018 with a corresponding impact on flooding. Similarly, the significant change in urban cover expected by 2060 in this study area implies that a future increase in urban surface is likely to have severe implications for flooding in the watershed.
The implications of the findings are profound, since urbanisation and loss of agricultural land were the dominant forces of land-use change in the GPH watershed. However, as in most regions, urbanisation in the study area was the main driving force in the past, and is likely to continue in the future, which implies that flood risk could increase dramatically due to likely drastic changes in the future. Even relatively small changes in the proportion of urban area can cause dramatic flooding effects [20,56]. If the catchment increases its percentage of impervious surface from 0-60% of its area due to urban growth, flow could increase two to five times of its pre-urbanisation state [70].
In 2012, Port Harcourt and other cities in the Niger River system suffered the effects of a disastrous flood event attributed to intense precipitation and developments upstream (GFDRR, 2013). The impact of the floods was reportedly very severe in Port Harcourt, leaving 363 people dead, 5851 injured, 3,891,314 others affected, and 3,871,530 displaced. Given the projected growth in urban area, urbanisation could increase the risk of flood in future if the development plan is poorly implemented.
Besides the flood risk implications, the planned expansion if poorly implemented could continue to significantly alter the natural landscape and disturb ecosystems as well as incur environmental, economic, and social costs [71][72][73]. For example, the rapid expansion of the city could also increase air pollution, traffic injury, and ecosystem services change. Moreover, the rapid land consumption may have other serious consequences, such as the loss of livelihoods due to a reduction of arable land, loss of biodiversity and critical habitats, and farmland displacement.
In addition to other residual driving forces such as population, land-use policy and planning could become the main driving factors of land-use change in the future due to the implementation of the master plan and associated land zoning. Although these could have positive impacts such as the provision of basic infrastructure and services, the construction of formal settlements, and the reduction of slums in some areas, this study argues that the proposed plan could amplify flood risk if poorly implemented. Moreover, urban expansion in the future could have negative consequences on environmental sustainability and quality of life if the policy is only applied in preferred areas.

What Was the Nature of Historical LULC Changes?
Regarding the nature of change, the data in this study revealed that three quarters (about 75%) of the watershed persisted, while as much as 25% of the watershed transitioned from one land-use category to another between 1986-2003 (Table 8). Unlike other LULC categories, this study found that the transition to urban land cover was the most dynamic in terms of gross gain and net change, which are the most important measures of urban expansion [33]. Moreover, urban area exhibited the largest gross gain (about 9% of the watershed) and the smallest gross loss, leading to a high net change of around 8.6%. Indeed, the most prominent transition was historically the conversion of agricultural land to urban land (Table 5), with 93.3% of all conversions to urban land resulting from agricultural land.
Beyond knowledge of the extent of change and land-use transitions, this study has extended knowledge of the watershed by providing new insights into the process of change that occurred historically between 1986-2003, demonstrating that urban transition was dominated by net changes (Tables 8-10). This implies that the transition from agricultural land to urban land experienced an actual change in extent rather than just a change in location (swap change).
It is expected that changes in the percentage of impermeable surfaces due to urbanisation would compound flood risk more than changes in a location. On the other hand, the prominent non-urban transitions such as agricultural land and forestland transitions were predominantly of swap type changes rather than net change, indicating that a considerable amount of agricultural and forest lands experienced relocation. While the effect is likely to be less dramatic than change to urban land use, the Urban Sci. 2018, 2, 108 20 of 24 relocation of forests could still have some flood implications because changes in surface roughness can have a significant impact on small local sub-basins [74].
The phenomena observed in the GPH watershed were comparable to findings in other watersheds such as the Mara River Basin in East Africa. For example, Mwangi et al. [75] found that swap changes accounted for more than 50% of the overall LULC changes. In our study, swap changes accounted for over 60% of the total changes in the GPH watershed, suggesting that the GPH watershed is very dynamic in terms of land cover relocation, resulting from accelerated farming activities in the area and high afforestation. Therefore, this study supports the argument that reporting only net changes can underestimate the total land-use change in landscapes [36].
The study also revealed that the prominent non-urban shifts were from forest to agricultural land (approximately 321.0 km 2 ) and from agricultural land to forest (about 229.5 km 2 ) between 1986-2003 (Table 8), which was not expected. These findings relate to the whole watershed and in part corroborate previous research findings, which covered a smaller part of the landscape. Although the impacts of such conversions vary at the watershed scale, the vulnerability to flooding may increase locally due to changes in soil infiltration capacity and increased surface runoff. Deforestation can also have severe implications for flooding because of reduced interception losses; however, the precise effect in large watersheds varies for different watersheds [76]. From the analysis, 73 km 2 of the forest was eventually lost between 1986-2003. However, the large fraction of forest land-use change was a swap type of change (568 km 2 ) rather than a change in a quantity. This observation implies forest disturbance, which can accelerate the rate that precipitation becomes streamflow depending on the catchment and storm size.

What Are the Dominant Forces and Tendencies of Land-Use Change?
This study has furthered the understanding of the landscape dynamics (in terms of the landscapes tendency to change) by revealing that urban land exhibited the highest gain-to-persistence ratio of 4.6 (far greater than one), as shown in Table 11. This means that urban land exhibited the greatest tendency to expand, while in contrast, agricultural land exhibited a low tendency to expand. Meanwhile, the observed changes in forest area was dynamic, showing a tendency to both expand and contract over time. Therefore, based on the findings of the extent of change, trend of change, process of change, and tendency of change presented, this study emphasises that urbanisation has been the most dominant force of land-use change in the watershed, which is likely to continue up to 2060.

Does the Application of the Multisource Data Affect Result Accuracy?
No change detection method is perfect, and similar to other studies, the results are presented with the caveat that uncertainties could arise due to possible errors from digitisation, input data quality, and the assumptions made. In particular, the use of multi-sourced data could have increased uncertainty in the results of this study due to the overclassification of some categories, as suggested by the mangrove class in the 1995 map ( Figure 4). Nevertheless, uncertainties as a result of urban land-cover change seemed insignificant. Results were generally deemed dependable because of the reliability of the classification method applied, the similarity in the trend of changes in the two classified TM maps used in this study, the similarity in the trend of changes of the classified maps in this study when compared to other studies of the area, the author's knowledge of the landscape, and the accuracy of assessment. This study was used to demonstrate that multi-sourced data could increase uncertainty in change detection analysis.

Conclusions
The Greater Port Harcourt watershed has experienced substantial changes in land use/land cover over the last few decades. The dominant force of land-use change in the watershed has been urbanisation. Prominent transitions include shifts from agricultural land to urban land, as well as agricultural land to forest land. The conversion to urban land was the most dominant land-use change in the watershed (about 415 km 2 ), and about 93% of this conversion was chiefly at the expense of agricultural land. The transition to urban land mainly experienced a net change process, which implies a change in actual quantity. The loss of agricultural land was another dominant force of land-use change; however, agricultural land mainly exhibited a swap process of change, i.e., mainly relocation. Reporting the net changes alone fails to capture the swap component, since it underestimates the total change, and could therefore be misleading. Beyond land-use expansion and transitions, this study provides new insights into the nature or process of change in the watershed that have not been previously reported. Generally, swap change explains up to 64% of the overall changes in this watershed. This means that more changes may have occurred in the watershed than those captured in previous studies.
The important lesson learnt from this study is that the scale of an urban landscape studied can have an impact on the examined extent of change compared to other similar studies. While the use of multi-source data could introduce error and increase uncertainty in estimating changes for data-sparse regions, the important messages are that: (1) urban expansion was rapid; (2) the expansion could rapidly increase by 2060 due to the GPH development; (3) the rapid urban expansion occurred mainly at the expense of agricultural land; and (4) changes to urban area were predominantly net rather than swap-type changes.