Rivers in steep mountainous landscapes control how, where, and when signals of base-level fall are transmitted to the surrounding topography. In doing so rivers play an important role in determining landscape evolution in response to external controls of tectonics and climate. However, tectonics and climate often covary and understanding how they influence landscape evolution remains a significant challenge. The Hawaiian Islands, where tectonics are minimized but climate signals are amplified, provide an opportunity to better understand how signals of climate are recorded by landscapes. Focusing on the Hawaiian Islands, I examine (1) how variability in rock mass properties and thresholds in sediment mobility determine where waterfalls form or stall along the Nāpali coast of Kauaʻi, (2) I then extend these findings to other volcanoes to test if observed physical limits in flood size, climate, and volcano gradient can determine where waterfalls form, and (3) I explore how thresholds in river incision below waterfalls limit information about the influence of climate on river incision rates. Findings from this analysis show that waterfalls form or stall where the maximum unit stream power is at or below a critical unit stream power for bedrock river incision. Climate appears to have little effect in determining where these conditions are met but where waterfalls stall or form does record information about discharge-area scaling for global maximum observed floods. Below waterfalls the maximum incision depth for rivers on the island of Kauaʻi (which formed ~ 4-5 million years ago) is approximately proportional to the inverse square root of mean annual rainfall. Though maximum river incision depths for some of the younger volcanoes do not exhibit the same dependency on mean annual rainfall rates they are comparable to the maximum incision depths observed on Kauaʻi even though they are a quarter to one-tenth the age of Kauaʻi. Importantly, these patterns of incision can be explained by thresholds in sediment mobility as recorded by river longitudinal profiles and indicate that the Hawaiian Islands are dominated by threshold conditions where signals of climate are recorded in the topography through controls on incision depth but not incision rates.

Geochronology and thermochronology are valuable tools for investigating the synergy between the deformational and erosional processes that shape mountainous terrains. Though numerous techniques have been developed to probe the rate and timing of events within these settings, the research presented here explores how scientists can use fewer samples to produce richer data products with broader contextual importance.

The beginning of this compilation focuses on establishing laboratory techniques to facilitate this goal. I developed a novel laser ablation ‘double dating’ (LADD) technique that rapidly yields paired U/Pb and (U-Th)/He dates for the accessory minerals zircon, titanite, and apatite. The technique obviates the need for geometric corrections typically applied during (U-Th)/He data reduction, enables the analysis of a broader spectrum of detrital crystals, and provides the opportunity for additional mapping and isotopic analyses that are traditionally challenging to procure and/or fraught with assumptions. Despite the technique’s promise, I also found it essential to weigh several considerations of relevance when attempting to date young (≤ Miocene) accessory minerals with low concentrations of U + Th. Consequently, I discuss the impact that such variables have on the magnitude of analytical imprecision and the data’s flexibility for geologic interpretation.

Beyond the lab, I collected a suite of bedrock and detrital samples from small catchments draining the southeastern Sierra Nevada mountains of California. Using the techniques described above as well as conventional methods for (U-Th)/He zircon dating, I compared the utility of both bedrock and detrital approaches for extrapolating local exhumation histories. I additionally tested the ability to employ detrital datasets to extrapolate cooling histories that span from mineral crystallization to rock exhumation through the upper crust. Employing principal mode dates from a combination of zircon and apatite LADD dates and detrital hornblende 40Ar/39Ar dates, I was able to derive thermal models that demonstrate the existence of significant variability in the cooling histories of various intrusive units along the eastern Sierra Nevada. While these results only scratch the surface of what’s possible within the realm of detrital-based research, this contribution demonstrates the utility of expanding the temporal and spatial scope of traditional detrital methodologies.

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.

The San Gabriel Mountains (SGM) of southern California provide the opportunity to study the topographic controls on erosion rate in a mountain range where climate and lithology are relatively constant. I use a combination of digital elevation model data, detailed channel survey data, decadal climate records, and catchment-averaged erosion rates quantified from 10Be concentrations in stream sands to investigate the style and rates of hillslope and channel processes across the transition from soil-mantled to rocky landscapes in the SGM. Specifically, I investigate (1) the interrelations among different topographic metrics and their variation with erosion rate, (2) how hillslopes respond to tectonic forcing in "threshold" landscapes, (3) the role of discharge variability and erosion thresholds in controlling the relationship between relief and erosion rate, and (4) the style and pace of transient adjustment in the western SGM to a recent increase in uplift rate. Millennial erosion rates in the SGM range from 0.03-1.1 mm/a, generally increasing from west to east. For low erosion rates (< 0.3 mm/a), hillslopes tend to be soil-mantled, and catchment-averaged erosion rates are positively correlated with catchment-averaged slope, channel steepness, and local relief. For erosion rates greater than 0.3 mm/a, hillslopes become increasingly rocky, catchment-mean hillslope angle becomes much less sensitive to erosion rate, and channels continue to steepen. I find that a non-linear relationship observed between channel steepness and erosion rate can be explained by a simple bedrock incision model that combines a threshold for erosion with a probability distribution of discharge events where large floods follow an inverse power-law. I also find that the timing of a two-staged increase in uplift rate in the western SGM based on stream profile analysis agrees with independent estimates. Field observations in the same region suggest that the relict topography that allows for this calculation has persisted for more than 7 Ma due to the stalling of migrating knickpoints by locally stronger bedrock and a lack of coarse sediment cover.