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

Global Urban Fractional Changes at a 1 Km Resolution Throughout 2100 under Eight Scenarios of Shared SocioeconomicPathways (ssps) and Representative Concentration Pathways (Rcps)

Global urban fractional changes at a 1 km resolution throughout 2100 under eight SSP-RCP scenarios

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

Global Urban Growth Between 1870 and 2100 from Integrated High Resolution Mapped Data and Urban Dynamic Modeling

A Cellular Automata Downscaling Based 1 Km Global Land Use Datasets (2010–2100)

Modeling Gridded Urban Fractional Change Using the Temporal Context Information in the Urban Cellular Automata Model

Predicting Land Use Changes under Shared Socioeconomic Pathway–Representative Concentration Pathway Scenarios to Support Sustainable Planning in High-Density Urban Areas: A Case Study of Hangzhou, Southeastern China

Projecting global urban area growth through 2100 based on historical time-series data and future Shared Socioeconomic Pathways

Simulation of urban land expansion in China at 30 m resolution through 2050 under shared socioeconomic pathways

Unveiling the Impacts of Climate Change and Human Activities on Land-Use Evolution in Ecologically Fragile Urbanizing Areas: A Case Study of China’s Central Plains Urban Agglomeration

Mapping the spatial heterogeneity of global land use and land cover from 2020 to 2100 at a 1 km resolution

1 Km Land Use/land Cover Change of China under Comprehensive Socioeconomic and Climate Scenarios for 2020–2100

Simulating the Dynamics of Urban Land Quantity in China from 2020 to 2070 under the Shared Socioeconomic Pathways

High-spatiotemporal-resolution Mapping of Global Urban Change from 1985 to 2015

Projecting LUCC dynamics and ecosystem services in an emerging urban agglomeration under SSP-RCP scenarios and their management implications

Global urbanization projections for the Shared Socioeconomic Pathways

Large Disagreements in Estimates of Urban Land Across Scales and Their Implications

An online tool with Google Earth Engine and cellular automata for seamlessly simulating global urban expansion at high resolutions

Simulation of Land-Use Spatiotemporal Changes under Ecological Quality Constraints: The Case of the Wuhan Urban Agglomeration Area, China, over 2020-2030

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