Abstract:The impacts of urbanization on terrestrial thermal and hydrological changes should be considered for sustainable development of urban environment and economy. However, previous studies mainly focused on the impacts at short time scales, which is insufficient to provide the needed information. This study augmented the high‐resolution Conjunctive Surface‐Subsurface Process version 2 (CSSPv2) land model with an urban canopy model that incorporates domestic wastewater recharge (DWWR) and dynamical impervious surface area (ISA) to develop a new CSSPv2_Urban model and conducted 38‐year simulations over the Beijing metropolitan area at 1 km resolution to quantify the urbanization impacts. Compared with CSSPv2, CSSPv2_Urban decreases the root‐mean‐square errors in surface latent/sensible heat fluxes, near‐surface air temperature and land surface temperature by 15%–66%, increases Kling‐Gupta efficiency for streamflow by 27%–34%, and reproduces the observed long‐term trends of land surface temperature and streamflow. It also reduces the corresponding errors in the modern reanalysis products by up to 70%. The model sensitivity analyses showed that the urbanization over Beijing area contributed to 31% of the surface warming during 1979–2017 and kept the warming trend over the new development zone during 2000–2013 when global warming slowed down. Although the large fraction of ISA in the city center increased streamflow by generating more surface runoff, 39%–62% of the increasing trend of streamflow was caused by the increased DWWR. Our study highlights the importance of the urban heat island effect from a climatological perspective and suggests that the influences of urbanization on streamflow may be underestimated if ignoring DWWR. Fast urbanization since last century has been modifying the terrestrial energy and water cycle significantly. However, quantifying urbanization's impacts on long‐term changes of terrestrial thermal and hydrological conditions at high spatial resolution, which is important for providing guidelines for urban sustainable development, is still in progress. Here, we added an urban canopy model that incorporates domestic wastewater recharge (DWWR) and dynamical impervious surface area (ISA) into the Conjunctive Surface–Subsurface Process version 2 (CSSPv2) land model and conducted high‐resolution (1 km) and long‐term (1979–2017) simulations over the Beijing metropolitan area to investigate the urbanization impacts. The newly developed CSSPv2_Urban model is superior to the CSSPv2 and modern reanalysis data sets in simulating the thermal and hydrological variables. Sensitivity analysis revealed that urbanization contributed to about 31% of the surface warming trend in Beijing and sustained the warming during 2000–2013 when the global warming slowed down. Moreover, increased ISA and DWWR during the urbanization were found to interpret the increasing streamflow trend by 30% and 51%, respectively. Thus, more attentions should be paid to the DWWR which is usually ignored in current urban models to better capture the hydrological influences of urbanization. CSSPv2_Urban was developed with urban energy model, dynamic impervious surface area (ISA), and domestic wastewater recharge (DWWR) scheme CSSPv2_Urban 1 km high‐resolution simulation outperformed modern reanalysis for energy and water cycle over Beijing metropolitan area Simulation suggests that ISA and DWWR contributed to 30% and 51% of the increasing streamflow trend during 1979–2017 near Beijing city center CSSPv2_Urban was developed with urban energy model, dynamic impervious surface area (ISA), and domestic wastewater recharge (DWWR) scheme CSSPv2_Urban 1 km high‐resolution simulation outperformed modern reanalysis for energy and water cycle over Beijing metropolitan area Simulation suggests that ISA and DWWR contributed to 30% and 51% of the increasing streamflow trend during 1979–2017 near Beijing city center

High‐Resolution Land Surface Modeling of the Effect of Long‐Term Urbanization on Hydrothermal Changes Over Beijing Metropolitan Area

Quantifying Changes in Water Use and Groundwater Availability in a Megacity Using Novel Integrated Systems Modeling

Simulation of the Urban Climate in A Chinese Megacity with Spatially Heterogeneous Anthropogenic Heat Data

Urbanization Effect on Spatiotemporal Thermal Patterns and Changes in Hangzhou (China)

A Long-term Simulation of Land Surface Conditions at High-Resolution Over Continental China

Combined Effects of Impervious Surface Change and Large-Scale Afforestation on the Surface Urban Heat Island Intensity of Beijing, China Based on Remote Sensing Analysis

Characterizing the spatiotemporal response of runoff to impervious surface dynamics across three highly urbanized cities in southern China from 2000 to 2017

Assessment of Regional Climate Effects of Urbanization around Subtropical City Wuhan in Summer Using Numerical Modeling

Assessment of the impact of climate change on urban flooding: a case study of Beijing, China

Land surface model CAS‐LSM: Model description and evaluation

Quantifying the Impacts of High-Resolution Urban Information on the Urban Thermal Environment

Revealing the Response of Evapotranspiration and Water Use Efficiency to Vegetation Changes in Different Types of Urban Green Spaces

The Impact of Urbanization-Induced Land Use Change on Land Surface Temperature

The Spatiotemporal Implications of Urbanization for Urban Heat Islands in Beijing: A Predictive Approach Based on CA–Markov Modeling (2004–2050)

Impacts of Water Layout on the Atmospheric Environment in Urban Areas

High‐Resolution Land Surface Modeling of Hydrological Changes over the Sanjiangyuan Region in the Eastern Tibetan Plateau: 1. Model Development and Evaluation

Global urban fractional changes at a 1 km resolution throughout 2100 under eight scenarios of Shared Socioeconomic Pathways (SSPs) and Representative Concentration Pathways (RCPs)

Hydrological Response to Urbanization at Different Spatio-Temporal Scales Simulated by Coupling of CLUE-S and the SWAT Model in the Yangtze River Delta Region

Impact of Urbanization-Driven Land Use Changes on Runoff in the Upstream Mountainous Basin of Baiyangdian, China: A Multi-Scenario Simulation Study

Urban water system theory and its model development and application

Regulation characteristics of underlying surface on runoff regime metrics and their spatial differences in typical urban communities across China