To enhance the capability of models in better characterizing the urban water cycle, physical parameterizations of urban hydrological processes have been implemented into the single-layer urban canopy model in the widely used Weather Research and Forecasting (WRF) Model. While the new model has been evaluated offline against field measurements at various cities, its performance in online settings via coupling to atmospheric dynamics requires further examination. In this study, the impact of urban hydrological processes on regional hydrometeorology of the fully integrated WRF–urban modeling system for two major cities in the United States, namely, Phoenix and Houston, is assessed. Results show that including hydrological processes improves prediction of the 2-m dewpoint temperature, an indicative measure of coupled thermal and hydrological processes. The implementation of green roof systems as an urban mitigation strategy is then tested at the annual scale. The reduction of environmental temperature and increase of humidity by green roofs indicate strong diurnal and seasonal variations and are significantly affected by geographical and climatic conditions. Comparison with offline studies reveals that land–atmosphere interactions play a crucial role in determining the effect of green roofs.
To meet the demand of a rapidly growing global population, urban areas have expanded considerably in recent decades (Seto et al. 2011). Through the modification of land surface energy and moisture balances, urbanization has determinant impacts on local and regional hydroclimates, leading to elevated temperature (Arnfield 2003; Yang et al. 2015a), reduced humidity (Unkašević et al. 2001), and change in precipitation patterns (Georgescu et al. 2012). To capture urban land–atmosphere interactions, numerous mesoscale atmosphere–urban modeling systems have been developed in the last decade (Best 2005; Chen et al. 2011; Martilli et al. 2002), among which the Weather Research and Forecasting (WRF) Model–urban system has been widely utilized and examined for major metropolitan regions around the world (Kusaka et al. 2012; Lin et al. 2008; Miao and Chen 2008).
The WRF Model includes several urban parameterization options, including the single-layer urban canopy model (SLUCM), which has been extensively studied and renders satisfactory capability in resolving urban land surface processes with a moderate requirement of input parameter space (Kusaka et al. 2001; Salamanca et al. 2011; Wang et al. 2011a,b). However, because of the oversimplified representation of urban hydrological processes, existing coupled atmosphere–urban modeling systems, including the WRF–SLUCM, have the least capacity in modeling the latent heat flux compared to other fluxes (Grimmond et al. 2010, 2011; Wang et al. 2013). The inadequate urban hydrological modeling consequently introduces errors into the atmospheric system via the lower boundary condition, impairing the reliability and accuracy of modeling land–atmosphere interactions (Song and Wang 2015a).
The importance of realistic representation of hydrological processes in the urban canopy model has long been recognized. Aiming at improving prediction of turbulent heat fluxes, Miao and Chen (2014) incorporated urban irrigation, oasis effect, and anthropogenic latent heat into the SLUCM. Evaluation against meteorological observations in Beijing showed that modeled latent heat flux was evidently improved. Inspired by this work and a state-of-the-art urban canopy model (Wang et al. 2011b, 2013), Yang et al. (2015a) further implemented physically based parameterization of evaporation over impervious surfaces and incorporated a green roof system into the SLUCM. Significant enhancement in prediction of latent heat flux was found for the four metropolitan areas investigated: Beijing (China), Vancouver (Canada), Phoenix (United States), and Montreal (Canada). However, tested in an offline (i.e., not dynamically coupled to the overlying atmosphere) setting, this study did not consider the interaction between the land and atmospheric system, and the potential omission of important feedback mechanisms can lead to significant uncertainty and potential errors in model results (Brubaker and Entekhabi 1996).
Using the online WRF–urban modeling system, Georgescu et al. (2011) studied the diurnal cycle of near-surface temperature over the urbanizing Phoenix metropolitan area, concluding that urban irrigation did not have a significant impact on near-surface temperatures. Vahmani and Hogue (2014) assessed the impact of irrigation over the Los Angeles metropolitan area by incorporating an irrigation module into the SLUCM. With irrigation, large biases in prediction of evapotranspiration from urban terrain were reduced. Nevertheless, other urban hydrological processes were neglected in their study, such as anthropogenic latent heat and evaporation from water-holding engineered pavements. It is also noteworthy that the aforementioned studies focused on turbulent heat fluxes, whereas the impact of urban hydrological processes on meteorological variables was rarely quantified. Though turbulent heat fluxes are closely related to temperature and humidity of the atmosphere, the interactions among them are complex with a variety of surface and meteorological conditions (Wang 2014a). Online studies that directly address the effect of hydrological processes on urban meteorology by coupling with the energy balance are still lacking in the literature.
One important function of the mesoscale atmosphere–urban modeling system is to evaluate potential strategies for sustainable cities. Urban areas, owing to impacts of global climate change, are projected to experience more frequent occurrences of climatic extremes in the future [e.g., heat wave (Meehl and Tebaldi 2004) and strong precipitation (Kripalani et al. 2007)], increasing the need for sustainable adaptation/mitigation strategies in areas where the majority of the globe’s inhabitants reside (IPCC 2012; Georgescu et al. 2014). Green (vegetated) roofs have significant potential and have been adopted in many cities (e.g., Chicago) to alleviate urban-induced heat stresses. The adoption of green roofs has been shown to reduce near-surface temperature (Yang and Wang 2014, 2015; Georgescu 2015), improve stormwater management (VanWoert et al. 2005; Carter and Jackson 2007), and enhance air quality (Yang et al. 2008; Rowe 2011). Although researches on meteorological impacts of green roofs are becoming increasingly widespread, most investigations have focused on building-resolving scales and explored with offline models where meteorological forcing is provided as boundary conditions (Sailor 2008; Sun et al. 2013). Upscaling the results of these studies for preparing guidance on green roof implementation for a city or at regional scales remains challenging because of the substantial influence of surface heterogeneity and land–atmosphere interactions (Ramamurthy et al. 2014; Song and Wang 2015b).
Consequently, there have been only a handful of studies investigating climatic (Georgescu et al. 2014; Georgescu 2015) and meteorological (Li et al. 2014) impacts of green roofs in a coupled atmosphere–urban modeling framework. Georgescu et al. (2014) explored the benefits of green roofs relative to highly reflective roofs and the potential to offset urban-induced warming at seasonal and annual time scales across the contiguous United States. Although impacts on near-surface temperature were less for green roofs relative to reflective roofs, a considerably reduced hydroclimatological trade-off was simulated for several regions via deployment of vegetated roofs. However, assuming green roofs were infinitely evaporating without water constraint, their results represented the maximum potential benefits of evaporating rooftop water pools rather than green roofs. Li et al. (2014) also compared the effectiveness of green roofs with white roofs by coupling the Princeton urban canopy model into the WRF system. They focused on a 3-day summer heat wave event, whereas the long-term performance of green roofs was not addressed. More importantly, urban hydrological processes were not adequately represented in these studies (urban irrigation, oasis effect, etc.), leading to potential uncertainties in the findings.
It is therefore imperative to implement urban hydrological processes into a coupled atmosphere–urban modeling system to investigate their impacts under a fully interacting environment. Enabled by the realistic resolution of urban hydrological processes in a recent study (Yang et al. 2015a), here we use the enhanced integrated WRF–urban modeling system to 1) evaluate the impact of hydrological processes on prediction of urban hydrometeorological variables and 2) assess the effect of green roofs at the regional scale with seasonal variability. To investigate model results under different geographical and climatic conditions, simulations are conducted for two major cities in the United States, namely, Phoenix and Houston.
a. WRF–urban modeling system
WRF is a fully compressible, nonhydrostatic modeling system that has been used for a variety of applications, ranging from local to global scale (Skamarock and Klemp 2008). Here we used WRF, version 3.4.1, to conduct simulations over study metropolitan areas. Initial meteorological conditions for the WRF simulations were obtained from the National Centers for Environmental Prediction Final Operational Global Analysis data, which were available on a 1° × 1° resolution with a 6-h temporal frequency (details can be found on http://rda.ucar.edu/datasets/ds083.2/). The Noah land surface model, coupled with the single-layer urban canopy model, was used to simulate land surface processes after initiation. Note that we adopted an enhanced version of SLUCM, which featured the integration of 1) anthropogenic latent heat, 2) urban irrigation, 3) evaporation from water-holding engineered pavements, 4) urban oasis effect, and 5) multilayer green roof system. Detailed information of individual processes can be found in previous work (Yang et al. 2015a). Other major physical parameterization schemes used in this study include: 1) the new Thompson scheme for microphysics (Thompson et al. 2008), 2) the Rapid Radiative Transfer Model for longwave radiation (Mlawer et al. 1997), 3) the Dudhia scheme for shortwave radiation (Dudhia 1989), 4) the MM5 similarity scheme for surface layer, and 5) the Yonsei University scheme for planetary boundary layer (Hong et al. 2006). Cumulus parameterization is turned on only for the outer and middle domain, using the Kain–Fritsch scheme (Kain 2004).
b. Experiment design
To compare the effect of urban hydrological processes under different geographical and climatic conditions, we selected Phoenix and Houston as our study sites. These two are among the top 10 most populous cities in the United States, whose urban heat island and hydroclimate has been extensively studied in the literature (Georgescu et al. 2012; Salamanca et al. 2011; Yang et al. 2015b). Distinct conditions in two regions (e.g., inland semiarid for Phoenix and coastal humid for Houston) facilitate a better understanding of urban hydrological processes under different geographical and climatic conditions.
For both areas, we used a two-way nested grid configuration with all three domains centered on the city (see Figs. 1a,b). Spatial resolution for the outer, middle, and inner domains was 32, 8, and 2 km, respectively. The outer domain covered a surface area of 1856 km × 1856 km, and the inner domain had a size of 212 km × 212 km. As the outer and middle domains cover portions of Mexico, MODIS global land-cover data were used (Friedl et al. 2002). For the inner domain, we used the National Land Cover Database (NLCD) 2006 (Fry et al. 2011) to represent the heterogeneous urban landscape that is subdivided into three categories (see Figs. 1c,d). We selected year 2006 for this study to represent a normal annual climatic condition for both cities. Simulations were initiated at 0000 UTC 1 November 2005 and concluded at 0000 UTC 1 December 2006. November 2005 was the spinup period and was not included in the subsequent analysis. Considering the time span of simulations and geographical locations, sea surface temperature was updated at an interval of 1 day. In this study, our analysis focused on the inner domain, and results of the other two domains are not discussed.
For each city, a total of three sets of simulation were conducted (see Table 1). The first case was a control run with the default SLUCM (hereafter “old SLUCM”) in WRF. The second case employed the recently enhanced SLUCM (Yang et al. 2015a; hereafter “new SLUCM”) with a more realistic representation of urban hydrological processes. The last case assumed a 100% areal fraction of green roof deployment over the study cities using the new SLUCM. With this experiment design, the impact of hydrological processes can be readily obtained by comparing results from the first and second cases. The difference in results between the second and last cases renders the regional impact of green roofs.
3. Impact of urban hydrological processes
Performance of the WRF simulations was evaluated against hourly meteorological observations from ground-based weather stations. Simulated 2-m air temperature T2 and 2-m dewpoint temperature at 1-h frequency were available for direct comparison to observed data. These two variables were selected because of their importance in fire weather prediction (Cheng and Steenburgh 2005). Besides, they are essential inputs to a variety of hydrological and ecological models for resolving evapotranspiration processes and plant productivity (Dodson and Marks 1997). For Phoenix, we utilized data from the Arizona Meteorological Network (AZMET) and NOAA’s National Centers for Environmental Information (NCEI). Details of the meteorological stations are summarized in Table 2. Based on NLCD 2006 land-use classification, four stations were identified as urban and the rest as rural. In Houston, only data from NCEI were used, among which three stations were urban and two were rural. Arithmetic average of ground-based measurements is compared against that of simulation results at corresponding model grids for evaluation.
Figure 2 compares the simulated annual average diurnal profiles of T2 and with the old and new SLUCM against the observations for Phoenix. Hydrological processes are expected to reduce air temperature and increase dewpoint temperature at the 2-m level of urban areas; nevertheless, Fig. 2 shows that the effect on 2-m air temperature is negligible. One important reason for this phenomenon is the limited effective area and time of hydrological processes. Among the implemented urban hydrological processes, urban irrigation and oasis effect are effective over vegetated area, which is only 20% of the urban land surface in Phoenix. Evaporation from water-holding engineered pavements functions during and shortly after rainfall; therefore, its long-term average impact is trivial. Another critical reason is the parameterization schemes used in the WRF Model. The WRF Model adopts a “tile” approach, where fluxes over built and vegetated surfaces are weighted by their respective areal fractions to calculate the total flux arise from the urban land surface. In this case, surface and air temperatures are largely determined by the built surface, whose areal fraction and temperature are significantly larger than those over the vegetated surface. Anthropogenic latent heat is directly added to the latent heat flux term so that it does not directly participate into the urban surface energy balance via the distribution of available solar radiation to dispersive and ground heat. On the other hand, moisture and humidity over the urban land surface are primarily controlled by vegetation, as there is no evaporation over the built surface most of the time.
In terms of , it is clear from Fig. 2 that modeled is significantly underestimated in WRF simulations with the old SLUCM. Diurnal minimum is found at 1600 local time (LT), which corresponds to maximum T2. Incorporating urban hydrological processes, model prediction agrees better with observations. Increase of can be up to about 1.5°C for the urban area across the diurnal cycle. Via urban–rural circulations, urban hydrological processes also have detectable effects on rural , though with a smaller magnitude. Results for Houston are plotted in Fig. 3. Located in a coastal area, Houston has a lower air temperature and a higher dewpoint temperature than Phoenix. Daytime onshore flow provides moisture for the urban area and therefore weakens the influence of urban hydrological processes. Increase in by hydrological processes is about 0.3°C for Houston.
Figures 2 and 3 illustrate that urban hydrological processes have limited effects on T2. Hence, we used the daily maximum, mean, and minimum 2-m dewpoint temperatures for statistical analysis. Evaluating these temperatures is very useful as they are indices for climate extremes (Alexander et al. 2006; Perkins et al. 2007). Seasonally averaged results for Phoenix and Houston are summarized in Table 3, which shows that with the old SLUCM, WRF simulations considerably underestimate daily maximum and mean . The new SLUCM with enhanced urban hydrological modeling enables improved predictions for the entire simulation period.
For Phoenix, increase in daily maximum and mean is about 1.2°, 0.4°, 0.3°, and 0.6°C for winter [December–February (DJF)], spring [March–May (MAM)], summer [June–August (JJA)], and fall [September–November (SON)], respectively. With respect to daily minimum , the improvement is less clear. Improvement is observed in winter while degradation is reported for fall. In Houston, the impact of urban hydrological processes is weak because of the presence of sea–land breezes. During spring and summer when temperature difference is distinct between the land and sea surface, strong onshore flow makes the effect of urban hydrological processes negligible. For fall and winter, average daily maximum, mean, and minimum is increased by about 0.6° and 0.3°C, respectively.
4. Regional hydroclimatic effect of green roofs
Using the model with enhanced urban hydrology, here we conducted simulations to investigate the regional effect of green roofs for both Phoenix and Houston. Our hypothetical scenario assumes that all rooftops of the two cities are replaced by green roofs, with results indicating the maximum possible effect. Here we select short grasses for green roof vegetation type with a 0.3-m deep loam soil layer. Sensitivity of green roof performance to parameters related to soil and vegetation type is referred to the previous study (Yang and Wang 2014).
Figure 4 shows the seasonal variability of impacts of green roofs on land surface temperature Ts at 1400 LT for Phoenix. We present the result at 1400 LT, as subsequent analysis finds the time corresponds to diurnal maximum effect (see Fig. 13, described in greater detail below). From Fig. 4, it is clear that green roofs can reduce Ts of the urban area by more than 4°C throughout the year (cf. Fig. 1c for the urban area in Phoenix). Compared to other seasons, fall (SON) has the smallest reduction in Ts, primarily because of the extensive amount of precipitation simulated in this season. Simulated accumulated precipitation depth for spring, summer, fall, and winter is about 47.9, 59.4, 100.5, and 4.8 mm, respectively. Seasonal variation of precipitation in the case with implemented green roofs is similar to that of the control case (see Table 4). Compared to in situ measurements, model prediction underestimates precipitation in summer and overestimates it in fall for Phoenix. The deviation in precipitation pattern can be caused by various physical parameterizations, such as microphysics, planetary boundary layer, and cumulus schemes. Closing the gap between simulated and observed precipitation requires a thorough sensitivity analysis in the future and is beyond the scope of this study.
In a previous offline study, Yang et al. (2015a) reported a green roof cooling of the Phoenix metropolitan area by about 8°C at 1400 LT in summer. This significant difference between offline and online simulation results indicates that feedback between the atmospheric system and land surface has notable influences on the performance of green roofs. Results in this study, derived from the fully coupled WRF–urban modeling system, are more representative of actual effects. To demonstrate impacts during nighttime hours, results at 0200 LT are shown in Fig. 5. With additional soil layers on top of buildings, green roofs are able to store extra solar energy during daytime as compared to conventional roofs. The energy is released and causes a considerable warming effect at night. Figure 5 demonstrates that increase in Ts is about 1°–2°C from spring to fall and is less than 1°C in winter. The magnitude of nighttime warming is much smaller than that of daytime cooling by green roofs for Phoenix. These results are consistent with recent high-resolution simulations for urbanizing regions in California, which similarly indicated an increased nighttime warming tendency for green roofs deployment that was considerably smaller in magnitude relative to daytime cooling (Georgescu 2015).
Cooling effect of green roofs on Ts at 1400 LT for Houston is shown in Fig. 6 (cf. Fig. 1d for the urban area in Houston). Compared to Phoenix, temperature reduction in fall and winter for Houston is much lower. Evaporative cooling of green roofs is mainly controlled by two factors: available energy and availability of water at the surface. As precipitation for Houston is abundant throughout the year, evapotranspiration arising from green roofs is largely determined by the available solar radiation. Houston is known to have a much cloudier weather and thus less total available solar radiation than the desert city Phoenix. In winter, when the sun angle is lower, solar radiation intensity decreases significantly and green roofs become relatively less effective. Precipitation also plays a role in determining the cooling effect. During the simulation period, Houston receives nearly double the amount of rainfall in fall as compared to spring (see Table 4), which indicates fewer clear days on average, leading to ineffectiveness of green roofs.
Difference in simulated 2-m air temperature between 0% and 100% green roof fraction cases at 1400 LT for Phoenix is shown in Fig. 7. Opposite to the trend of surface temperature, it is found that the strongest cooling effect on T2 of more than 1.2°C occurs in winter, while the smallest reduction is less than 0.8°C in summer. A reason for this phenomenon is that nonlinear relation exists between surface temperature and 2-m air temperature. When green roofs reduce Ts, buoyancy effect is also reduced such that the reduction of T2 is smaller than the reduction of Ts. Another critical factor contributing to the phenomenon is the warming effect caused by green roofs at night, as demonstrated in Fig. 8. Compared to winter, the urban land surface in summer receives a considerably enhanced solar radiative flux, which is stored via a large thermal mass of manmade structures and is subsequently released at night. In the absence of incoming solar radiation, vertical mixing over urban terrain in nighttime is weak so that the evolution of air temperature is steady (Poulos et al. 2002). As a consequence, increase in T2 by heat released from green roofs dissipates slowly until sunrise when surface heating modifies the stability condition of the boundary layer. Figure 8 clearly illustrates that increase in T2 in summer is more significant than that in winter, in terms of both the influence area and the magnitude. This nighttime warming impedes cooling of air temperature in daytime and results in the stronger cooling of T2 in winter as compared to summer.
Figure 9 demonstrates the regional effect of green roofs on T2 at 1400 LT for Houston. It is noteworthy that unlike in Phoenix, the order of reduction in T2 among different seasons generally follows that of Ts for Houston. This is primarily due to the negligible nighttime warming of air temperature in Houston throughout the year (results not shown here). In a coastal area, different surface cooling over land and sea results in a temperature gap in overlying air layers and consequently leads to nighttime advection of marine air. Simulated results of the 10-m wind speed at 2100 LT (sunset around 2000 LT) for Houston during summer are presented in Fig. 10. Advection of marine air toward the land tends to reduce T2 over sea and increase T2 over land. As illustrated in Fig. 8, green roofs tend to increase T2 over urban areas at night. The increase in urban T2 reduces the land–sea air temperature difference, weakens the nocturnal advection, and eventually offsets the warming effect on T2 over urban areas. Figure 10 shows that green roofs decrease 10-m wind speed by about 1 m s−1 in the bay area. The combined effect of green roofs on nocturnal T2 is therefore insignificant. With an insignificant nighttime warming, daytime cooling of T2 follows the trend of reduction in Ts. Reduction in T2 at 1400 LT for Houston is less than 0.8°C in winter and can be up to more than 1.2°C in summer. It is worth mentioning that the cooling effect on T2 has a larger spatial coverage in Houston because of the existence of land–sea circulation, especially in spring and summer when there is a considerable gap between land and sea surface temperature.
Impact of green roofs on 2-m dewpoint temperature for Phoenix at 1400 LT is shown in Fig. 11. Through evaporative cooling, green roofs are able to increase moisture and decrease temperature of the near-surface air layer, thus leading to a substantial rise in for the entire simulation period. Sunwoo et al. (2006a,b) suggested relative humidity should be maintained at greater than 30% to avoid dryness of the eyes and skin. Therefore, increased air humidity can enhance the thermal comfort of pedestrians in a dry environment, such as the premonsoon season in Phoenix (relative humidity ≈ 12%). However, extra moisture in the monsoon season can aggravate the thermal discomfort of residents, as illustrated in a recent study on urban irrigation (Yang and Wang 2015). This two-sided effect of green roofs needs special attention, especially in humid regions like Houston. Figure 11 demonstrates that increase of can be up to more than 1.2°C across the year. It is easy to recognize that seasonal variation of the increase in is similar to that of the decrease in T2 for Phoenix because of the nonlinear relationship between saturated vapor pressure with air temperature. According to the Clausius–Clapeyron equation, a same amount of increase in absolute humidity of air will cause a larger increase of dewpoint temperature at a lower air temperature. At night, the evapotranspiration rate becomes much slower as the driving force (solar radiation) disappears, and thus the influence of green roofs on becomes insignificant (results not shown here). Figure 12 presents the results at 1400 LT for Houston. The relation between green roofs’ effects on and T2 is consistent in both cities.
Besides spatial variation, temporal variation of the impact of green roofs is investigated. Realizing the maximum and minimum effects of green roofs in a temporal cycle has important implications for urban planning. In fact, the time at which spatial effect of green roofs was presented (e.g., 1400 and 0200 LT in above context) is selected based on diurnal results in Fig. 13. Figure 13 demonstrates the diurnal impact of green roofs averaged over the entire Phoenix urban area. As expected, latent heat flux (LE) from green roofs increases with intensity of solar radiation at the surface, and the largest increment of more than 70 W m−2 is found in summer. Additionally, daytime sunshine duration controls the effective period of green roofs. It is indicated from Fig. 13a that green roofs function about 4 h more in summer than in winter.
While the increase in LE is the largest in summer, it does not necessarily lead to the greatest reduction in Ts. As shown in Fig. 13b, the strongest cooling of the urban land surface by green roofs occurs in spring instead of summer, owing to the monsoon period from July to September in Arizona. The extensive amount of precipitation in fall also results in a smaller reduction of Ts than that in winter. With respect to the nighttime warming, the increase of Ts by green roofs from the largest to the smallest is in summer, fall, spring, and winter. The order is the same for the increase in nighttime T2. The average increment of nighttime T2 is about 1.1°C in summer and about 0.3°C in winter. As mentioned, the difference in nighttime warming has significant implications for the daytime cooling process. Consequently, the largest reduction of daytime T2 and the largest increase of by green roofs occur in winter. In this particular period from December 2015 to February 2016, Phoenix was abnormally dry, leading to low evapotranspiration from green roofs, such that a few rainfall events in this period cause the spike in latent heat flux in Fig. 13d.
The average impact of green roofs on studied variables for Houston is qualitatively similar to that for Phoenix; however, the seasonal variation of the impact differs considerably. With sufficient supply of water from precipitation, effectiveness of green roofs in Houston largely depends on the duration and strength of incoming solar radiation. Figure 14a shows that increased LE by green roofs can be up to more than 130 W m−2 in spring and summer, which is remarkably larger than the increase of about 80 W m−2 in fall. With respect to Ts, Fig. 14b demonstrates that daytime cooling effect is the strongest in summer and the weakest in winter, while nighttime warming is almost negligible. Diurnal impact of green roofs on T2 across various seasons is similar to that on Ts. The peak cooling effect is found to be about 1°C in spring and summer. As land–sea circulation mixes the air layer of coastal area, increased by green roofs has a relatively limited seasonal variation.
5. Concluding remarks
In this study, we applied the WRF Model framework with an enhanced single-layer urban canopy model to assess the effect of hydrological processes on urban hydrometeorology. Evaluation against field measurements illustrates that more realistic representation of urban hydrological processes improves the prediction of the 2-m dewpoint temperature. Results of regional hydroclimate simulations indicate that green roofs are effective in reducing daytime air temperature and increasing dewpoint temperature over urban areas. The impact of green roofs exhibits strong diurnal and seasonal variability and depends on geographical and climatic conditions. It is noteworthy that urban vegetation is largely represented as grasses and short corps in the WRF–urban modeling system, whereas physical resolution of more diverse urban vegetation types, for example, tall trees, and their hydrometeorological effect, such as on radiative energy exchange, remains an open challenge in regional simulations (Krayenhoff et al. 2014; Wang 2014b; Wang et al. 2016).
Numerical experiments clearly demonstrate the effect of urban hydrological processes for Phoenix and Houston in this study. However, with a limited computational resource, uncertainty of simulation results is not adequately addressed after fulfilling the fine spatial resolution and long simulation period in experimental setup. The uncertainty consists mainly of two parts: 1) sensitivity of model results to initial and boundary conditions and 2) sensitivity to physical parameterizations in the WRF Model. To reduce the uncertainty and provide a better quantitative estimation of the effect, an ensemble approach will be needed in future work.
Comparing results from this study and a previous offline study, it is indicated that land–atmosphere interactions cannot be ignored in quantifying the influence of surface hydrological process. In the coastal area, land–sea circulation mixes the near-surface air layer, leading to a weaker effect of hydrological processes on the meteorological field than that of the inland area. To accurately evaluate sustainable adaptation/mitigation strategies for urban areas, numerical experiments should be carried out with a fully interacting land–atmosphere modeling system. In addition, modification of urban landscape will have strong implications for hydrometeorology of surrounding rural areas, necessitating serious consideration and planning prior to large-scale implementation. Finally, this study explores high-resolution numerical simulation of green roofs at the annual scale; we expect the findings to provide a useful guidance for sustainable development of other cities.
This work is supported by the National Science Foundation (NSF) under Grants CBET-1435881 and CBET-1444758 and the USDA-NIFA Agriculture and Food Research Initiative (Awards 2015-67003-23508 and 2015-67003-23460).