Rapid urbanization and population growth occurring in the cities of South Western

United States have led to significant modifications in its environment at local and

regional scales. Both local and regional climate changes are expected to have massive

impacts on the hydrology of Colorado River Basin (CRB), thereby accentuating the need

of study of hydro-climatic impacts on water resource management in this region. This

thesis is devoted to understanding the impact of land use and land cover (LULC) changes

on the local and regional hydroclimate, with the goal to address urban planning issues

and provide guidance for sustainable development.

In this study, three densely populated urban areas, viz. Phoenix, Las Vegas and

Denver in the CRB are selected to capture the various dimensions of the impacts of land

use changes on the regional hydroclimate in the entire CRB. Weather Research and

Forecast (WRF) model, incorporating the latest urban modeling system, is adopted for

regional climate modeling. Two major types of urban LULC changes are studied in this

Thesis: (1) incorporation of urban trees with their radiative cooling effect, tested in

Phoenix metropolitan, and (2) projected urban expansion in 2100 obtained from

Integrated Climate and Land Use Scenarios (ICLUS) developed by the US

Environmental Protection Agency for all three cities.

The results demonstrated prominent nocturnal cooling effect of due to radiative

shading effect of the urban trees for Phoenix reducing urban surface and air temperature

by about 2~9 °C and 1~5 °C respectively and increasing relative humidity by 10~20%

during an mean diurnal cycle. The simulations of urban growth in CRB demonstratedii

nocturnal warming of about 0.36 °C, 1.07 °C, and 0.94 °C 2m-air temperature and

comparatively insignificant change in daytime temperature, with the thermal environment

of Denver being the most sensitive the urban growth. The urban hydroclimatic study

carried out in the thesis assists in identifying both context specific and generalizable

relationships, patterns among the cities, and is expected to facilitate urban planning and

management in local (cities) and regional scales.

The combination of rapid urban growth and climate change places stringent constraints on multisector sustainability of cities. Green infrastructure provides a great potential for mitigating anthropogenic-induced urban environmental problems; nevertheless, studies at city and regional scales are inhibited by the deficiency in modelling the complex transport coupled water and energy inside urban canopies. This dissertation is devoted to incorporating hydrological processes and urban green infrastructure into an integrated atmosphere-urban modelling system, with the goal to improve the reliability and predictability of existing numerical tools. Based on the enhanced numerical tool, the effects of urban green infrastructure on environmental sustainability of cities are examined.

Findings indicate that the deployment of green roofs will cool the urban environment in daytime and warm it at night, via evapotranspiration and soil insulation. At the annual scale, green roofs are effective in decreasing building energy demands for both summer cooling and winter heating. For cities in arid and semiarid environments, an optimal trade-off between water and energy resources can be achieved via innovative design of smart urban irrigation schemes, enabled by meticulous analysis of the water-energy nexus. Using water-saving plants alleviates water shortage induced by population growth, but comes at the price of an exacerbated urban thermal environment. Realizing the potential water buffering capacity of urban green infrastructure is crucial for the long-term water sustainability and subsequently multisector sustainability of cities. Environmental performance of urban green infrastructure is determined by land-atmosphere interactions, geographic and meteorological conditions, and hence it is recommended that analysis should be conducted on a city-by-city basis before actual implementation of green infrastructure.

Rapid urban expansion and the associated landscape modifications have led to significant changes of surface processes in built environments. These changes further interact with the overlying atmospheric boundary layer and strongly modulate urban microclimate. To capture the impacts of urban land surface processes on urban boundary layer dynamics, a coupled urban land-atmospheric modeling framework has been developed. The urban land surface is parameterized by an advanced single-layer urban canopy model (SLUCM) with realistic representations of urban green infrastructures such as lawn, tree, and green roof, etc. The urban atmospheric boundary layer is simulated by a single column model (SCM) with both convective and stable schemes. This coupled SLUCM-SCM framework can simulate the time evolution and vertical profile of different meteorological variables such as virtual potential temperature, specific humidity and carbon dioxide concentration. The coupled framework has been calibrated and validated in the metropolitan Phoenix area, Arizona. To quantify the model sensitivity, an advanced stochastic approach based on Markov-Chain Monte Carlo procedure has been applied. It is found that the development of urban boundary layer is highly sensitive to surface characteristics of built terrains, including urban land use, geometry, roughness of momentum, and vegetation fraction. In particular, different types of urban vegetation (mesic/xeric) affect the boundary layer dynamics through different mechanisms. Furthermore, this framework can be implanted into large-scale models such as Weather Research and Forecasting model to assess the impact of urbanization on regional climate.

Land surface fluxes of energy and mass developed over heterogeneous mountain landscapes are fundamental to atmospheric processes. However, due to their high complexity and the lack of spatial observations, land surface processes and land-atmosphere interactions are not fully understood in mountain regions. This thesis investigates land surface processes and their impact on convective precipitation by conducting numerical modeling experiments at multiple scales over the North American Monsoon (NAM) region. Specifically, the following scientific questions are addressed: (1) how do land surface conditions evolve during the monsoon season, and what are their main controls?, (2) how do the diurnal cycles of surface energy fluxes vary during the monsoon season for the major ecosystems?, and (3) what are the impacts of surface soil moisture and vegetation condition on convective precipitation?

Hydrologic simulation using the TIN-based Real-time Integrated Basin Simulator (tRIBS) is firstly carried out to examine the seasonal evolution of land surface conditions. Results reveal that the spatial heterogeneity of land surface temperature and soil moisture increases dramatically with the onset of monsoon, which is related to seasonal changes in topographic and vegetation controls. Similar results are found at regional basin scale using the uncoupled WRF-Hydro model. Meanwhile, the diurnal cycles of surface energy fluxes show large variation between the major ecosystems. Differences in both the peak magnitude and peak timing of plant transpiration induce mesoscale heterogeneity in land surface conditions. Lastly, this dissertation examines the upscale effect of land surface heterogeneity on atmospheric condition through fully-coupled WRF-Hydro simulations. A series of process-based experiments were conducted to identify the pathways of soil moisture-rainfall feedback mechanism over the NAM region. While modeling experiments confirm the existence of positive soil moisture/vegetation-rainfall feedback, their exact pathways are slightly different. Interactions between soil moisture, vegetation cover, and rainfall through a series of land surface and atmospheric boundary layer processes highlight the strong land-atmosphere coupling in the NAM region, and have important implications on convective rainfall prediction. Overall, this dissertation advances the study of complex land surface processes over the NAM region, and made important contributions in linking complex hydrologic, ecologic and atmospheric processes through numerical modeling.

This research work uses the Weather Research and Forecasting Model to study the effect of large wind farms with an area of 900 square kilometers and a high power density of 7.58 W/m2 on regional climate. Simulations were performed with a wind farm parameterization scheme turned on in south Oregon. Control cases were also run with the parameterization scheme turned off. The primary emphasis was on offshore wind farms. Some analysis on onshore wind farms was also performed. The effects of these wind farms were studied on the vertical profiles of temperature, wind speed, and moisture as well as on temperature and on wind speed near the surface and at hub height. The effects during the day and at night were compared. Seasonal variations were also studied by performing simulations in January and in July. It was seen that wind farms produce a reduction in wind speed at hub height and that the downward propagation of this reduction in wind speed lessens as the atmosphere becomes more stable. In all the cases studied, the wind farms produced a warming effect near the surface, with greater atmospheric stability leading to higher near-surface temperatures. It was also observed that wind farms caused a drying effect below the hub height and a moistening effect above it, because they had facilitated vertical transport of moisture in the air from the lower layers of the atmosphere to the layers of the atmosphere above the wind farm.

Bioretention basins are a common stormwater best management practice (BMP) used to mitigate the hydrologic consequences of urbanization. Dry wells, also known as vadose-zone wells, have been used extensively in bioretention basins in Maricopa County, Arizona to decrease total drain time and recharge groundwater. A mixed integer nonlinear programming (MINLP) model has been developed for the minimum cost design of bioretention basins with dry wells.

The model developed simultaneously determines the peak stormwater inflow from watershed parameters and optimizes the size of the basin and the number and depth of dry wells based on infiltration, evapotranspiration (ET), and dry well characteristics and cost inputs. The modified rational method is used for the design storm hydrograph, and the Green-Ampt method is used for infiltration. ET rates are calculated using the Penman Monteith method or the Hargreaves-Samani method. The dry well flow rate is determined using an equation developed for reverse auger-hole flow.

The first phase of development of the model is to expand a nonlinear programming (NLP) for the optimal design of infiltration basins for use with bioretention basins. Next a single dry well is added to the NLP bioretention basin optimization model. Finally the number of dry wells in the basin is modeled as an integer variable creating a MINLP problem. The NLP models and MINLP model are solved using the General Algebraic Modeling System (GAMS). Two example applications demonstrate the efficiency and practicality of the model.

The increasingly recurrent extraordinary flood events in the metropolitan area of Monterrey, Mexico have led to significant stakeholder interest in understanding the hydrologic response of the Santa Catarina watershed to extreme events. This study analyzes a flood mitigation strategy proposed by stakeholders through a participatory workshop and are assessed using two hydrological models: The Hydrological Modeling System (HEC-HMS) and the Triangulated Irregular Network (TIN)-based Real-time Integrated Basin Simulator (tRIBS).

The stakeholder-derived flood mitigation strategy consists of placing new hydraulic infrastructure in addition to the current flood controls in the basin. This is done by simulating three scenarios: (1) evaluate the impact of the current structure, (2) implementing a large dam similar to the Rompepicos dam and (3) the inclusion of three small detention dams. These mitigation strategies are assessed in the context of a major flood event caused by the landfall of Hurricane Alex in July 2010 through a consistent application of the two modeling tools. To do so, spatial information on topography, soil, land cover and meteorological forcing were assembled, quality-controlled and input into each model. Calibration was performed for each model based on streamflow observations and maximum observed reservoir levels from the National Water Commission in Mexico.

Simulation analyses focuses on the differential capability of the two models in capturing the spatial variability in rainfall, topographic conditions, soil hydraulic properties and its effect on the flood response in the presence of the different flood mitigation structures. The implementation of new hydraulic infrastructure is shown to have a positive impact on mitigating the flood peak with a more favorable reduction in the peak at the outlet from the larger dam (16.5% in tRIBS and 23% in HEC-HMS) than the collective effect from the small structures (12% in tRIBS and 10% in HEC-HMS). Furthermore, flood peak mitigation depends strongly on the number and locations of the new dam sites in relation to the spatial distribution of rainfall and flood generation. Comparison of the two modeling approaches complements the analysis of available observations for the flood event and provides a framework within which to derive a multi-model approach for stakeholder-driven solutions.

Vegetative filter strips (VFS) are an effective methodology used for storm water management particularly for large urban parking lots. An optimization model for the design of vegetative filter strips that minimizes the amount of land required for stormwater management using the VFS is developed in this study. The resulting optimization model is based upon the kinematic wave equation for overland sheet flow along with equations defining the cumulative infiltration and infiltration rate.

In addition to the stormwater management function, Vegetative filter strips (VFS) are effective mechanisms for control of sediment flow and soil erosion from agricultural and urban lands. Erosion is a major problem associated with areas subjected to high runoffs or steep slopes across the globe. In order to effect economy in the design of grass filter strips as a mechanism for sediment control & stormwater management, an optimization model is required that minimizes the land requirements for the VFS. The optimization model presented in this study includes an intricate system of equations including the equations defining the sheet flow on the paved and grassed area combined with the equations defining the sediment transport over the vegetative filter strip using a non-linear programming optimization model. In this study, the optimization model has been applied using a sensitivity analysis of parameters such as different soil types, rainfall characteristics etc., performed to validate the model

A model is presented for real-time, river-reservoir operation systems. It epitomizes forward-thinking and efficient approaches to reservoir operations during flooding events. The optimization/simulation includes five major components. The components are a mix of hydrologic and hydraulic modeling, short-term rainfall forecasting, and optimization and reservoir operation models. The optimization/simulation model is designed for ultimate accessibility and efficiency. The optimization model uses the meta-heuristic approach, which has the capability to simultaneously search for multiple optimal solutions. The dynamics of the river are simulated by applying an unsteady flow-routing method. The rainfall-runoff simulation uses the National Weather Service NexRad gridded rainfall data, since it provides critical information regarding real storm events. The short-term rainfall-forecasting model utilizes a stochastic method. The reservoir-operation is simulated by a mass-balance approach. The optimization/simulation model offers more possible optimal solutions by using the Genetic Algorithm approach as opposed to traditional gradient methods that can only compute one optimal solution at a time. The optimization/simulation was developed for the 2010 flood event that occurred in the Cumberland River basin in Nashville, Tennessee. It revealed that the reservoir upstream of Nashville was more contained and that an optimal gate release schedule could have significantly decreased the floodwater levels in downtown Nashville. The model is for demonstrative purposes only but is perfectly suitable for real-world application.

Land-atmosphere interactions of semiarid shrublands have garnered significant scientific interest. One of the main tools used for this research is the eddy covariance (EC) method, which measures fluxes of energy, water vapor, and carbon dioxide. EC fluxes can be difficult to interpret due to complexities within the EC footprint (i.e. the surface conditions that contribute to the flux measurements). Most EC studies use a small number of soil probes to estimate the land surface states underlying the measured fluxes, which likely undersamples the footprint-scale conditions, especially in semiarid shrublands which are characterized by high spatial and temporal variability. In this study, I installed a dense network of soil moisture and temperature probe profiles in the footprint region of an EC tower at two semiarid sites: a woody savanna in southern Arizona and a mixed shrubland in southern New Mexico. For data from May to September 2013, I link land surface states to EC fluxes through daily footprints estimated using an analytical model. Novel approaches are utilized to partition evapotranspiration, estimate EC footprint soil states, connect differences in fluxes to footprint composition, and assess key drivers behind soil state variability. I verify the hypothesis that a small number of soil probes poorly estimates the footprint conditions for soil moisture, due to its high spatial variability. Soil temperature, however, behaves more consistently in time and space. As such, distributed surface measurements within the EC footprint allow for stronger ties between evapotranspiration and moisture, but demonstrate no significant improvement in connecting sensible heat flux and temperature. I also find that in these systems vegetation cover appears to have stronger controls on soil moisture and temperature than does soil texture. Further, I explore the influence of footprint vegetation composition on the measured fluxes, which reveals that during the monsoon season evaporative fraction tends to increase with footprint bare soil coverage for the New Mexico site and that the ratio of daily transpiration to evapotranspiration increases with grass coverage at the Arizona site. The thesis results are useful for understanding the land-atmosphere interactions of these ecosystems and for guiding future EC studies in heterogeneous landscapes.