In this work, we develop a fast hierarchical solver for solving large, sparse least squares problems. We build upon the algorithm, spaQR (sparsified QR Gnanasekaran and Darve in SIAM J Matrix Anal Appl 43(1):94–123, 2022), that was developed by authors to solve large linear systems. Our algorithm is built on top of Nested Dissection based multifrontal approach. use low-rank approximations front...