
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39237653
71689
10.1038/s41598-024-71689-5
Article
Remote sensing analysis of spatiotemporal impacts of anthropogenic influence on mountain landscape ecology in Pir Chinasi national park
Farooq Muhammad Akhlaq 1
Ghufran Muhammad Asad 1
Ahmed Naeem mnaeem@numl.edu.pk

2
Attia Kotb A. 3
Mohammed Arif Ahmed 3
Hafeez Yaser M. 4
Amanat Aamir 2
Farooq Muhammad Shahbaz 5
Uzair Muhammad uzairbreeder@gmail.com

6
Naz Saima 7
1 https://ror.org/047w75g40 grid.411727.6 0000 0001 2201 6036 Department of Environmental Science, International Islamic University, Islamabad, Pakistan
2 https://ror.org/008dh2426 grid.444798.2 0000 0004 0607 5732 Department of Economics, National University of Modern Languages, Islamabad, Pakistan
3 https://ror.org/02f81g417 grid.56302.32 0000 0004 1773 5396 Department of Biochemistry, College of Science, King Saud University, P.O. Box 2455, 11451 Riyadh, Saudi Arabia
4 grid.411978.2 0000 0004 0578 3577 Kafr-El-Sheikh University, Kafr el-Sheikh, 33516 Egypt
5 grid.419165.e 0000 0001 0775 7565 Rice Research Program, Crop Sciences Institute (CSI), National Agricultural Research Centre (NARC), Park Road 44000, Islamabad, Pakistan
6 grid.419165.e 0000 0001 0775 7565 National Institute for Genomics and Advanced Biotechnology (NIGAB), National Agricultural Research Centre (NARC), Park Road 45500, Islamabad, Pakistan
7 grid.418920.6 0000 0004 0607 0704 COMSATS University, Islamabad, Pakistan
5 9 2024
5 9 2024
2024
14 206955 12 2023
30 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Mountain landscapes can be fragmented due to various human activities such as tourism, road construction, urbanization, and agriculture. It can also be due to natural factors such as flash floods, glacial lake outbursts, land sliding, and climate change such as rising temperatures, heavy rains, or drought.The study’s objective was to analyze the mountain landscape ecology of Pir Chinasi National Park under anthropogenic influence and investigate the impact of anthropogenic activities on the vegetation. This study observed spatiotemporal changes in vegetation due to human activities and associated climate change for the past 25 years (1995–2020) around Pir Chinasi National Park, Muzaffrabad, Pakistan. A structured questionnaire was distributed to 200 residents to evaluate their perceptions of land use and its effects on local vegetation. The findings reveal that 60% of respondents perceived spatiotemporal pressure on the park. On the other hand, the Landsat-oriented Normalized Difference Vegetation Index (NDVI) was utilized for the less than 10% cloud-covered images of Landsat 5, 7, and 8 to investigate the vegetation degradation trends of the study area. During the entire study period, the mean maximum NDVI was approximately 0.28 in 1995, whereas the mean minimum NDVI was − 2.8 in 2010. QGIS 3.8.2 was used for the data presentation. The impact of temperature on vegetation was also investigated for the study period and increasing temperature trends were observed. The study found that 10.81% (1469.08 km2) of the area experienced substantial deterioration, while 23.57% (3202.39 km2) experienced minor degradation. The total area of degraded lands was 34.38% (or 4671.47 km2). A marginal improvement in plant cover was observed in 24.88% of the regions, while 9.69% of the regions experienced a major improvement. According to the NDVI-Rainfall relationships, the area was found to be significantly impacted by human pressures and activities (r ≤ 0.50) driving vegetation changes covering 24.67% of the total area (3352.03 km2). The area under the influence of climatic variability and change (r ≥ 0.50 ≥ 0.90) accounted for 55.84% (7587.26 km2), and the area under both climatic and human stressors (r ≥ 0.50 < 0.70) was 64%. Sustainable land management practices of conservation tillage, integrated pest management, and agroforestry help preserve soil health, water quality, and biodiversity while reducing erosion, pollution, and the degradation of natural resources. landscape restoration projects of reforestation, wetland restoration, soil erosion control, and the removal of invasive species are essential to achieve land degradation neutrality at the watershed scale.

Keywords

Landsat NDVI
Spatiotemporal analysis
Land surface temperature
Anthropogenic
Overgrazing
Pir-Chinasi national park
Subject terms

Climate sciences
Ecology
Environmental sciences
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Vegetation is an essential component of the biological cycle. Vegetation trends highlight changes in ecological systems and help determine the best strategies for mitigating climate change. Highlighting the changes in vegetation due to anthropogenic activities as well as climate change is essential for improvement. An important topic in the study of the local ecological environment is the components of the ecosystem represented by vegetation under the effect of anthropogenic activities1. The dynamics of regional vegetation have a substantial impact on ecological security2, ecosystem services and, are frequently used as important indicators of ecological variations in the environment3. As human activity has risen in recent decades, alterations in vegetation has deeply captured the trails of human activity, which have been made worse by climate change4,5. Anthropogenic stresses are thought to have significant effects on vegetation and ecosystem services as we go through the Anthropocene6,7. The global ecosystem is now changing as a result of climate change and land cover variations, which are identified as key influences on the dynamics of vegetation under universal change. Numerous anthropogenic stresses along with climate change are causing vegetation to alter and degrade in highland locations, which further compromises the ecosystem services provided by mountains and the way of life of a small number of mountain people. In the realm of research on global change, understanding the linkages between natural vegetation and cultivated vegetation has been a crucial problem that is receiving more and more attention from the scientific community8. To investigate the change in vegetation, Normalized Difference Vegetation Index (NDVI) is mostly recommended1,9,10. Utilizing the NDVI, a surrogate technique for the greenness of landscape and biological dynamics of change has made satellite Remote Sensing (RS) of vegetation simple1,9,10. The investigation of vegetation status under changing climatic patterns and the tracking of the ecological environment’s quality relies heavily on the monitoring and assessment of NDVI changes in vegetation11. Limited research examined how vegetation dynamics relate to both human activity and climate change1,12–21. Moreover, residents of these mountains can better highlight the changes in vegetation that occurred due to human intervention as well as climate change22–24. The objective of this study was to understand the mechanisms behind grassland degradation of Pir Chinasi National Park and managing damaged grasslands, which needs quantitative assessments of the relative effects of anthropogenic activity and climate change on grasslands25. As for the significance of the mountain, the UN 2030 Sustainable Development Agenda has included targets for protecting and sustainably developing mountain ecosystems26. The cover index (SDG 15.4.2) of the mountain forests of Pir Chinasi National Park has been utilized as a first approximation for measuring the sustainability and protection of mountain habitats as well as the progress made towards the aforementioned aim.

Study area

As the heights vary from roughly 700 m Above Sea Level (ASL), the Pir Chinasi National Park in the Muzaffarabad (Fig. 1) neighborhood exhibits a topographic variety27. The study area’s height is 2900 m (9,500 ft)28. The area has natural beauty with its dense trees, swiftly flowing rivers, and twisting streams. The mean extreme temperature for January and June respectively, ranged from − 2.6 °C to 45.2 °C. The average rainfall is found to be between 1000 and 1300 mm or around 680 mm and falls over the four months from May to August. During the day, the wind blows from west to east, while at night, it blows from southeast to north and the wind is brisker in the afternoon. The area of Pir Chinasi is found to be dry Subtropical type in which Acacia modesta, Olea ferruginea, and Chir Pine are the most dominant tree species in the area29. The vegetation of the area is comprised of a spacious variety of herbs, trees, climbers, and shrubs. The ground cover is comprised of a wide variety of angiosperms, mosses, and ferns.Fig. 1 Location map of the Pir Chinasi National Park, Muzaffarabad (District of Azad Jammu & Kashmir). The map was created using ArcGIS Desktop software (version 10.8.2, Esri lnc., Redlands,CA, USA).

The western Himalayan region, a biodiversity hotspot, is impacted in several ways by anthropogenic climate change. Due to the great seasonality of the western Himalayan region, including the Muzaffarabad district, investigations on the interactions between the timing of phenological periodic occurrences and climatic seasonality are important. The different land users of Pir Chinasi National Park benefit from the area’s natural resources and environment. During dry seasons, herders who live on the mountain frequently light fires to encourage the growth of forage for cattle grazing. This study assesses the dynamics of vegetation degradation during five years (1995–2020). This serves as the baseline information for tracking the watershed-level progress towards the objectives and pathways of SDGs 15.3.1 and 15.4.2 related to land degradation neutrality.

Materials and methods

Data collection

This study was conducted in Pir Chinasi National Park, Muzaffarabad, Kashmir, Pakistan. Pir Chinasi National Park is a tourist spot, and many tourists visit this area annually. A one-shot Household Survey (HHS) was conducted from 200 residents using convenient sampling in this study. It was a hilly area with a scattered population, so it was difficult to collect a larger sample size. A team of enumerators and the chief investigator spoke with the heads of the families. Family members were contacted in case of the unavailability of the head of household. The purpose of each question of the survey was explained to the respondents. Participants were provided with informed consent, and they were informed of the study’s advantages, goals, and financing. All surveys were conducted following the shared research principles and ethics. The survey asked about households’ occupation, age, gender, education, environmental harm to vegetation, and government support. Moreover, plants/seeds were handled under the direct supervision of Dr. Naeem Ahmed, National University of Modern Languages, Islamabad-Pakistan following the proper national and international strategies.

Data processing

The authors confirmed that the collection and execution of the experiment complied with the IUCN statement on research, involving species at risk of extinction and by the convention on the trade in endangered species of wild fauna and flora. For data collection, Google Forms were utilized and for data processing and analysis, SPSS was used. All methods were performed following the relevant guidelines and regulations.

To investigate the effects of human activities on the flora and ecology of the area, Remote Sensing (RS) and Geographic Information System (GIS) techniques and datasets have been utilized in Table 1. RS and GIS techniques datasets were proven excellent for large areas and difficult terrain30,31. It was difficult to collect data manually in such an area; hence, remote detection and data collection are more appropriate and easier.Table 1 Landsat satellite imagery of the years 1995 to 2020.

Scene ID	Date of acquisition	Sensor	Resolution	
LT05_L2SP_150036_19950508_20200913_02_T1	08. 05. 1995	Landsat 5	30 m	
LE07_L2SP_150036_20000411_20200918_02_T1	11. 05. 2000	Landsat 7	30 m	
LE07_L2SP_150036_20050527_20200914_02_T1	27. 05. 2005	Landsat 7	30 m	
LE07_L2SP_150036_20100509_20200911_02_T1	09. 05. 2010	Landsat 7	30 m	
LC08_L2SP_150036_20150515_20200909_02_T1	15.05.2015	Landsat 8	30 m	
LC08_L2SP_149035_20201015_20211025_02_T1	15. 05. 2020	Landsat 8	30 m	
LC08_L2SP_149035_20211015_20211025_02_T1	15. 05. 2021	Landsat 8	30 m	
LC08_L2SP_150036_20200512_20200820_02_T1	12.05.2020	Landsat 9	30 m	
aTropical Livestock Units (TLU).

Landsat data

Landsat imagery is an efficient data source for analyzing mountain landscapes, particularly difficult terrain because of its wide and remote coverage. In this study, Landsat 4/5 Thematic Mapper (TM), 7 Enhanced Thematic Mapper Plus (ETM +) & 8 Operational Land Imager (OLI) data were utilized for the temporal analysis and vegetation indices analysis. Multiple Landsat images were downloaded from the official website of USGS (www.usgs.gov).

The spatial resolution of Landsat 7 is 30 × 30 m, the revisit time of the satellite is 16 days whereas, the swath width of the satellite is 185 km. Researchers used this data in different kinds of studies due to its large area coverage, fine resolution, and result accuracy, therefore in this study, Landsat data was selected to be used.

Digital elevation model

The Digital Elevation Model (DEM) is an effective mode of 3D representation of the earth’s surface. DEM with a spatial resolution of 30m is utilized for ground elevation estimation. DEM of 30-m resolution have been utilized in the study to check the elevation and Slope of the area. Elevation and slope affect the LULC of the area, therefore 30-m resolution was sufficient for the required analysis. DEM data is freely available on the official website. Data is freely available and can be downloaded from the official website of Earth Explorer.

After data collection from the various websites, data was segregated and compiled for analysis. The study area was analyzed for the twenty-five years from 1995 to 2020 and time series data was collected with a five-year gap e.g. (1995, 2000, 2005, 2010, 2015, and 2020). Scan Line error (SLE) was present in Landsat 7 (2005 & 2010 Imagery). SLE was removed from the Landsat 7 imagery via the Landsat toolbox extension.

Error handling

Due to the missing data of the Scan Line Corrector (SLC) of the Landsat instrument in 2003, approximately 20 to 22% of data went missing, and gaps were generated in the imagery. This error made the image difficult to detect. Therefore, the technique of Gap-Fill proposed by the reference32 had been applied and error was removed from the 2007 and 2010 imagery.

Index calculation

Raw satellite data was collected freely from the website for further processing and calculation of the vegetation indices. The satellite data was converted from Digital Numbers (DN) to the values of the Top-of Atmospheric Reflectance (TOA). It was a two-step method (1) radiance of DN and (2) radiance to TOA33,34.

Due to the presence of a Scan Line Error (SLE) in Landsat 7 DN, conversion was done manually. The conversion of Landsat 7 data to DN was presented in Eq. 1.1 Lλ=gainλ×DN7+biasλ

In this equation where Lλ was found to be the calculated radiance [in Watts per square meter ∗ μm ∗ ster)], DN7 was the Landsat 7 ETM + DN data, and the gain and bias were band-specific numbers21. Radiance to TOA conversion was also done for Landsat 7. The reflectance could be thought of as a “planetary albedo”, or a fraction of the sun’s energy that was reflected by the surface (22; Eq. (2).2 Rλ=Π×Lλ×d2Esun,λ×sinθSE

Rλ = π ∗ Lλ ∗ d2 Esun,λ ∗ sin (θSE) (3) where Rλ was the reflectance (unitless ratio), Lλ was the radiance calculated in Eq. (2), d was the earth-sun distance (in astronomical units), Esun, λ was the band-specific radiance emitted by the sun, π was a constant value. The process of (TOA) was done for Landsat 8 via band 10. Thermal Infra-Red DN could be converted into TOA spectral radiance utilizing rescaling factors. A Set of equations had been used for LST calculation. DN to TOA conversion was presented in equation Eq. (3).3 Lλ=ML×Qcal+AL-Oi

where, Lλ = TOA spectral radiance (in Watts/(sq. meter ∗ μm), ML = Radiance multiplicative Band No, AL = Radiance add band, Qcal = DN, Qi band correction value. The temperature was converted from Kelvin to Celsius via Eq. (4).4 BT=K2Ink1Lλ+1-273.15

Finally, NDVI was calculated. The reflectance known as NDVI, a measure of greenness and proxy for vegetation degradation, was captured (Eq. 5).5 NDVI=NIR - REDNIR + RED

In this equation RED = DN value of RED band and NIR = DN value of Near-Infrared band, after that land Surface Emissivity was calculated via square of NDVI. Lastly, LST was calculated for the selected study area using Eq. 66 LST=BT1+λ×BTc2×InE

In this equation BT = TOA, whereas, brightness in 0C, λ = wavelength, similarly E = Land surface Emissivity, and C = velocity of light.

Results

Demographic and socioeconomic information

The data in Table 2 shows that the surveyed population is primarily composed of males (56.7%) and females (43.3%). The percentages were calculated based on the total number of respondents.Table 2 Gender distribution of the surveyed population.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	Male	89	56.7	56.7	56.7	
Female	68	43.3	43.3	100.0	
Total	157	100.0	100.0		

The data in Table 3 indicates that the majority of the sample respondents fell within the 18–24 years age range (60.5%), followed by smaller proportions in older age ranges. This age distribution provides insight into the demographic composition of the surveyed population and can be useful for understanding how different age groups perceive and respond to the issues related to protecting natural resources.Table 3 Age distribution of the surveyed population.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	18–24 Years	95	60.5	60.5	60.5	
25–34 Years	20	12.7	12.7	73.2	
35–44 Years	23	14.6	14.6	87.9	
45–54 Years	9	5.7	5.7	93.6	
55–64 Years	7	4.5	4.5	98.1	
65 and above	3	1.9	1.9	100.0	
Total	157	100.0	100.0		

The data in Table 4 provides insight into the different occupational roles present within the surveyed population. The majority of respondents were students (53.5%), followed by government/private/NGO employees (18.5%) and individuals with other occupations such as housewives, the unemployed, and those involved in agriculture, livestock farming, daily wage labor, and business. This distribution can offer context for understanding how various segments of the population view and engage with natural resource protection measures.Table 4 The occupational distribution of the sample respondents.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	Government/Private/NGO employee	29	18.5	18.5	18.5	
Student	84	53.5	53.5	72.0	
Housewife	8	5.1	5.1	77.1	
Unemployed	8	5.1	5.1	82.2	
Agriculture	6	3.8	3.8	86.0	
Livestock farming	10	6.4	6.4	92.4	
Daily wage laborer	5	3.2	3.2	95.5	
Business	7	4.5	4.5	100.0	
Total	157	100.0	100.0		

The data in Table 5 provides insight into the educational backgrounds of the surveyed individuals. The majority had a college/university degree (59.2%), followed by postgraduate degree holders (14.0%). Smaller proportions had secondary education, primary education, or no formal education. This distribution offers context for understanding how education levels might influence perceptions and attitudes toward natural resource protection measures.Table 5 The educational distribution of the surveyed population.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	Postgraduate degree or Higher	22	14.0	14.0	14.0	
College/University degree	93	59.2	59.2	73.2	
Secondary education	13	8.3	8.3	81.5	
Primary education	12	7.6	7.6	89.2	
No formal education	17	10.8	10.8	100.0	
Total	157	100.0	100.0		

Table 6 provides insight into the household size of the surveyed individuals. The most common household sizes were those with 5–6 people (35.0%) and 7–8 people (25.5%). It is observed that the educated respondents had more awareness about environmental degradation and the impact of anthropogenic activities. This distribution can offer context for understanding the living arrangements of the surveyed population, which might impact their resource consumption patterns and attitudes toward environmental protection.Table 6 The distribution of the size of households of the surveyed individuals.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	1–2	3	1.9	1.9	1.9	
3–4	30	19.1	19.1	21.0	
5–6	55	35.0	35.0	56.1	
7–8	40	25.5	25.5	81.5	
9–10	15	9.6	9.6	91.1	
More than 10	14	8.9	8.9	100.0	
Total	157	100.0	100.0		

In each group of Table 7, the frequency and percent values indicate the number and proportion of respondents involved in the specified type of agricultural activity. This data provided insights into the agricultural practices and activities of the surveyed population, showing the prevalence of different crops and livestock within their households. This information can be useful for understanding the composition of their agricultural practices and their potential impact on the environment and natural resources.Table 7 Distribution of different crops and livestock groups within the surveyed population.

Group	Crop or livestock	Frequency (N = 157)	Percent	
Group 1: Crops (Wheat; Maize; Rice)	Wheat	78	49.7	
Rice	34	21.7	
Maize	37	23.6	
Group 2: Vegetables and Fruits (Vegetables; Fruits)	Vegetables	56	35.7	
Fruits	51	32.5	
Group 3: Livestock (Cattle; Sheep; Goats; Poultry)	Cattle	51	32.5	
Sheep	41	26.1	
Goats	54	34.4	
Poultry	40	25.5	

Table 8 provides insights into the purpose of agricultural activities within the sample respondents. The majority of respondents engaged in subsistence farming (49.0%), followed by those involved in both subsistence and commercial farming (38.9%), and a smaller proportion focused on commercial farming alone (12.1%). Although subsistence farming contains a major share of the area, but still commercial farming has a second major share. Commercial farming uses fertilizers and pesticides which harm human health as well as deteriorate the environment. This distribution sheds light on the nature of agricultural practices within the surveyed community and how they might impact natural resource utilization and conservation efforts.Table 8 Distribution of agricultural practices by the respondents.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	Subsistence	77	49.0	49.0	49.0	
Commercial	19	12.1	12.1	61.1	
Both	61	38.9	38.9	100.0	
Total	157	100.0	100.0		

Table 9 revealed that more than half of the surveyed individuals were aware of the environmental issues in the Pir Chinasi National Park (52.9%), while the remaining respondents were not aware (47.1%). This awareness information was crucial for understanding the level of knowledge within the surveyed population about the environmental challenges in the region and can guide efforts to improve awareness and engagement in addressing those issues.Table 9 Awareness of the surveyed population regarding environmental issues in the Pir Chinasi National Park.

	Frequency	Percent	Valid Percent	Cumulative Percent	
Valid	Yes	83	52.9	52.9	52.9	
No	74	47.1	47.1	100.0	
Total	157	100.0	100.0		

These percentages in Table 10 indicate the proportion of respondents who perceive each anthropogenic activity as a significant contributor to the depletion of natural resources in the area. The respondents’ opinions pointed to deforestation, road construction, livestock grazing, and tourism as the main activities causing resource depletion, with varying levels of consensus among the surveyed population. This information can be valuable for understanding the perceived drivers of environmental degradation and a reference35 can help guide strategies for mitigating these activities to ensure sustainable resource management36.Table 10 Respondents’ opinions about anthropogenic activities causing the depletion of natural resources in Pir-Chinasi.

Anthropogenic activities	Frequency	Percent	
(N = 157)	
Deforestation	114	72.6	
Road construction	68	43.3	
Agriculture	30	19.1	
Livestock grazing	43	27.4	
Tourism	62	39.5	

Table 11 indicates the varied perceptions within the surveyed population regarding the impact of climate change on the natural resources of the Pir Chinasi National Park. Responses ranged from highly negative to highly positive, with a significant portion of respondents expressing a somewhat negative perception. These responses provide insight into the diverse perspectives on the potential consequences of climate change on the local environment and resources31.Table 11 Perceptions of respondents about the impact of climate change on the natural resources of the Pir Chinasi National Park area.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	Highly negative	29	18.5	18.5	18.5	
Somewhat negative	38	24.2	24.2	42.7	
Somewhat positive	44	28.0	28.0	70.7	
Somewhat positive	11	7.0	7.0	77.7	
Highly positive	3	1.9	1.9	79.6	
Total	157	100.0	100.0		

The percentages in Table 12 indicate the proportion of respondents who utilize each type of fuel for cooking and heating purposes. The data provides insights into the fuel preferences within the surveyed community, which have significant implications for energy consumption patterns and their potential environmental impact. The use of firewood, LPG, natural gas, and electric heaters appears to be relatively common among the surveyed population37.Table 12 Types of fuel used for cooking and heating within the surveyed population.

Fuel for cooking and heating	Frequency (N = 157)	Percent	
Firewood	77	49	
Charcoal	16	10.2	
LPG	47	29.9	
Natural Gas	54	34.4	
Kerosene	08	5.1	
Electric Heater	17	10.8	

Table 13 reveals that a significant portion of the surveyed population did not have access to alternative sources of energy, such as solar power (65.0%), while a minority did have access (35.0%). This information provides insights into the availability and adoption of renewable energy solutions within the surveyed community36,37, which can have significant implications for energy security, environmental sustainability, and quality of life.Table 13 Accessibility of alternative sources of energy, such as solar power, within the surveyed population.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	Yes	55	35.0	35.0	35.0	
No	102	65.0	65.0	100.0	
Total	157	100.0	100.0		

Table 14 provides insights into the livestock ownership patterns within the surveyed community. The majority owned 1–5 animals (43.9%), with smaller proportions owning 6–10 animals (17.2%) and 11–20 animals (7.6%). A small percentage owned more than 50 animals (1.3%), and a significant proportion did not own any livestock (29.9%). This distribution offers an understanding of the diversity in livestock ownership and its potential implications for resource use and management.Table 14 Livestock animals within the surveyed population.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	1–5 (TLU1)	69	43.9	43.9	43.9	
6–10 (TLU)	27	17.2	17.2	61.1	
11–20 (TLU)	12	7.6	7.6	68.8	
More than 50 (TLU)	2	1.3	1.3	70.1	
N/A	47	29.9	29.9	100.0	
Total	157	100.0	100.0		
TLU tropical livestock units

Table 15 provides insight into where the surveyed individuals’ livestock are allowed to graze. The responses indicated that animals were allowed to graze on private land (26.1%), common land (27.4%), and forests (17.2%), while a notable percentage did not have animals that graze (29.3%). Understanding the locations where livestock is grazed can help assess the potential impact on these areas and inform resource management strategies.Table 15 Locations where the surveyed population’s animals are allowed to graze.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	Private land	41	26.1	26.1	26.1	
Common land	43	27.4	27.4	53.5	
Forests	27	17.2	17.2	70.7	
N/A	46	29.3	29.3	100.0	
Total	157	100.0	100.0		

Table 16 reveals that a majority of the surveyed population had observed changes in the natural resources of the Pir Chinasi National Park area over the past 5–10 years (59.9%), whereas a smaller proportion (40.1%) has not perceived such changes. This information provides insights into the perceived dynamics of natural resource changes in the region and can help in understanding the evolving environmental conditions and potential factors contributing to these changes.Table 16 Perceptions of the surveyed population about any changes in the natural resources (forests, soil, water, and biodiversity) in the Pir Chinasi National Park area over the past 5–10 years.

	Frequency	Percent	Valid percent	Cumulative percent	
Valid	Yes	94	59.9	59.9	59.9	
No	63	40.1	40.1	100.0	
Total	157	100.0	100.0		

The percentages in Table 17 indicate the proportion of respondents who believed that each measure should be taken to protect the natural resources in the area. It appears that reforestation programs had the highest agreement (58.0%), followed by strict enforcement of environmental protection laws (53.2%). Other measures such as public awareness campaigns, collaboration with local communities, sustainable tourism practices, and sustainable agriculture practices also received notable support from the surveyed population. These opinions provide insight into the potential strategies that could be pursued to safeguard the natural resources in the Pir Chinasi National Park area according to the surveyed community.Table 17 Opinions of the surveyed population on additional measures that should be taken to protect the natural resources in the Pir Chinasi National Park area.

Additional measures should be taken to protect the natural resources	Frequency (N = 157)	Percent	
Strict enforcement of environmental protection laws	84	53.2	
Reforestation programs	91	58	
Sustainable agriculture practices	44	28	
Sustainable tourism practices	56	35.7	
Public awareness campaigns	76	48.7	
Collaboration with local communities	57	36.3	

NDVI

Time series analysis of vegetation change dynamics of the selected area was evaluated via NDVI. Spatial and temporal dynamics of vegetation degradation, depending on the time scale were analyzed. The status of the landscape was reflected in vegetation indices (VIs), and their interpretation across time explained patterns in vegetative greening (land improvement) and browning (land deterioration). Based on a long-term investigation of vegetation dynamics, Figs. 2–4 illustrate swings between periods of gradual deterioration and minor improvements (greening). A slightly declining trend of NDVI was observed during the study period (1995–2020). During the entire study period mean NDVI (maximum = 0.28 (1995)/minimum − 2.8 (2010) was identified see Fig. 2 for instance. In 2000, mean NDVI was monitored at 0.24 (Fig. 2b). Similarly, the NDVI of 2005 was also calculated via the set of formulas mentioned in the materials and method Section38–40.Fig. 2 shows the satellite-based NDVI variation temporally from 1995 to 2020.

Fig. 3 Slope patterns of the study area.

Fig. 4 Land Surface Temperature (LST) ° C trends of the study area during the study period. The map was created using ArcGIS Desktop software (version 10.8.2, Esri lnc., Redlands,CA, USA). ArcGIS is widely used for creating slope maps from Digital Elevation Models (DEMs).

The mean value of NDVI was 0.22 to 0.76. Pir Chinasi National Park was known for tourism, this decline can be attributed to the expansion of built-up areas to accommodate growing tourist demands, leading to habitat loss and ecological changes. This decline can be attributed to the expansion of built-up areas to accommodate growing tourist demands, leading to habitat loss and ecological changes. NDVI for the year 2015 and 2020 was observed 0.07–0.53 to 0.08–0.57 (Fig. 2e and f).

Land surface temperature (LST)

Surface temperature had a significant impact on the local ecology and biodiversity. Therefore, in this study, LST for the study period was also calculated to investigate the impacts of temperature on vegetation. The temperature was estimated via Landsat imagery for the selected images from 1995 to 2020 (Fig. 4). The temperature range for 1995 was 19 minimum to 29o C maximum (Fig. 4a), 23° to 36° for 2000 (Fig. 4b), whereas 21° to 35° temperature range was observed for 2005 (Fig. 4c). 29 °C to 39 °C for 2010 see Fig. 4d. For the year 2015 and 2020 (9°–26 °C and 17–37 °C) were observed respectively.

It has tools specifically for slope calculation, where the slope can be expressed in degrees or as a percentage. Slope analysis, spatial analysis, and 3D visualization using extensions like Spatial Analyst or 3D Analyst. User-friendly interface, extensive documentation, and strong support for various geospatial analyses.

Ground truthing

Data collection from the field was time-consuming due to the monetary constraints23,24,33, therefore, Google Earth data was utilized for the investigation of urbanization patterns. Deforestation was observed where recreational activities were occurring, see Fig. 5.Fig. 5 Google Earth image of Pir Chinasi National Park.

Figure 5 highlighted the ground truthing based on detailed information about vegetation to validate the remote sensing. The figure shows urbanization patterns. Ground truthing with physical sampling strengthened the validation of the NDVI analysis. Map was generated using the open-source software (QGIS 3.8) along with that an open-source software Google Earth Pro (GEP) was utilized for image data collection. KML files were generated using GEP and converted into shapefiles in the QGIS environment. The geographical coordinates (34°23′22"N 73°32′57"E) were used. https://earth.google.com/web/search/Pir+Chinasi/@34.38987286,73.55007956,2825.08631358a,2608.57706282d,35y,0h,0t,0r/data=CigiJgokCXkkdsjBB0FAEWoFBBCW-EBAGZ7DNELhfFJAIah7yxvnaFJAOgMKATA

Discussion

Numerous reports have been published relating NDVI to rainfall patterns. In a similar vein, by reference41 discovered that precipitation in Ethiopia’s semi-arid areas was strongly correlated with NDVI levels throughout the growing season. According to the reference42, 80% of the Caatinga vegetation productivity anomaly may be accounted for by the rainfall anomaly using the vegetation index as a proxy taken by reference43 in the Horn of Africa reported similar findings.

In China, a significant link between NDVI and precipitation was discovered by researchers17,44. According to reference45, NDVI across Northern Mongolia showed a positive association with yearly precipitation and a negative correlation with temperature16. Researchers13,46,47, found that NDVI was more closely connected with rainfall than temperature on the Tibetan Plateau and the Loess Plateau, respectively.

In actuality, when evaluating vegetation response to climate change the perceptions of the residents can also be taken48. In this study, the respondents perceived a major change in vegetation due to climate change and anthropogenic activities. A weak or negative (−) NDVI-Rainfall link implies that vegetation dynamics are influenced by the climate; a robust and positive ( +) connection indicates changes in vegetation caused by humans. The researcher49 also found comparable results in its worldwide survey.

According to researcher50, precipitation is the key environmental element influencing NDVI changes in the Heihe River Basin, China. According to another researcher51, precipitation and land cover in Kenya’s Mara River Basin revealed a highly positive association. According to reports, the primary human activities responsible for the deterioration of the local vegetation and the disappearance of the forest in this region include grazing, agriculture-related settlement growth, and wood collecting52,53. According to researcher51, the highlands’ growing pace of forest degradation is a sign of how much the high population density, 350inch/km22,24,39 value and exploit this forest on an economic, social, ethnological, and cultural level. It has been studied how some of the aforementioned factors relate to NDVI. The NDVI of the Wei and Jing River Basins, China, was shown to be primarily influenced by three variables: temperature, soil moisture, and precipitation. The researcher54 demonstrated a tight link between the dynamics of grassland and temperature using NDVI-max at an annual time scale. Future studies could focus on exploring the connection between the environmental elements and the dynamics of the vegetation on the Pir Chinasi National Park.

Conclusions

The study described the dynamics of vegetation and how it responded to changing rainfall patterns and human pressures. It was proven that vegetation change dynamics and degradation are impacted by climate change as indicated by long-term rainfall anomaly (RAI). The study found that 10.81% of the area experienced substantial deterioration, while 23.57% experienced minor degradation. The total area of degraded lands was 34.38% suggesting a major change. To overcome these phenomena, steps of adaptation by the public–private sector are necessary. Improvement in plant cover was observed which is a positive sign. According to the NDVI-Rainfall relationships, the area was found to be significantly impacted by human pressures which require control by public–private partnership. The area under the influence of climatic variability and change is also large enough which requires immediate attention. Sustainable land management practices of conservation tillage, integrated pest management, and agroforestry help preserve soil health, water quality, and biodiversity while reducing erosion, pollution, and the degradation of natural resources. landscape restoration projects of reforestation, wetland restoration, soil erosion control, and the removal of invasive species are essential to achieve land degradation neutrality at the watershed scale. The quality and health of the vegetation were further harmed by unsustainable human activities such as agricultural development, overgrazing, settlement growth, and wood exploitation. However, climate change proved to be a driver of these phenomena. Even though these were the major causes, the degree to which the flora on the plateau tends to green up or brown out, could also be influenced by environmental variables such as soil moisture, solar radiation intensity and duration, drainage density, and topographic restrictions. This study does not conclude that the sole variables affecting vegetation dynamics are rainfall and human influences.

In areas of the park where vegetation patterns indicate deterioration or loss of forest cover, afforestation, and replanting efforts could be done. These initiatives have the potential to improve ecosystem resilience to climate change, restore biodiversity, and sequester carbon. Implement sustainable land management techniques that boost agricultural production, and preserve and improve plant cover, such as conservation agriculture, rotational grazing, and agroforestry. These methods could lessen soil erosion, increase water retention, and protect the benefits provided by ecosystem services. In response to climate change and human intervention establish buffer zones, wildlife corridors, and protected areas to preserve important ecosystems and enhance landscape connectivity.

Acknowledgements

The authors extend their appreciation to Researchers Supporting Project number (RSP-2024 R369), King Saud University, Riyadh, Saudi Arabia.

Author contributions

Conceptualization, Muhammad Akhlaq Farooq and Naeem Ahmed; Data curation, Muhammad Akhlaq Farooq and Muhammad Shahbaz Farooq; Formal analysis, Muhammad Akhlaq Farooq, Muhammad Ghufran, Aamir Amanat, and Saima Naz; Funding acquisition, Naeem Ahmed, Muhammad Shahbaz Farooq, and Muhammad Uzair; Investigation, Muhammad Ghufran, Yasser M Hafez, and Aamir Amanat; Methodology, Aamir Amanat and Saima Naz; Project administration, Naeem Ahmed; Resources, Naeem Ahmed, Muhammad Uzair and Saima Naz; Software, Aamir Amanat; Supervision, Naeem Ahmed; Validation, Aamir Amanat; Visualization, Muhammad Akhlaq Farooq, Kotb A. Attia, Arif Ahmed Mohammed; Writing – original draft, Muhammad Akhlaq Farooq; Writing – review & editing, Muhammad Akhlaq Farooq, Muhammad Ghufran, Naeem Ahmed, Muhammad Shahbaz Farooq, Muhammad Uzair, Kotb A. Attia, Arif Ahmed Mohammed and Saima Naz. All authors have read and agreed to the published version of the manuscript.

Data availability

Data is provided within the manuscript or supplementary information files.

Competing interests

The authors declare no competing interests.

Research Ethics

Experimental research and field studies on plants (either cultivated or wild), including the collection of plant material, complied with relevant institutional, national, and international guidelines and legislation. Prior approval was undertaken from the Offices of Research, Innovation and Commercialization, National University of Modern Languages, Islamabad-Pakistan. We also took appropriate permission from the farm or field owner during specimens’ collection and experimentation. We confirm that during the collection and execution of the experiment, the authors have complied with the IUCN Statement on Research Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora. All methods were performed in accordance with the relevant guidelines and regulations. All experimental protocols were approved by Offices of Research, Innovation and Commercialization, National University of Modern Languages, Islamabad-Pakistan.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Guan X Liu M Meng Y A comprehensive ecological compensation indicator based on pollution damage—Protection bidirectional model for river basin Ecol. Indic. 2021 126 107708 10.1016/j.ecolind.2021.107708
Guan, X., Liu, M. & Meng, Y. A comprehensive ecological compensation indicator based on pollution damage—Protection bidirectional model for river basin. Ecol. Indic. 126, 107708 (2021).10.1016/j.ecolind.2021.107708
2. Fu B Current condition and protection strategies of Qinghai-Tibet Plateau ecological security barrier Bull. Chinese Acad. Sci. 2021 36 1298 1306
Fu, B. et al. Current condition and protection strategies of Qinghai-Tibet Plateau ecological security barrier. Bull. Chinese Acad. Sci. 36, 1298–1306 (2021).
3. Wen Z Wu S Chen J Lü M NDVI indicated long-term interannual changes in vegetation activities and their responses to climatic and anthropogenic factors in the three gorges reservoir region, China Sci. Total Environ. 2017 574 947 959 10.1016/j.scitotenv.2016.09.049 27665454
Wen, Z., Wu, S., Chen, J. & Lü, M. NDVI indicated long-term interannual changes in vegetation activities and their responses to climatic and anthropogenic factors in the three gorges reservoir region, China. Sci. Total Environ. 574, 947–959 (2017).27665454 10.1016/j.scitotenv.2016.09.049
4. Shiweda M Shivute F Sales AR Pereira MJ Climate change and anthropogenic factors are influencing the loss of habitats and emerging human-elephant conflict in the Namib desert Sustainability 2023 15 12400 10.3390/su151612400
Shiweda, M., Shivute, F., Sales, A. R. & Pereira, M. J. Climate change and anthropogenic factors are influencing the loss of habitats and emerging human-elephant conflict in the Namib desert. Sustainability 15, 12400 (2023).10.3390/su151612400
5. Du Y Human activity coupled with climate change strengthens the role of lakes as an active pipe of dissolved organic matter Earth’s Futur. 2023 10.1029/2022EF003412
Du, Y. et al. Human activity coupled with climate change strengthens the role of lakes as an active pipe of dissolved organic matter. Earth’s Futur.10.1029/2022EF003412 (2023).10.1029/2022EF003412
6. Gabriele M Brumana R Previtali M Cazzani A A combined GIS and remote sensing approach for monitoring climate change-related land degradation to support landscape preservation and planning tools: The Basilicata case study Appl. Geomat. 2023 15 497 532 10.1007/s12518-022-00437-z
Gabriele, M., Brumana, R., Previtali, M. & Cazzani, A. A combined GIS and remote sensing approach for monitoring climate change-related land degradation to support landscape preservation and planning tools: The Basilicata case study. Appl. Geomat. 15, 497–532 (2023).10.1007/s12518-022-00437-z
7. Subramanian A Long-term impacts of climate change on coastal and transitional eco-systems in India: An overview of its current status, future projections, solutions, and policies RSC Adv. 2023 13 12204 12228 10.1039/D2RA07448F 37091602
Subramanian, A. et al. Long-term impacts of climate change on coastal and transitional eco-systems in India: An overview of its current status, future projections, solutions, and policies. RSC Adv. 13, 12204–12228 (2023).37091602 10.1039/D2RA07448F
8. Zheng K Impacts of climate change and anthropogenic activities on vegetation change: Evidence from typical areas in China Ecol. Indic. 2021 126 107648 10.1016/j.ecolind.2021.107648
Zheng, K. et al. Impacts of climate change and anthropogenic activities on vegetation change: Evidence from typical areas in China. Ecol. Indic. 126, 107648 (2021).10.1016/j.ecolind.2021.107648
9. Du J Spatiotemporal changes in vegetation coverage and its causes in China since the Chinese economic reform Environ. Sci. Pollut. Res. 2020 27 1144 1159 10.1007/s11356-019-06609-6
Du, J. et al. Spatiotemporal changes in vegetation coverage and its causes in China since the Chinese economic reform. Environ. Sci. Pollut. Res. 27, 1144–1159 (2020).10.1007/s11356-019-06609-6
10. Lamchin M Long-term trend of and correlation between vegetation greenness and climate variables in Asia based on satellite data MethodsX 2018 5 803 807 10.1016/j.mex.2018.07.006 30105213
Lamchin, M. et al. Long-term trend of and correlation between vegetation greenness and climate variables in Asia based on satellite data. MethodsX 5, 803–807 (2018).30105213 10.1016/j.mex.2018.07.006
11. Pettorelli N Using the satellite-derived NDVI to assess ecological responses to environmental change Trends Ecol. Evol. 2005 20 503 510 10.1016/j.tree.2005.05.011 16701427
Pettorelli, N. et al. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends Ecol. Evol. 20, 503–510 (2005).16701427 10.1016/j.tree.2005.05.011
12. Sharon E Nicholson MLD A comparison of the vegetation response to rainfall in the sahel and east africa, using normalized difference vegetation from NOAA AVHRR Clim. Change 1990 21 99 104
Sharon, E. & Nicholson, M. L. D. A comparison of the vegetation response to rainfall in the sahel and east africa, using normalized difference vegetation from NOAA AVHRR. Clim. Change 21, 99–104 (1990).
13. Chu D Lu L Zhang T Sensitivity of normalized difference vegetation index (NDVI) to seasonal and interannual climate conditions in the Lhasa area, Tibetan Plateau, China Arctic Antarct. Alp. Res. 2007 39 635 641 10.1657/1523-0430(07-501)[CHU]2.0.CO;2
Chu, D., Lu, L. & Zhang, T. Sensitivity of normalized difference vegetation index (NDVI) to seasonal and interannual climate conditions in the Lhasa area, Tibetan Plateau, China. Arctic Antarct. Alp. Res. 39, 635–641 (2007).10.1657/1523-0430(07-501)[CHU]2.0.CO;2
14. Zawude Bakure B Hundera K Abara M Review on the effect of climate change on ecosystem services IOP Conf. Ser. Earth Environ. Sci. 2022 1016 012055 10.1088/1755-1315/1016/1/012055
Zawude Bakure, B., Hundera, K. & Abara, M. Review on the effect of climate change on ecosystem services. IOP Conf. Ser. Earth Environ. Sci. 1016, 012055 (2022).10.1088/1755-1315/1016/1/012055
15. Salimon C Anderson L How strong is the relationship between rainfall variability and caatinga productivity? A case study under a changing climate An. Acad. Bras. Cienc. 2018 90 2121 2127 10.1590/0001-3765201720170143 28538819
Salimon, C. & Anderson, L. How strong is the relationship between rainfall variability and caatinga productivity? A case study under a changing climate. An. Acad. Bras. Cienc. 90, 2121–2127 (2018).28538819 10.1590/0001-3765201720170143
16. Nanzad L NDVI anomaly for drought monitoring and its correlation with climate factors over Mongolia from 2000 to 2016 J. Arid Environ. 2019 164 69 77 10.1016/j.jaridenv.2019.01.019
Nanzad, L. et al. NDVI anomaly for drought monitoring and its correlation with climate factors over Mongolia from 2000 to 2016. J. Arid Environ. 164, 69–77 (2019).10.1016/j.jaridenv.2019.01.019
17. Chen H The impacts of climate change and human activities on biogeochemical cycles on the Qinghai-Tibetan Plateau Glob. Chang. Biol. 2013 19 2940 2955 10.1111/gcb.12277 23744573
Chen, H. et al. The impacts of climate change and human activities on biogeochemical cycles on the Qinghai-Tibetan Plateau. Glob. Chang. Biol. 19, 2940–2955 (2013).23744573 10.1111/gcb.12277
18. Liu Z Liu Y Li Y Anthropogenic contributions dominate trends of vegetation cover change over the farming-pastoral ecotone of northern China Ecol. Indic. 2018 95 370 378 10.1016/j.ecolind.2018.07.063
Liu, Z., Liu, Y. & Li, Y. Anthropogenic contributions dominate trends of vegetation cover change over the farming-pastoral ecotone of northern China. Ecol. Indic. 95, 370–378 (2018).10.1016/j.ecolind.2018.07.063
19. Zhang B Spatial distribution and suitability evaluation of nighttime tourism in Kunming utilizing multi-source data Heliyon 2023 9 e16826 10.1016/j.heliyon.2023.e16826 37313168
Zhang, B. et al. Spatial distribution and suitability evaluation of nighttime tourism in Kunming utilizing multi-source data. Heliyon 9, e16826 (2023).37313168 10.1016/j.heliyon.2023.e16826
20. Liu Y Zhang M Liu R The impact of income inequality on carbon emissions in china: A household-level analysis Sustain 2020 12 1 22
Liu, Y., Zhang, M. & Liu, R. The impact of income inequality on carbon emissions in china: A household-level analysis. Sustain 12, 1–22 (2020).
21. XU J The melting Himalayas: Cascading effects of climate change on water, biodiversity, and livelihoods Conserv. Biol. 2009 23 520 530 10.1111/j.1523-1739.2009.01237.x 22748090
XU, J. et al. The melting Himalayas: Cascading effects of climate change on water, biodiversity, and livelihoods. Conserv. Biol. 23, 520–530 (2009).22748090 10.1111/j.1523-1739.2009.01237.x
22. Hussein AJ Salih SA Quantifying the reliability of volumetric and areal calculation with UAV-generated DEMs: A comparative study with ground truth data Tikrit J. Eng. Sci. 2024 31 43 55 10.25130/tjes.31.1.5
Hussein, A. J. & Salih, S. A. Quantifying the reliability of volumetric and areal calculation with UAV-generated DEMs: A comparative study with ground truth data. Tikrit J. Eng. Sci. 31, 43–55 (2024).10.25130/tjes.31.1.5
23. Sirianni H Richter J Sirianni MJ Pettyjohn S Shoreline classification maps and ground truth data for the Neuse River Estuary North Carolina Sci. Data 2024 11 103 10.1038/s41597-024-02954-5 38253576
Sirianni, H., Richter, J., Sirianni, M. J. & Pettyjohn, S. Shoreline classification maps and ground truth data for the Neuse River Estuary North Carolina. Sci. Data 11, 103 (2024).38253576 10.1038/s41597-024-02954-5
24. Kraatz S Comparing NISAR (Using Sentinel-1), USDA/NASS CDL, and ground truth crop/Non-crop areas in an urban agricultural region Sensors 2023 23 8595 10.3390/s23208595 37896688
Kraatz, S. et al. Comparing NISAR (Using Sentinel-1), USDA/NASS CDL, and ground truth crop/Non-crop areas in an urban agricultural region. Sensors 23, 8595 (2023).37896688 10.3390/s23208595
25. Zoungrana, B. J. B. Vegetation dynamics in the southwest of Burkina Faso in response to rainfall variability and land use. Thesis 162 (2016).
26. Musau J Patil S Sheffield J Marshall M Vegetation dynamics and responses to climate anomalies in East Africa Earth Syst. Dyn. Discuss. 2018 2018 1 27
Musau, J., Patil, S., Sheffield, J. & Marshall, M. Vegetation dynamics and responses to climate anomalies in East Africa. Earth Syst. Dyn. Discuss. 2018, 1–27 (2018).
27. Chinasi, P. Blog. 2024 https://eagleeye.com.pk/pttl/pir-chinasi/.
28. Ministry of Finance. Economic Survey of Pakistan 2018–19. Chapter 14- Energy 183–196 (2019).
29. Sharma N Behera MD Das AP Panda RM Plant richness pattern in an elevation gradient in the Eastern Himalaya Biodivers. Conserv. 2019 28 2085 2104 10.1007/s10531-019-01699-7
Sharma, N., Behera, M. D., Das, A. P. & Panda, R. M. Plant richness pattern in an elevation gradient in the Eastern Himalaya. Biodivers. Conserv. 28, 2085–2104 (2019).10.1007/s10531-019-01699-7
30. Ahmed N Spatial vulnerability assessment of urbanites of Rawalpindi to the combined effects of climate change and pandemics using GIS mapping Local Environ. 2024 10.1080/13549839.2024.2353048
Ahmed, N. Spatial vulnerability assessment of urbanites of Rawalpindi to the combined effects of climate change and pandemics using GIS mapping. Local Environ.10.1080/13549839.2024.2353048 (2024).10.1080/13549839.2024.2353048
31. Ahmed N Padda IUH Khan A Mapping adaptive capacity of urban residents: a case study of Rawalpindi Pakistan NUST J. Soc. Sci. Humanit. 2023 8 151 169
Ahmed, N., Padda, I. U. H. & Khan, A. Mapping adaptive capacity of urban residents: a case study of Rawalpindi Pakistan. NUST J. Soc. Sci. Humanit. 8, 151–169 (2023).
32. Muir C Southworth J Khatami R Herrero H Akyapı B Vegetation dynamics and climatological drivers in ethiopia at the turn of the century Remote Sens. 2021 13 3267 10.3390/rs13163267
Muir, C., Southworth, J., Khatami, R., Herrero, H. & Akyapı, B. Vegetation dynamics and climatological drivers in ethiopia at the turn of the century. Remote Sens. 13, 3267 (2021).10.3390/rs13163267
33. Puente M Vélez VN Ground-truthing as critical race Feminista methodology: Toward an embodied and community-centered GIS in educational inquiry Int. J. Qual. Stud. Educ. 2024 37 1287 1306 10.1080/09518398.2024.2336479
Puente, M. & Vélez, V. N. Ground-truthing as critical race Feminista methodology: Toward an embodied and community-centered GIS in educational inquiry. Int. J. Qual. Stud. Educ. 37, 1287–1306 (2024).10.1080/09518398.2024.2336479
34. Baltagi BH Random effects and spatial autocorrelation with equal weights Econ. Theory 2006 22 973 984 10.1017/S0266466606060464
Baltagi, B. H. Random effects and spatial autocorrelation with equal weights. Econ. Theory 22, 973–984 (2006).10.1017/S0266466606060464
35. Ahmed N Climate change adaption strategies in urban communities: New evidence from Islamabad Pakistan Environ. Sci. Pollut. Res. 2023 10.1007/s11356-023-25316-x
Ahmed, N. et al. Climate change adaption strategies in urban communities: New evidence from Islamabad Pakistan. Environ. Sci. Pollut. Res.10.1007/s11356-023-25316-x (2023).10.1007/s11356-023-25316-x
36. Ahmed N Luqman M Explaining urban communities’ adaptation strategies for climate change risk: Novel evidence from Rawalpindi Pakistan Nat. Hazards 2024 10.1007/s11069-024-06501-8
Ahmed, N. & Luqman, M. Explaining urban communities’ adaptation strategies for climate change risk: Novel evidence from Rawalpindi Pakistan. Nat. Hazards10.1007/s11069-024-06501-8 (2024).10.1007/s11069-024-06501-8
37. Batool Z Ahmed N Luqman M Examining the role of ICT, transportation energy consumption, and urbanization in CO2 emissions in Asia: A threshold analysis Environ. Sci. Pollut. Res. 2023 30 78482 78494 10.1007/s11356-023-27995-y
Batool, Z., Ahmed, N. & Luqman, M. Examining the role of ICT, transportation energy consumption, and urbanization in CO2 emissions in Asia: A threshold analysis. Environ. Sci. Pollut. Res. 30, 78482–78494 (2023).10.1007/s11356-023-27995-y
38. Bedair H Global warming status in the african continent: Sources, challenges, policies, and future direction Int. J. Environ. Res. 2023 17 45 10.1007/s41742-023-00534-w
Bedair, H. et al. Global warming status in the african continent: Sources, challenges, policies, and future direction. Int. J. Environ. Res. 17, 45 (2023).10.1007/s41742-023-00534-w
39. Gajendiran K Kandasamy S Narayanan M Influences of wildfire on the forest ecosystem and climate change: A comprehensive study Environ. Res. 2024 240 117537 10.1016/j.envres.2023.117537 37914016
Gajendiran, K., Kandasamy, S. & Narayanan, M. Influences of wildfire on the forest ecosystem and climate change: A comprehensive study. Environ. Res. 240, 117537 (2024).37914016 10.1016/j.envres.2023.117537
40. Garoon EM Mule MB A multiscale analysis and classification of normalized difference vegetation index value method in Kannad Taluka, Aurangabad District, using remote sensing and GIS J. Environ. Biol. 2023 44 335 344 10.22438/jeb/44/3/MRN-5055
Garoon, E. M. & Mule, M. B. A multiscale analysis and classification of normalized difference vegetation index value method in Kannad Taluka, Aurangabad District, using remote sensing and GIS. J. Environ. Biol. 44, 335–344 (2023).10.22438/jeb/44/3/MRN-5055
41. Zewdie, W. & Csaplovics, E. Evaluation of rainfall and NDVI anomalies using distributed lag models. in Algorithms and Technologies for Multispectral, Hyperspectral, and Ultraspectral Imagery XXI (eds. Velez-Reyes, M. & Kruse, F. A.) vol. 9472 94721O (2015).
42. Andersen AN Vasconcelos HL Historical biogeography shapes functional ecology: Inter-continental contrasts in responses of savanna ant communities to stress and disturbance J. Biogeogr. 2022 49 590 599 10.1111/jbi.14343
Andersen, A. N. & Vasconcelos, H. L. Historical biogeography shapes functional ecology: Inter-continental contrasts in responses of savanna ant communities to stress and disturbance. J. Biogeogr. 49, 590–599 (2022).10.1111/jbi.14343
43. Ghebrezgabher MG Yang T Yang X Eyassu Sereke T Assessment of NDVI variations in responses to climate change in the Horn of Africa Egypt J. Remote Sens. Sp. Sci. 2020 23 249 261
Ghebrezgabher, M. G., Yang, T., Yang, X. & Eyassu Sereke, T. Assessment of NDVI variations in responses to climate change in the Horn of Africa Egypt. J. Remote Sens. Sp. Sci. 23, 249–261 (2020).
44. Chen D How visitors perceive heritage value—A quantitative study on visitors’ perceived value and satisfaction of architectural heritage through SEM Sustainability 2023 15 9002 10.3390/su15119002
Chen, D. How visitors perceive heritage value—A quantitative study on visitors’ perceived value and satisfaction of architectural heritage through SEM. Sustainability 15, 9002 (2023).10.3390/su15119002
45. Enebish B Dashkhuu D Renchin M Russell M Singh P Impact of climate on the NDVI of northern Mongolia J. Indian Soc. Remote Sens. 2020 48 333 340 10.1007/s12524-019-01080-9
Enebish, B., Dashkhuu, D., Renchin, M., Russell, M. & Singh, P. Impact of climate on the NDVI of northern Mongolia. J. Indian Soc. Remote Sens. 48, 333–340 (2020).10.1007/s12524-019-01080-9
46. Hao Y States, trends, and future of aquaponics research Sustainability 2020 12 7783 10.3390/su12187783
Hao, Y. et al. States, trends, and future of aquaponics research. Sustainability 12, 7783 (2020).10.3390/su12187783
47. Ning T Liu W Lin W Song X NDVI variation and its responses to climate change on the Northern loess plateau of China from 1998 to 2012 Adv. Meteorol. 2015 2015 1 10 10.1155/2015/725427
Ning, T., Liu, W., Lin, W. & Song, X. NDVI variation and its responses to climate change on the Northern loess plateau of China from 1998 to 2012. Adv. Meteorol. 2015, 1–10 (2015).10.1155/2015/725427
48. Sapkale JB Mane MM Susware NK Sapkale SJ Dynamic changes in coastal geomorphology of Shiroda coasts, using remote sensing and GIS: An approach to climate change and coastal disaster risk Disaster Adv. 2023 16 20 32 10.25303/1612da020032
Sapkale, J. B., Mane, M. M., Susware, N. K. & Sapkale, S. J. Dynamic changes in coastal geomorphology of Shiroda coasts, using remote sensing and GIS: An approach to climate change and coastal disaster risk. Disaster Adv. 16, 20–32 (2023).10.25303/1612da020032
49. Bai ZG Dent DL Olsson L Schaepman ME Proxy global assessment of land degradation Soil Use Manag. 2008 24 223 234 10.1111/j.1475-2743.2008.00169.x
Bai, Z. G., Dent, D. L., Olsson, L. & Schaepman, M. E. Proxy global assessment of land degradation. Soil Use Manag. 24, 223–234 (2008).10.1111/j.1475-2743.2008.00169.x
50. Yuan L Chen X Wang X Xiong Z Song C Spatial associations between NDVI and environmental factors in the Heihe River Basin J. Geogr. Sci. 2019 29 1548 1564 10.1007/s11442-019-1676-0
Yuan, L., Chen, X., Wang, X., Xiong, Z. & Song, C. Spatial associations between NDVI and environmental factors in the Heihe River Basin. J. Geogr. Sci. 29, 1548–1564 (2019).10.1007/s11442-019-1676-0
51. Mngube FM Kapiyo R Aboum P Anyona D Dida GO Subtle impacts of temperature and rainfall patterns on land cover change overtime and future projections in the Mara river Basin, Kenya Open J. Soil Sci. 2020 10 327 358 10.4236/ojss.2020.109018
Mngube, F. M., Kapiyo, R., Aboum, P., Anyona, D. & Dida, G. O. Subtle impacts of temperature and rainfall patterns on land cover change overtime and future projections in the Mara river Basin, Kenya. Open J. Soil Sci. 10, 327–358 (2020).10.4236/ojss.2020.109018
52. Rawat A Kumar D Khati BS A review on climate change impacts, models, and its consequences on different sectors: A systematic approach J. Water Clim. Chang. 2024 15 104 126 10.2166/wcc.2023.536
Rawat, A., Kumar, D. & Khati, B. S. A review on climate change impacts, models, and its consequences on different sectors: A systematic approach. J. Water Clim. Chang. 15, 104–126 (2024).10.2166/wcc.2023.536
53. Vinet L Zhedanov A A ‘missing’ family of classical orthogonal polynomials J. Phys. A Math. Theor. 2011 44 1 22 10.1088/1751-8113/44/8/085201
Vinet, L. & Zhedanov, A. A ‘missing’ family of classical orthogonal polynomials. J. Phys. A Math. Theor. 44, 1–22 (2011).10.1088/1751-8113/44/8/085201
54. Huang S A novel fusion method for generating surface soil moisture data with high accuracy, high spatial resolution, and high Spatio-temporal continuity Water Resour. Res. 2022 58 1 26 10.1029/2021WR030827
Huang, S. et al. A novel fusion method for generating surface soil moisture data with high accuracy, high spatial resolution, and high Spatio-temporal continuity. Water Resour. Res. 58, 1–26 (2022).10.1029/2021WR030827
