Data Analysis of Land Use Change and Urban and Rural Impacts in Lagos State, Nigeria

This study examines land use change and impacts on urban and rural activity in Lagos State, Nigeria. To achieve this, multi-temporal land use and land cover (LULC) datasets derived from the GlobeLand30 product of years 2000 and 2010 for urban and rural areas of Lagos State were imported into ArcMap 10.6 and converted to raster files (raster thematic maps) for spatial analysis in the FRAGSTATS situated in the Patch Analyst. Thus, different landscape metrics were computed to generate statistical results. The results have shown that fragmentation of cultivated lands increased in the rural areas but decreased in the urban areas. Also, the findings display that land-use change resulted in incremental fragmentation of forest in the urban areas, and reduction in the rural areas. The fragmentation measure of diversity increased in the urban areas, while it decreased in the rural areas during the period of study. These results suggest that cultivated land fragmentation is a complex process connected with socio-economic trends at regional and local levels. In addition, this study has shown that landscape metrics can be used to understand the spatial pattern of LULC change in an urban-rural context. Finally, the outcomes of this study will help the policymakers at the three levels of governments in Nigeria to make crucial informed decisions about sustainable land use.

Landscapes comprise of patches that are natural or altered by few human influences. These patches differ in their sizes, shapes, and arrays [22]. Rural landscapes are the effect of a dynamic process motivated by environmental and anthropogenic factors [23,24]. Landscape metrics are tools extensively used for analyzing [25], monitoring, and planning of landscape patterns [26]. These metrics are often used to analyze the spatial changes in the composition and configuration of landscapes [27]. Also, analyses such as fragmentation, diversity, changes in landscape patterns, and descriptions of the landscape structure are offered by landscape metrics. Besides, other techniques such as Remote Sensing (RS), Geographic Information System (GIS), land-use models, and statistical approaches have been used to examine changes in landscape patterns [28][29][30][31].
Urban growth in developing countries is diverse and alarming. As well, the residential and industrial incursions of growing urban areas lead to a constant decrease of high-quality agricultural land and forests in the marginal areas of cities and in rural areas [32,33]. Research shows an overall trend of declines in farmland during urbanization in many regions around the world [34][35][36][37][38][39]. From a regional context, semi-natural land reveals a decreasing trend, principally because of agricultural progression [40][41][42]. In Nigeria, studies that concern the mapping of urban growth, land use change and incursion into the rural lands include the following: A study was conducted by Dekolo, Oduwaye, and Nwokoro [43] using a multi-temporal RS and GIS analysis to find out urban sprawl and loss of agricultural land over a 40-year period in Ikorodu municipality, Lagos State. Findings display that agricultural lands were reduced by a significant amount of 25%, and this correlates with rapid urban changes. Olayiwola and Igbavboa [44] applied remotely sensed data to examine the trend of urban growth and its consequence on the land use pattern in Benin City, Nigeria. Research findings showed that the expansion of Benin City resulted in an inverse spillover effect on the vegetation resources of the adjoining rural settlements.
Landscape fragmentation is a process that brings about habitat loss and size reduction [45]; its main causes include agricultural growth, infrastructural projects, e.g., road construction, and urban sprawl [46]. Unrestrained urban expansion resulted in fragmentation of landscapes, destruction of wildlife habitat, and a decrease in biodiversity [47]. Several landscape metrics have been reported for quantifying the spatial patterns of landscapes at regional or local scales over the past years [48][49][50][51][52][53][54][55]. For instance, Crews-Meyer [56,57] used landscape metrics to evaluate the stability of agricultural land through the adoption of farmland patches as the elementary analytical element. Also, Southworth et al. [58] applied landscape metrics to pixel-based change trajectories to evaluate the degree of fragmentation of forest land cover over three periods. As well, Vaz et al. [12] employed landscape metrics to analyze the impacts of agricultural land changes and connection of tourism growth on these changes in Algarve, Portugal. Research findings showed an increase in urbanization while agriculture decreased significantly in Algarve between 1990 and 2006, which subsequently disturb the equilibrium of the environment [12]. The study by Luck and Wu [59] on gradient analysis of the urban landscape pattern in Phoenix, USA showed that landscape metrics are suitable to identify changes in land use from urban to rural environments. Thus, the benefit of landscape metrics for analyzing agricultural and socio-economic changes is extremely closely connected to the possibility of realizing the changing aspects of land use [60].
Much work has been accomplished using landscape metrics and remote sensing; however, little research has been concentrated on land use change and impacts on urban and rural activity in Lagos State using these tools. Therefore, this study seeks to understand if urban sprawl is altering rural landscape of the Lagos State, Nigeria. Likewise, the primary objective is to investigate land use change and impacts on urban and rural areas. The study specifically intends to (1) characterize landscape pattern changes across Lagos State from 2000 to 2010; and (2) quantify the impacts of land use change on the landscape structure of urban and rural areas using landscape metrics.

Study Area
The study area is Lagos State, which is situated in the low-lying coastal region of South West Nigeria. It is geographically situated on longitudes 2 • 42 -4 • 20 E and latitudes 6 • 22 -6 • 42 N. The State is bordered on the west by the Republic of Benin and on the north and east by Ogun State of Nigeria, as well as the Atlantic coast as its southern boundary. Lagos State is the most populous state in Nigeria with 9,113,605 residents in 2006 (see Table 1), which represented about 6.44% of the total population of Nigeria [61]. Through land area, Lagos State is the smallest in size in the country with 3577 km 2 or 357,700 hectares of land and 779 km 2 waterbodies.  [62]. Also, the four rural LGAs comprise (7) Ikorodu; (8) Epe; (11) Badagry; and (20) Ibeju-Lekki (see Figure 1); they are recognized as peri-urban settlements because they are located outside the metropolitan region of Lagos State. Furthermore, Lagos State experiences the typical tropical rainforest climate of southern Nigeria with very high temperatures and humidity throughout the year. There are two main seasons: a wet season which occurs from April to November, and a dry season which happens from December to March. The wet and dry seasons are also accompanied by Harmattan dry winds between December and early February of every year.

Study Area
The study area is Lagos State, which is situated in the low-lying coastal region of South West Nigeria. It is geographically situated on longitudes 2°42'-4°20' East and latitudes 6°22'-6°42' North. The State is bordered on the west by the Republic of Benin and on the north and east by Ogun State of Nigeria, as well as the Atlantic coast as its southern boundary. Lagos State is the most populous state in Nigeria with 9,113,605 residents in 2006 (see Table 1), which represented about 6.44% of the total population of Nigeria [61]. Through land area, Lagos State is the smallest in size in the country with 3577 km 2 or 357,700 hectares of land and 779 km 2 waterbodies.
There are 20 Local Government Areas (LGAs) in Lagos State. For the purpose of this study, these LGAs were divided into two sections: the urban (metropolitan); and the rural. The urban segment contains 16 LGAs  (19) Eti-Osa which covers approximately 90% of the state's total population [62]. Also, the four rural LGAs comprise (7) Ikorodu; (8) Epe; (11) Badagry; and (20) Ibeju-Lekki (see Figure 1); they are recognized as peri-urban settlements because they are located outside the metropolitan region of Lagos State. Furthermore, Lagos State experiences the typical tropical rainforest climate of southern Nigeria with very high temperatures and humidity throughout the year. There are two main seasons: a wet season which occurs from April to November, and a dry season which happens from December to March. The wet and dry seasons are also accompanied by Harmattan dry winds between December and early February of every year. In the past few years, the Lagos region has witnessed a fast-growing population and urban development which has altered its landscape. The Lagos region metamorphosed from a less important agricultural community into its current megacity status [63]. In the past few years, the Lagos region has witnessed a fast-growing population and urban development which has altered its landscape. The Lagos region metamorphosed from a less important agricultural community into its current megacity status [63].
However, the effect of population change as a driving force on land use has attracted much attention lately. Thus, deductions could be drawn about the impacts of population pressure after different types of comparative empirical work. But this empirical evidence is still lacking in developing countries [64]. In Nigeria for instance, it is known that the changing population-land affiliation and the related land use change have noticeable implications for the socio-economic conditions of the people and the environment [65].

Data Description
For this study, the description of the remotely sensed data used is illustrated in Table 2. The GlobeLand30, which is a land cover dataset produced by the National Geomatics Center of China (NGCC), was utilized. This dataset has a 30 m high resolution imagery that covers the period from 2000 to 2010 [66]; they are produced based on over 10,000 data images obtained mainly from Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper plus (ETM+). GlobeLand30 was a product of integration of a pixel-object knowledge-based method [67,68], and the Landsat images has minimal cloud cover and noiseless [69], which makes it appropriate for this study. A global accuracy of over 80% (i.e., >80%) were reported for the dataset [68,70,71]. As well, GlobeLand30 contain 10 land cover types [68], 8 (cultivated land, forest, grassland, shrubland, wetland, water bodies, artificial surfaces, and bareland) of which are applicable to Nigeria as itemized in Table 3.
Accuracy assessment of the dataset was conducted using high-resolution Google Earth images for study on land use land cover (LULC) change and statistical modeling of cultivated land change drivers in Nigeria [72]. The overall accuracy of >73% in 2000 and >75% in 2010 was achieved; thus, the suitability of this dataset for the analysis of this study was found to be high. Table 3. GlobeLand30 land-cover types used for the study (Chen et al. [68]).

Land Cover Code
Land Cover Type Description

Cultivated land
Lands used for agriculture, horticulture and gardens, including paddy fields, irrigated and dry farmland, vegetation and fruit gardens, etc.

Forest
Lands covered with trees, with vegetation cover over 30%, including deciduous and coniferous forests, and sparse woodland with 10%-30% cover, etc.

Grassland
Lands covered by natural grass with a cover over 10% including typical grassland, meadow grassland, alpine grassland, desert grassland.

Shrubland
Lands covered with shrubs with a cover over 30% including deciduous and evergreen shrubs and desert steppe with a cover over 10%, etc.

Wetland
Lands covered with wetland plants and water bodies including inland marsh, lake marsh, sea marsh, river floodplain wetland, forest/shrub wetland, peat bogs, mangrove and salt marsh, etc.

60
Water bodies Water bodies in the land area, including rivers, lakes, reservoirs, fish ponds, etc.

Methods
This section explains the methods used for the land use change and impacts in urban and rural areas of Lagos State. Figure 2 illustrates the flowchart of the study, which consists of: (a) land use data-GlobeLand30 product; (b) creation of suitability maps of land use change using Esri ArcGIS 10.6; and (c) integration of a FRAGSTATS (Patch Analyst raster version) to compute the landscape-level and class-level metrics for each land cover type. This implies that the methods of this study were based on the evaluation of LULC changes and landscape metrics analysis.

Methods
This section explains the methods used for the land use change and impacts in urban and rural areas of Lagos State. Figure 2 illustrates the flowchart of the study, which consists of: (a) land use data -GlobeLand30 product; (b) creation of suitability maps of land use change using Esri ArcGIS 10.6; and (c) integration of a FRAGSTATS (Patch Analyst raster version) to compute the landscapelevel and class-level metrics for each land cover type. This implies that the methods of this study were based on the evaluation of LULC changes and landscape metrics analysis.

LULC Maps Used
In this study, the Esri ArcGIS 10.6 software was used to perform the analysis of the GlobeLand30 product of years 2000 and 2010 in Lagos State. The Universal Transverse Mercator (UTM) projection system within zone 31 North with a World Geodetic System (WGS) 1984 at 30 m spatial resolution was used as the reference system. Figure 3 shows the classified LULC maps of Lagos State used during the period of study.
Data 2020, 5, x FOR PEER REVIEW 6 of 19

LULC Maps Used
In this study, the Esri ArcGIS 10.6 software was used to perform the analysis of the GlobeLand30 product of years 2000 and 2010 in Lagos State. The Universal Transverse Mercator (UTM) projection system within zone 31 North with a World Geodetic System (WGS) 1984 at 30 m spatial resolution was used as the reference system. Figure 3 shows the classified LULC maps of Lagos State used during the period of study.

Landscape Metrics Analysis
Landscape metrics (LMs) are quantitative tools and environmental indicators used to define the structures and spatial pattern of a landscape [48,73,74]. They have been recognized to be appropriate for landscape analysis, although they are restricted because of up and downscaling of the generated results [18,31,46,58,75]. Since the 1990s, the introduction and advancement of landscape metrics has resulted in speedy progress of mapping landscape ecology [74].
Landscape metrics can be calculated on three different levels: patch, class, and landscape levels [76]. In this study, landscape-level and class-level metrics were employed. Class-level metrics are used to qualify the characteristics of the same LULC type and return a unique value for each class in the landscape [77]. Quite a lot of different metrics [73,78] were developed and improved based on the work of O'Neill et al. [48]. Thus, software was developed to update computation of landscape metrics [52], which helped to improve their use. In this study, LMs were calculated using the Patch Analyst which is an extension to the ArcGIS (ArcMap) 10.6 software. The Patch Analyst extension calculates spatial statistics on both polygon files (e.g., shape files) and raster files (e.g., Arc grids). Thus, landscape metrics at the landscape-level and class-level for individual land cover types were computed with the raster version of the Patch Analyst. In this case, the input data included GlobeLand30 LULC Maps of 2000 and 2010 (the generated maps were divided into categories-urban and rural) of Lagos State (see Figure 2), which were imported into ArcMap and converted into grid (raster) files for analysis in the FRAGSTATS situated in the Patch Analyst.

Landscape Metrics Analysis
Landscape metrics (LMs) are quantitative tools and environmental indicators used to define the structures and spatial pattern of a landscape [48,73,74]. They have been recognized to be appropriate for landscape analysis, although they are restricted because of up and downscaling of the generated results [18,31,46,58,75]. Since the 1990s, the introduction and advancement of landscape metrics has resulted in speedy progress of mapping landscape ecology [74].
Landscape metrics can be calculated on three different levels: patch, class, and landscape levels [76]. In this study, landscape-level and class-level metrics were employed. Class-level metrics are used to qualify the characteristics of the same LULC type and return a unique value for each class in the landscape [77]. Quite a lot of different metrics [73,78] were developed and improved based on the work of O'Neill et al. [48]. Thus, software was developed to update computation of landscape metrics [52], which helped to improve their use. In this study, LMs were calculated using the Patch Analyst which is an extension to the ArcGIS (ArcMap) 10.6 software. The Patch Analyst extension calculates spatial statistics on both polygon files (e.g., shape files) and raster files (e.g., Arc grids). Thus, landscape metrics at the landscape-level and class-level for individual land cover types were computed with the raster version of the Patch Analyst. In this case, the input data included GlobeLand30 LULC Maps of 2000 and 2010 (the generated maps were divided into categories-urban and rural) of Lagos State (see Figure 2), which were imported into ArcMap and converted into grid (raster) files for analysis in the FRAGSTATS situated in the Patch Analyst.
Moreover, the following 12 landscape metrics in the Patch Analyst extension program were included in the analysis: AWMSI (area weighted mean shape index), MSI (mean shape index), MPFD (mean patch fractal dimension), AWMPFD (area weighted mean patch fractal dimension), TE (total edge), ED (edge density), MPS (mean patch size), NUMP (number of patches), PSCOV (patch size coefficient of variance), PSSD (patch size standard deviation), CA (class area), and SDI (Shannon's diversity index). The descriptions of the landscape metrics that were used in this study are obtainable in the Patch Analyst User's Manual [79]. Also, the ecological denotations of these landscape metrics were detailed in Herold et al. [27,80,81] (Visit this Website for information about the denotations of LMs: http://www.umass.edu/landeco/research/fragstats/documents/fragstats.help.4.2.pdf). Table 4 shows the 12-landscape metrics adapted from the above-mentioned scholars. The table displays the names, acronyms, and units of measure of the metrics used. The Shannon's diversity index, which is only available at the landscape level, is a relative measure of patch diversity. The SDI is equal to zero when there is only one patch in the landscape, and it increases as the number of patch types increases [73]. In this study, SDI was used to analyze the richness of the diversity of different land use types within the administrative boundaries of Lagos State LGAs. Finally, the selection and combination of numerous metrics for this study corroborate the suggestions mentioned by scholars indicating that a mixture of several indices allow a better understanding of the fragmentation process [74]. After the computation of the landscape fragmentation indices for all the classified LULC maps using FRAGSTATS, research results are shown in Tables 5-7.

Results and Discussions
The

Patch Analysis: Overall Landscape Metrics and Class-Specific Landscape Metrics for Urban Areas in Lagos State
The FRAGSTATS generated statistical results for urban areas of Lagos State are explained subsequently. The overall landscape fragmentation for Lagos State urban areas display slight increases for AWMSI, MSI, MPFD, AWMPFD, and MPS, while the PSSD showed a substantial increase from 2000 to 2010 (Table 5). Also, TE, ED, NUMP, and PSCOV decreased during the period of study. Total Edge exhibited the highest decrease, next is the NUMP while ED slightly decreased, and CA generated zero numbers during the time of study ( Table 5). The zero numbers of CA denote that the patch types became increasing scarce in the Lagos State landscape between 2000 and 2010. The land-cover class data for urban areas of Lagos State in 2000 and 2010 is presented in Table 6. During the time of study, the landscape of each land cover type showed different trends of change. It is important to mention that the FRAGSTATS statistical results for class-specific landscape metrics generated numbers for bareland in year 2000, but not any numbers were generated for this land cover type in 2010. As a result, bareland was not itemized in the class-specific landscape metric Table 6 for discussion. Firstly, the AWMSI metric for grassland and artificial surfaces displayed increments with grassland presented a higher increment amount of 1.88 than artificial surfaces of 0.16 from 2000 to Thus, the following analysis focuses on the spatiotemporal dynamics of land use change and impacts on urban and rural areas of study.

Patch Analysis: Overall Landscape Metrics and Class-Specific Landscape Metrics for Urban Areas in Lagos State
The FRAGSTATS generated statistical results for urban areas of Lagos State are explained subsequently. The overall landscape fragmentation for Lagos State urban areas display slight increases for AWMSI, MSI, MPFD, AWMPFD, and MPS, while the PSSD showed a substantial increase from 2000 to 2010 (Table 5). Also, TE, ED, NUMP, and PSCOV decreased during the period of study. Total Edge exhibited the highest decrease, next is the NUMP while ED slightly decreased, and CA generated zero numbers during the time of study ( Table 5). The zero numbers of CA denote that the patch types became increasing scarce in the Lagos State landscape between 2000 and 2010.
The land-cover class data for urban areas of Lagos State in 2000 and 2010 is presented in Table 6. During the time of study, the landscape of each land cover type showed different trends of change. It is important to mention that the FRAGSTATS statistical results for class-specific landscape metrics generated numbers for bareland in year 2000, but not any numbers were generated for this land cover type in 2010. As a result, bareland was not itemized in the class-specific landscape metric Table 6 for discussion. Firstly, the AWMSI metric for grassland and artificial surfaces displayed increments with grassland presented a higher increment amount of 1.88 than artificial surfaces of 0.16 from 2000 to 2010. Also, there were declines in cultivated land, forest, shrubland, wetland, and waterbodies during the period of study. The most noticeable decline occurred with forest with an amount of 1.23, and wetland showed the lowest decline amount of 0.20. In conclusion, the AWMSI of cultivated land, forest, shrubland, wetland, and waterbodies showed that shape irregularity in 2000 had decreased in 2010. Besides, the increase in artificial surfaces for AWMSI metric could suggest an urban sprawl that possibly was not well-planned in Lagos State.
The MSI of forest, shrubland and wetland increases in the urban areas; wetland showed the highest increase at the rate of 0.19, and shrubland displayed the lowest increase amount of 0.04 from 2000 to 2010. In addition, cultivated land, grassland, waterbodies, and artificial surfaces showed reductions in urban areas of Lagos State. The greatest reduction of 0.23 was displayed by waterbodies, while artificial surfaces showed the lowest reduction with the sum of 0.03.
The MPFD of forest and wetland increases in the urban areas of Lagos State during the period of study. It is worth mentioning that both forest and wetland exhibited the same increment of 0.02 from 2000 to 2010. Moreover, cultivated land, grassland, and waterbodies decreases during the study period, and the rate of decreases, i.e., 0.01 are the same for these three land cover types. Furthermore, shrubland and artificial surfaces presented no variations, thus, their amounts remain unchanged from 2000 to 2010. The MSI and MPFD metrics quantify the complexity of a patch shape [80]. As shown in Table 6, both metrics show a similar trend in change, and the values of both show growing trends for forest and wetland as well as decreasing trends for cultivated land, grassland, and waterbodies. The shape index signifies the regularity of a land cover type; a higher value represents higher irregularity while a decrease value implies a more regular shape of the patch. Results of this study in urban areas of Lagos State suggest that the shape of the forest and wetland became more irregular, but the shapes of the cultivated land, grassland, and waterbodies changed into more compact and regular shapes.
The AWMPFD metric for grassland increased marginally from 1.19 in 2000 to 1.20 in 2010 in Lagos State urban areas. Also, cultivated land, forest, shrubland, and waterbodies showed slight decreases; the greatest decrease of 0.04 was recorded for cultivated land and waterbodies exhibited the smallest reduction of the sum of 0.01 during the study period. As well, wetland and artificial surfaces presented no variations, thus, their amounts were unchanged from 2000 to 2010. The AWMPFD can be used to calculate the complexity of patch shape at a given observation scale. The lower value of AWMPFD lead to a more regular patch shape in the landscape and greater human disturbance on the patches. Finally, the decreases recorded for cultivated land, forest, shrubland, and waterbodies imply an increment of the intensity of human activities on these land cover types in urban areas of Lagos State; while the slight increase of AWMPFD for grassland meant the lessening of human disturbances on grassland in the study area.
Total edge (TE) equals the sum of the lengths (m) of all edge segments involving the corresponding patch type. It is an absolute measure of total edge length of a patch type (class level) or of all patch types (landscape level) [81]. The TE of forest and artificial surfaces increased in the urban areas for the period of 10 years studied. Forest exhibited a higher increase of 852,180 m in contrast to an increment from artificial surfaces. Also, the TE of cultivated land, grassland, shrubland, wetland and waterbodies declined; the TE of grassland exhibited the highest decrease amount of 1,031,040 m while waterbodies exhibited the least decrease amount of 17,940 m from 2000 to 2010, respectively.
Edge density equals the sum of the lengths (m) of all edge segments in the landscape, divided by the total land area (TLA) (m 2 ), multiplied by 10,000 (to convert to hectares) [81]. From 2000 to 2010, the ED of forest and artificial surfaces increased in Lagos State urban areas with forest showing the higher increase of 6.76 m/ha against 0.28 m/ha for artificial surfaces. For the period of study, the ED of cultivated land, grassland, shrubland, wetland, and waterbodies decreased in the urban areas. Grassland presented the highest decrease of 8.13 m/ha, while waterbodies showed the lowest decrease of 0.13 m/ha from 2000 to 2010.
Mean patch size (MPS) is the amount of central tendency in the patch characteristics of the entire landscape [81]. In this study, the MPS metric shows that forest, grassland, shrubland, and wetland in urban areas displayed increased amounts between 2000 and 2010. Wetland showed the highest increase of 15.92 ha and shrubland displayed the lowest increase of 0.21 ha. Moreover, the MPS metric for cultivated land, waterbodies, and artificial surfaces diminished; water bodies displayed the greatest decrease with the sum of 390.28 ha, and cultivated land showed the lowest decrease with the sum of 0.16 ha from 2000 to 2010. The fragmentation of cultivated land can be observed by the decrease of MPS. Thus, in this study, the reduction of MPS of cultivated land implies that cultivated lands had become more fragmented because of the loss of agricultural land.
Number of patches (NUMP) is an important metric which measures the degree of division or fragmentation of a specific patch type. The subdivision deals with the rate to which patch types are broken up into distinct patches [80]. The NUMP of forest, water bodies, and artificial surfaces increased in the urban areas; the highest metric increase of 18 occurred to water bodies while the lowest increase amount of 3 was displayed by artificial surfaces from 2000 to 2010. Between the period of 2000 and 2010, the NUMP of cultivated land, grassland, shrubland, and wetland displayed decreasing amounts in urban areas with the greatest decrease occurring to grassland at the rate of 1643, and the lowest decrease at the rate of 47 was shown by wetland. This decrease in patches indicate a lesser number of disturbances to the cultivated landscape in Lagos State during this ten-year duration. This reduction of the size of cultivated land patches substantiates the views of scholars that cultivated lands show a process of reduction in size to lesser exploitations [82].
Regarding discussion on the NUMP, MPS, and ED metrics for cultivated land in urban areas; decreases in values of NUMP, MPS and ED metrics for cultivated land indicate that cultivated land decreased in amounts, size, and edge in urban areas of Lagos State (see the urban segment of Table 6).
As itemized in Table 6, the PSCOV of grassland, waterbodies, and artificial surfaces showed increments in the Lagos State urban areas between 2000 and 2010. Grassland showed the highest increase of 747.21, while artificial surfaces displayed the lowest increase of 15.25. Results of spatial analysis also showed declines for cultivated land, forest, shrubland, and wetland. The highest decrease with the sum of 474.96 was recorded for forest, and wetland showed the lowest decrease with an amount of 26.25.
The PSSD of forest and grassland increased in the urban areas, and the higher increase was documented for forest at the quantity of 101.70 ha. In addition, the PSSD of cultivated land, shrubland, wetland, waterbodies, and artificial surfaces reduces during the study period in the urban areas of Lagos State. The greatest reduction with the sum of 1197.93 ha occurred with waterbodies, while cultivated land presented the lowest reduction with the quantity of 1.24 ha.
The CA of forest increased in the urban areas from 16,286.04 ha in 2000 to 31,116.42 ha in 2010. The CA of other land cover types-cultivated land, grassland, shrubland, wetland, waterbodies, and artificial surfaces showed declines in urban areas during the period of this study, but the highest decrease among these six land-cover types occurred with wetland at the amount of 10,458.81 ha while cultivated land showed the lowest decrease with the sum of 148.23 ha from the 2000 to 2010 study period.
Subsequently, cultivated land and artificial surfaces classes presented some differences. Results of spatial analysis of the landscape metrics (except MSI, MPS, PSSD, and CA) showed decreases in cultivated lands and increments in artificial surfaces in the urban LGAs of Lagos State between 2000 and 2010 (see Table 6). Thus, the growth of artificial surfaces was responsible for cultivated land decline. In 2010, the NUMP of cultivated lands had less intensive fragmentation than in 2000. The consistent decline in other fragmentation metrics, such as AWMSI, MSI, MPFD, AWMPFD, TE, ED, MPS, PSCOV, PSSD, and CA suggests that cultivated landscapes in 2010 were less fragmented. During the period of study, the NUMP of forest, water bodies, and artificial surfaces increased in the urban areas of Lagos State. Also, the NUMP of grassland, shrubland, and wetland decreased in urban areas (see Table 6). Relating the declines in cultivated lands to the increment in forest, it can be deduced that the lost cultivated land was being replaced by forest, which was due to the 2008 urban renewal and tree planting program carried out by the Lagos State government. This program was envisioned to help in alleviation of the harmful effects of climate change, thereby promoting a greener, healthier, and sustainable environment [83].

Patch Analysis: Overall Landscape Metrics and Land-Cover Class Level for Rural Areas in Lagos State
The spatiotemporal dynamics of landscapes for Lagos State rural areas were contrasted from the urban areas. The degree of landscape fragmentation in rural areas increased for the MSI, MPS, and PSSD metrics (Table 5). AWMSI, AWMPFD, TE, ED, NUMP, and PSCOV exhibited declines from 2000 to 2010. Substantial decreases were presented for TE and NUMP, while the AWMPFD metric showed a slight decrease. Besides, the overall land scape metric for CA in rural areas generated zero figures alongside urban areas during the time of study. Although AWMSI seemed to decrease, the extent of this decrease was higher than for MSI ( Table 5), implying that larger patches have more complex shapes; the increment of MSI portrayed the increased irregularity in patch shapes.
Results of the FRAGSTATS statistical analysis for rural areas of Lagos State (see Table 6) showed that the AWMSI metric for all the land-cover types (cultivated land, forest, grassland, shrubland, wetland, water bodies, and artificial surfaces) decreased between 2000 and 2010. Grassland exhibited the highest decrease with the quantity of 8.01, and water bodies recorded the lowest decrease with the quantity of 0.10.
The MSI metrics for forest, shrubland, and artificial surfaces presented increments in the rural areas of Lagos State during the period of this study. Artificial surfaces displayed the highest increment of the sum of 0.21 while forest showed the lowest increment of the sum of 0.02 from 2000 to 2010. The MSI of cultivated land, grassland, wetland, and water bodies showed decreases; cultivated land displayed the greatest decrease with an amount of 0.42, while grassland exhibited the lowest decrease with an amount of 0.05.
Results showed that the MPFD metric for artificial surfaces increased from 1.08 in 2000 to 1.10 in 2010. Also, the MPFD of cultivated land and water bodies exhibited slight declines; the greater decrease occurred from cultivated land with an amount of 0.04. It is worth mentioning that the quantities of MPFD metric for forest, grassland, shrubland, and wetland were unchanged from 2000 to 2010 (see Table 6). For example, the metric of forest, grassland, and shrubland generated the same amount, i.e., 1.04 in 2000 and remained the same at 1.04 in 2010. Also, the amount of wetland was 1.12 in 2000, and it was unchanged and remained at 1.12 in 2010 (see Table 6).
The AWMPFD of cultivated land, forest, grassland, shrubland, and artificial surfaces decreased in the rural areas of Lagos State from 2000 to 2010. Grassland showed the greatest decrease with the quantity of 0.12, and both shrubland and artificial surfaces displayed the lowest decrease with the quantity of 0.01 between 2000 and 2010, respectively. The AWMPFD of wetland and water bodies generated the same quantity of 1.18 in 2000 and remained at 1.18 in 2010 apiece without changing ( Table 6).
The TE of water bodies and artificial surfaces in rural areas of Lagos State showed increments for the period of 10 years studied. Artificial surfaces recorded the larger increment with the sum of 142,380 m. Also, the TE of cultivated land, forest, grassland, shrubland, and wetland decreases during the study period. The TE of grassland exhibited the highest decrease with the sum of 7,140,660 m and wetland presented the lowest decrease with the sum of 154,560 m from 2000 to 2010.
From 2000 to 2010, the ED of water bodies and artificial surfaces in rural areas of Lagos State showed increases. Artificial surfaces recorded the larger increase with the sum of 0.57 m/ha. Moreover, the TE of cultivated land, forest, grassland, shrubland, and wetland reduced during the study period. The TE of grassland exhibited the highest reduction with the sum of 28.40 m/ha and wetland presented the lowest reduction with an amount of 0.61 m/ha between 2000 and 2010.
The MPS of forest, shrubland, and artificial surfaces increased in the Lagos State rural areas between 2000 and 2010. Forest showed the greatest increase with the quantity of 48.92 ha while shrubland displayed the lowest increase with a quantity of 2.08 ha. Thus, the increment of MPS for artificial surfaces could imply there was rapid economic development, which resulted to an increasing demand for land use for construction and housing with consequent loss of agricultural land in rural areas of Lagos State. Also, cultivated land, grassland, wetland, and waterbodies showed decreases; the greatest decrease occurred with waterbodies with the sum of 538.26 ha, and grassland exhibited the lowest decrease with the sum of 0.93 ha in the rural areas during the study period.
Moreover, the NUMP of cultivated land, waterbodies, and artificial surfaces increased from 2000 to 2010. The greatest increase was documented for cultivated land with the sum of 72 while artificial surfaces demonstrated the lowest increase by a sum of 8 during the period of study. Also, the amounts of NUMP for forest, grassland, shrubland, and wetland showed reductions. Grassland presented the greatest reduction with a quantity of 8951, while wetland showed the lowest reduction of 6 between 2000 and 2010 in the rural areas.
The PSCOV of cultivated land, shrubland, wetland, and water bodies presented increments in the rural areas between 2000 and 2010 study period. The highest increment was recorded for water bodies with the amount of 120.37 and the lowest increment amount of 0.97 was documented for wetland. Besides, the PSCOV of forest, grassland, and artificial surfaces reduced; grassland showed the greatest reduction at the rate of 1648.99 while artificial surfaces displayed the lowest reduction at the rate of 21.41 in the Lagos State rural areas.
The PSSD of forest and shrubland rises in the rural areas during the period of study with the higher rise of 802.76 ha displayed by the forest. However, the PSSD of cultivated land, grassland, wetland, water bodies, and artificial surfaces decreased. The metric for water bodies displayed the greatest decrease at the rate of 1510.44 ha, while artificial surfaces exhibited the lowest decrease at the rate of 10.05 ha between 2000 and 2010.
Lastly, the CA of forest, shrubland, and artificial surfaces increased in the rural areas, and the greatest increase was displayed by forest at the rate of 25,696.89 ha while artificial surfaces showed the lowest increase at the rate of 989.37 ha between 2000 and 2010. But, the CA of cultivated land, grassland, wetland, and water bodies all showed reductions, and the greatest reduction was presented by grassland at the quantity of 22,600.62 ha while water bodies showed the lowest reduction at the quantity of 854.01 ha during the period of the study.
It is worth mentioning that the increasing fragmentation of cultivated lands in the rural areas of Lagos State may be attributed to rising anthropogenic incursion into cultivated landscapes that occurred as the outcome of urban expansion. The NUMP of artificial surfaces in the rural areas presented an increase from 96 in 2000 to 104 in 2010 (see rural segment of Table 6), and this concurs with the possibility that urban sprawl from Lagos Metropolitan Region or expansion of artificial surfaces in the rural parts of Lagos State might have resulted in an increase in cultivated lands fragmentation. These research findings correspond with Dekolo et al.'s [43] study in Ikorodu municipality, which found a significant reduction in agricultural lands due to rapid urban changes. Moreover, the NUMP of forest, grassland, shrubland, and wetland (except water bodies) decreased in the rural areas of Lagos State during the study period. Concerning the land-cover types (forest, grassland, shrubland, and wetland) and their relationships with the landscape metrics in rural areas, the grassland showed declines in all the metrics, while CA recorded the greatest and most substantial decline. Similarly, wetland diminished for almost all the metrics except for PSCOV. However, forest decreased for most metrics with increments in MSI, MPS, PSSD, and CA (see rural segment of Table 6).

Change of Landscape Diversity in Lagos State
Shannon's diversity index is extensively used to detect the change of landscape diversity. A higher value of SDI signifies more landscape elements or more diverse landscape patterns. Table 7 below illustrates the variation in SDI between 2000 and 2010 in Lagos State. During the period of this study, the richness of the extant patches at the local government level enables us to comprehend whether significant changes occur in the land use and land cover. Thus, the SDI calculation made it possible to evaluate the richness of the diversity of the total land cover types for each LGA.
As depicted in Table 7, findings from the Patch Analyst Software-FRAGSTATS concerning this landscape analysis showed a predominant decrease in diversity in ten LGAs, an increase in diversity occurring in eight LGAs, and there was no change in diversity in two LGAs throughout the years. Also, all the rural LGAs (Badagry, Epe, Ibeju-Lekki, and Ikorodu) mostly exhibited a decline in the diversity of patches-A negative SDI, while the more urbanized areas of Lagos (e.g., Agege, Ajeromi-Ifelodun, Ifako-Ijaiye, Ikeja, Lagos Island, Lagos Mainland, Shomolu, and Surulere) showed more change in the SDI.
Results substantiate Vaz et al.'s [12] study that applied SDI to quantify changes in agricultural areas in Algarve region, Portugal; the study found a lower SDI in rural areas while the urban areas reported a larger SDI. Thus, the diverse urban region of Lagos State could be associated with significant LULC changes. It is noteworthy that SDI variation for Eti-Osa and Mushin urban LGAs were equal to zero because no changes in landscape diversity were recorded in these two areas from 2000 to 2010. Finally, from the spatial concept, the findings demonstrate a signal of the loss of definite cultivated land cover types in the study area.
In summary, the major differences in spatiotemporal dynamics in urban and rural areas of Lagos State were: (a) overall, landscape fragmentation decreased in both areas, nevertheless, the decrease was larger in the rural areas. Cultivated land for all the metrics decreased in both areas but the proportion of the decrease was larger in the urban areas. Studies display a general trend of decreases in agricultural land during urbanization in several regions around the world [34][35][36][37][38][39]. Similarly, the rising agricultural land loss could have resulted from the constant propensity to use this land for urban activities, and may have been caused by the uncommonly rapid dynamics in the agricultural sector [12]; (b) fragmentation of artificial surfaces increased in both areas but the increment was higher in rural areas; (c) fragmentation of grassland, shrubland, and wetland reduces in both areas, but reductions of grassland and shrubland were larger in rural areas, and the reduction of wetland was larger in urban areas; (d) in comparison, fragmentation of cultivated land increased in the rural areas, but decreased in the urban areas; (e) fragmentation of forest showed an increase in the urban areas and a reduction in the rural areas; (f) in the urban areas, the fragmentation measure of diversity increased, in contrast to a decline in the rural areas.

Conclusions
Land use patterns are quickly changing within the urban and surrounding area of African cities [84]. Also, LULC is widely known as a local environmental issue worldwide [85]. Landscape metrics (LMs) have mainly been utilized in landscape ecology [12]; and they are effective instrument in LULC change detection and other environmental research [86]. Much work has been accomplished using landscape metrics and remote sensing; however, little research has been concentrated on land use change and impacts on urban and rural activity in Lagos State using these tools. This study was carried out to characterize landscape pattern changes and to quantify the impacts of land use changes in urban and rural areas of Lagos State using landscape metrics. The input data included GlobeLand30 dataset of 2000 and 2010 for Lagos State which were imported into ArcMap 10.6 and converted to grid (raster) files for statistical analysis in the FRAGSTATS algorithm situated in the Patch Analyst Software. Thus, landscape metrics at the landscape level and class-level for individual land cover types were computed with the raster version of the Patch Analyst to generate statistical results. The study of landscape metrics indicated loss of cultivated lands and changes in fragmentation between 2000 and 2010 in Lagos State. Firstly, it was found that fragmentation in cultivated lands in the rural areas increases during the time of study with the growth of urban region and diversification of land uses. Findings concur with the study which found that land use changes occur at the periphery of big city concentration where urbanization and industrialization pressures often result in loss of major agricultural lands and tree cover [87].
Next, land-use diversification leads to larger fragmentation of forest apart from near cultivated lands. Findings about increment in cultivated land fragmentation in rural areas, resulting from other land-use activities in urban region could lead to more deterioration of agricultural activities in Lagos State. Fragmentation that occurs through the conversion of agricultural lands to residential or manufacturing land uses may also have negative effects on the wildlife and on sustainable ecosystems [88,89]. The spatial variability of rural landscape attributes and pattern changes necessitate the development of land use management policies at the local, state, or federal levels [90]. In rural areas (such as Badagry, Epe, Ibeju-Lekki, and Ikorodu) around metropolitan region of Lagos State, environmental policy should focus on the enhancement of land use sustainability. However, the management of urbanized land requires decentralization of urban patches/urban land-use planning. Decentralization of urban patches specifies that urban areas such as Lagos metropolis with large physical size and high population density should reduce their growth cycles and enable peri-urban centers to grow. These growing peri-urban centers are characterized by lower levels of ecological footprints [91]. Because increasing the number of patches is an indication of ecosystem degradation [92], this strategy of decentralization of urban patches is a way of reducing the number of patches and addressing ecological problems on a regional scale.
This research has shown that LMs can be used to understand the spatial pattern of LULC change in an urban-rural context [11,93,94]. Moreover, the combination of GIS and LMs methods at the regional level with the spatiotemporal and statistical analysis can help to solve the problems of loss of agricultural lands and loss of the natural resource areas [12]. Also, the application of LMs to analyze landscape pattern changes reveal functional changes in the landscape, especially in a rapidly developing region [95]. Furthermore, LMs are universal, potent, and consistent; this suggests that landscape structures in different regions (e.g., Lagos State, Nigeria) have similar attributes including size, area similarity, and fragmentation, which can be studied by similar set of the metrics [77]. As well, these analytical tools may help decision makers to design more appropriate actions for land use management and sustainable environment. In conclusion, there is a need to adapt to the environmental policies and programs at the federal, state, and local levels to ensure the conservation and proper use of the Lagos State landscape.