Population growth within drylands is occurring faster than growth in any other ecologic zone, putting pressure on already stressed water resources. Because the availability of surface water supplies in drylands tends to be highly variable, many of these populations rely on groundwater. A critical process contributing to groundwater recharge is the interaction between ephemeral channels and groundwater aquifers. Generally, it has been found that ephemeral channels contribute to groundwater recharge when streamflow infiltrates into the sandy bottoms of channels. This process has traditionally been studied in channels that drain large areas (10s to 100s km2). In this dissertation, I study the interactions between surface water and groundwater via ephemeral channels in a first-order watershed located on an arid piedmont slope within the Jornada Experimental Range (JER) in the Chihuahuan Desert. To achieve this, I utilize a combination of high-resolution observations and computer simulations using a modified hydrologic model to quantify groundwater recharge and shed light on the geomorphic and ecologic processes that affect the rate of recharge. Observational results indicate that runoff generated within the piedmont slope contributes significantly to deep percolation. During the short-term (6 yr) study period, we estimated 385 mm of total percolation, 62 mm/year, or a ratio of percolation to rainfall of 0.25. Based on the instrument network, we identified that percolation occurs inside channel areas when these receive overland sheetflow from hillslopes. By utilizing a modified version of the hydrologic model, TIN-based Real-time Integrated Basin Simulator (tRIBS), that was calibrated and validated using the observational dataset, I quantified the effects of changing watershed properties on groundwater recharge. Distributed model simulations quantify how deep percolation is produced during the streamflow generation process, and indicate that it plays a significant role in moderating the production of streamflow. Sensitivity analyses reveal that hillslope properties control the amount of rainfall necessary to initiate percolation while channel properties control the partitioning of hillslope runoff into streamflow and deep percolation. Synthetic vegetation experiments show that woody plant encroachment leads to increases in both deep percolation and streamflow. Further woody plant encroachment may result in the unexpected enhancement of dryland aquifer sustainability.

Ponderosa pine forests are a dominant land cover type in semiarid montane areas. Water supplies in major rivers of the southwestern United States depend on ponderosa pine forests since these ecosystems: (1) receive a significant amount of rainfall and snowfall, (2) intercept precipitation and transpire water, and (3) indirectly influence runoff by impacting the infiltration rate. However, the hydrologic patterns in these ecosystems with strong seasonality are poorly understood. In this study, we used a distributed hydrologic model evaluated against field observations to improve our understandings on spatial controls of hydrologic patterns, appropriate model resolution to simulate ponderosa pine ecosystems and hydrologic responses in the context of contrasting winter to summer transitions. Our modeling effort is focused on the hydrologic responses during the North American Monsoon (NAM), winter and spring periods. In Chapter 2, we utilized a distributed model explore the spatial controls on simulated soil moisture and temporal evolution of these spatial controls as a function of seasonal wetness. Our findings indicate that vegetation and topographic curvature are spatial controls. Vegetation controlled patterns during dry summer period switch to fine-scale terrain curvature controlled patterns during persistently wet NAM period. Thus, a climatic threshold involving rainfall and weather conditions during the NAM is identified when high rainfall amount (such as 146 mm rain in August, 1997) activates lateral flux of soil moisture and frequent cloudy cover (such as 42% cloud cover during daytime of August, 1997) lowers evapotranspiration. In Chapter 3, we investigate the impacts of model coarsening on simulated soil moisture patterns during the NAM. Results indicate that model aggregation quickly eradicates curvature features and its spatial control on hydrologic patterns. A threshold resolution of ~10% of the original terrain is identified through analyses of homogeneity indices, correlation coefficients and spatial errors beyond which the fidelity of simulated soil moisture is no longer reliable. Based on spatial error analyses, we detected that the concave areas (~28% of hillslope) are very sensitive to model coarsening and root mean square error (RMSE) is higher than residual soil moisture content (~0.07 m3/m3 soil moisture) for concave areas. Thus, concave areas need to be sampled for capturing appropriate hillslope response for this hillslope. In Chapter 4, we investigate the impacts of contrasting winter to summer transitions on hillslope hydrologic responses. We use a distributed hydrologic model to generate a consistent set of high-resolution hydrologic estimates. Our model is evaluated against the snow depth, soil moisture and runoff observations over two water years yielding reliable spatial distributions during the winter to summer transitions. We find that a wet winter followed by a dry summer promotes evapotranspiration losses (spatial averaged ~193 mm spring ET and ~ 600 mm summer ET) that dry the soil and disconnect lateral fluxes in the forested hillslope, leading to soil moisture patterns resembling vegetation patches. Conversely, a dry winter prior to a wet summer results in soil moisture increases due to high rainfall and low ET during the spring (spatially averaged 78 mm ET and 232 mm rainfall) and summer period (spatially averaged 147 mm ET and 247 mm rainfall) which promote lateral connectivity and soil moisture patterns with the signature of terrain curvature. An opposing temporal switch between infiltration and saturation excess runoff is also identified. These contrasting responses indicate that the inverse relation has significant consequences on hillslope water availability and its spatial distribution with implications on other ecohydrological processes including vegetation phenology, groundwater recharge and geomorphic development. Results from this work have implications on the design of hillslope experiments, the resolution of hillslope scale models, and the prediction of hydrologic conditions in ponderosa pine ecosystems. In addition, our findings can be used to select future hillslope sites for detailed ecohydrological investigations. Further, the proposed methodology can be useful for predicting responses to climate and land cover changes that are anticipated for the southwestern United States.

Numerous studies have examined the interplay of climate, tectonics, biota and erosion and found that these variables are intertwined in a complicated system of feedbacks and as a result, some of these factors are often oversimplified or simply neglected. To understand the interplay of these factors one must understand the processes that transport or inhibit transport of soil. This study uses the short-lived, fallout-derived, radionuclides 137Cs and 210Pb to identify soil transport processes and to quantify soil transport using the profile distribution model for 137Cs. Using five field sites in the San Gabriel Mountains of California, I address four questions: (1) Is there a process transition between high and low gradient slopes observable with short-lived isotopes? (2) Do convex hilltops reflect short-term equilibrium erosion rates? (3) Do linear transects of pits accurately characterize hillslope averaged erosion rates? and (4) What role does fire play on short-term soil transport and isotope distribution? I find no evidence supporting a process transition from low gradient to high gradient slopes but also find that significant spatial variability of erosion rates exist. This spatial variability is the result of sensitivity of the method to small scale variations in isotopes and indicates that small scale processes may dominate broader scale trends. I find that short-term erosion rates are not at equilibrium on a convex hilltop and suggest the possibility of a headward incision signal. Data from a post-fire landscape indicates that fires may create complications in 137Cs and 210Pb distribution that current models for erosion calculation do not account for. I also find that across all my field sites soil transport processes can be identified and quantified using short-lived isotopes and I suggest high resolution grid sampling be used instead of linear transects so that small scale variability can be averaged out.