Remote Sensing-Based Spatiotemporal Assessment of Heat Risk in the Guangdong–Hong Kong–Macao Greater Bay Area

Under the dual pressures of climate change and rapid urbanization, extreme heat events pose growing risks to densely populated megaregions. The Guangdong–Hong Kong–Macao Greater Bay Area (GBA), a densely populated and economically vital region, serves as a critical hotspot for heat risk aggregation. This study develops a high-resolution multi-dimensional framework to assess the spatiotemporal evolution of its heat risk profile from 2000 to 2020. A Heat Risk Index (HRI) integrating heat hazard and vulnerability components to measure potential heat-related impacts is calculated as the product of the Heat Hazard Index (HHI) and Heat Vulnerability Index (HVI) for 1 km grids in GBA. The HHI integrates the frequency of hot days and hot nights. HVI incorporates population density, GDP, remote-sensing nighttime light data, and MODIS-based landscape indicators (e.g., NDVI, NDWI, and NDBI), with weights determined objectively using the static Entropy Weight Method to ensure spatiotemporal comparability. The findings reveal an escalation of heat risk, expanding at an average rate of 342 km2 per year (p = 0.008), with the proportion of areas classified as high-risk or above increasing from 21.8% in 2000 to 33.3% in 2020. This trend was characterized by (a) a pronounced asymmetric warming pattern, with nighttime temperatures rising more rapidly than daytime temperatures; (b) high vulnerability dominated by the concentration of population and economic assets, as indicated by high EWM-based weights; and (c) isolated high-risk hotspots (Guangzhou and Hong Kong) in 2000, which have expanded into a high-risk belt across the Pearl River Delta’s industrial heartland, like Foshan seeing their high-risk area expand from 3.4% to 27.0%. By combining remote sensing and socioeconomic data, this study provides a transferable framework that moves beyond coarse-scale assessments to identify specific intra-regional risk hotspots. The resulting high-resolution risk maps offer a quantitative foundation for developing spatially explicit climate adaptation strategies in the GBA and other rapidly urbanizing megaregions.