This research evaluates the climatic summertime representation of the diurnal cycle of near-surface temperature using the Weather Research and Forecasting System (WRF) over the rapidly urbanizing and water-vulnerable Phoenix metropolitan area. A suite of monthly, high-resolution (2 km grid spacing) simulations are conducted during the month of July with both a contemporary landscape and a hypothetical presettlement scenario. WRF demonstrates excellent agreement in the representation of the daily to monthly diurnal cycle of near-surface temperatures, including the accurate simulation of maximum daytime temperature timing. Thermal sensitivity to anthropogenic land use and land cover change (LULCC), assessed via replacement of the modern-day landscape with natural shrubland, is small on the regional scale. The WRF-simulated characterization of the diurnal cycle, supported by previous observational analyses, illustrates two distinct and opposing impacts on the urbanized diurnal cycle of the Phoenix metro area, with evening and nighttime warming partially offset by daytime cooling. The simulated nighttime urban heat island (UHI) over this semiarid urban complex is explained by well-known mechanisms (slow release of heat from within the urban fabric stored during daytime and increased emission of longwave radiation from the urban canopy toward the surface). During daylight hours, the limited vegetation and dry semidesert region surrounding metro Phoenix warms at greater rates than the urban complex. Although prior work has suggested that daytime temperatures are lower within the urban complex owing to the addition of residential and agricultural irrigation (i.e., “oasis effect”) we show that modification of Phoenix's surrounding environment to a biome more representative of temperate regions eliminates the daytime urban cooling. Our results indicate that surrounding environmental conditions, including land cover and availability of soil moisture, play a principal role in establishing the nature and evolution of the diurnal cycle of near-surface temperature for the greater Phoenix, Arizona, metropolitan area relative to its rural and undeveloped counterpart.

Evaluation of built environment energy demand is necessary in light of global projections of urban expansion. Of particular concern are rapidly expanding urban areas in environments where consumption requirements for cooling are excessive. Here, we simulate urban air conditioning (AC) electric consumption for several extreme heat events during summertime over a semiarid metropolitan area with the Weather Research and Forecasting (WRF) model coupled to a multilayer building energy scheme. Observed total load values obtained from an electric utility company were split into two parts, one linked to meteorology (i.e., AC consumption) which was compared to WRF simulations, and another to human behavior. WRF-simulated non-dimensional AC consumption profiles compared favorably to diurnal observations in terms of both amplitude and timing. The hourly ratio of AC to total electricity consumption accounted for ~53% of diurnally averaged total electric demand, ranging from ~35% during early morning to ~65% during evening hours. Our work highlights the importance of modeling AC electricity consumption and its role for the sustainable planning of future urban energy needs. Finally, the methodology presented in this article establishes a new energy consumption-modeling framework that can be applied to any urban environment where the use of AC systems is prevalent.

Maricopa County, Arizona, anchor to the fastest growing megapolitan area in the United States, is located in a hot desert climate where extreme temperatures are associated with elevated risk of mortality. Continued urbanization in the region will impact atmospheric temperatures and, as a result, potentially affect human health. We aimed to quantify the number of excess deaths attributable to heat in Maricopa County based on three future urbanization and adaptation scenarios and multiple exposure variables.

Two scenarios (low and high growth projections) represent the maximum possible uncertainty range associated with urbanization in central Arizona, and a third represents the adaptation of high-albedo cool roof technology. Using a Poisson regression model, we related temperature to mortality using data spanning 1983–2007. Regional climate model simulations based on 2050-projected urbanization scenarios for Maricopa County generated distributions of temperature change, and from these predicted changes future excess heat-related mortality was estimated. Subject to urbanization scenario and exposure variable utilized, projections of heat-related mortality ranged from a decrease of 46 deaths per year (− 95%) to an increase of 339 deaths per year (+ 359%).

Projections based on minimum temperature showed the greatest increase for all expansion and adaptation scenarios and were substantially higher than those for daily mean temperature. Projections based on maximum temperature were largely associated with declining mortality. Low-growth and adaptation scenarios led to the smallest increase in predicted heat-related mortality based on mean temperature projections. Use of only one exposure variable to project future heat-related deaths may therefore be misrepresentative in terms of direction of change and magnitude of effects. Because urbanization-induced impacts can vary across the diurnal cycle, projections of heat-related health outcomes that do not consider place-based, time-varying urban heat island effects are neglecting essential elements for policy relevant decision-making.

In this study, WRF-Chem is utilized at high resolution (1.333 km grid spacing for the innermost domain) to investigate impacts of southern California anthropogenic emissions (SoCal) on Phoenix ground-level ozone concentrations ([O³]) for a pair of recent exceedance episodes. First, WRF-Chem control simulations, based on the US Environmental Protection Agency (EPA) 2005 National Emissions Inventories (NEI05), are conducted to evaluate model performance. Compared with surface observations of hourly ozone, CO, NO^X, and wind fields, the control simulations reproduce observed variability well. Simulated [O³] are comparable with the previous studies in this region. Next, the relative contribution of SoCal and Arizona local anthropogenic emissions (AZ) to ozone exceedances within the Phoenix metropolitan area is investigated via a trio of sensitivity simulations: (1) SoCal emissions are excluded, with all other emissions as in Control; (2) AZ emissions are excluded with all other emissions as in Control; and (3) SoCal and AZ emissions are excluded (i.e., all anthropogenic emissions are eliminated) to account only for Biogenic emissions and lateral boundary inflow (BILB). Based on the USEPA NEI05, results for the selected events indicate the impacts of AZ emissions are dominant on daily maximum 8 h average (DMA8) [O³] in Phoenix. SoCal contributions to DMA8 [O³] for the Phoenix metropolitan area range from a few ppbv to over 30 ppbv (10–30 % relative to Control experiments). [O³] from SoCal and AZ emissions exhibit the expected diurnal characteristics that are determined by physical and photochemical processes, while BILB contributions to DMA8 [O³] in Phoenix also play a key role.

Physical mechanisms of incongruency between observations and Weather Research and Forecasting (WRF) Model predictions are examined. Limitations of evaluation are constrained by (i) parameterizations of model physics, (ii) parameterizations of input data, (iii) model resolution, and (iv) flux observation resolution. Observations from a new 22.1-m flux tower situated within a residential neighborhood in Phoenix, Arizona, are utilized to evaluate the ability of the urbanized WRF to resolve finescale surface energy balance (SEB) when using the urban classes derived from the 30-m-resolution National Land Cover Database. Modeled SEB response to a large seasonal variation of net radiation forcing was tested during synoptically quiescent periods of high pressure in winter 2011 and premonsoon summer 2012. Results are presented from simulations employing five nested domains down to 333-m horizontal resolution. A comparative analysis of model cases testing parameterization of physical processes was done using four configurations of urban parameterization for the bulk urban scheme versus three representations with the Urban Canopy Model (UCM) scheme, and also for two types of planetary boundary layer parameterization: the local Mellor–Yamada–Janjić scheme and the nonlocal Yonsei University scheme. Diurnal variation in SEB constituent fluxes is examined in relation to surface-layer stability and modeled diagnostic variables. Improvement is found when adapting UCM for Phoenix with reduced errors in the SEB components. Finer model resolution is seen to have insignificant (<1 standard deviation) influence on mean absolute percent difference of 30-min diurnal mean SEB terms.

The effects of urbanization on ozone levels have been widely investigated over cities primarily located in temperate and/or humid regions. In this study, nested WRF-Chem simulations with a finest grid resolution of 1 km are conducted to investigate ozone concentrations O₃ due to urbanization within cities in arid/semi-arid environments. First, a method based on a shape preserving Monotonic Cubic Interpolation (MCI) is developed and used to downscale anthropogenic emissions from the 4 km resolution 2005 National Emissions Inventory (NEI05) to the finest model resolution of 1 km. Using the rapidly expanding Phoenix metropolitan region as the area of focus, we demonstrate the proposed MCI method achieves ozone simulation results with appreciably improved correspondence to observations relative to the default interpolation method of the WRF-Chem system. Next, two additional sets of experiments are conducted, with the recommended MCI approach, to examine impacts of urbanization on ozone production: (1) the urban land cover is included (i.e., urbanization experiments) and, (2) the urban land cover is replaced with the region's native shrubland. Impacts due to the presence of the built environment on O₃ are highly heterogeneous across the metropolitan area. Increased near surface O₃ due to urbanization of 10–20 ppb is predominantly a nighttime phenomenon while simulated impacts during daytime are negligible. Urbanization narrows the daily O₃ range (by virtue of increasing nighttime minima), an impact largely due to the region's urban heat island. Our results demonstrate the importance of the MCI method for accurate representation of the diurnal profile of ozone, and highlight its utility for high-resolution air quality simulations for urban areas.