Properties of atomic van der Waals heterostructures are profoundly influenced by interlayer coupling, which critically depends on stacking of the proximal layers. Rotational misalignment or lattice mismatch of the layers gives rise to a periodic modulation of the stacking, the moiré superlattice. Provided the superlattice period extends over many unit cells, the coupled layers undergo lattice relaxation, leading to the concentration of strain at line defects - solitons - separating large area commensurate domains. We visualize such long-range periodic superstructures in thin crystals of hexagonal boron nitride using atomic-force microscopy and nano-infrared spectroscopy. The solitons form sub-surface hexagonal networks with periods of a few hundred nanometers. We analyze the topography and infrared contrast of these networks to obtain spatial distribution of local strain and its effect on the infrared-active phonons of hBN.