Recovery Ice Stream has a substantial number of active subglacial lakes that are observed, with satellite altimetry, to grow and drain over multiple years. These lakes store and release water that could be important for controlling the velocity of the ice stream. We apply a subglacial hydrology model to analyze lake growth and drainage characteristics together with the simultaneous development of the ice stream hydrological network. Our outputs produce a good match between modeled lake location and those identified using satellite altimetry for many of the lakes. The modeled subglacial system demonstrates development of pressure waves that initiate at the ice stream neck and transit to within 100 km of the terminus. These waves alter the hydraulic potential of the ice stream and encourage growth and drainage of the subglacial lakes. Lake drainage can cause large R-channels to develop between basal overdeepenings that persist for multiple years. The pressure waves, along with lake growth and drainage rates, do not identically repeat over multiple years, due to basal network development. This suggests that the subglacial hydrology of Recovery Ice Stream is influenced by regional drainage development on the scale of hundreds of kilometers rather than local conditions over tens of kilometers.