We find a systematic reappearance of massive Dirac features at the edges of consecutive minibands formed at magnetic fields B p/q = p q φ0/S providing rational magnetic flux through a unit cell of the moiré superlattice created by a hexagonal substrate for electrons in graphene. The Diractype features in the minibands at B = B p/q determine a hierarchy of gaps in the surrounding fractal spectrum, and show that these minibands have topological insulator properties. Using the additional q-fold degeneracy of magnetic minibands at B p/q , we trace the hierarchy of the gaps to their manifestation in the form of incompressible states upon variation of the carrier density and magnetic field.PACS numbers: 73.22.Pr, 73.21.Cd, The fractal spectrum of electron waves in crystals subjected to a strong magnetic field [1, 2] is a fundamental result in the quantum theory of solids [3]. In particular, the bandstructure of electrons on a two-dimensional (2D) lattice is fractured into multiple bands which occur at the values B p q = p q φ 0 /S of the field providing a rational fraction of magnetic flux quantum φ 0 = h/e per unit cell area S [1, 3]. Its image [4], obtained for a square lattice hopping model, known as the Hofstadter butterfly, has stimulated numerous attempts to observe the fractal spectrum of electrons via quantum transport measurements. Since the sparsity of the spectrum increases for larger values of the denominator q in p q (hence smaller gaps), the observation of fractal magnetic bands in real crystals would require unsustainably strong magnetic fields. Hence, early efforts were focused on 2D electrons in periodically patterned GaAs/AlGaAS heterostructures [5], where the superimposed superlattice period was made large enough to obtain low denominator fractions within the experimentally available steady magnetic field range. The more recent observation of moiré superlattices, both for twisted bilayer graphene [6] and graphene residing on substrates with hexagonal facets [7], has shown an alternative way to create a long-range periodic potential for electrons: by making lattice-aligned graphene heterostructures using a hexagonal crystal with an almost commensurate period. For this, hexagonal boron nitride (hBN, with the lattice constant of l hBN = 2.50Å vs l = 2.46Å in graphene) provides a perfect match. Several observations of moiré superlattice effects in graphene-hBN heterostructures have already been reported [8][9][10].In this article, we study magnetic minibands caused by a moiré superlattice in graphene. When the surface layer of the substrate is inversion symmetric, the zeromagnetic-field spectra displays clear secondary Dirac points at the first miniband edge [11][12][13]. We find that generations of massive Dirac electrons systematically reappear at the edges of magnetic minibands at2 ), unit cell area S = √ 3 2 a 2 ) and effective perturbations for electrons in graphene on an hexagonal substrate [θ = 0, and V + = 0.063vb,q , and the surrounding fractal spectrum groups around q-fold degenerate Landau le...