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

Simulation model for urban water dissipation

Dissipation of Water in Urban Area, Mechanism and Modelling with the Consideration of Anthropogenic Impacts: A Case Study in Xiamen

Prediction of Water Demand and Pollution Pressure of Urban Households in China

Simulation of Water Dissipation in Urban Buildings

How Does Urban Form Explain Water Cooling Effect Heterogeneity: A Case Study in Wuhan, China

Water Dissipation Mechanism of Residential and Office Buildings in Urban Areas

Estimating Irrigation Water Demand for Green Spaces in Humid Areas: Seeking a Sustainable Water Management Strategy

Urban water system theory and its model development and application

Using Nighttime Light Data to Estimate Water Evaporation Inside Buildings in China's Urban Areas

Simulation of evapotranspiration in urban areas based on URMOD and validation with remote sensing data

Simulation and Analysis on Mechanisms of Urban Households Water Use

Impacts of Water Layout on the Atmospheric Environment in Urban Areas

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

Mechanisms and Empirical Modeling of Evaporation from Hardened Surfaces in Urban Areas.

Principles of urban hydrological model

Revealing the Response of Urban Heat Island Effect to Water Body Evaporation from Main Urban and Suburb Areas

Hydrologically Enhanced Distributed Urban Drainage Model and Its Application in Beijing City

A dynamic water resources management approach in Beijing using system dynamics model

Evapotranspiration Estimation Considering Anthropogenic Heat Based on Remote Sensing in Urban Area

Statistical Interpretation of the Daily Variation of Urban Water Consumption in Beijing, China

Assessing the Impact of Urbanization on Urban Evapotranspiration and Its Components Using a Novel Four-Source Energy Balance Model