Morocco’s Témara Plain relies heavily on its aquifer system as a critical resource for drinking water, irrigation, and industrial activities. However, this essential groundwater reserve is increasingly threatened by over-extraction, seawater intrusion, and complex hydrogeochemical processes driven by the region’s geological characteristics and anthropogenic pressures. This study aims to assess groundwater quality and its vulnerability to pollution risks and map the spatial distribution of key hydrochemical processes through an integrated approach combining Geographic Information System (GIS) techniques and multivariate statistical analysis, as well as applying the DRASTIC model to evaluate water vulnerability. A total of fifty-eight groundwater samples were collected across the plain and analyzed for major ions to identify dominant hydrochemical facies. Spatial interpolation using Inverse Distance Weighting (IDW) within GIS revealed distinct patterns of sodium chloride (Na-Cl) facies near the coastal areas with chloride concentrations exceeding the World Health Organization (WHO) drinking water guideline of 250 mg/L—indicative of seawater intrusion. In addition to marine intrusion, agricultural pollution constitutes a major diffuse pressure across the aquifer. Shallow groundwater zones in agricultural areas show heightened vulnerability to salinization and nitrate contamination, with nitrate concentrations reaching up to 152.3 mg/L, far surpassing the WHO limit of 45 mg/L. Furthermore, other anthropogenic pollution sources—such as wastewater discharges from septic tanks in peri-urban zones lacking proper sanitation infrastructure and potential leachate infiltration from informal waste disposal sites—intensify stress on the aquifer. Principal Component Analysis (PCA) identified three key factors influencing groundwater quality: natural mineralization due to carbonate rock dissolution, agricultural inputs, and salinization driven by seawater intrusion. Additionally, The DRASTIC model was used within the GIS environment to create a vulnerability map based on seven key parameters. The map revealed that low-lying coastal areas are most vulnerable to contamination.

Rapid urbanization and climate change pose significant challenges to air quality in arid metropolitan areas, with critical implications for public health and sustainable development. This study projects the evolution of air pollution in Riyadh, Saudi Arabia, through 2070 using an integrated modeling approach that combines CMIP6 climate projections with localized air quality data. We analyzed daily concentrations of major pollutants (SO₂, NO₂) across 15 strategically selected monitoring stations representing diverse urban environments, including traffic corridors, residential areas, healthcare facilities, and semi-natural zones. Climate data from two Earth System Models (CNRM-ESM2-1 and MPI-ESM1.2) were bias-corrected and integrated with historical pollution measurements (2000–2015) using hierarchical Bayesian statistical modeling under SSP2-4.5 and SSP5-8.5 emission scenarios. Our results revealed substantial deterioration in air quality, with projected increases of 80–130% for SO₂ and 45–55% for NO₂ concentrations by 2070 under high-emission scenarios. Spatial analysis demonstrated pronounced pollution gradients, with traffic corridors (Eastern Ring Road, Northern Ring Road, Southern Ring Road) and densely urbanized areas (King Fahad Road, Makkah Road) experiencing the most severe increases, exceeding WHO guidelines by factors of 2–3. Even semi-natural areas showed significant increases in pollution due to regional transport effects. The hierarchical Bayesian framework effectively quantified uncertainties while revealing consistent degradation trends across both climate models, with the MPI-ESM1.2 model showing a greater sensitivity to anthropogenic forcing. Future concentrations are projected to reach up to 70 μg m⁻³ for SO₂ and exceed 100 μg m⁻³ for NO₂ in heavily trafficked areas by 2070, representing 2–3 times the Traffic corridors showed concentration increases of 21–24% compared to historical baselines, with some stations (R5, R13, and R14) recording projected levels above 4.0 ppb for SO₂ under the SSP5-8.5 scenario. These findings highlight the urgent need for comprehensive emission reduction strategies, accelerated renewable energy transition, and reformed urban planning approaches in rapidly developing arid cities.

Vegetation degradation in arid and semi-arid regions is intensifying due to rising temperatures, declining rainfall, soil exposure, and persistent human pressures. Drylands cover over 41% of the global land surface and support nearly two billion people, making their degradation a major environmental and socio-economic concern. However, many remote sensing and GIS-based assessment approaches remain fragmented and difficult to reproduce. This study proposes a Computational Sustainability Framework for vegetation degradation assessment that integrates multi-source satellite data, biophysical indicators, automated geospatial preprocessing, and the Analytical Hierarchy Process (AHP) within a transparent and reproducible workflow. The framework comprises four phases: data preprocessing, indicator extraction and normalization, AHP-based modeling, and spatial classification with qualitative validation. The framework was applied to the Al-Khunfah and Harrat al-Harrah Protected Areas in northern Saudi Arabia using multi-source datasets for the January–April 2023 period, including Sentinel-2, Landsat-8, CHIRPS precipitation, ESA-CCI land cover, FAO soil data, and SRTM DEM. High degradation zones were associated with low NDVI (0.276), and elevated LST (>49 °C), whereas low degradation areas were concentrated near wadis and relatively more fertile soils. Overall, the proposed framework provides a scalable and interpretable tool for early-stage vegetation degradation screening in arid environments, supporting the prioritization of areas for ecological investigation and restoration planning.

This paper proposes a Computational Sustainability Framework for assessing vegetation degradation and desertification in arid ecosystems. The framework integrates remote-sensing data, GIS modeling, and multi-criteria decision algorithms to analyze environmental conditions at scale. It consists of four computational phases: (1) data acquisition and preprocessing, (2) feature extraction and normalization, (3) multi-criteria decision modeling using the Analytical Hierarchy Process (AHP), and (4) visualization and validation. The framework was applied to two Saudi protected areas, Al Khunfah and Ḩorrah al-Ḥorra, using Sentinel-2, Landsat-8, and ancillary datasets. Results show that high degradation zones correspond to high BSI (>0.27) and low NDVI (<0.09) with surface temperatures exceeding 48 °C. The AHP model achieved a consistency ratio of 0.0957, confirming reliability, and identified three suitability classes covering approximately 4 000 km² (high), 26 000 km² (moderate), and 6 000 km² (low). The study demonstrates how computational integration of spatial data and decision algorithms can support large-scale environmental monitoring and restoration planning.