Divergent Sensitivity of Gross Primary Productivity to Compound Drought and Heatwaves Across China’s Three Major Urban Agglomerations

Compound Drought and Heatwave (CDH) events increasingly threaten terrestrial carbon uptake, yet the spatiotemporal heterogeneity of Gross Primary Productivity (GPP) responses in urban agglomerations remains unclear. This study analyzed CDH impacts in China’s three major urban agglomerations, namely the Beijing–Tianjin–Hebei (BTH), Yangtze River Delta (YRD), and Pearl River Delta (PRD) regions, using ERA5 and satellite GPP data (GOSIF and FluxSat) for representative CDH years (2007 for BTH; 2022 for YRD and PRD). CDH conditions exhibited a coherent hot–dry coupling, with temperature anomalies of 0.46–1.26 K and soil moisture deficits of −0.042 to −0.169 m3 m−3, accompanied by enhanced atmospheric dryness. Pronounced spatial heterogeneity in GPP responses aligned with regional climatic regimes and ecosystem types. The water-limited BTH region exhibited significant GPP deficits, with anomalies of −1.13 Standard Deviations (STD) and −0.96 STD for GPPFluxSat and GPPGOSIF, respectively. Conversely, the energy-limited regions showed positive anomalies: the YRD recorded +0.32 and +1.79 STD, while the PRD reached +1.86 and +1.06 STD for GPPFluxSat and GPPGOSIF, respectively. Mechanistically, the north–south contrast suggests a transition from water-limited vulnerability to energy-limited resilience, with vegetation traits and management (e.g., potential irrigation buffering in croplands and deeper water access in forests) modulating sensitivity to atmospheric dryness. These findings provide quantitative benchmarks for improving regional carbon-cycle assessments and adaptation planning under increasing compound extremes.