Proglacial lakes, formed during retreat of the Laurentide ice sheet, evolved quickly as outlets became ice-free and the earth deformed through glacial isostatic adjustment. With high-resolution digital elevation models (DEMs) and GIS methods, it is possible to reconstruct the evolution of surface hydrology. When a DEM deforms through time as predicted by our model of viscoelastic earth relaxation, the entire surface hydrologic system with its lakes, outlets, shorelines and rivers also evolves without requiring assumptions of outlet position. The method is applied to proglacial Lake Oshkosh in Wisconsin (13,600 to 12,900 cal yr BP). Comparison of predicted to observed shoreline tilt indicates the ice sheet was about 400 m thick over the Great Lakes region. During ice sheet recession, each of the five outlets are predicted to uplift more than 100 m and then subside approximately 30 m. At its maximum extent, Lake Oshkosh covered 6600 km2 with a volume of 111 km3. Using the Hydrologic Engineering Center-River Analysis System model, flow velocities during glacial outburst floods up to 9 m/s and peak discharge of 140,000 m3/s are predicted, which could drain 33.5 km3 of lake water in 10 days and transport boulders up to 3 m in diameter.