Pakistan is a flood-prone country and almost every year, it is hit by floods of varying magnitudes. This study was conducted to generate a flash flood map using analytical hierarchy process (AHP) and frequency ratio (FR) models in the ArcGIS 10.6 environment. Eight flash-flood-causing physical parameters were considered for this study. Five parameters were based on the digital elevation model (DEM), Advanced Land Observation Satellite (ALOS), and Sentinel-2 satellite, including distance from the river and drainage density slope, elevation, and land cover, respectively. Two other parameters were geology and soil, consisting of different rock and soil formations, respectively, where both layers were classified based on their resistance against water percolation. One parameter was rainfall. Rainfall observation data obtained from five meteorological stations exist close to the Chitral District, Pakistan. According to its significant importance in the occurrence of a flash flood, each criterion was allotted an estimated weight with the help of AHP and FR. In the end, all the parameters were integrated using weighted overlay analysis in which the influence value of the drainage density was given the highest value. This gave the output in terms of five flood risk zones: very high risk, high risk, moderate risk, low risk, and very low risk. According to the results, 1168 km², that is, 8% of the total area, showed a very high risk of flood occurrence. Reshun, Mastuj, Booni, Colony, and some other villages were identified as high-risk zones of the study area, which have been drastically damaged many times by flash floods. This study is pioneering in its field and provides policy guidelines for risk managers, emergency and disaster response services, urban and infrastructure planners, hydrologists, and climate scientists.

Understanding the spatial growth of cities is crucial for proactive planning and sustainable urbanization. The largest and most densely inhabited megapolis of Pakistan, Karachi, has experienced massive spatial growth not only in the core areas of the city, but also in the city’s suburbs and outskirts over the past decades. In this study, the land use/land cover (LULC) in Karachi was classified using Landsat data and the random forest algorithm from the Google Earth Engine cloud platform for the years 1990, 2000, 2010, and 2020. Land use/land cover classification maps as well as an urban sprawl matrix technique were used to analyze the geographical patterns and trends of urban sprawl. Six urban classes, namely, the primary urban core, secondary urban core, sub-urban fringe, scatter settlement, urban open space, and non-urban area, were determined for the exploration of urban landscape changes. Future scenarios of LULC for 2030 were predicted using a CA–Markov model. The study found that the built-up area had expanded in a considerably unpredictable manner, primarily at the expense of agricultural land. The increase in mangroves and grassland and shrub land proved the effectiveness of afforestation programs in improving vegetation coverage in the study area. The investigation of urban landscape alteration revealed that the primary urban core expanded from the core districts, namely, the Central, South, and East districts, and a new urban secondary core emerged in Malir in 2020. The CA–Markov model showed that the total urban built-up area could potentially increase from 584.78 km² in 2020 to 652.59 km² in 2030. The integrated method combining remote sensing, GIS, and an urban sprawl matrix has proven invaluable for the investigation of urban sprawl in a rapidly growing city.

Environmental managers and policymakers increasingly discuss trade-offs between ecosystem services (ESs). However, few studies have used nonlinear models to provide scenario-specific land-use planning. This study determined the effects of different future land use/land cover (LULC) scenarios on ESs in the Yili River Valley, China, and analyzed the trade-offs and synergistic response characteristics. We simulated land-use changes in the Yili River Valley during 2020–2030 under three different scenarios using a patch-generating land-use simulation (PLUS) model—business as usual (BAU), economic development (ED), and ecological conservation (EC). Subsequently, we evaluated the water yield (WY), carbon storage (CS), soil retention (SR), and nutrient export (NE) ESs by combining the PLUS and integrated valuation of ecosystem services and trade-offs (InVEST) models, thus exploring multiple trade-offs among these four ESs at a regional scale. For the BAU scenario, there are some synergistic effects between WY and SR in the Yili River Valley, in addition to significant trade-off effects between CS and NE. For the ED scenario, the rapid expansion of cropland and constructed land is at the expense of forested grassland, leading to a significant decline in ESs. For the EC scenario, the model predicted that the cumulative regional net future carbon storage, cumulative water retention, and cumulative soil conservation would all increase due to ecological engineering and the revegetation of riparian zones and that formerly steep agricultural land can be effective in improving ESs. Meanwhile, the trade-off effect would be significantly weakened between CS and NE. These results can inform decision makers on specific sites where ecological engineering is implemented. Our findings can enhance stakeholders’ understanding of the interactions between ESs indicators in different scenarios.

Understanding the spatiotemporal patterns of urban heat islands and the factors that influence this phenomenon can help to alleviate the heat stress exacerbated by urban warming and strengthen heat-related urban resilience, thereby contributing to the achievement of the United Nations Sustainable Development Goals. The association between surface urban heat island (SUHI) effects and land use/land cover features has been studied extensively, but the situation in tropical cities is not well-understood due to the lack of consistent data. This study aimed to explore land use/land cover (LULC) changes and their impact on the urban thermal environment in a tropical megacity—Karachi, Pakistan. Land cover maps were produced, and the land surface temperature (LST) was estimated using Landsat images from five different years over the period 2000–2020. The surface urban heat island intensity (SUHII) was then quantified based on the LST data. Statistical analyses, including geographically weighted regression (GWR) and correlation analyses, were performed in order to analyze the relationship between the land cover composition and LST. The results indicated that the built-up area of Karachi increased from 97.6 km² to 325.33 km² during the period 2000–2020. Among the different land cover types, the areas classified as built-up or bare land exhibited the highest LST, and a change from vegetation to bare land led to an increase in LST. The correlation analysis indicated that the correlation coefficients between the normalized difference built-up index (NDBI) and LST ranged from 0.14 to 0.18 between 2000 and 2020 and that NDBI plays a dominant role in influencing the LST. The GWR analysis revealed the spatial variation in the association between the land cover composition and the SUHII. Parks with large areas of medium- and high-density vegetation play a significant role in regulating the thermal environment, whereas the scattered vegetation patches in the urban core do not have a significant relationship with the LST. These findings can be used to inform adaptive land use planning that aims to mitigate the effects of the UHI and aid efforts to achieve sustainable urban growth.

The objective of this study is to adopt a methodology for analysing spatial patterns of danger of forest fire at Margalla Hills, Islamabad, Pakistan. The work is concentrated on burnt areas using Landsat data and to classify forest fire severity with different parameters (climatic, vegetation, topography and human activities). In addition to these four variables, the extent of the burned areas was measured. Statistical analysis at each fire scene was used to measure the effect on the variables. To calculate the fire severity ratio correlated to each variable, logistic and stepwise regressions were used. The results showed that the burned areas have increased at a rate of 25.848 ha/day (R² = 0.98) if the number of total days since the start of fire has increased. As a result, forest density, distance to roads, average quarterly maximum temperature and average quarterly mean wind speed were highly correlated with the fire severity. Only average quarterly maximum temperature and forest density affected the size of the burnt areas. Prediction maps indicate that 53% of forests are in the very low severity level (0.25–0.45), 25% in the low level (0.45–0.65) and 22% in high and very high levels (>0.65).

Prescribed burning is a common strategy for minimizing forest fire risk. Fire is introduced under specific environmental conditions, with explicit duration, intensity, and rate of spread. Such conditions deviate from those encountered during the fire season. Prescribed burns mostly affect surface fuels and understory vegetation, an outcome markedly different when compared to wildfires. Data on prescribed burning are crucial for evaluating whether land management targets have been reached. This research developed a methodology to quantify the effects of prescribed burns using multi-temporal Sentinel-1 Synthetic Aperture Radar (SAR) imagery in the forests of southeastern Australia. C-band SAR datasets were specifically used to statistically explore changes in radar backscatter coefficients with the intensity of prescribed burns. Two modeling approaches based on pre- and post-fire ratios were applied for evaluating prescribed burn impacts. The effects of prescribed burns were documented with an overall accuracy of 82.3% using cross-polarized backscatter (VH) SAR data under dry conditions. The VV polarization indicated some potential to detect burned areas under wet conditions. The findings in this study indicate that the C-band SAR backscatter coefficient has the potential to evaluate the effectiveness of prescribed burns due to its sensitivity to changes in vegetation structure.