Our understanding of the full range of natural variability in streamflow, including how modern flow compares to the past, is poorly understood for the Upper Indus Basin (UIB) because of short instrumental gauge records. To help address this challenge, we use Hierarchical Bayesian Regression (HBR) with partial pooling to develop six centuries long (1394-2008 C.E.) streamflow reconstructions at three UIB gauges (Doyian, Gilgit, and Kachora), concurrently demonstrating that HBR can be used to reconstruct short records with interspersed missing data. At one gauge (Partab Bridge), with a longer instrumental record (47 years), we develop reconstructions using both Bayesian Regression (BR) and the more conventionally used Principal Components Regression (PCR). The reconstructions produced by PCR and BR at Partab Bridge are nearly identical and yield comparable reconstruction skill statistics, highlighting that the resulting tree-ring reconstruction of streamflow is not dependent on the choice of statistical method. Reconstructions at all four reconstructions indicate flow levels in the 1990s were higher than mean flow for the past six centuries. While streamflow appears most sensitive to accumulated winter (January-March) precipitation and summer (MJJAS) temperature, with warm summers contributing to high flow through increased melt of snow and glaciers, shifts in winter precipitation and summer temperatures cannot explain the anomalously high flow during the 1990s. Regardless, the sensitivity of streamflow to summer temperatures suggests that projected warming may increase streamflow in coming decades, though long-term water risk will additionally depend on changes in snowfall and glacial mass balance.