We study the spreading of initially-local operators under unitary time evolution in a onedimensional random quantum circuit model which is constrained to conserve a U (1) charge and also the dipole moment of this charge. These constraints are motivated by the quantum dynamics of fracton phases. We discover that charge remains localized at its initial position, providing a crisp example of a non-ergodic dynamical phase of random circuit dynamics. This localization can be understood as a consequence of the return properties of low dimensional random walks, through a mechanism reminiscent of weak localization, but insensitive to dephasing. The charge dynamics is well-described by a system of coupled hydrodynamic equations, which makes several non-trivial predictions which are all in good agreement with numerics in one dimension. Importantly, these equations also predict localization in two-dimensional fractonic random circuits. We further find that the immobile fractonic charge emits non-conserved operators, whose spreading is governed by exponents distinct from those observed in non-fractonic circuits. These non-standard exponents are also explained by our coupled hydrodynamic equations. Where entanglement properties are concerned, we find that fractonic operators exhibit a short time linear growth of observable entanglement with saturation to an area law, as well as a subthermal volume law for operator entanglement. The entanglement spectrum is found to follow semi-Poisson statistics, similar to eigenstates of many-body localized systems. The non-ergodic phenomenology is found to persist to initial conditions containing non-zero density of dipolar or fractonic charge, including states near the sector of maximal charge. Our work implies that low-dimensional fracton systems should preserve forever a memory of their initial conditions in local observables under noisy quantum dynamics, thereby constituting ideal memories. It also implies that one-and two-dimensional fracton systems should realize true many-body localization (MBL) under Hamiltonian dynamics, even in the absence of disorder, with the obstructions to MBL in translation invariant systems and in spatial dimensions greater than one being evaded by the nature of the mechanism responsible for localization. We also suggest a possible route to new non-ergodic phases in high dimensions.