<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JGIS</journal-id><journal-title-group><journal-title>Journal of Geographic Information System</journal-title></journal-title-group><issn pub-type="epub">2151-1950</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jgis.2019.115034</article-id><article-id pub-id-type="publisher-id">JGIS-95802</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Spatio-Temporal Dynamics and Evolution of Land Use Land Cover Using Remote Sensing and GIS in Sebou Estuary, Morocco
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ayoub</surname><given-names>El Adnani</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Adnane</surname><given-names>Habib</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Khalid</surname><given-names>El khalidi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Bendahhou</surname><given-names>Zourarah</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Laboratoire Géosciences Marines et Sciences des sols—LGMSS—URAC 45, Faculté des sciences El Jadida, El Jadida, Morocco</addr-line></aff><aff id="aff2"><addr-line>DBAR_ICoE-EOAWARE, (LGMSS-URAC45), Chouaib Doukkali University, BD Jabran Khalil Jabran, EL Jadida, Morocco</addr-line></aff><pub-date pub-type="epub"><day>27</day><month>09</month><year>2019</year></pub-date><volume>11</volume><issue>05</issue><fpage>551</fpage><lpage>566</lpage><history><date date-type="received"><day>14,</day>	<month>September</month>	<year>2019</year></date><date date-type="rev-recd"><day>15,</day>	<month>October</month>	<year>2019</year>	</date><date date-type="accepted"><day>18,</day>	<month>October</month>	<year>2019</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Land use and land cover (LULC) represent the ongoing challenge of environmental variation. The understanding of the level and process of its change is the basis for any environmental planning and management. In Morocco, as everywhere in the world, human population densities are constantly increasing on the coastal zones. This results in a continuous and rapid acceleration of the use of coastal space and an increase in pressures on ecosystems and the different species they contain. The purpose of this study is the analysis of the changes in LULC from 1985 to 2017 in the coastal area of Sebou estuary, situated in the Northwest of the Moroccan Atlantic coast. The changes were identified and assessed after classifying a series of Landsat images taken during 1985, 2002 and 2017. The algorithm used for the classification is the Support Vector Machine (SVM), which yielded results with accuracy higher than 85%. The results of the land use land cover change describe phenomenal urbanization and deforestation, as well as an evolution of the agricultural sector, indicating the impact of anthropization in this vulnerable environment.
 
</p></abstract><kwd-group><kwd>LULC</kwd><kwd> Change Detection</kwd><kwd> Accuracy Assessment</kwd><kwd> Landsat</kwd><kwd> Sebou Estuary</kwd><kwd> Anthropization</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Land use and land cover (LULC) is one of the fundamental topics in the study of environmental changes. During the 20<sup>th</sup> century, the world experienced a phenomenal and rapid population growth, which goes from 1.6 Billion in 1900 to 6.1 Billion in 2006 [<xref ref-type="bibr" rid="scirp.95802-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref3">3</xref>]. The vast majority of this growth has occurred in developing countries [<xref ref-type="bibr" rid="scirp.95802-ref3">3</xref>]. In Mediterranean countries, more than a third of this growth takes place in coastal regions, especially in coastal cities [<xref ref-type="bibr" rid="scirp.95802-ref4">4</xref>].</p><p>According to the World Development Indicators (WDI), the urbanization rate in Morocco has increased from 29% in 1960 to 60.3% in 2014, where 20.4 million of Moroccan’s population live in urban areas, while 13.4 million live in rural areas. This evolution leads to an over-exploitation of natural resources, rural exodus, uncontrolled growth of cities and also the creation of new urban centers [<xref ref-type="bibr" rid="scirp.95802-ref5">5</xref>]. Therefore, these changes boost the modification of land use causing a real danger to the ecosystems and biodiversity of the affected areas.</p><p>The evolution of urban planning and agriculture in an estuarine area is inevitable given the multiplicity of resources that this area can offer. Over the past three decades, the world has undergone very strong change and massive migration to coastal areas. The study of land use change is a key factor in environmental changes [<xref ref-type="bibr" rid="scirp.95802-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref10">10</xref>]. Those changes have a significant impact on natural resources through their impacts on soil and water quality, ecosystem processes and functions, and global climate systems [<xref ref-type="bibr" rid="scirp.95802-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref14">14</xref>]. Recently, remote sensing with multi-temporal satellite data has become a strong tool for monitoring aspects such as vegetation cover, soil degradation, urban expansion and more generally for most types of LULC changes [<xref ref-type="bibr" rid="scirp.95802-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref19">19</xref>].</p><p>The city of Kenitra, one of the biggest cities in Morocco, has experienced rapid population growth due to the creation and concentration of new and large industrial units making a large number of job opportunities and consequently leading to massive immigration of rural populations. All these parameters have led to environmental degradation in the absence of safeguard action plans integrated into the socio-economic development of the region [<xref ref-type="bibr" rid="scirp.95802-ref20">20</xref>]. In literature, no study has been conducted on the evolution and impact of LULC in the Sebou estuary. Therefore, this original work will try to better understand the evolution of the LULC in the study area and to serve as a basis for further environmental studies as well as it will give an overall assessment of anthropogenic development (agriculture, urban ...).</p><p>The objective of this research is to study the temporal changes of LULC for the years 1985, 2002, 2017 and to evaluate its transformation between 1985 and 2017. This analysis will make possible the identification of urban sprawl, as well as changes in forest and agricultural cover, and to quantify the impact of anthropogenic pressure on coastal environments. The produced results will make an essential contribution to decision-making in the fields of environmental management and future planning.</p></sec><sec id="s2"><title>2. Study Area</title><p>The Sebou Estuary situated in Northwest of the Moroccan Atlantic coast belonging to the El Gharb plain, located between 34˚08'21&quot;N and 34˚26'22&quot;N latitude, and 6˚44'06&quot;W and 6˚16'22&quot;W longitude (<xref ref-type="fig" rid="fig1">Figure 1</xref>). It is characterized by a</p><p>Mediterranean climate with mild, moderate and rainy Atlantic influence during winter, humid and temperate weather during summer with hot wind coming from the East. The rainfall average between 1973 and 2011 is about 537 mm/year, and the average temperature is about 15˚C. However, the temperature is characterized by an apparent variability (minimum temperature of 4˚C during winter and maximum 40˚C during summer) [<xref ref-type="bibr" rid="scirp.95802-ref21">21</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref22">22</xref>].</p><p>The Kenitra city being the capital of the El Gharb region, with a concentration of 90% of the region’s industrial factories, has become in a few decades a territory with strong anthropogenic attraction, causing a subsequent expansion of the urban area. The Kenitra and Mehdia agglomerations have been established on an exceptional natural site, represented by two wide loops of the Estuary Sebou’s meanders. The left bank is occupied by unsuitable constructions for this vulnerable environment, namely two port areas, a military air base, a penitentiary complex and an industrial area, while the right bank is occupied specially by agricultural lands.</p></sec><sec id="s3"><title>3. Material and Methods</title><p>In this research, three Landsat images were selected to map and assess the LULC changes over the last thirty years. The scenes are available with a spatial resolution of 30 m and acquired during the summer period for years 1985, 2002 and 2017. After the processing of Landsat images, LULC maps and LULC changes maps were produced and assessed. The flowchart in <xref ref-type="fig" rid="fig2">Figure 2</xref> summarizes the methodology used for this study.</p><sec id="s3_1"><title>3.1. Preprocessing</title><p>Images acquired by Landsat sensors are subject to the perturbation due to the effects of sensors, sun, atmosphere and topography. The step of image preprocessing attempts to minimize these effects to the extent required for a particular application [<xref ref-type="bibr" rid="scirp.95802-ref23">23</xref>]. First, Digital Number (DN) values were converted to radiance values. Later, the atmospheric correction was performed using FLAASH model (Fast Line of sight Atmospheric Analysis of Spectral Hypercubes) which incorporates a radiative transfer code based on MODerate resolution atmospheric TRANs mission (MODTRAN4) [<xref ref-type="bibr" rid="scirp.95802-ref24">24</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref25">25</xref>].</p></sec><sec id="s3_2"><title>3.2. LULC Mapping</title><p>LULC maps of the Sebou Estuary were produced for the years 1985, 2002 and 2017, using a supervised classification algorithm SVM (Support Vector Machine), following a meticulous selection of pixel samples, based on the spectral variation of each class. The generalized images were reclassified into 7 categories (<xref ref-type="table" rid="table1">Table 1</xref>).</p><p>Generally, a supervised classification requires learning samples as well as the definition of the size and number of learning samples to achieve a specific result,</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The identified LULC types in the study area</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Class</th><th align="center" valign="middle" >Description</th></tr></thead><tr><td align="center" valign="middle" >Built-up</td><td align="center" valign="middle" >Residential, industrial, commercial structures and highways</td></tr><tr><td align="center" valign="middle" >Agricultural land</td><td align="center" valign="middle" >Areas of land prepared for agricultural activity, including areas currently under cultivation and land under preparation, and greenhouse</td></tr><tr><td align="center" valign="middle" >Shrubs and forest land</td><td align="center" valign="middle" >Areas dominated by forest and covered by herbaceous plants such as grass and shrubs</td></tr><tr><td align="center" valign="middle" >Bare land</td><td align="center" valign="middle" >Unvegetated land</td></tr><tr><td align="center" valign="middle" >Water body</td><td align="center" valign="middle" >Sea, river, lake and marshes, ponds, and reservoirs</td></tr><tr><td align="center" valign="middle" >Beaches dunes sand</td><td align="center" valign="middle" >Beach sands and coastal dune sands</td></tr><tr><td align="center" valign="middle" >Spears vegetation</td><td align="center" valign="middle" >Areas consisting of arid lands with short and disperse vegetation</td></tr></tbody></table></table-wrap><p>which is one of the most critical problems of the supervised classification [<xref ref-type="bibr" rid="scirp.95802-ref26">26</xref>]. Previous studies show that SVM classifier (Support Vector Machine) is not related to sample’s size, and it has been improved to give better results with quality and limited quantity learning samples [<xref ref-type="bibr" rid="scirp.95802-ref27">27</xref>]. Although SVM is less known than other classifiers for LULC mapping, it has been shown to achieve very good performance and has the advantage of dealing well with small ground training samples, which is of particular interest for this case study [<xref ref-type="bibr" rid="scirp.95802-ref27">27</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref28">28</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref29">29</xref>].</p></sec><sec id="s3_3"><title>3.3. Accuracy Assessment</title><p>The accuracy assessment was based on confusion matrices between classification maps and training samples. For this purpose, two indices were evaluated which are the overall accuracy and Kappa index.</p><p>The overall accuracy is calculated by summing the number of correctly classified values and dividing by the total number of values according to the equation below [<xref ref-type="bibr" rid="scirp.95802-ref30">30</xref>]:</p><p>O A = C A ∗ 100 (1)</p><p>where:</p><p>OA is the overall classification accuracy.</p><p>C is the number of correct points.</p><p>A is the total number of reference points.</p><p>The kappa coefficient measures the agreement between classification and real values. A kappa value of 1 represents perfect agreement, while a value of 0 represents no agreement [<xref ref-type="bibr" rid="scirp.95802-ref30">30</xref>]. The kappa coefficient is computed as follows:</p><p>k = N ∑ i = 1 n m i , i − ∑ i = 1 n ( G i C i ) N 2 − ∑ i = 1 n ( G i C i ) (2)</p><p>where:</p><p>i is the class number.</p><p>N is the total number of classified values compared to truth values.</p><p>m<sub>i</sub><sub>,</sub><sub>i</sub> is the number of values belonging to the truth class i that have also been classified as class i (i.e., values found along the diagonal of the confusion matrix).</p><p>C<sub>i</sub> is the total number of predicted values belonging to class i.</p><p>G<sub>i</sub> is the total number of truth values belonging to class i.</p></sec><sec id="s3_4"><title>3.4. Change Detection Analysis</title><p>Change detection analysis describes and quantifies the differences between images of the same scene at different times. This analysis was performed with ENVI software, providing a simple approach to measure changes between a pair of images that represent an initial state and final state.</p></sec></sec><sec id="s4"><title>4. Results and Discussion</title><sec id="s4_1"><title>4.1. LULC Classification and Accuracy Assessment</title><p>Based on the accuracy assessment, it can be shown that the accuracy of the three produced maps is very similar (higher than 84%). The overall accuracies and kappa indices are, respectively, 88% and 84% for 1985, whereas years 2002 and 2017 show an overall accuracy of 93% and 88.4%, while the kappa indices are around 90.6% and 85%, respectively. <xref ref-type="table" rid="table2">Table 2</xref> shows the producer and user accuracy for the related years. The choice of the SVM classifier made this task easier, especially for 1985 and 2002 images, where there is a lack of the data needed to establish the different samples classes ROIs (Regions of Interest). As explained above, this classification is based on the quality, not on the number of samples. Therefore, we have been able to use topographic maps, field data and data collected from the High Commission for Water and Forests and the Kenitra Urban Agency. This approach allows for a reliable classification and assessment of the obtained results using a limited number of quality samples. The classified images of LUCL for the years 1985 and 2002 and 2017 are presented in <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p><p>The Landsat images used in this study were taken in the same month of the year. All images belong to the summer season, the choice of season is based on</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Accuracy assessment of classification results based on the confusion matrix for the years 1985, 2002 and 2017</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle"  colspan="2"  >1985</th><th align="center" valign="middle"  colspan="2"  >2002</th><th align="center" valign="middle"  colspan="2"  >2017</th></tr></thead><tr><td align="center" valign="middle" >Acc. %</td><td align="center" valign="middle" >Prod</td><td align="center" valign="middle" >User</td><td align="center" valign="middle" >Prod</td><td align="center" valign="middle" >User</td><td align="center" valign="middle" >Prod</td><td align="center" valign="middle" >User</td></tr><tr><td align="center" valign="middle" >Agricultural land</td><td align="center" valign="middle" >92.65</td><td align="center" valign="middle" >85.94</td><td align="center" valign="middle" >93.34</td><td align="center" valign="middle" >91.89</td><td align="center" valign="middle" >89.47</td><td align="center" valign="middle" >81.46</td></tr><tr><td align="center" valign="middle" >Shrub and Forest land</td><td align="center" valign="middle" >85.76</td><td align="center" valign="middle" >77.32</td><td align="center" valign="middle" >87.50</td><td align="center" valign="middle" >93.84</td><td align="center" valign="middle" >55.14</td><td align="center" valign="middle" >81.19</td></tr><tr><td align="center" valign="middle" >Bare land</td><td align="center" valign="middle" >72.13</td><td align="center" valign="middle" >58.96</td><td align="center" valign="middle" >80.21</td><td align="center" valign="middle" >77.75</td><td align="center" valign="middle" >87.37</td><td align="center" valign="middle" >82.37</td></tr><tr><td align="center" valign="middle" >Water body</td><td align="center" valign="middle" >98.80</td><td align="center" valign="middle" >99.64</td><td align="center" valign="middle" >98.98</td><td align="center" valign="middle" >99.77</td><td align="center" valign="middle" >98.94</td><td align="center" valign="middle" >99.70</td></tr><tr><td align="center" valign="middle" >Beaches dunes sand</td><td align="center" valign="middle" >78.87</td><td align="center" valign="middle" >95.79</td><td align="center" valign="middle" >94.89</td><td align="center" valign="middle" >91.24</td><td align="center" valign="middle" >99.46</td><td align="center" valign="middle" >90.20</td></tr><tr><td align="center" valign="middle" >Sparse vegetation</td><td align="center" valign="middle" >62.88</td><td align="center" valign="middle" >84.83</td><td align="center" valign="middle" >75.20</td><td align="center" valign="middle" >80.00</td><td align="center" valign="middle" >71.39</td><td align="center" valign="middle" >75.74</td></tr><tr><td align="center" valign="middle" >Built-up</td><td align="center" valign="middle" >85.83</td><td align="center" valign="middle" >80.05</td><td align="center" valign="middle" >82.43</td><td align="center" valign="middle" >92.74</td><td align="center" valign="middle" >89.31</td><td align="center" valign="middle" >93.99</td></tr></tbody></table></table-wrap><p>our needs to avoid seasonal plants and is to highlight the degradation or evolution of shrubs and forests in the study area. The corresponding classifications have been presented in the form of maps and graphs, which are illustrated in <xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref>. <xref ref-type="fig" rid="fig4">Figure 4</xref> shows the area corresponding to each LULC category and its percentage of the total area. Accordingly, a significant increase in built-up areas from 6.7 km<sup>2</sup> (0.43%) in 1985 to 49.97 km<sup>2</sup> (3.45%) in 2017 can be noted. On the other hand, it can be observed that the class occupying the majority of the study area since 1985 is the agricultural land which has been known a general increase from 333.1 km<sup>2</sup> (23%) to 468.6 (32.37%) with a peak in 2002, expressed by a surface of 496.3 km<sup>2</sup> (34.31%). Since 1985 Shrubs and Forest land class has known a phenomenal drop from 310.6 km<sup>2</sup> (21.45%) to 74.2 km<sup>2</sup> (5.42%) in 2017. The bare land, water body, beaches dunes sand and sparse vegetation classes do not show significant changes representing a variation of 1 to 3% over the last 32 years.</p></sec><sec id="s4_2"><title>4.2. Change Detection</title><p>The analysis of the magnitude of changes for the periods 1985-2002, 2002-2017 and 1985-2017 shows that the study area has undergone very significant changes, especially in three categories, which are shrubs and forest land, agricultural land and built-up (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p><sec id="s4_2_1"><title>4.2.1. Change in Built-Up</title><p>At a rapid rate, LULC has changed in the cities of Kenitra and Mehdia over the last 32 years, as a result of a sevenfold increase in built-up space between 1985 and 2017. During the periods 1985-2002 and 2002-2017, the built-up area increased by 14.6 km<sup>2</sup> and 29.2 km<sup>2</sup> respectively (<xref ref-type="fig" rid="fig5">Figure 5</xref>). Due to the decline in mortality, the agriculture development, the establishment of infrastructure for economic and social purposes and the establishment of too industrial parks (e.g. the free zone of Kenitra); the city has become the target of the rural exodus. Population census data between 1982 and 2014 shows a higher growth rate in this area than the national recorded rate. Indeed, the growth percentage reached 2.3% between 1982 and 1994, 1.5% between 1994 and 2004 then a rate of 1.3% between 2004 and 2014 against 2.06%, 1.3% and 1.25% respectively for the national growth rate [<xref ref-type="bibr" rid="scirp.95802-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref31">31</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref32">32</xref>]. This population increase directly involves an expansion of the urban area.</p><p>A built-up evolution map shows that excessive growth has mainly occurred south of the Sebou river, following three main directions (<xref ref-type="fig" rid="fig6">Figure 6</xref>), which are the south-east, south-west and west directions. South-east direction has expanded starting from the old medina till the west and south boundaries of Si Ali Boucheta marshes. On the other hand, the south-west direction starts from the old medina’s west boundary till the highway in the south. An important gap is located between the first and second directions. It can be explained by two reasons: the forest protection status (Kenitra Urban forest), as well as the elevation of the area (80 m), compared to the south-east and south-west extensions (20 - 30 m). The creation and expansion of industrial zones, which can have a direct effect on directions of expansion, such as the creation of the Bir Rami industrial zone in 1984 in the extreme southwest of the city. the expansion in 1982 of the former industrial zone created in colonial times since 1912, northeast of marshes Si Ali Boucheta. The third evolution direction, which occurs along the military airport till Mehdia city, can be explained by the touristic and economic attractiveness of the Atlantic coast.</p></sec><sec id="s4_2_2"><title>4.2.2. Change in Shrub and Forest Land</title><p>The forest domain of the study area shows enormous damage, with a loss of 75% of its area between 1985 and 2017 illustrated in the map below (<xref ref-type="fig" rid="fig7">Figure 7</xref>). The results obtained from the change detection statistics show that the greatest damage was recorded in the period 1985-2004 with a loss of 192.6 km<sup>2</sup> while during the period between 2004 and 2017 the losses were not significant (43.8 km<sup>2</sup>) (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p><p>According to the map of change of the shrub and forest land category, it can be seen that the major part of the class is located in the south of the Sebou River. This area is the northwestern limit of the Mamora forest considered to be the largest continuous plain cork oak forests in the world. It covered 133,000 ha at the beginning of the 20th century [<xref ref-type="bibr" rid="scirp.95802-ref33">33</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref34">34</xref>]. The main forest species are cork oaks, eucalyptus, acacias and pines according to the High Commissariat of Water and Forests. Morocco experienced between 1993 and 2004 the most severe periods of drought [<xref ref-type="bibr" rid="scirp.95802-ref35">35</xref>], hence, cork exploitation exploded during this period as a result of the increased demand for wood for socio-economic needs (the delimbing). Consequently, a loss of more than two thirds of the shrubs and forest land classes was reported between 1985 and 2017.</p><p>Despite the efforts made by the state to preserve the forest, it still suffers from a severe dysfunction. According to the technical report of the Food and Agriculture Organization of the United Nations FAO carried out in 2015 [<xref ref-type="bibr" rid="scirp.95802-ref36">36</xref>], the agents and causes of deforestation in the Mamora forest are: livestock farming, the creation of agricultural land, domestic and commercial firewood harvesting, coal production, urbanization, industrial and road infrastructure, and forest fires.</p></sec><sec id="s4_2_3"><title>4.2.3. Change in Agricultural Land</title><p>The agricultural sector is one of the major assets of the economic activity in the study area. It benefits from different local factors such as a diverse range of high-quality soils, abundant water resources, availability of labor and proximity to Europe and major consumption centers. All these factors make agriculture the leading activity in the region.</p><p>The agricultural lands experienced a general evolution between 1985 and 2017 with a spread of 135.6 km<sup>2</sup>, knowing that between 1985 and 2002 this category reached 163.7 km<sup>2</sup> before it decreased by 23.7 km<sup>2</sup> between 2002 and 2017 (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p><p>Since 1985, the agricultural land has occupied the north part of Kenitra city and the Sebou river lowlands as shown on the map of agricultural land change (<xref ref-type="fig" rid="fig8">Figure 8</xref>). From 1985 to 2002 the expansion followed a south-eastern direction and occupies especially the valleys of the dry wadis. Between 2004 and 2017 the agricultural land reached the Atlantic coast occupying the largest part of the coastal band.</p></sec><sec id="s4_2_4"><title>4.2.4. Change Detection Map</title><p>In order to identify anthropogenic action in the study area, we have clustered the LULC change trajectories together, into six clusters (<xref ref-type="fig" rid="fig9">Figure 9</xref>), which are urbanization, cultivation, deforestation, afforestation, abandonment and no change. The urbanization cluster includes all classes transformed into built-up. The cultivation cluster includes all classes transformed into agricultural land. The deforestation cluster includes changes from the shrub and forest land class to the other classes, while the afforestation cluster includes all classes transformed into shrub and forest land. The abandonment cluster includes all classes transformed</p><p>into spears vegetation and bare land. On the other hand, the no change cluster shows the areas that have not undergone any transformation (e.g. built-up to built-up, agricultural land to agricultural land, etc.).</p><p>It is clear that the study area has undergone enormous changes in LULC, those changes due mainly to human activity given the increase in the population, which implies land needs for the expansion of the city in wood and building materials as well as an expansion and search for new agricultural lands [<xref ref-type="bibr" rid="scirp.95802-ref34">34</xref>]. This has perfectly led to an imbalance in the fragile ecosystems of the Sebou estuary. This evolution implies, in addition to changes in LULC, chemical pollution of the estuary’s water and soil [<xref ref-type="bibr" rid="scirp.95802-ref37">37</xref>] [<xref ref-type="bibr" rid="scirp.95802-ref38">38</xref>].</p></sec></sec></sec><sec id="s5"><title>5. Conclusions</title><p>The understanding of the LULC changes and the environmental assessment are the basis for developing an appropriate urban planning policy with a focus on sustainable development. In this study, we produced LULC maps based on multi-temporal Landsat images from 1985, 2002 and 2017. They were used to produce LULC maps and assess the impact of land use change over the last three decades in the coastal area of Sebou Estuary, Northwest of the Moroccan Atlantic coast. The present study provided a piece of evidence of significant changes in LULC, especially at the level of agricultural land, shrub and forest land and built-up categories. The shrub and forest have decreased dramatically from 310.6 km<sup>2</sup> (21.5%) to 74.2 km<sup>2</sup> (5.4%). The built-up area increased from 6.7 km<sup>2</sup> (0.4%) in 1985 to 50 km<sup>2</sup> (3.5%) in 2017. The agricultural land increased from 333.9 km<sup>2</sup> (23%) to 468.6 (32.4%).</p><p>The phenomenal evolution of the built-up areas, the degradation of the forest as well as the green spaces in the study area, the expansion and creation of agricultural land are leading to a degradation of environmental quality throughout the study area.</p><p>Over the last 32 years, Kenitra city and its surroundings have become more and more a large urban region. We recognize that urban development may be beneficial for a variety of social and economic reasons in a region, such as Kenitra. However, it seems necessary to take into consideration the environmental and landscape concerns more seriously in order to achieve sustainable development.</p><p>LULC mapping and documentation may not provide the ultimate explanation for all the problems associated with environmental degradation. They can provide a descriptive overview of the evolution of the causes of anthropogenic expansion, but the impact of this evolution cannot be assessed without several biogeochemical studies on sediments and estuarine waters. However, it is one of the most important steps for better understanding of trends and possible causes of this degradation.</p></sec><sec id="s6"><title>Acknowledgements</title><p>This research was supported by LGMSS-URAC45 Laboratory in the Faculty of Sciences El Jadida. The authors would like to thank the anonymous reviewers for their helpful and precious suggestions that largely improved the original version of the manuscript.</p></sec><sec id="s7"><title>Author Contributions</title><p>All the authors contributed extensively to this paper. Ayoub El Adnani conceived and designed the research, processed and analyzed the data, and wrote the manuscript draft. Adnane Habib helped conceive the research and review the manuscript. Additional reviews and edits were done by Khalid El Khalidi. Bendahhou Zourarah supervised the research and validated the manuscript.</p></sec><sec id="s8"><title>Conflicts of Interest</title><p>The authors declare no conflict of interest.</p></sec><sec id="s9"><title>Cite this paper</title><p>El Adnani, A., Habib, A., El Khalidi, K. and Zourarah, B. (2019) Spatio-Temporal Dynamics and Evolution of Land Use Land Cover Using Remote Sensing and GIS in Sebou Estuary, Morocco. Journal of Geographic Information System, 11, 551-566. https://doi.org/10.4236/jgis.2019.115034</p></sec></body><back><ref-list><title>References</title><ref id="scirp.95802-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Livi-Bacci, M. (2017) A Concise History of World Population. 
https://doi.org/10.1002/9781119406822</mixed-citation></ref><ref id="scirp.95802-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">World Population Prospect (2007) World Population Prospects: The 2006 Revision Volume II: Sex and Age Distribution of the World Population.</mixed-citation></ref><ref id="scirp.95802-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Cohen, B. (2006) Urbanization in Developing Countries: Current Trends, Future Projections, and Key Challenges for Sustainability. Technology in Society, 28, 63-80. 
https://doi.org/10.1016/j.techsoc.2005.10.005</mixed-citation></ref><ref id="scirp.95802-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Houpin, S. (2011) Urban Mobility and Sustainable Development in the Mediterranean: Regional Diagnostic Outlook, Blue Plan.</mixed-citation></ref><ref id="scirp.95802-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">CERED (2005) Démographie Marocaine: Tendances passées et perspectives d ’avenir.</mixed-citation></ref><ref id="scirp.95802-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Shi, Z.H., Chen, L.D., Hao, J.P., Wang, T.W. and Cai, C.F. (2009) The Effects of Land Use Change on Environmental Quality in the Red Soil Hilly Region, China: A Case Study in Xianning County. Environmental Monitoring and Assessment, 150, 295-306. https://doi.org/10.1007/s10661-008-0231-8</mixed-citation></ref><ref id="scirp.95802-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Dale, V.H. (1997) The Relationship between Land-Use Change and Climate Change. Ecological Applications, 7, 753-769.  
https://doi.org/10.1890/1051-0761(1997)007[0753:TRBLUC]2.0.CO;2</mixed-citation></ref><ref id="scirp.95802-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Imbernon, J. (1999) Pattern and Development of Land-Use Changes in the Kenyan Highlands Since the 1950s. Agriculture, Ecosystems &amp; Environment, 76, 67-73. 
https://doi.org/10.1016/S0167-8809(99)00061-4</mixed-citation></ref><ref id="scirp.95802-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Meyer, W.B. and Turner, B.L. (1994) Changes in Land Use and Land Cover: A Global Perspective.</mixed-citation></ref><ref id="scirp.95802-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Jain, M., Dawa, D., Mehta, R., Dimri, A.P. and Pandit, M.K. (2016) Monitoring Land Use Change and Its Drivers in Delhi, India Using Multi-Temporal Satellite Data. Modeling Earth Systems and Environment, 2, 19.  
https://doi.org/10.1007/s40808-016-0075-0</mixed-citation></ref><ref id="scirp.95802-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Chen, L., Wang, J., Fu, B. and Qiu, Y. (2001) Land-Use Change in a Small Catchment of Northern Loess Plateau, China. Agriculture, Ecosystems &amp; Environment, 86, 163-172. https://doi.org/10.1016/S0167-8809(00)00271-1</mixed-citation></ref><ref id="scirp.95802-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Kilic, S., Evrendilek, F., Berberoglu, S. and Demirkesen, A.C. (2006) Environmental Monitoring of Land-Use and Land-Cover Changes in a Mediterranean Region of Turkey. Environmental Monitoring and Assessment, 114, 157-168.  
https://doi.org/10.1007/s10661-006-2525-z</mixed-citation></ref><ref id="scirp.95802-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Turner, A.B.L., Meyer, W.B. and Skole, D.L. (1994) Global Land-Use/Land-Cover Change: Towards an Integrated Study. Ambio, 23, 91-95.</mixed-citation></ref><ref id="scirp.95802-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Choudhury, D., Das, K. and Das, A. (2018) Assessment of Land Use Land Cover Changes and Its Impact on Variations of Land Surface Temperature in Asansol-Durgapur Development Region. The Egyptian Journal of Remote Sensing and Space Science, 22, 203-218. https://doi.org/10.1016/j.ejrs.2018.05.004</mixed-citation></ref><ref id="scirp.95802-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Güler, M., Yomraliogu, T. and Reis, S. (2007) Using Landsat Data to Determine Land Use/Land Cover Changes in Samsun, Turkey. Environmental Monitoring and Assessment, 127, 155-167. https://doi.org/10.1007/s10661-006-9270-1</mixed-citation></ref><ref id="scirp.95802-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Mundia, C.N. andAniya, M. (2005) Analysis of Land Use/Cover Changes and Urban Expansion of Nairobi City Using Remote Sensing and GIS. International Journal of Remote Sensing, 26, 2831-2849. https://doi.org/10.1080/01431160500117865</mixed-citation></ref><ref id="scirp.95802-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Yuan, F., Sawaya, K.E., Loeffelholz, B.C. and Bauer, M.E. (2005) Land Cover Classification and Change Analysis of the Twin Cities (Minnesota) Metropolitan Area by Multitemporal Landsat Remote Sensing. Remote Sensing of Environment, 98, 317-328. https://doi.org/10.1016/j.rse.2005.08.006</mixed-citation></ref><ref id="scirp.95802-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Yildirim, H., Alparslan, E. and Ozel, M.E. (1995) Temporal Change Detection by Principal Component Transformation on Satellite Imagery. 1995 International Geoscience and Remote Sensing Symposium, Firenze, Italy, 10-14 July 1995, 1227-1229.</mixed-citation></ref><ref id="scirp.95802-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Gidey, E., Dikinya, O., Sebego, R., Segosebe, E. and Zenebe, A. (2017) Modeling the Spatio-Temporal Dynamics and Evolution of Land Use and Land Cover (1984-2015) Using Remote Sensing and GIS in Raya, Northern Ethiopia. Modeling Earth Systems and Environment, 3, 1285-1301.  
https://doi.org/10.1007/s40808-017-0375-z</mixed-citation></ref><ref id="scirp.95802-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Mergaoui, L., Fekhaoui, M., Bouya, D., Gheit, A. and Stambouli, A. (2003) Qualité des eaux et macrofaune benthique d’un milieu estuarien du Maroc: Cas de l’estuaire de Sebou. Bulletin de l’Institut Scientifique, section Sciences de la Vie, No. 25, 67-75.</mixed-citation></ref><ref id="scirp.95802-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">ORMVAG (2014) Données des stations météorologiques de la zone de l’Office Régional de Mise en Valeur Agricole du Gharb (ORMVAG). Kénitra.</mixed-citation></ref><ref id="scirp.95802-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">El Jihad, M.-D., Peyrusaubes, D. and El Bouzidi, A. (2014) Sécheresses saisonnières et changement climatique dans le Gharb (Maroc). Rur @lités, 4, 14-25.</mixed-citation></ref><ref id="scirp.95802-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Young, N.E., Anderson, R.S., Chignell, S.M., Vorster, A.G., Lawrence, R. and Evangelista, P.H. (2017) A Survival Guide to Landsat Preprocessing. Ecology, 98, 920-932. https://doi.org/10.1002/ecy.1730</mixed-citation></ref><ref id="scirp.95802-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Abreu, L.W. and Anderson, G.P. (1996) The MODTRAN 2/3 Report and LOWTRAN 7 MODEL.</mixed-citation></ref><ref id="scirp.95802-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Berk, A., Feldeb, A.J.R.G.W., Anderson, G.P., Hoke, M.L. and Cooky, T. (2002) Exploiting Atmospheric Correction: The FLAASH Algorithm. Proceedings of 5th International Conference on Information Fusion, FUSION 2002, Annapolis, MD, 8-11 July 2002, 798-803.</mixed-citation></ref><ref id="scirp.95802-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Chi, M., Feng, R. and Bruzzone, L. (2008) Classification of Hyperspectral Remote-Sensing Data with Primal SVM for Small-Sized Training Dataset Problem. Advances in Space Research, 41, 1793-1799.  
https://doi.org/10.1016/j.asr.2008.02.012</mixed-citation></ref><ref id="scirp.95802-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Mountrakis, G., Im, J. and Ogole, C. (2010) Support Vector Machines in Remote Sensing: A Review. I ISPRS Journal of Photogrammetry and Remote Sensing, 1, 1-13.</mixed-citation></ref><ref id="scirp.95802-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Fernandez-Delgado, M., Cernadas, E. and Barro, S. (2014) Transboundary Flows of Western Entertainment across the Iron Curtain. Journal of Machine Learning Research, 15, 3133-3181.</mixed-citation></ref><ref id="scirp.95802-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Melgani, F. and Bruzzone, L. (2004) Classification of Hyperspectral Remote Sensing Images with Support Vector Machines. IEEE Transactions on Geoscience and Remote Sensing, 42, 1778-1790. https://doi.org/10.1109/TGRS.2004.831865</mixed-citation></ref><ref id="scirp.95802-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Conglaton, R.G. (1991) A Review of Assessing the Accuracy of Classifications of Remotely Sensed Data. Remote Sensing of Environment, 37, 35-46. 
https://doi.org/10.1016/0034-4257(91)90048-B</mixed-citation></ref><ref id="scirp.95802-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">RGPH and HCP (2014) Caracteristiques socio-économique et démographique de la population de la province de Kénitra.</mixed-citation></ref><ref id="scirp.95802-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">RGPH and HCP (2004) Caracteristiques demographiques et socio-economiques: Region de Gharb-Chrarda-Beni Hssen.</mixed-citation></ref><ref id="scirp.95802-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">FFEM, FAO, ON-FInternational, and HCEFLCD (2015) Caractérisation des agents et causes la déforestation et de la dégradation forestière dans le site de la Maamora au Maroc.</mixed-citation></ref><ref id="scirp.95802-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Noumonvi, K.D., Mounir, F. and Belghazi, B. (2017) Spatial Multi-Criteria Based Analysis to Assess Dynamics and Vulnerability of Forest Ecosystems to Global Changes: Case of Maamora Forest-Morocco. OALib: A New Open Access Journal, 4, 1-16.</mixed-citation></ref><ref id="scirp.95802-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Stour, L. and Agoumi, A. (2008) Sécheresse climatique au Maroc durant les dernières décennies (Climatic Drought in Morocco during the Last Decades) (Translated by Author). Hydroécologie Appliquée, 16, 215-232.  
https://doi.org/10.1051/hydro/2009003</mixed-citation></ref><ref id="scirp.95802-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">FAO (2015) Agents et causes de la déforestation et dégradation dans les sites pilotes du projet.</mixed-citation></ref><ref id="scirp.95802-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Cheggour, M., Chafik, A., Fisher, N.S. and Benbrahim, S. (2005) Metal Concentrations in Sediments and Clams in Four Moroccan Estuaries. Marine Environmental Research, 59, 119-137. https://doi.org/10.1016/j.marenvres.2004.04.002</mixed-citation></ref><ref id="scirp.95802-ref38"><label>38</label><mixed-citation publication-type="other" xlink:type="simple">Loukili, A. and Belghyti, D. (2007) The Dynamics of the Nematode Anguillicola crassus, Kuvahara 1974 in Eel Anguilla anguilla (L. 1758) in the Sebou estuary (Morocco). Parasitology Research, 100, 683-686.  
https://doi.org/10.1007/s00436-006-0349-y</mixed-citation></ref></ref-list></back></article>