Urban land–atmosphere interactions can be captured by numerical modeling framework with coupled land surface and atmospheric processes, while the model performance depends largely on accurate input parameters. In this study, we use an advanced stochastic approach to quantify parameter uncertainty and model sensitivity of a coupled numerical framework for urban land–atmosphere interactions. It is found that the development of urban boundary layer is highly sensitive to surface characteristics of built terrains. Changes of both urban land use and geometry impose significant impact on the overlying urban boundary layer dynamics through modification on bottom boundary conditions, i.e., by altering surface energy partitioning and surface aerodynamic resistance, respectively. Hydrothermal properties of conventional and green roofs have different impacts on atmospheric dynamics due to different surface energy partitioning mechanisms. Urban geometry (represented by the canyon aspect ratio), however, has a significant nonlinear impact on boundary layer structure and temperature. Besides, managing rooftop roughness provides an alternative option to change the boundary layer thermal state through modification of the vertical turbulent transport. The sensitivity analysis deepens our insight into the fundamental physics of urban land–atmosphere interactions and provides useful guidance for urban planning under challenges of changing climate and continuous global urbanization.

This study examines the spatial and temporal patterns of the surface urban heat island (SUHI) intensity in the Phoenix metropolitan area and the relationship with land use land cover (LULC) change between 2000 and 2014. The objective is to identify specific regions in Phoenix that have been increasingly heated and cooled to further understand how LULC change influences the SUHI intensity. The data employed include MODerate-resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) 8-day composite June imagery, and classified LULC maps generated using 2000 and 2014 Landsat imagery. Results show that the regions that experienced the most significant LST changes during the study period are primarily on the outskirts of the Phoenix metropolitan area for both daytime and nighttime. The conversion to urban, residential, and impervious surfaces from all other LULC types has been identified as the primary cause of the UHI effect in Phoenix. Vegetation cover has been shown to significantly lower LST for both daytime and nighttime due to its strong cooling effect by producing more latent heat flux and less sensible heat flux. We suggest that urban planners, decision-makers, and city managers formulate new policies and regulations that encourage residential, commercial, and industrial developers to include more vegetation when planning new construction.