In many partially glacierized watersheds glacier recession driven by a warming climate could lead to complex patterns of streamflow response over time, often marked with rapid increases followed by sharp declines, depending on initial glacier ice cover and rate of climate change. Capturing such "phases'' of hydrologic response is critical in regions where communities rely on glacier meltwater, particularly during low flows. In this paper, we investigate glacio-hydrologic response in the headwaters of the Zongo River, Bolivia, under climate change using a distributed glacio-hydrological model over the period of 1987-2100. Model predictions are evaluated through comparisons with satellite-derived glacier extent estimates, glacier surface velocity, in situ glacier mass balance, surface energy flux, and stream discharge measurements. Historically (1987-2010) modeled glacier melt accounts for 27% of annual runoff, and 61% of dry season (JJA) runoff on average. During this period the relative glacier cover was observed to decline from 35 to 21% of the watershed. In the future, annual and dry season discharge is projected to decrease by 4% and 27% by midcentury and 25% and 57% by the end of the century, respectively, following the loss of 81% of the ice in the watershed. Modeled runoff patterns evolve through the interplay of positive and negative trends in glacier melt and increased evapotranspiration as the climate warms. Sensitivity analyses demonstrate that the selection of model surface energy balance parameters greatly influences the trajectory of hydrological change projected during the first half of the 21st century. These model results underscore the importance of coupled glacio-hydrology modeling.