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

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

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

Role of Satellite and Reanalysis Precipitation Products in Streamflow and Sediment Modeling over a Typical Alpine and Gorge Region in Southwest China.

Effects of meteorological forcings and land surface model on soil moisture simulation over China

On the Reliability of 12 High-Resolution Precipitation Products for Process-Based Hydrological Modeling in China

The Evaluation of Snow Depth Simulated by Different Land Surface Models in China Based on Station Observations

Elucidating Dominant Factors Affecting Land Surface Hydrological Simulations of the Community Land Model over China

Insight into Historical and Future Spring Snow Cover from Satellite Observation and Model Simulations over the Northern Hemisphere

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

Land surface model CAS‐LSM: Model description and evaluation

Long-term Reconstruction of Satellite-Based Precipitation, Soil Moisture, and Snow Water Equivalent in China

Validation of Remotely Sensed Estimates of Snow Water Equivalent Using Multiple Reference Datasets from the Middle and High Latitudes of China

Assessment of the Performance of Three Dynamical Climate Downscaling Methods Using Different Land Surface Information over China

Investigating Spatial And Temporal Variation Of Hydrological Processes In Western China Driven By Cmads

Coupled Snow and Frozen Ground Physics Improves Cold Region Hydrological Simulations: An Evaluation at the upper Yangtze River Basin (Tibetan Plateau)

Enhancement of land surface information and its impact on atmospheric modeling in the Heihe River Basin, northwest China

Development of A Land Surface Model with Coupled Snow and Frozen Soil Physics

Application and Evaluation of the China Meteorological Assimilation Driving Datasets for the SWAT Model (CMADS) in Poorly Gauged Regions in Western China

Synergistic Effects of High‐Resolution Factors for Improving Soil Moisture Simulations Over China

Evaluation of Ten Fresh Snow Density Parameterization Schemes for Simulating Snow Depth and Surface Energy Fluxes on the Eastern Tibetan Plateau

Evaluating the Performance of BCC-CSM2-MR Model in Simulating the Land Surface Processes in China