<jats:p>The present study addresses the reaction zone structure and burning mechanism of unstable detonations. Experiments investigated mainly two-dimensional methane–oxygen cellular detonations in a thin channel geometry. The sufficiently high temporal resolution permitted the determination of the probability density function of the shock distribution, a power law with an exponent of<jats:inline-formula><jats:alternatives><jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="gif" xlink:type="simple" xlink:href="S0022112017001458_inline1" /><jats:tex-math>$-3$</jats:tex-math></jats:alternatives></jats:inline-formula>, and the burning rate of unreacted pockets from their edges – through surface turbulent flames with a speed approximately 3–7 times larger than the laminar one at the local conditions. Numerical simulations were performed using a novel large-eddy simulation method where the reactions due to both autoignition and turbulent transport were treated exactly at the subgrid scale in a reaction–diffusion formulation. The model is an extension of Kerstein and Menon’s linear eddy model for large-eddy simulation to treat flows with shock waves and rapid gas-dynamic transients. The two-dimensional simulations recovered well the amplification of the laminar flame speed due to the turbulence generated mainly by the shear layers originating from the triple points and subsequent Richtmyer–Meshkov instability associated with the internal pressure waves. The simulations clarified how the level of turbulence generated controlled the burning rate of the pockets, the hydrodynamic thickness of the wave, the cellular structure and its distribution. Three-dimensional simulations were found to be in general good agreement with the two-dimensional ones, in that the subgrid-scale model captured the ensuing turbulent burning once the scales associated with the cellular dynamics, where turbulent kinetic energy is injected, are well resolved.</jats:p>