Belowground urban stormwater networks (BUSNs) are critical for removing excess rainfall from impervious urban areas and preventing or mitigating urban flooding. However, available BUSN data are sparse, preventing the modeling and analysis of urban hydrologic processes at regional and larger scales. We propose a novel algorithm for estimating BUSNs by drawing on concepts from graph theory and existing, extensively available land surface data, such as street network, topography, and land use/land cover. First, we derive the causal relationships between the topology of BUSNs and urban surface features based on graph theory concepts. We then apply the causal relationships and estimate BUSNs using web-service data retrieval, spatial analysis, and high-performance computing techniques. Finally, we validate the derived BUSNs in the metropolitan areas of Los Angeles, Seattle, Houston, and Baltimore in the US, where real BUSN data are partly available to the public. Results show that our algorithm can effectively capture 59 %–76 % of the topology of real BUSN data, depending on the supporting data quality. This algorithm has promising potential to support large-scale urban hydrologic modeling and future urban drainage system planning.