ABSTRACT In densely populated cities, people confront multiple impacts of climate change, including stormwater runoff, flooding, heat waves, and health issues. To address these impacts, nature-based solutions (NbS) have gained significant attention in recent years due to their potential to contribute to sustainability, resilience, and ecosystem services. Nevertheless, more research is required to understand the wider potential of NbS and to identify priority areas for their deployment. This study developed an analytical approach and implemented it in the Tehran Metropolitan Area (TMA) by considering factors related to NbS co-benefits. The approach included spatial prioritization, spatial correlation analysis, and suitability analysis. It combines a Geographic Information System, Integrated Valuation of Ecosystem Services and Trade-offs (InVEST) model, Zonation, Best-Worst Method, Bivariate Local Indicator of Spatial Autocorrelation analysis (BiLISA), a Fuzzy Inference System, and Boolean overlay analysis. The results include maps showing priority areas for different types of NbS and spatial co-benefits and mismatches. High-priority areas are located in disadvantaged areas. In total, 105.73 km2 (17% of the TMA) was designated as highly prioritized. We show that various NbS interventions are applicable, with green roofs and rainwater harvesting being the most feasible. This study can assist decision-makers in optimizing NbS to deliver maximum co-benefits.