We compute the topological susceptibility χ t of lattice QCD with 2 + 1 dynamical quark flavors described by the Möbius domain wall fermion. Violation of chiral symmetry as measured by the residual mass is kept at ∼1 MeV or smaller. We measure the fluctuation of the topological charge density in a "slab" sub-volume of the simulated lattice using the method proposed by Bietenholz et al. The quark mass dependence of χ t is consistent with the prediction of chiral perturbation theory, from which the chiral condensate is extracted as Σ MS (2GeV) = [274(13)(29)MeV] 3 , where the first error is statistical and the second one is systematic. Combining the results for the pion mass M π and decay constant F π , we obtain χ t = 0.229(03)(13)M 2 π F 2 π at the physical point.