Stable isotopes of water are important tracers in hydrologic research for understanding water partitioning between vegetation, groundwater, and runoff, but are rarely applied to large watersheds with persistent snowpack and complex topography. We combined an extensive isotope dataset with a coupled hydrologic and snow isotope fractionation model to assess mechanisms of isotopic inputs into the soil zone and implications on recharge dynamics within a large, snow-dominated watershed of the Upper Colorado River Basin. Results indicate seasonal isotopic variability and isotope lapse rates of net precipitation are the dominant control on isotopic inputs to the basin. Snowpack fractionation processes account for <5% annual isotope influx variability. Isotopic fractionation processes are most important in the shrub-dominated upper montane. Effects of isotopic fractionation are less important in the low-density conifer forests of the upper subalpine due to vegetative shading, low aridity, and a deep, persistent snowpack that buffers small sublimation losses. Melt fractionation can have sub-seasonal effects on snowmelt isotope ratios with initial snowmelt depleted but later snowmelt relatively enriched in heavy isotopes through the isotopic mass balance of the remaining snowpack, with the efficiency of isotopic exchange between ice and liquid water declining as snow ablation progresses. Hydrologic analysis indicates maximum recharge in the upper subalpine with wet years producing more isotopically depleted snowmelt (1-2‰ reduction in d18O) through reduced aridity when energy-limited. The five-year volume-weighted d18O in this zone (18.2±0.4‰) matches groundwater observations from multiple deep wells, providing evidence that the upper subalpine is a preferential recharge zone in mountain systems.
