<?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">ARS</journal-id><journal-title-group><journal-title>Advances in Remote Sensing</journal-title></journal-title-group><issn pub-type="epub">2169-267X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ars.2021.103006</article-id><article-id pub-id-type="publisher-id">ARS-112019</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Forest Dynamics with Sentinel 2 in Antanambe between 2005 and 2016 with the Snap Tool
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Tsiorinantenaina</surname><given-names>René Rakotoarison</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>Aimé</surname><given-names>Richard Hajalalaina</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>Elysa</surname><given-names>Nambinintsoa Safidinirina</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib></contrib-group><aff id="aff3"><addr-line>School of Management and Technological Innovation, University of Fianarantsoa, Fianarantsoa, Madagascar</addr-line></aff><aff id="aff2"><addr-line>Computer Science, Geomatics, Mathematics and Applications, Hosting Team Fianarantsoa, Fianarantsoa, Madagascar</addr-line></aff><aff id="aff1"><addr-line>Laboratory of Computer Science and Mathematics Applied to Development, University of Fianarantsoa, Fianarantsoa, Madagas-car</addr-line></aff><pub-date pub-type="epub"><day>28</day><month>07</month><year>2021</year></pub-date><volume>10</volume><issue>03</issue><fpage>92</fpage><lpage>101</lpage><history><date date-type="received"><day>30,</day>	<month>July</month>	<year>2021</year></date><date date-type="rev-recd"><day>15,</day>	<month>September</month>	<year>2021</year>	</date><date date-type="accepted"><day>18,</day>	<month>September</month>	<year>2021</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>
 
 
  In order to protect and sustainably manage the forest in Madagascar, which is 
  currently one of the countries still covered by forests, it is essential to use
   technological advances, particularly with regard to remote sensing. It provides valuable data, and sometimes free with a wide range of spatial, spectral and temporal resolutions to meet the demands for information on forest resources that are increasingly numerous and requires ever increasing levels of accuracy. The present work presents a methodology for the analysis of forest dynamics in the Antanambe area for the period 2005-2016 for monitoring forest degradation in this forest area to be conserved. The Random Forest algorithm was used to classify a Sentinel 2 image collected on November 07, 2016 and compare with a classification result with LandSat 5 in 2005 to detect change. The per-pixel change detection of both results captured the change map to better interpret the situation.
 
</p></abstract><kwd-group><kwd>Random Forest</kwd><kwd> Detection Change</kwd><kwd> Remote Sensing</kwd><kwd> Forest</kwd><kwd> Madagascar</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The demands for information on forest ecosystems are increasing and require ever increasing levels of precision. This calls for technological advances, especially in mapping and remote sensing means (diversity and specificity of sensors and platforms) provide valuable, and sometimes free, data for a wide range of spatial, spectral and temporal resolutions [<xref ref-type="bibr" rid="scirp.112019-ref1">1</xref>]. There are several research works that are concerned with forest dynamics in Madagascar such as the works of [<xref ref-type="bibr" rid="scirp.112019-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.112019-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.112019-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.112019-ref5">5</xref>]. In 2007, the analysis of forest cover evolution in Madagascar is done by processing Landsat satellite images that cover the period 1990-2000-2005. This analysis gives for 1990, an area significantly less than 10.7 million hectares (RPP, 2014). This difference can relatively be identified by the fact that the 2007 analysis used a more strict determination of forest. This study estimates the total forest area to be about 9,990,000 hectares in 2000 and about 9,725,000 hectares in 2005, which means an annual deforestation rate of 0.83% for the period 1990-2000 and 0.53% between 2000 and 2005 [<xref ref-type="bibr" rid="scirp.112019-ref6">6</xref>]. However, this study presents an analysis of forest dynamics in the Mananara Nord Madagascar area for the period 2005-2016 with the Random Forest (RF) classification algorithm. Several studies have used and highlighted the effectiveness of the RF algorithm among existing algorithms. Based on the accuracy evaluation and visual interpretation of the resulting maps, Schneider (2012) claimed that RF is more effective than the maximum likelihood classifier and support vector machine classification. In this study, a Sentinel-2 image acquired on November 07, 2016 was used and compared to a classification result of the LandSat image in 2005 [<xref ref-type="bibr" rid="scirp.112019-ref7">7</xref>]. The results represent the statistics and mapping of the dynamics of the forests of Mananara Nord between 2005 and 2016.</p></sec><sec id="s2"><title>2. Materials and Methods</title><sec id="s2_1"><title>2.1. Study Area</title><p>The Mananara Nord National Park is located in the northeastern part of Madagascar. More precisely, this area is located in the district of Mananara Nord and the region of Analanjorofo. Its geographical coordinates are oriented between latitude 16˚10'00&quot; and longitude 49˚46'00&quot;. This part presents tropical rainforests. We simulated a sample on a part of this zone which the commune of Antanambe (<xref ref-type="fig" rid="fig1">Figure 1</xref>).</p></sec><sec id="s2_2"><title>2.2. Presentation of the Data</title><p>A Sentinel-2 image is used that was acquired on November 07, 2016. The Sentinel-2 image comes directly from the European Space Agency (ESA) download site or more precisely in https://scihub.copernicus.eu.</p><p>The Sentinel-2 satellite is launched as part of the European Union’s Copernicus program which is equipped with the MSI sensor, it consists of two satellites, Sentinel-2A and Sentinel-2B: the first was launched in June 2015 and the second in March 2016. Each of these satellites provides all-land observation every 5 days, 10 m to 60 m resolution. The MSI sensor will provide 13 spectral bands from visible to mid-infrared illustrated in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>In 2014, inventory work was carried out in Madagascar by [<xref ref-type="bibr" rid="scirp.112019-ref7">7</xref>] on forest stratification. Landsat 5 images are used for the year 2005. A classification result from this work was used in the study.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Sentinel-2 bands</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Bands</th><th align="center" valign="middle" >Resolution (m)</th><th align="center" valign="middle" >wavelength (nm)</th><th align="center" valign="middle" >Width of the bands (nm)</th></tr></thead><tr><td align="center" valign="middle" >B1</td><td align="center" valign="middle" >60</td><td align="center" valign="middle" >443</td><td align="center" valign="middle" >20</td></tr><tr><td align="center" valign="middle" >B2</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >490</td><td align="center" valign="middle" >65</td></tr><tr><td align="center" valign="middle" >B3</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >560</td><td align="center" valign="middle" >35</td></tr><tr><td align="center" valign="middle" >B4</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >665</td><td align="center" valign="middle" >30</td></tr><tr><td align="center" valign="middle" >B5</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >705</td><td align="center" valign="middle" >15</td></tr><tr><td align="center" valign="middle" >B6</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >740</td><td align="center" valign="middle" >15</td></tr><tr><td align="center" valign="middle" >B7</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >783</td><td align="center" valign="middle" >20</td></tr><tr><td align="center" valign="middle" >B8</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >842</td><td align="center" valign="middle" >115</td></tr><tr><td align="center" valign="middle" >B8A</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >865</td><td align="center" valign="middle" >20</td></tr><tr><td align="center" valign="middle" >B9</td><td align="center" valign="middle" >60</td><td align="center" valign="middle" >945</td><td align="center" valign="middle" >20</td></tr><tr><td align="center" valign="middle" >B10</td><td align="center" valign="middle" >60</td><td align="center" valign="middle" >1375</td><td align="center" valign="middle" >30</td></tr><tr><td align="center" valign="middle" >B11</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >1610</td><td align="center" valign="middle" >90</td></tr><tr><td align="center" valign="middle" >B12</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >2190</td><td align="center" valign="middle" >180</td></tr></tbody></table></table-wrap></sec><sec id="s2_3"><title>2.3. Tool</title><p>SNAP (Sentinel Application Platform) is an ESA product specifically for processing Sentinel images. SNAP contains Sentinel toolkits with common architectures developed by Brockmann Consult, SkyWatch and CS. SNAP is an open-source solution with a number of toolboxes and plugins.</p></sec><sec id="s2_4"><title>2.4. Methodology</title><p>The Sentinel-2 image that we acquired is at pre-processing level 2a, i.e. it is no longer necessary to make the radiometric and geometric correction. We continue the pre-processing by combining the mono-spectral bands into a single multi-spectral band because when downloading Sentinel-2 images to the ESA site, the images are still separated into channels, so we need a transformation to obtain a multi-spectral image. Then, we cut the study area in order to minimize the processing time.</p><p>The NDVI thresholding technique is used for a pre-classification to discriminate vegetation to other land cover types [<xref ref-type="bibr" rid="scirp.112019-ref8">8</xref>]. NDVI values can range from −1 and 1. Higher values of the index correspond to the presence of dense and healthy vegetation cover. They are generally between 0.1 and 0.7 while clouds and snow result in NDVI values close to 0 [<xref ref-type="bibr" rid="scirp.112019-ref9">9</xref>].</p><p>The supervised classification [<xref ref-type="bibr" rid="scirp.112019-ref10">10</xref>] was used to determine the forest cover of the study area. We implemented the existing Random forest algorithm in the SNAP tool which is particularly effective in identifying relationships between a variable to be explained and explanatory variables [<xref ref-type="bibr" rid="scirp.112019-ref11">11</xref>].</p><p>To wait for the objective, we performed change detection between the 2016 classified Sentinel-2 image and the 2005 classification result from LandSat 5 [<xref ref-type="bibr" rid="scirp.112019-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.112019-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.112019-ref14">14</xref>]. The methodology is illustrated in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p></sec></sec><sec id="s3"><title>3. Results and Discussions</title><sec id="s3_1"><title>3.1. Vegetation Index NDVI</title><p>We have here a pre-classification result after the NDVI calculation. This index uses the near infrared band (PIR) and the red band (R) of the satellite image. In the case of Sentinel-2 images, the PIR is the eighth band and the R is the fourth band.</p><p>NDVI = ( P I R − R ) / ( P I R + R ) (1)</p><p>By applying Equation (1), we have the result presented in <xref ref-type="fig" rid="fig3">Figure 3</xref> where the higher values of the index correspond to the presence of a dense and healthy vegetation cover.</p></sec><sec id="s3_2"><title>3.2. Random Forest</title><p>Supervised classification is used using the Random forest algorithm. It is a classification algorithm based on parallel learning on multiple decision trees that improves prediction accuracy [<xref ref-type="bibr" rid="scirp.112019-ref11">11</xref>]. After integrating the field validation points, the SNAP tool gave the confusion matrix presented in the <xref ref-type="table" rid="table2">Table 2</xref> below, which deduces an acceptable Kappa index of 94, 8863%. The result is presented in two classes (forest and non-forest) in <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Confusion matrix</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Classes</th><th align="center" valign="middle" >Forest</th><th align="center" valign="middle" >No forest</th><th align="center" valign="middle" >Production accuracy</th></tr></thead><tr><td align="center" valign="middle" >Forest</td><td align="center" valign="middle" >99.31</td><td align="center" valign="middle" >1.76</td><td align="center" valign="middle" >60.38</td></tr><tr><td align="center" valign="middle" >No forest</td><td align="center" valign="middle" >0.69</td><td align="center" valign="middle" >98.24</td><td align="center" valign="middle" >39.62</td></tr></tbody></table></table-wrap></sec><sec id="s3_3"><title>3.3. Change Detection</title><p>Change detection consists in comparing the classification in 2005 with Landsat 5 and the classification in 2016 with Sentinel 2 (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p><p>The transition matrix (<xref ref-type="table" rid="table3">Table 3</xref>) represents the percentage change in pixels for the two classes (Forest and No-Forest) between the two dates of 2005 and 2016.</p><p>According to the table, 70.044% of the forests classes and 75.652% of the no-forests classes remain stable the two dates of 2005 and 2016. During this 11-year period, the analysis shows a significant deforestation (Forest → Non-Forest) of 29.956% but also regeneration (Non-Forest → Forest) of 24.348%. The result is presented in <xref ref-type="fig" rid="fig6">Figure 6</xref>.</p><p>In terms of area, 16,687 ha of forest has been lost in 11 years. This means that there is a significant deforestation between two dates 2005 and 2016 (<xref ref-type="table" rid="table4">Table 4</xref>).</p><p>According to the Strategic Environmental and Social Assessment of Madagascar’s REDD+ program, the annual deforestation rate for the rainforest ecoregion</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Transition matrix</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="2"   rowspan="2"  >Percentage of pixels per class</th><th align="center" valign="middle"  colspan="2"  >Initial state (2005)</th></tr></thead><tr><td align="center" valign="middle" >Forest</td><td align="center" valign="middle" >No forest</td></tr><tr><td align="center" valign="middle" >Final state</td><td align="center" valign="middle" >Forest</td><td align="center" valign="middle" >73.044%</td><td align="center" valign="middle" >24.348%</td></tr><tr><td align="center" valign="middle" >(2016)</td><td align="center" valign="middle" >No forest</td><td align="center" valign="middle" >26.956%</td><td align="center" valign="middle" >75.652%</td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Surface changed</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Areas (Ha)</th><th align="center" valign="middle" ></th><th align="center" valign="middle" ></th></tr></thead><tr><td align="center" valign="middle" >Forest</td><td align="center" valign="middle" >63,071.2</td><td align="center" valign="middle" >46,383.3694</td></tr><tr><td align="center" valign="middle" >No forest</td><td align="center" valign="middle" >1288.1636</td><td align="center" valign="middle" >17,975.9942</td></tr></tbody></table></table-wrap><p>in Madagascar is 0.94% per year for the period 2010-2013. The Analanjorofo region is found among the regions with the highest annual deforestation rate. In terms of area, the Analanjorofo region lost 24,000 ha of forest for the period 2010-2013. Surveys conducted by ONE in 2016 were able to group four root causes of deforestation and forest degradation: The weakness of governance and institutions, the non-rational use of forest resources and forest areas, the weak coordination of the development of rural areas to the detriment of forests and the rural poverty and the lack of alternatives development alternatives.</p></sec></sec><sec id="s4"><title>4. Conclusion and Perspectives</title><p>In this experience, we retain the technique of monitoring deforestation by remote sensing. This computerized solution is a facilitating means for experts in the field of environmental management that offers interesting results generally without direct contact to all the studied covers. We performed a supervised classification of a part of the National Park of Mananara Nord, in the commune of Antanambe by the satellite image Sentinel 2 acquired in 2016. This result is compared to a LandSat 5 classification in 2005 to detect the change between these two dates. The analysis shows the existence of more deforestation than generation in the study area. Deforestation monitoring mobilizes many entities in Madagascar. The results play a very important role for decision making in the field of environmental conservation in Madagascar. For this, the discussion is open for the capitalization of the experience in the form of knowledge that can be reused in another study area or other input data. The technique is applicable to a wider coverage such as the whole of Madagascar by inserting task automation (workflow) and parallel processing in the distributed system (cloud computing, grid computing) in order to optimize processing time.</p></sec><sec id="s5"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s6"><title>Cite this paper</title><p>Rakotoarison, T.R., Hajalalaina, A.R. and Safidinirina, E.N. (2021) Forest Dynamics with Sentinel 2 in Antanambe between 2005 and 2016 with the Snap Tool. Advances in Remote Sensing, 10, 92-101. https://doi.org/10.4236/ars.2021.103006</p></sec></body><back><ref-list><title>References</title><ref id="scirp.112019-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Bouvier, M., Durrieu, S., Fournier, R.A. and Renaud, J.P. (2015) Generalizing Predictive Models of Forest Inventory Attributes Using an Area-Based Approach with Airborne LiDAR Data. Remote Sensing of Environment, 156, 322-334.  
https://doi.org/10.1016/j.rse.2014.10.004</mixed-citation></ref><ref id="scirp.112019-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Giri, C. and Muhlhausen, J. (2008) Mangrove Forest Distributions and Dynamics in Madagascar (1975-2005). Sensors, 8, 2104-2117. https://doi.org/10.3390/s8042104</mixed-citation></ref><ref id="scirp.112019-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Grinand, C., Rakotomalala, F., Gond, V., Vaudry, R., Bernoux, M. and Vieilledent, G. (2013) Estimating Deforestation in Tropical Humid and Dry Forests in Madagascar from 2000 to 2010 Using Multi-Date Landsat Satellite Images and the Random Forests Classifier. Remote Sensing of Environment, 139, 68-80.  
https://doi.org/10.1016/j.rse.2013.07.008</mixed-citation></ref><ref id="scirp.112019-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Hajalalaina, A.R., Grizonnet, M., Delaitre, E., Rakotondraompiana, S. and Herve, D. (2013) Discrimination des zones humides en foret malgache, proposition d’une methodologie multiresolution et multisource utilisant orfeo toolbox. Revue Franaise de Photogrammétrie et de Télédétection, 201, 37-48.  
https://doi.org/10.52638/rfpt.2013.44</mixed-citation></ref><ref id="scirp.112019-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Razafinimaro, A., Hajalalaina, A.R., Reziky, Z.T., Delaitre, E. and Andrianarivo, A. (2021) Landsat8 Satellite Image Classification with ERDAS for Mapping the Kalambatritra Special Reserve. American Journal of Remote Sensing, 9, 16.</mixed-citation></ref><ref id="scirp.112019-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Waeber, P.O., Wilmé, L., Ramamonjisoa, B., Garcia, C., Rakotomalala, D., Rabemananjara, Z.H., Kull, C.A., Ganzhorn, J.U. and Sorg, J.P. (2015) Dry Forests in Madagascar: Neglected and under Pressure. International Forestry Review, 17, 127-148. https://doi.org/10.1505/146554815815834822</mixed-citation></ref><ref id="scirp.112019-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">WCS, ONE, MNP, ETC (2014) Analyse Historique de la déforestation: 2005-2010-2013.  
https://bnc-redd.mg/images/documents/rapports/20170822/141210-FCC_051013_PERR-FH_2014.pdf</mixed-citation></ref><ref id="scirp.112019-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Gandhi, G.M., Parthiban, S., Thummalu, N. and Christy, A. (2015) Ndvi: Vegetation Change Detection Using Remote Sensing and Gis—A Case Study of Vellore District. Procedia Computer Science, 57, 1199-1210.  
https://doi.org/10.1016/j.procs.2015.07.415</mixed-citation></ref><ref id="scirp.112019-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Ngom, R. and Gosselin, P. (2013) Development of a Remote Sensing-Based Method to Map Likelihood of Common Ragweed (Ambrosia artemisiifolia) Presence in Urban Areas. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 7, 126-139. https://doi.org/10.1109/JSTARS.2013.2254469</mixed-citation></ref><ref id="scirp.112019-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Perumal, K. and Bhaskaran, R. (2010) Supervised Classification Performance of Multispectral Images. arXiv:1002.4046.</mixed-citation></ref><ref id="scirp.112019-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Breiman, L. (2001) Random Forests. Machine Learning, 45, 5-32.  
https://doi.org/10.1023/A:1010933404324</mixed-citation></ref><ref id="scirp.112019-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Dalmiya, C.P., Santhi, N. and Sathyabama, B. (2019) A Novel Feature Descriptor for Automatic Change Detection in Remote Sensing Images. The Egyptian Journal of Remote Sensing and Space Science, 22, 183-192.  
https://doi.org/10.1016/j.ejrs.2018.03.005</mixed-citation></ref><ref id="scirp.112019-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Stehman, S.V. and Foody, G.M. (2019) Key Issues in Rigorous Accuracy Assessment of Land Cover Products. Remote Sensing of Environment, 231, 111199.  
https://doi.org/10.1016/j.rse.2019.05.018</mixed-citation></ref><ref id="scirp.112019-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Negassa, M.D., Mallie, D.T. and Gemeda, D.O. (2020) Forest Cover Change Detection Using Geographic Information Systems and Remote Sensing Techniques: A Spatio-Temporal Study on Komto Protected Forest Priority Area, East Wollega Zone, Ethiopia. Environmental Systems Research, 9, 1.</mixed-citation></ref></ref-list></back></article>