In this work, we present our DLPNO-CCSDT(Q) algorithm as the first, to the best of our knowledge, linear-scaling treatment of coupled-cluster quadruples through the use of local pair natural orbitals. Our algorithm offers the accuracy of frozen natural orbital-based approaches (in the kJ mol -1 range) at a fraction of the cost, due to our quadruples natural orbitals (QNOs) being computed from DLPNO-CCSDT doubles amplitudes, as opposed to MP2 frozen natural orbitals. We will show that the target accuracy can be reached even at a relatively loose QNO occupation number cutoff of 10 -6 . Our algorithm also inherits the linear scaling properties of DLPNO approaches at lower levels of theory, and we show this through timings run on a series of growing linear alkanes (up to 10 carbons and 608 basis functions) and water clusters (up to 49 molecules and 2842 basis functions) in a triple-zeta basis set.