<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "https://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">

<article article-type="research-article" dtd-version="1.1" xml:lang="en" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">

	<front>
		<journal-meta>
			<journal-id journal-id-type="publisher-id">Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA)</journal-id>
			<journal-title-group>
				<journal-title>Forest System</journal-title>
				<abbrev-journal-title>For Syst</abbrev-journal-title>
			</journal-title-group>
			<issn pub-type="epub">2171-9845</issn>
			<publisher>
				<publisher-name>Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA)</publisher-name>
			</publisher>
		</journal-meta>
		<article-meta>
			<article-id pub-id-type="publisher-id">20348</article-id>
			<article-id pub-id-type="doi">10.5424/fs/2023322-20348</article-id>
			<article-categories>
				<subj-group subj-group-type="heading">
					<subject>RESEARCH ARTICLE</subject>
				</subj-group>
			</article-categories>

			<title-group>
				<article-title>Monitoring the dynamic changes in vegetation cover and driving factors from 2000 to 2020 in the Maoershan Forest Farm region, China, using satellite remote sensing data</article-title>
			</title-group>

			<contrib-group>
				<contrib contrib-type="author">
					<contrib-id contrib-id-type="orcid">http://orcid.org/0009-0001-5210-5512</contrib-id>
					<name>
						<surname>Li</surname>
						<given-names>Teng </given-names>
					</name>
					<aff id="aff1"><institution>Liaoning Vocational College of Ecological Engineering, </institution><addr-line>Shenyang 110101, </addr-line><country>China.</country></aff>
				</contrib>

				<contrib contrib-type="author">
					<contrib-id contrib-id-type="orcid">http://orcid.org/0009-0007-2250-2596</contrib-id>
					<name>
						<surname>Gao</surname>
						<given-names>Yuanke</given-names>
					</name>
					<aff id="aff2"><institution>Department of Teaching and Research Management, Maoershan Teaching Area, Northeast Forestry University, </institution><addr-line>Shangzhi 150611, </addr-line><country>China.</country></aff>
				</contrib>

			</contrib-group>

			<pub-date pub-type="epub">
				<day>06</day>
				<month>07</month>
				<year>2023</year>
			</pub-date>			
			<pub-date pub-type="collection">
				<month>08</month>
				<year>2023</year>
			</pub-date>
			<volume>32</volume>
			<issue>2</issue>
			<elocation-id>e015</elocation-id>
			<history>
				<date date-type="received">
					<day>23</day>
					<month>03</month>
					<year>2023</year>
				</date>
				<date date-type="accepted">
					<day>06</day>
					<month>07</month>
					<year>2023</year>
				</date>
				<date date-type="pub">
					<day>06</day>
					<month>07</month>
					<year>2023</year>
				</date>
			</history>			
			<permissions>
				<copyright-statement>&#xa9;2023 CSIC</copyright-statement>
				<copyright-year>2023</copyright-year>
				<license license-type="open-access" xlink:href="https://creativecommons.org/licenses/by/4.0/">
					<license-p>This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC BY 4.0) License.</license-p>
				</license>
			</permissions>

			<self-uri xlink:href="https://doi.org/10.5424/fs/2023322-20348"/>
			<abstract>
				<title>Aim of study:</title>
				<p>Natural climate change is a central driver of global ecosystem and forest change. Climate change and topographical factors have had the greatest impact on different types of forests around the world. We used remote sensing technology to detect and analyze the temporal and spatial changes of forest vegetation to provide reference for regional management.</p>
				<title>Area of study:</title>
				<p>Maoershan Forest Farm, China.</p>
				<title>Material and methods:</title>
				<p>The Landsat images were preprocessed using ArcGIS and ENVI software. The normalized difference vegetation index (NDVI) was calculated to identify vegetation changes from 2000 to 2020. In addition, the vegetation fraction cover (VFC) was calculated using the pixel binary model. The driving factors and their influences on vegetation changes in this region were determined using the random forest algorithm and Pearson correlation analysis method.</p>
				<title>Main results:</title>
				<p>From 2000 to 2020, the NDVI showed an overall increasing trend. The results indicated that compared with the climatic factors, topographic factors were more important to vegetation growth in the study area. Among the topographic factors, elevation was the most important factor affecting vegetation growth and both showed a negative correlation. Among the climatic factors, relative humidity was the primary driving factor affecting vegetation growth and both showed a positive correlation.</p>
				<title>Research highlights:</title>
				<p>Accurate and timely assessment of vegetation change and its relationship to climate and topographical changes can provide very useful information for policy makers, governments and planners in formulating management policies.</p>
			</abstract>

			<kwd-group>
				<kwd>forest management;</kwd>
				<kwd>climate change;</kwd>
				<kwd>random forest;</kwd>
				<kwd>NDVI;</kwd>
			</kwd-group>

			<abbrev>DEM
				<def>(Digital Elevation Model)</def>
			</abbrev>
			<abbrev>GIMMS
				<def>(Global Inventory Modelling and Mapping Studies)</def>
			</abbrev>
			<abbrev>MODIS
				<def>(Moderate Resolution Imaging Spectroradiometer)</def>
			</abbrev>
			<abbrev>NDVI
				<def>(Normalized Difference Vegetation Index)</def>
			</abbrev>
			<abbrev>VFC
				<def>(Vegetation Fraction Coverage)</def>
			</abbrev>

			<supplementary-material>
				<label>Supplementary material</label>
				<caption>
					<p>(Figs. S1, S2 and S3) accompanies the paper on <em>Forest System</em>´s website.</p>
				</caption>
			</supplementary-material>			
							
		</article-meta>

		<funding-group id="fw-01">
			<award-group id="aw1">
				<funding-source>The authors received no specific funding for this work.</funding-source>
			</award-group>				
		</funding-group>	

	</front>
	
	<body>
		<sec id="sec1" sec-type="intro">
			<title>Introduction</title>

			<p>Vegetation cover is an important component of the terrestrial ecosystem. It is not only a linkage between the natural elements such as the atmosphere, soil, and water, but also an important indicator of global climate and ecosystem change (Seidl et al., 2017). In the past few years, with the intensification of global warming and human activities, monitoring and understanding the dynamics of vegetation cover in response to global change has become the core focus of earth system science. Understanding the vegetation cover development and changes contributes to a comprehensive understanding of the regional ecological environment as well as helps in evaluating the quality of the regional ecological environment system, which can provide a guide for formulating sustainable environmental protection strategies (Sun W et al., 2015). Therefore, studies on the long-term dynamic changes in the vegetation cover have become one of the hot spots in the field of ecology.</p>

			<p>At present, field-measured and remote sensing inversion are two basic approaches to monitoring changes in the vegetation cover (Ge et al., 2018). The surface measurement method primarily obtains the vegetation growth and cover status through a comprehensive investigation of the vegetation status and digital camera photographs. This method usually involves a detailed survey of the height, density, greenness and distribution of the vegetation. The ground survey method can provide accurate information on the vegetation structure and distribution at a small scale (Talukdar et al., 2020). This method is a way to obtain vegetation structure and distribution information through field survey, visual and manual measurement of vegetation height, greenness and distribution, and record these data. However, as they involve much manpower and material resources and the scope of the survey is small, they could not cover a large area, so it is difficult to effectively estimate the changes in the vegetation cover at a regional scale. With the development of remote sensing technology, the method based on remote sensing provides a unique opportunity for the inversion of vegetation cover in large regions and even the entire country.</p>

			<p>The normalized difference vegetation index (NDVI) is very sensitive to the biophysical characteristics of vegetation and is an effective indicator of vegetation growth status. It is the most commonly used remote sensing monitoring index for evaluating the changes in vegetation cover and one of the most effective parameters for evaluating vegetation status (Jiang et al., 2006). The vegetation fraction coverage (VFC) is a comprehensive quantitative index of the terrestrial plant community dynamics, and an important parameter of ecological, climate, and global change. It helps to provide basic and important details about the ecosystems (Amiri et al., 2009). Quantifying VFC is important for assessing the vegetation growth status and its changes at the regional and global scales.</p>

			<p>NDVI is the most widely used vegetation index with several studies being carried out to monitor regional long-term vegetation changes based on NDVI data provided by different sensors. The results obtained from these studies have been helpful for qualitative and quantitative evaluation of the growth and vigour of surface vegetation cover (Martínez &amp; Gilabert, 2009; Huang et al., 2021). Studies have also been conducted to monitor the dynamic changes in vegetation based on Satellite Pour l’ Observation de la Terre (SPOT) NDVI data or Moderate Resolution Imaging Spectroradiometer (MODIS) data to analyse the changes in vegetation cover and understand the vegetation growth dynamics (Klisch &amp; Atzberger, 2016). Pixel dichotomy assumes that each pixel of information received by the sensor is composed of green vegetation and bare soil. The pixel dichotomy model is developed to estimate the fraction of green vegetation and bare soil in the pixel, based on which the VFC of the pixel can be obtained (Li et al., 2014).</p>

			<p>The dynamic characteristics of the changes in vegetation cover are influenced by several external factors (Campbell et al., 2005). Vegetation growth in the natural state is mainly affected by climate change and terrain. In forests, the growth of vegetation is closely related to meteorological factors. All these factors, including temperature, humidity, precipitation, and wind affect the moisture adequacy of vegetation, thereby directly affecting the growth of forest vegetation (Sun W et al., 2015; Wang et al., 2015). For topographic factors, elevation and slope can affect the loss of water from vegetation, while the aspect factor can affect the amount of solar radiation received by the vegetation, which directly affects the degree of vegetation dryness (Tian et al., 2022). Jia et al. (2020) investigated the spatio-temporal variations in characteristics of vegetation cover in the Xijiang River Basin from 2000 to 2013; analysed the correlation between precipitation, slope, and vegetation index; and discussed the influence of climatic and topographic factors on the changes in vegetation cover in Xijiang River Basin. Sun et al. (2022) concluded that the variational trend of vegetation cover and its relationship with climatic factors in the source area of the three rivers in China were mainly affected by altitude and temperature; the vegetation cover of the Gannan Plateau in the upper reaches of the Yellow River had a significant decreasing trend.</p>

			<p>With a vast territory and complex natural conditions as well as with the rise in global climate change and human activities in recent years, especially due to global warming, rapid urbanization and large-scale ecological restoration projects since the end of the 20th century, vegetation coverage in China has experienced complex changes. However, there are still some gaps in the study on the long-term dynamic changes in the vegetation cover (Liu et al., 2019). First of all, there are gaps in understanding the fine-scale macro pattern of the changes in the vegetation cover, especially the spatio-temporal variational characteristics at long-time scales. Dynamic monitoring of the changes in the vegetation cover and determining its relationship with climatic and topographic factors are of great significance for preserving vegetation.</p>

			<p>Previous studies have primarily relied on time series analysis of moderate to low spatial resolution remote sensing images such as MODIS and GIMMS products to assess changes in vegetation cover at a large scale. However, for analysing regional changes in vegetation cover, it is crucial to use higher spatial resolution imagery. Therefore, in this study, we utilized Landsat imagery to study vegetation changes detection in Maoershan (Landsat 5 TM, Landsat 7 ETM, Landsat 8 OLI). We used Landsat images from 2000 to 2020 as our data source to calculate NDVI, and the pixel dichotomy method was used to calculate VFC to estimate the multi-period vegetation coverage in the study area (image time was 2000, 2005, 2010, 2015, 2020, five years interval). Furthermore, we analyzed changes in vegetation growth and cover during the aforementioned period, as well as the spatial distribution characteristics of vegetation cover and the correlation between vegetation growth and climatic and topographic factors in Maoershan. The key objective of this study was to monitor vegetation changes in Maoershan from 2000 to 2020 using Landsat images and to investigate the relationship between NDVI and climate and topographical changes.</p>

		</sec> <!--/sec 1-->

		<sec id="sec2" sec-type="materials|methods">
			<title>Material and methods</title>
			
			<sec id="sec2.1">
				<title>Study area</title>

				<p>The study area is located in Maoershan Experimental Forest Farm (45°14’-29’ N, 127°29’-44’ E) of Northeast Forestry University in Shangzhi City, Heilongjiang Province, which is the remnant of Xiaoling on the western slope of Zhangguangcai Mountain in Changbai Mountain system, with a total area of 26,496 hectares (Fig. 1). The landform of the experimental forest farm is a low hilly region as a whole. The terrain rises from south to north, with an average elevation of 300 m. The climate in this region is mild and humid, belonging to the continental monsoon climate type. The average annual temperature is about 3℃, and the average annual precipitation is about 723 mm. The major vegetation in this region is typical natural secondary forest, including hard broad-leaved mixed forest dominated by valuable broad-leaved tree species such as <em>Fraxinus mandshurica</em>, <em>Junglus mandshurica</em>, and <em>Phellodendron amurense</em>; soft broad-leaved mixed forest dominated by <em>Betula platyphalla</em>, <em>Populus davidiana</em>, and <em>Betula costata</em>; coniferous forest dominated by <em>Korean pine</em>, <em>Larix gmelina</em>, and <em>Pinus sylvestris</em> var. <em>mongolica</em>. The plantation primarily consists of <em>L. gmelini</em> and <em>P. sylvestris</em> var. mongolica, forming a mosaic of natural secondary forests and plantations (Chang et al., 2020).</p>

				<fig id="f1">
					<label>Figure 1</label>
					<caption>
						<title>Study area: Maoershan Experimental Forest Farm, Heilongjiang Province, China.</title>
					</caption>
					<graphic id="gra-1" xlink:href="img/e015-fig1.jpg"/>
				</fig>				

			</sec> <!--/sec 2.1-->

			<sec id="sec2.2">
				<title>Data sources</title>

				<p>The Landsat images used in this study were downloaded from the Geospatial Data Cloud (https://www.gscloud.cn/) with a spatial resolution of 30 m. In this study, remote sensing images collected by Landsat-5 TM, Landsat-7 ETM and Landsat-8 OLI were used to monitor vegetation status (Table 1). In the first step, radiometric calibration, atmospheric correction and geometric correction are performed using ENVI software based on the Landsat image of the study area. These data pre-processing steps must be done to obtain high-quality results. A total of 10 images were obtained in this study, and the specific data information is shown in Table 1.</p>

				<table-wrap id="t1">
					<label>Table 1</label>
					<caption>
						<title>Collected Landsat images. Source: Geospatial Data Cloud (www.gscloud.cn/).</title>
					</caption>
					<table>
						<thead>
							<tr>
								<th align="center">Data type</th>
								<th align="center">Acquisition date</th>
								<th align="center">Spatial resolution</th>
								<th align="center">Path/Row</th>
							</tr>
						</thead>
						<tbody>
						    <tr>
						        <td rowspan="2" valign="top">Landsat-5 TM</td>
						        <td align="center" rowspan="2" valign="top">2000-8-18</td>
						        <td align="center" rowspan="2" valign="top">30 m</td>
						        <td align="center">117/028</td>
						    </tr>
						    <tr>
						        <td align="center">117/029</td>
						    </tr>
						    <tr>
						        <td rowspan="2" valign="top">Landsat-7 ETM</td>
						        <td align="center" rowspan="2" valign="top">2005-8-08</td>
						        <td align="center" rowspan="2" valign="top">30 m</td>
						        <td align="center">117/028</td>
						    </tr>
						    <tr>
						        <td align="center">117/029</td>
						    </tr>
						    <tr>
						        <td rowspan="2" valign="top">Landsat-5 TM</td>
						        <td align="center" rowspan="2" valign="top">2010-9-15</td>
						        <td align="center" rowspan="2" valign="top">30 m</td>
						        <td align="center">117/028</td>
						    </tr>
						    <tr>
						        <td align="center">117/029</td>
						    </tr>
						    <tr>
						        <td rowspan="2" valign="top">Landsat-7 ETM</td>
						        <td align="center" rowspan="2" valign="top">2015-9-05</td>
						        <td align="center" rowspan="2" valign="top">30 m</td>
						        <td align="center">117/028</td>
						    </tr>
						    <tr>
						        <td align="center">117/029</td>
						    </tr>
						    <tr>
						        <td rowspan="2" valign="top">Landsat-8 OLI</td>
						        <td align="center" rowspan="2" valign="top">2020-7-08</td>
						        <td align="center" rowspan="2" valign="top">30 m</td>
						        <td align="center">117/028</td>
						    </tr>
						    <tr>
						        <td align="center">117/029</td>
						    </tr>
						</tbody>
					</table>
				</table-wrap>

			</sec> <!--sec2.2-->

			<sec id="sec2.3">
				<title>Vegetation cover change</title>

				<p>The NDVI is an important index for characterizing the vegetation coverage of the land surface, and its value ranges from -1 to 1. In this study, NDVI was used to indicate the vegetation growth status of the current year. VFC, an important index reflecting vegetation coverage and ecological environment monitoring, represents the ratio of the vertical projection area of the vegetation in a region to the total area of the region (Gao et al., 2020); it ranges from 0 to 1. VFC is usually calculated using the NDVI-based pixel binary model. In this study, each pixel of NDVI was decomposed into pure vegetation and pure bare soil. VFC was calculated as follows:</p>

				<disp-formula id="e1">
				<math id="e1">
					<mi>NDVI</mi>
				  	<mo>=</mo>
				    <mo>(</mo>
						<mi>NIR</mi>
						<mo>-</mo>
						<mi>R</mi>
					<mo>)</mo>
					<mo>/</mo>
				    <mo>(</mo>
				    	<mi>NIR</mi>
				    	<mo>+</mo>
				    	<mi>R</mi>
				    <mo>)</mo>
				</math>
				<label>(1)</label>
				</disp-formula>


				<disp-formula id="e2">
				<math id="e2">
					<mi>VFC</mi>
				  	<mo>=</mo>
				    <mo>(</mo>
						<mi>NDVI</mi>
						<mo>-</mo>
						<msub>
							<mi>NDVI</mi>
							<mi>soil</mi>
						</msub>
					<mo>)</mo>
					<mo>/</mo>
				    <mo>(</mo>
				    	<msub>
				    		<mi>NDVI</mi>
				    		<mi>veg</mi>
				    	</msub>
				    	<mo>-</mo>
						<msub>
				    		<mi>NDVI</mi>
				    		<mi>soil</mi>
				    	</msub>
				    <mo>)</mo>
				</math>
				<label>(2)</label>
				</disp-formula>

				<p>where <em>NIR</em> represents near infrared band, <em>R</em> represents red band, <em>VFC</em> denotes vegetation coverage, <em>NDVI<sub>soil</sub></em> denotes NDVI value of completely bare soil or non-vegetation covered area, <em>NDVI<sub>veg</sub></em> represents NDVI value of pixels completely covered by vegetation, namely, NDVI value of pure vegetation pixels (Li et al., 2021).</p>

			</sec> <!--sec2.3-->

			<sec id="sec2.4">
				<title>Variable importance analysis</title>

				<p><em>Driving factors</em></p>

				<p>Considering the environmental factors affecting vegetation growth, we decided to collect meteorological and topographic factors as influencing factors. Therefore, the data of four factors (daily average temperature, relative humidity, precipitation, and wind speed) were selected as meteorological factors and collected from the China Meteorological Data Network. For topographic factors, the influencing factors were divided into two categories, climate and terrain, including seven variables (Table 2). The altitude, slope, and aspect factors of the study area were extracted from digital elevation model (DEM) data, which were acquired from the geospatial data cloud.</p>

				<table-wrap id="t2">
					<label>Table 2</label>
					<caption>
						<title>Driving factors affecting vegetation change.</title>
					</caption>
					<table>
						<thead>
							<tr>
								<th align="center">Variable factors</th>
								<th align="center">Variable name</th>
								<th align="center">Code</th>
								<th align="center">Source</th>
							</tr>
						</thead>
						<tbody>
						    <tr>
						        <td rowspan="4" valign="top">Climatic</td>
						        <td>Daily average temperature</td>
						        <td align="center">Temp_avg</td>
						        <td align="center" rospan="4" valign="top">China Meteorological Data Network (<a href="http://data.cma.cn/" target="_blank">http://data.cma.cn/</a>)</td>
						    </tr>
						    <tr>
						        <td>Daily average relative humidity</td>
						        <td align="center">Rh_avg</td>
						    </tr>
						    <tr>
						        <td>Daily average precipitation</td>
						        <td align="center">Pre_avg</td>
						    </tr>
						    <tr>
						        <td>Daily average wind speed</td>
						        <td align="center">Ws_avg</td>
						    </tr>
						    <tr>
						        <td rowspan="3" valign="top">Terrain</td>
						        <td>Altitude</td>
						        <td align="center">ALT</td>
						        <td align="center" rospan="2" valign="top">Geospatial Data Cloud (<a href="www.gscloud.cn/" target="_blank">www.gscloud.cn/</a>)</td>
						    </tr>
						    <tr>
						        <td>Slope</td>
						        <td align="center">Slope</td>
						    </tr>
						    <tr>
						        <td>Aspect</td>
						        <td align="center">Aspect</td>
						    </tr>
						</tbody>
					</table>
				</table-wrap>				

				<p><em>Random forest</em></p>

				<p>The importance of the identified driving factors was ranked by random forest method to determine the key factors affecting the changes in vegetation cover in Maoershan. This method directly measures the influence of each factor on the prediction accuracy of the model. The fundamental idea is to rearrange the eigenvalues and observe the reduction in the accuracy of the model. For unimportant features, this method has little impact on the accuracy of the model, but for important features, it largely reduces the accuracy of the model (Ma et al., 2020). In this study, the “Random Forest” package in R was used for calculating and sorting the importance of each variable. The importance of a certain feature X in the random forest was calculated as follows: (1) for each decision tree in the random forest, the corresponding out-of-pocket data (OOB) was used to calculate its out-of-pocket data error, denoted as errOOB1; (2) noise interference was randomly added to feature X of all samples of the data outside the bag (the sample value at feature X can be randomly changed), and the error of its data outside the bag was calculated again, denoted as errOOB2.</p>

				<p>If there are N variables in the random forest, then the importance equation of feature X is:</p>

				<disp-formula id="e3">
				<math id="e3">
					<mtext>Importance</mtext>
				  	<mo>=</mo>
				  	<mfrac>
				  		<mn>1</mn>
				  		<mi>N</mi>
				  	</mfrac>
				    <mo>&#8721;</mo>
				    <mo>(</mo>
						<mi>errOOB</mi>
						<mn>2</mn>
						<mo>-</mo>
						<mi>errOOB</mi>
						<mn>1</mn>						
					<mo>)</mo>
				</math>
				<label>(3)</label>
				</disp-formula>

				<p>Spearman correlation analysis was also used in this study to analyze the correlation between vegetation change and several factors, test the results of random forest, and further determine the influence of different factors on vegetation change.</p>

			</sec> <!--sec2.4-->

		</sec> <!--/sec 2-->

		<sec id="sec3" sec-type="results">
			<title>Results</title>

			<sec id="3.1">
				<title>Variations in vegetation growth status</title>

				<p>From 2000 to 2020, the average NDVI of vegetation in this region was 0.84, the lowest being 0.73 in 2010, and the highest being 0.90 in 2005 (Fig. S1 [suppl]). Further analysis demonstrated an increasing trend of NDVI from 2000 to 2005, a sharp reduction in 2010, and then an accelerating trend until 2020.</p>

				<p>The spatial distribution of NDVI (Fig. 2) indicates that the high value was mainly distributed in the northern and central parts of the region. From 2000 to 2020, vegetation NDVI in most areas of Maoershan showed an increasing trend. The areas with a rapid increase rate were mainly located in the north and east of the region; the areas with a slow increase were mainly located in the middle of the region. In addition, vegetation NDVI showed a decreasing trend in the southern and western regions (Fig. 2f).</p>

				<fig id="f2">
					<label>Figure 2</label>
					<caption>
						<title>Spatial distribution of normalized difference vegetation index (NDVI) in Maoershan from 2000 to 2020. (a)-(e) represent spatial distribution of NDVI from 2000 to 2020, respectively, and (f) represent changes of NDVI in the past 20 years.</title>
					</caption>
					<graphic id="gra-2" xlink:href="img/e015-fig2.jpg"/>
				</fig>					

				<p>From 2000 to 2020, VFC in most areas of Maoershan had a decreasing trend with a smaller amplitude (Fig. 3). The spatial distribution of VFC changed greatly over time. Especially in 2010, the part and overall VFC showed a sharp decline (Fig. 3c). A rapid rate of decline was primarily observed in the central part of the region. The areas with a slower rate of decline were mainly located in the northern and eastern parts of the region. In addition, VFC showed an increasing trend in the southwest and other local areas (Fig. 3f). Figs. 2 and 3 show that the period 2005-2010 marked a change due to the loss of vegetation and the poor climatic conditions of 2010, which then drags on the NDVIs and VFCs of the following years, giving negative balances for all those areas that suffered a loss of vegetation.</p>

				<fig id="f3">
					<label>Figure 3</label>
					<caption>
						<title>Spatial distribution of vegetation fraction cover (VFC) in Maoershan from 2000 to 2020. (a)-(e) represent the spatial distribution of VFC from 2000 to 2020, respectively, and (f) represent the changes of VFC in the past 20 years.</title>
					</caption>
					<graphic id="gra-3" xlink:href="img/e015-fig3.jpg"/>
				</fig>					

			</sec> <!--/sec 3.1-->


			<sec id="3.2">
				<title>Analysis of driving factors of vegetation growth</title>

				<p>From 2000 to 2020, the average annual temperature in the Maoershan area shows an increasing trend (Fig. S2 [suppl]). The highest average temperature was 5.16℃ in 2005, the lowest average temperature was 1.99℃ in 2010, and the average temperature in the past 20 years was 3.96℃. The annual precipitation also had an increasing trend. The maximum precipitation was 1596 mm in 2020, the lowest was 454 mm in 2000, and the average precipitation was 738 mm (Fig. S2a). Only a slight variation in the annual average relative humidity was observed, with the values initially increasing followed by a decreasing, the average value being 73.3%. The annual mean wind speed was also stable, ranging from 0.02 to 0.49, and the mean value in the 20 years was 2.64 m/s (Fig. S2b).</p>
				
				<p>The topographic conditions (elevation, slope, and aspect) in this region are shown in Fig. S3 [suppl]. The high elevations are located mainly in the north and east of the region, with gradients mostly below 25°.</p>

				<p>By utilizing random forest and spearman correlation methods, the importance of climate and terrain factors on NDVI and their specific relationships with NDVI were analyzed. As shown in Fig. 4, terrain factors were found to be more important than climate factors. Among all variables, altitude was the most important driving factor for vegetation growth, followed by relative humidity, temperature, wind speed, slope, precipitation, and aspect.</p>

				<fig id="f4">
					<label>Figure 4</label>
					<caption>
						<title>Importance and correlation of the seven variables. The length of the bar indicates the importance: the longer the bar, the higher the importance. The color of the bar indicates correlation: red, positive correlation; green, negative correlation. ALT, altitude; Rh_avg, daily average relative humidity; Temp_avg, daily average temperature; Ws_avg, daily average wind speed; Pre_avg, daily average precipitation.</title>
					</caption>
					<graphic id="gra-4" xlink:href="img/e015-fig4.jpg"/>
				</fig>				

				<p>The Spearman results can accurately determine the influence of changes in each factor on NDVI values. For climate factors, all of them showed a positive correlation with NDVI, including relative humidity, temperature, wind speed, and precipitation (Fig. 4). This indicates that with the increase of relative humidity, temperature, wind speed and precipitation values, NDVI values increased. In other words, the greater the value of these four meteorological factors were related with a greater positive effect on vegetation growth. For the topographic factors, NDVI was positively correlated with slope, and negatively correlated with altitude and aspect (Fig. 4). This indicated that, in the study area, greater slope values were related with higher NDVI values; conversely, the greater the altitude and aspect values (aspect value: northwest > west > southwest > south > southeast > east > northeast > north > flat) were related with lower NDVI values.</p>

			</sec> <!--/sec 3.2-->

		</sec> <!--/sec 3-->

		<sec id="sec4" sec-type="discussion">
			<title>Discussion</title>

			<p>This study evaluated the changes in NDVI and VFC as well as their distribution in the study area from 2000 to 2020 using Landsat time series image data and analysed the spatio-temporal differences in the vegetation cover in Maoershan. The topographic and climatic factors in this region were identified using the random forest algorithm and the correlation analysis method, and the driving factors and their influences on vegetation change in this region were analysed.</p>

			<p>The variations in the NDVI time series in the study area from 2000 to 2020 indicated an initial increase in the NDVI value in Maoershan followed by a decrease, and then again, an increase. However, the entire study area demonstrated an increasing trend of NDVI. The high values were mainly distributed in the northern and central parts of the region. The areas with a faster rate of increase in NDVI were located in the north and east of the region, while the areas with a slower rate of increase were mainly located in the central part of the region. These results were similar to those from previous studies on the changes in the vegetation cover in Northeast China (Mao et al., 2012). In addition, although vegetation NDVI demonstrated an overall increasing trend, it followed a decreasing trend in the southern and western regions. This contrasting phenomenon of regional changes in vegetation may be caused due to the differences in terrain and different vegetation types.</p>

			<p>In this study, we used NDVI values acquired every five years to monitor vegetation changes in the study area. The choice of five-year intervals was based on the availability of long-term NDVI data and the need to capture significant changes in vegetation cover and productivity. While using 5-year NDVI values may not capture short-term changes in vegetation, it can provide a broader perspective on long-term trends and changes in vegetation cover. It is also a practical approach, especially in areas with limited data availability or for large-scale assessments. However, we acknowledge that using 5-year intervals may not capture all changes in vegetation, such as those due to natural disturbances or human activities. Therefore, future studies could consider using finer temporal scales or complementary data sources, such as high-resolution imagery or field observations, to improve the accuracy of vegetation monitoring.</p>

			<p>In recent years, numerous researchers (e.g., Avetisyan et al., 2022; Yunus &amp; Polat, 2023) have endeavoured to employ alternative vegetation indices in addition to NDVI for the purpose of monitoring and assessing vegetation growth. These vegetation indexes include enhanced vegetation index (EVI), soil adjusted vegetation index (SAVI), normalized difference greenness index (NDGI. They take a more comprehensive approach to considering the various factors that impact vegetation. EVI enhances the stability and sensitivity of vegetation indices by mitigating the impact of atmospheric scattering and shading. SAVI, through modifying the soil reflection, enables a more comprehensive consideration of the interaction between bare soil and vegetation. By contrast, the NDGI not only captures the presence of sparse vegetation well, but also compensates for the effect of dry grass on green-up date extraction. Although these indexes have advantages over NDVI, they also have some limitations. Some studies have suggested that these indices have better applicability for certain types of vegetation (Dai et al., 2023). In addition, differences between these indices can also lead to some difficult to interpret results when comparing vegetation growth in the same area. Therefore, despite its limitations, NDVI remains one of the most widely used vegetation indices due to its high repeatability and broad applicability, providing valuable insights into the study of vegetation growth. Meanwhile, researchers should also prioritize the advancement of novel vegetation indices and select the most appropriate index based on specific analytical requirements.</p>

			<p>When using NDVI for vegetation studies, it provides information about the presence of vegetation, but not about the structure of vegetation. This is because NDVI only reflects the chlorophyll content of vegetation, but it cannot distinguish different types of vegetation or different levels of vegetation. Therefore, in order to obtain more comprehensive and accurate structural measurement information, many studies use LiDAR technology, which uses laser scanning, to obtain 3D vegetation structure information (Tang et al., 2022). Compared with NDVI, which can only provide vegetation surface information, LiDAR can obtain vegetation height, coverage, volume, density and other multi-dimensional information. This information cannot only help determine if vegetation is changing, but also provide researchers with more detailed information about the composition, distribution, and changes in vegetation. However, there are some limitations to using the LiDAR technique. First, it requires a lot of data processing and analysis, which requires a high level of technology and computing power. Secondly, the application scope of LiDAR technology is relatively small, mainly limited to some relatively flat terrain and easily accessible vegetation (Dhargay et al., 2022). Overall, combining LiDAR technology can provide more vegetation information, which will help to assess vegetation change more accurately. However, due to the limitation of data collection and study area, the selection of LiDAR technology should be comprehensively considered. Of course, other data and indices can be used to obtain more vegetation information to obtain a more comprehensive assessment of vegetation change.</p>

			<p>The driving factors and their influences on vegetation changes in Maoershan were analyzed using a random forest algorithm and correlation analysis. The topographic and climatic factors were the major factors affecting vegetation growth in Maoershan.</p>

			<p>In this study, altitude, slope and aspect were selected as three terrain factors affecting vegetation growth. Compared with climate factors, terrain factors had a greater impact on vegetation growth. Among all the variables, altitude was the most important factor affecting vegetation growth, and slope was more important than aspect. Among them, NDVI was negatively correlated with altitude, indicating that the higher altitude values were related with lower NDVI values, which meant that the negative impact on vegetation growth was higher. At very high altitudes, vegetation growth can be inhibited to a certain extent, particularly due to low-pressure and low-oxygen environments in high-altitude areas, and insufficient oxygen can retard the growth of vegetation (Lou et al., 2021). The aspect was less important to vegetation change than elevation. This may be because the aspect variations in the study area led to little differences in the solar radiation received, resulting in a minimal impact of the aspect on vegetation growth (Xiong &amp; Wang, 2022). Previous studies in northern China also reached similar conclusions (Sun Y et al., 2015).</p>

			<p>Among the climate factors, the daily average relative humidity was found to be the most important factor affecting the growth of forest vegetation in the study area followed by temperature, and wind speed was a more important factor affecting vegetation than precipitation. Relative humidity can affect the water content of forest vegetation. As high relative humidity increases the water absorbed by vegetation, the growth of vegetation is accelerated (Guang-Lei et al., 2012). Sufficient precipitation can improve soil conditions and increase the water content of vegetation, thereby promoting the growth of vegetation (Wang et al., 2003). The results indicated that temperature had a positive effect on NDVI, indicating that an increase in temperature can promote the growth of vegetation. The major reason for this effect is the direct effect of temperature on the photosynthesis of vegetation. As the increase in temperature accelerates the photosynthesis of vegetation, it contributes to the growth of forest vegetation (Chuai et al., 2013). The variational trend of NDVI from 2000 to 2020 indicated an increase from 2000 to 2005, a sharp decrease in 2010, and then an accelerating trend until 2020 (Fig. S1). NDVI decreases may be due to loss of vegetation or lower local photosynthetic activity. The period 2005-2010 marked a change due to the loss of vegetation and the poor climatic conditions of 2010 (Figs. 2 and 3), which then drags on the NDVIs and VFCs of the following years, giving negative balances for all those areas that have suffered a loss of vegetation. It is worth noting that studies have demonstrated a threshold for the influence of temperature on vegetation growth, and negative effects may occur when the positive effects of temperature on vegetation reach a threshold (Julien et al., 2006). Extremely high temperatures can accelerate the evaporation of vegetation and reduce the water content, which is not conducive to the growth of vegetation. Interestingly, the region demonstrated a sharp increase in relative humidity in 2010, despite a rapid decline in NDVI. The possible reason for this phenomenon is that a rapid decline in temperature might have affected the photosynthesis of vegetation. At this time, the increase in relative humidity factor did not completely improve the growth of vegetation, thereby leading to a sudden decline in NDVI in this region in 2010.</p>

			<p>In this study, Landsat images were selected to calculate the time series of NDVI. Compared with previous time series analysis of large-scale which mainly used medium and low-resolution remote sensing images (MODIS, GIMMS products, etc.) to study vegetation change, there was a large improvement in the spatial resolution of the images in this study, with these images more accurately reflecting the changes in the regional vegetation. Moreover, the image selected in this study was easy to obtain and unrestricted, and the time series was long, which increases its use in research and practical value (Ke et al., 2015).</p>

			<p>Although topographic and climatic factors are the major factors affecting vegetation growth, and they were analysed in this study, the influence of human activities (including forest management, experimental logging, etc.) should not be ignored, and this is a limitation of this paper. The influence of human activities on vegetation change has both positive and negative aspects. With the increase in population, there is an increase in the negative impact of human activities on the ecological environment, such as coal mining, overgrazing, urban expansion, etc. However, humans can improve the regional ecological environment quality by reducing environmental pollution, returning farmland to forest, improving energy structure, and implementing greening projects. As this research area is located in an experimental forest farm, there is little impact from sudden logging and other management activities. Therefore, the factor ‘forest management activities’ was not considered a driving factor in this study. However, subsequent studies should take this factor into account. The results of changes in the vegetation cover highlighted in this study are important for the development of a sustainable forest management strategy considering the requirements of all the stakeholders. The results of this study can support long-term monitoring and management of vegetation in the future and contribute to sustainable environmental development.</p>

			<p>In summary, by using remote sensing and GIS tools, we analysed the spatiotemporal characteristics and driving factors of vegetation cover on Maoershan from 2000 to 2020. According to our findings, altitude had the greatest impact on vegetation growth, and it was negatively correlated with NDVI. These results suggest that vegetation cover in this region is heavily influenced by environmental factors, which can provide scientific reference for the analysis of vegetation status and evaluation of regional ecological environment. Exploring the impact of driving factors on vegetation, further revealing the reasons for vegetation cover changes, and constructing long-term and high-spatial-temporal resolution vegetation cover change data are the focus of the next step of research.</p>

						
		</sec> <!--/sec 4-->
		
	</body>
	<back>
		
		<author-notes>
			<title>Authors’ contributions</title>
			<fn>Conceptualization: T. Li.</fn>
			<fn>Data curation: T. Li.</fn>
			<fn>Formal analysis: Y. Gao.</fn>
			<fn>Funding acquisition: T. Li.</fn>
			<fn>Investigation: Y. Gao.</fn>
			<fn>Methodology: T. Li.</fn>
			<fn>Project administration: T. Li.</fn>
			<fn>Resources: T. Li.</fn>
			<fn>Software: Y. Gao.</fn>
			<fn>Supervision: Y. Gao.</fn>
			<fn>Validation: T. Li.</fn>
			<fn>Visualization: Y. Gao.</fn>
			<fn>Writing – original draft: T. Li.</fn>
			<fn>Writing – review &amp; editing: Y. Gao.</fn>
		</author-notes>

		<ref-list>
			<title>References</title>
		
			<ref id="R1">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Amiri</surname>
			                <given-names>R</given-names>
			            </string-name>
			            <string-name>
			                <surname>Weng</surname>
			                <given-names>Q</given-names>
			            </string-name>
			            <string-name>
			                <surname>Alimohammadi</surname>
			                <given-names>A</given-names>
			            </string-name>
			            <string-name>
			                <surname>Alavipanah</surname>
			                <given-names>SK</given-names>
			            </string-name>
			        </person-group>
			        <year>2009</year>
			        <article-title>Spatial-temporal dynamics of land surface temperature in relation to fractional vegetation cover and land use/cover in the Tabriz urban area, Iran.</article-title>
			        <source>Remote Sens Environ</source>
			        <issue>113</issue>
			        <fpage>2606</fpage>
			        <lpage>2617</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.rse.2009.07.021"/>
			    </mixed-citation>
			</ref>

			<ref id="R2">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Avetisyan</surname>
			                <given-names>D</given-names>
			            </string-name>
			            <string-name>
			                <surname>Velizarova</surname>
			                <given-names>E</given-names>
			            </string-name>
			            <string-name>
			                <surname>Filchev</surname>
			                <given-names>L</given-names>
			            </string-name>
			        </person-group>
			        <year>2022</year>
			        <article-title>Post-fire forest vegetation state monitoring through satellite remote sensing and in situ data.</article-title>
			        <source>Remote Sens</source>
			        <issue>14</issue>
			        <fpage>6266</fpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/rs14246266"/>
			    </mixed-citation>
			</ref>

			<ref id="R3">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Campbell</surname>
			                <given-names>DJ</given-names>
			            </string-name>
			            <string-name>
			                <surname>Lusch</surname>
			                <given-names>DP</given-names>
			            </string-name>
			            <string-name>
			                <surname>Smucker</surname>
			                <given-names>TA</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wangui</surname>
			                <given-names>EE</given-names>
			            </string-name>
			        </person-group>
			        <year>2005</year>
			        <article-title>Multiple methods in the study of driving forces of land use and land cover change: a case study of SE Kajiado District, Kenya.</article-title>
			        <source>Human Ecol</source>
			        <issue>33</issue>
			        <fpage>763</fpage>
			        <lpage>794</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1007/s10745-005-8210-y"/>
			    </mixed-citation>
			</ref>

			<ref id="R4">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Chang</surname>
			                <given-names>L</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhang</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Li</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Mao</surname>
			                <given-names>X</given-names>
			            </string-name>
			        </person-group>
			        <year>2020</year>
			        <article-title>Multilevel extraction of vegetation type based on airborne LiDAR data.</article-title>
			        <source>Can J Remote Sens</source>
			        <issue>46</issue>
			        <fpage>681</fpage>
			        <lpage>694</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1080/07038992.2020.1850248"/>
			    </mixed-citation>
			</ref>

			<ref id="R5">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Chuai</surname>
			                <given-names>X</given-names>
			            </string-name>
			            <string-name>
			                <surname>Huang</surname>
			                <given-names>X</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>W</given-names>
			            </string-name>
			            <string-name>
			                <surname>Bao</surname>
			                <given-names>G</given-names>
			            </string-name>
			        </person-group>
			        <year>2013</year>
			        <article-title>NDVI, temperature and precipitation changes and their relationships with different vegetation types during 1998-2007 in Inner Mongolia, China.</article-title>
			        <source>Int J Climatol</source>
			        <issue>33</issue>
			        <fpage>1696</fpage>
			        <lpage>1706</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1002/joc.3543"/>
			    </mixed-citation>
			</ref>

			<ref id="R6">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Dai</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Luo</surname>
			                <given-names>H</given-names>
			            </string-name>
			            <string-name>
			                <surname>Hu</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zheng</surname>
			                <given-names>Q</given-names>
			            </string-name>
			            <string-name>
			                <surname>Li</surname>
			                <given-names>H</given-names>
			            </string-name>
			            <string-name>
			                <surname>Li</surname>
			                <given-names>M</given-names>
			            </string-name>
			        </person-group>
			        <year>2023</year>
			        <article-title>Retrieving leaf area index of rubber plantation in Hainan Island using empirical and neural network models with Landsat images.</article-title>
			        <source>J Appl Remote Sens</source>
			        <issue>17</issue>
			        <elocation-id>014503</elocation-id>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1117/1.JRS.17.014503"/>
			    </mixed-citation>
			</ref>

			<ref id="R7">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Dhargay</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Lyell</surname>
			                <given-names>CS</given-names>
			            </string-name>
			            <string-name>
			                <surname>Brown</surname>
			                <given-names>TP</given-names>
			            </string-name>
			            <string-name>
			                <surname>Inbar</surname>
			                <given-names>A</given-names>
			            </string-name>
			            <string-name>
			                <surname>Sheridan</surname>
			                <given-names>GJ</given-names>
			            </string-name>
			            <string-name>
			                <surname>Lane</surname>
			                <given-names>PN</given-names>
			            </string-name>
			        </person-group>
			        <year>2022</year>
			        <article-title>Performance of GEDI space-borne LiDAR for quantifying structural variation in the temperate forests of South-Eastern Australia.</article-title>
			        <source>Remote Sens</source>
			        <issue>14</issue>
			        <fpage>3615</fpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/rs14153615"/>
			    </mixed-citation>
			</ref>

			<ref id="R8">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Gao</surname>
			                <given-names>L</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>X</given-names>
			            </string-name>
			            <string-name>
			                <surname>Johnson</surname>
			                <given-names>BA</given-names>
			            </string-name>
			            <string-name>
			                <surname>Tian</surname>
			                <given-names>Q</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Verrelst</surname>
			                <given-names>J</given-names>
			            </string-name>
			        </person-group>
			        <year>2020</year>
			        <article-title>Remote sensing algorithms for estimation of fractional vegetation cover using pure vegetation index values: A review.</article-title>
			        <source>ISPRS J Photogram Remote Sens</source>
			        <issue>159</issue>
			        <fpage>364</fpage>
			        <lpage>377</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.isprsjprs.2019.11.018"/>
			    </mixed-citation>
			</ref>

			<ref id="R9">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Ge</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Meng</surname>
			                <given-names>B</given-names>
			            </string-name>
			            <string-name>
			                <surname>Liang</surname>
			                <given-names>T</given-names>
			            </string-name>
			            <string-name>
			                <surname>Feng</surname>
			                <given-names>Q</given-names>
			            </string-name>
			            <string-name>
			                <surname>Gao</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Yang</surname>
			                <given-names>S</given-names>
			            </string-name>
			        </person-group>
			        <year>2018</year>
			        <article-title>Modeling alpine grassland cover based on MODIS data and support vector machine regression in the headwater region of the Huanghe River, China.</article-title>
			        <source>Remote Sens Environ</source>
			        <issue>218</issue>
			        <fpage>162</fpage>
			        <lpage>173</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.rse.2018.09.019"/>
			    </mixed-citation>
			</ref>

			<ref id="R10">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Guang-Lei</surname>
			                <given-names>H</given-names>
			            </string-name>
			            <string-name>
			                <surname>De-Ying</surname>
			                <given-names>L</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zheng-Xiang</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Dan</surname>
			                <given-names>G</given-names>
			            </string-name>
			            <string-name>
			                <surname>Hong-Yan</surname>
			                <given-names>Z</given-names>
			            </string-name>
			        </person-group>
			        <year>2012</year>
			        <article-title>Response of NDVI in different climatic zone to climate change in Songnen plain.</article-title>
			        <source>Chin J Agrometeorol</source>
			        <issue>33</issue>
			        <fpage>271</fpage>
			    </mixed-citation>
			</ref>

			<ref id="R11">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Huang</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Tang</surname>
			                <given-names>L</given-names>
			            </string-name>
			            <string-name>
			                <surname>Hupy</surname>
			                <given-names>JP</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Shao</surname>
			                <given-names>G</given-names>
			            </string-name>
			        </person-group>
			        <year>2021</year>
			        <article-title>A commentary review on the use of normalized difference vegetation index (NDVI) in the era of popular remote sensing.</article-title>
			        <source>J Forest Res</source>
			        <issue>32</issue>
			        <fpage>1</fpage>
			        <lpage>6</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1007/s11676-020-01155-1"/>
			    </mixed-citation>
			</ref>

			<ref id="R12">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Jia</surname>
			                <given-names>L</given-names>
			            </string-name>
			            <string-name>
			                <surname>Li</surname>
			                <given-names>ZB</given-names>
			            </string-name>
			            <string-name>
			                <surname>Xu</surname>
			                <given-names>GC</given-names>
			            </string-name>
			            <string-name>
			                <surname>Ren</surname>
			                <given-names>ZP</given-names>
			            </string-name>
			            <string-name>
			                <surname>Li</surname>
			                <given-names>P</given-names>
			            </string-name>
			            <string-name>
			                <surname>Cheng</surname>
			                <given-names>YT</given-names>
			            </string-name>
			        </person-group>
			        <year>2020</year>
			        <article-title>Dynamic change of vegetation and its response to climate and topographic factors in the Xijiang River basin, China.</article-title>
			        <source>Environ Sci Pollut Res</source>
			        <issue>27</issue>
			        <fpage>11637</fpage>
			        <lpage>11648</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1007/s11356-020-07692-w"/>
			    </mixed-citation>
			</ref>

			<ref id="R13">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Jiang</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Huete</surname>
			                <given-names>AR</given-names>
			            </string-name>
			            <string-name>
			                <surname>Chen</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Chen</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Li</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Yan</surname>
			                <given-names>G</given-names>
			            </string-name>
			        </person-group>
			        <year>2006</year>
			        <article-title>Analysis of NDVI and scaled difference vegetation index retrievals of vegetation fraction.</article-title>
			        <source>Remote Sens Environ</source>
			        <issue>101</issue>
			        <fpage>366</fpage>
			        <lpage>378</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.rse.2006.01.003"/>
			    </mixed-citation>
			</ref>

			<ref id="R14">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Julien</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Sobrino</surname>
			                <given-names>JA</given-names>
			            </string-name>
			            <string-name>
			                <surname>Verhoef</surname>
			                <given-names>W</given-names>
			            </string-name>
			        </person-group>
			        <year>2006</year>
			        <article-title>Changes in land surface temperatures and NDVI values over Europe between 1982 and 1999.</article-title>
			        <source>Remote Sens Environ</source>
			        <issue>103</issue>
			        <fpage>43</fpage>
			        <lpage>55</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.rse.2006.03.011"/>
			    </mixed-citation>
			</ref>

			<ref id="R15">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Ke</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Im</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Lee</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Gong</surname>
			                <given-names>H</given-names>
			            </string-name>
			            <string-name>
			                <surname>Ryu</surname>
			                <given-names>Y</given-names>
			            </string-name>
			        </person-group>
			        <year>2015</year>
			        <article-title>Characteristics of Landsat 8 OLI-derived NDVI by comparison with multiple satellite sensors and in-situ observations.</article-title>
			        <source>Remote Sens Environ</source>
			        <issue>164</issue>
			        <fpage>298</fpage>
			        <lpage>313</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.rse.2015.04.004"/>
			    </mixed-citation>
			</ref>

			<ref id="R16">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Klisch</surname>
			                <given-names>A</given-names>
			            </string-name>
			            <string-name>
			                <surname>Atzberger</surname>
			                <given-names>C</given-names>
			            </string-name>
			        </person-group>
			        <year>2016</year>
			        <article-title>Operational drought monitoring in Kenya using MODIS NDVI time series.</article-title>
			        <source>Remote Sens</source>
			        <issue>8</issue>
			        <fpage>267</fpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/rs8040267"/>
			    </mixed-citation>
			</ref>

			<ref id="R17">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Li</surname>
			                <given-names>F</given-names>
			            </string-name>
			            <string-name>
			                <surname>Chen</surname>
			                <given-names>W</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zeng</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhao</surname>
			                <given-names>Q</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wu</surname>
			                <given-names>B</given-names>
			            </string-name>
			        </person-group>
			        <year>2014</year>
			        <article-title>Improving estimates of grassland fractional vegetation cover based on a pixel dichotomy model: A case study in Inner Mongolia, China.</article-title>
			        <source>Remote Sens</source>
			        <issue>6</issue>
			        <fpage>4705</fpage>
			        <lpage>4722</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/rs6064705"/>
			    </mixed-citation>
			</ref>

			<ref id="R18">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Li</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhang</surname>
			                <given-names>M</given-names>
			            </string-name>
			            <string-name>
			                <surname>Tang</surname>
			                <given-names>Q</given-names>
			            </string-name>
			        </person-group>
			        <year>2021</year>
			        <article-title>Characterizing and attributing the vegetation coverage changes in North Shanxi coal base of China from 1987 to 2020.</article-title>
			        <source>Resour Policy</source>
			        <issue>74</issue>
			        <elocation-id>102331</elocation-id>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.resourpol.2021.102331"/>
			    </mixed-citation>
			</ref>

			<ref id="R19">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Liu</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Huang</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Xie</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>H</given-names>
			            </string-name>
			            <string-name>
			                <surname>Huang</surname>
			                <given-names>Q</given-names>
			            </string-name>
			            <string-name>
			                <surname>Leng</surname>
			                <given-names>G</given-names>
			            </string-name>
			        </person-group>
			        <year>2019</year>
			        <article-title>Spatial-temporal changes in vegetation cover in a typical semi-humid and semi-arid region in China: Changing patterns, causes and implications.</article-title>
			        <source>Ecol Ind</source>
			        <issue>98</issue>
			        <fpage>462</fpage>
			        <lpage>475</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.ecolind.2018.11.037"/>
			    </mixed-citation>
			</ref>

			<ref id="R20">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Lou</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Xu</surname>
			                <given-names>G</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Yang</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Ni</surname>
			                <given-names>S</given-names>
			            </string-name>
			        </person-group>
			        <year>2021</year>
			        <article-title>Multi-year NDVI values as indicator of the relationship between spatiotemporal vegetation dynamics and environmental factors in the Qaidam Basin, China.</article-title>
			        <source>Remote Sens</source>
			        <issue>13</issue>
			        <fpage>1240</fpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/rs13071240"/>
			    </mixed-citation>
			</ref>

			<ref id="R21">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Ma</surname>
			                <given-names>W</given-names>
			            </string-name>
			            <string-name>
			                <surname>Feng</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Cheng</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Chen</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>F</given-names>
			            </string-name>
			        </person-group>
			        <year>2020</year>
			        <article-title>Identifying forest fire driving factors and related impacts in China using random forest algorithm.</article-title>
			        <source>Forests</source>
			        <issue>11</issue>
			        <fpage>507</fpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/f11050507"/>
			    </mixed-citation>
			</ref>

			<ref id="R22">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Mao</surname>
			                <given-names>D</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Luo</surname>
			                <given-names>L</given-names>
			            </string-name>
			            <string-name>
			                <surname>Ren</surname>
			                <given-names>C</given-names>
			            </string-name>
			        </person-group>
			        <year>2012</year>
			        <article-title>Integrating AVHRR and MODIS data to monitor NDVI changes and their relationships with climatic parameters in Northeast China.</article-title>
			        <source>Int J Appl Earth Observ Geoinform</source>
			        <issue>18</issue>
			        <fpage>528</fpage>
			        <lpage>536</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.jag.2011.10.007"/>
			    </mixed-citation>
			</ref>

			<ref id="R23">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Martínez</surname>
			                <given-names>B</given-names>
			            </string-name>
			            <string-name>
			                <surname>Gilabert</surname>
			                <given-names>MA</given-names>
			            </string-name>
			        </person-group>
			        <year>2009</year>
			        <article-title>Vegetation dynamics from NDVI time series analysis using the wavelet transform.</article-title>
			        <source>Remote Sens Environ</source>
			        <issue>113</issue>
			        <fpage>1823</fpage>
			        <lpage>1842</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.rse.2009.04.016"/>
			    </mixed-citation>
			</ref>

			<ref id="R24">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Seidl</surname>
			                <given-names>R</given-names>
			            </string-name>
			            <string-name>
			                <surname>Thom</surname>
			                <given-names>D</given-names>
			            </string-name>
			            <string-name>
			                <surname>Kautz</surname>
			                <given-names>M</given-names>
			            </string-name>
			            <string-name>
			                <surname>Martin-Benito</surname>
			                <given-names>D</given-names>
			            </string-name>
			            <string-name>
			                <surname>Peltoniemi</surname>
			                <given-names>M</given-names>
			            </string-name>
			            <string-name>
			                <surname>Vacchiano</surname>
			                <given-names>G</given-names>
			            </string-name>
			        </person-group>
			        <year>2017</year>
			        <article-title>Forest disturbances under climate change.</article-title>
			        <source>Nature Climate Change</source>
			        <issue>7</issue>
			        <fpage>395</fpage>
			        <lpage>402</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1038/nclimate3303"/>
			    </mixed-citation>
			</ref>

			<ref id="R25">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Sun</surname>
			                <given-names>L</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhao</surname>
			                <given-names>D</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhang</surname>
			                <given-names>G</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wu</surname>
			                <given-names>X</given-names>
			            </string-name>
			            <string-name>
			                <surname>Yang</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>Z</given-names>
			            </string-name>
			        </person-group>
			        <year>2022</year>
			        <article-title>Using SPOT VEGETATION for analyzing dynamic changes and influencing factors on vegetation restoration in the Three-River Headwaters Region in the last 20 years (2000-2019), China.</article-title>
			        <source>Ecology Engineering</source>
			        <issue>183</issue>
			        <elocation-id>106742</elocation-id>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.ecoleng.2022.106742"/>
			    </mixed-citation>
			</ref>

			<ref id="R26">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Sun</surname>
			                <given-names>W</given-names>
			            </string-name>
			            <string-name>
			                <surname>Song</surname>
			                <given-names>X</given-names>
			            </string-name>
			            <string-name>
			                <surname>Mu</surname>
			                <given-names>X</given-names>
			            </string-name>
			            <string-name>
			                <surname>Gao</surname>
			                <given-names>P</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>F</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhao</surname>
			                <given-names>G</given-names>
			            </string-name>
			        </person-group>
			        <year>2015</year>
			        <article-title>Spatiotemporal vegetation cover variations associated with climate change and ecological restoration in the Loess Plateau.</article-title>
			        <source>Agricultural and Forest Meteorology</source>
			        <issue>209</issue>
			        <fpage>87</fpage>
			        <lpage>99</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.agrformet.2015.05.002"/>
			    </mixed-citation>
			</ref>

			<ref id="R27">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Sun</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Yang</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhang</surname>
			                <given-names>L</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>Z</given-names>
			            </string-name>
			        </person-group>
			        <year>2015</year>
			        <article-title>The relative roles of climate variations and human activities in vegetation change in North China.</article-title>
			        <source>Physics and Chemistry of the Earth</source>
			        <issue>87</issue>
			        <fpage>67</fpage>
			        <lpage>78</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.pce.2015.09.017"/>
			    </mixed-citation>
			</ref>

			<ref id="R28">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Talukdar</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Singha</surname>
			                <given-names>P</given-names>
			            </string-name>
			            <string-name>
			                <surname>Mahato</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Pal</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Liou</surname>
			                <given-names>YA</given-names>
			            </string-name>
			            <string-name>
			                <surname>Rahman</surname>
			                <given-names>A</given-names>
			            </string-name>
			        </person-group>
			        <year>2020</year>
			        <article-title>Land-use land-cover classification by machine learning classifiers for satellite observations-A review.</article-title>
			        <source>Remote Sensing</source>
			        <issue>12</issue>
			        <fpage>1135</fpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/rs12071135"/>
			    </mixed-citation>
			</ref>

			<ref id="R29">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Tang</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Liang</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Yang</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhang</surname>
			                <given-names>S</given-names>
			            </string-name>
			            <string-name>
			                <surname>Hou</surname>
			                <given-names>H</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhu</surname>
			                <given-names>X</given-names>
			            </string-name>
			        </person-group>
			        <year>2022</year>
			        <article-title>Revealing the structure and composition of the restored vegetation cover in semi-arid mine dumps based on LiDAR and hyperspectral images.</article-title>
			        <source>Remote Sensing</source>
			        <issue>14</issue>
			        <fpage>978</fpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/rs14040978"/>
			    </mixed-citation>
			</ref>

			<ref id="R30">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Tian</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wu</surname>
			                <given-names>Z</given-names>
			            </string-name>
			            <string-name>
			                <surname>Li</surname>
			                <given-names>M</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>B</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhang</surname>
			                <given-names>X</given-names>
			            </string-name>
			        </person-group>
			        <year>2022</year>
			        <article-title>Forest fire spread monitoring and vegetation dynamics detection based on multi-source remote sensing images.</article-title>
			        <source>Remote Sensing</source>
			        <issue>14</issue>
			        <fpage>4431</fpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.3390/rs14184431"/>
			    </mixed-citation>
			</ref>

			<ref id="R31">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Rich</surname>
			                <given-names>PM</given-names>
			            </string-name>
			            <string-name>
			                <surname>Price</surname>
			                <given-names>KP</given-names>
			            </string-name>
			        </person-group>
			        <year>2003</year>
			        <article-title>Temporal responses of NDVI to precipitation and temperature in the central Great Plains, USA.</article-title>
			        <source>International Journal of Remote Sensing</source>
			        <issue>24</issue>
			        <fpage>2345</fpage>
			        <lpage>2364</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1080/01431160210154812"/>
			    </mixed-citation>
			</ref>

			<ref id="R32">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>J</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>K</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhang</surname>
			                <given-names>M</given-names>
			            </string-name>
			            <string-name>
			                <surname>Zhang</surname>
			                <given-names>C</given-names>
			            </string-name>
			        </person-group>
			        <year>2015</year>
			        <article-title>Impacts of climate change and human activities on vegetation cover in hilly southern China.</article-title>
			        <source>Ecological Engineering</source>
			        <issue>81</issue>
			        <fpage>451</fpage>
			        <lpage>461</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.ecoleng.2015.04.022"/>
			    </mixed-citation>
			</ref>

			<ref id="R33">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Xiong</surname>
			                <given-names>Y</given-names>
			            </string-name>
			            <string-name>
			                <surname>Wang</surname>
			                <given-names>H</given-names>
			            </string-name>
			        </person-group>
			        <year>2022</year>
			        <article-title>Spatial relationships between NDVI and topographic factors at multiple scales in a watershed of the Minjiang River, China.</article-title>
			        <source>Ecological Informatics</source>
			        <issue>69</issue>
			        <elocation-id>101617</elocation-id>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.ecoinf.2022.101617"/>
			    </mixed-citation>
			</ref>

			<ref id="R34">
			    <mixed-citation publication-type="journal">
			        <person-group person-group-type="author">
			            <string-name>
			                <surname>Yunus</surname>
			                <given-names>K</given-names>
			            </string-name>
			            <string-name>
			                <surname>Polat</surname>
			                <given-names>N</given-names>
			            </string-name>
			        </person-group>
			        <year>2023</year>
			        <article-title>A linear approach for wheat yield prediction by using different spectral vegetation indices.</article-title>
			        <source>International Journal of Engineering and Geosciences</source>
			        <issue>8</issue>
			        <fpage>52</fpage>
			        <lpage>62</lpage>
			        <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.26833/ijeg.1035037"/>
			    </mixed-citation>
			</ref>
		</ref-list>	
		</back>
	</article>