Maxima of the scalar dissipation rate in turbulence appear in form of sheets and correspond to the potentially most intensive scalar mixing events. Their cross-section extension determines a locally varying diffusion scale of the mixing process and extends the classical Batchelor picture of one mean diffusion scale. The distribution of the local diffusion scales is analysed for different Reynolds and Schmidt numbers with a fast multiscale technique applied to very high-resolution simulation data. The scales take always values across the whole Batchelor range and beyond. Furthermore, their distribution is traced back to the distribution of the contractive short-time Lyapunov exponent of the flow. When a scalar concentration field θ(x, t) is transported in a turbulent flow very large scalar gradients are generated which can be associated with potentially intensive mixing [1]. Frequently, such large-amplitude gradient regions exist across scales that are finer than the smallest turbulent eddies, a case which is known as the Batchelor regime of scalar mixing [2]. Their cross-section is usually estimated by a single mean diffusion scale, the Batchelor scale η B , that equilibrates advection by flow and scalar diffusion. However, scalar gradients are known to fluctuate strongly in turbulent mixing. These fluctuations are caused by the fluctuating scalar amplitudes and by the varying spatial sections across which the scalar differences are built up. Both aspects will cause the strong spatial variability of potentially intensive mixing regions. Thus, a whole range of local diffusion scales l d , which quantifies exactly this variability, will exist around the mean diffusion scale η B . Their distribution is interesting for several reasons. It can enter mixing efficiency measures [3]. The scales prescribe the extension of chemically reactive layers in which the combustion of fuel takes place [4] or the variations of the fluorescence signal as used for the measurement of zooplankton patchiness in the ocean [5].In this Letter, we want to calculate this local diffusion scale distribution p(l d ). It arises from the competition of two dynamic processes. On one hand, the scale distribution will be affected by molecular diffusion that causes a diminishing of existing steep gradients as well as the completion of their formation by reconnection [6,7]. On the other hand, the scales will be determined by the statistics of the local advection flow patterns that pile up scalar differences. While the strongest scalar gradients were related to instantaneous velocity gradients in Ref.[8], we will go one step further in the second part of the Letter and relate the local diffusion scale distribution to the distribution of Lagrangian contraction rates of the flow, as given by the smallest of the three finite-time Lyapunov exponents along Lagrangian trajectories. The Lagrangian approach incorporates the temporal evolution of the persistent flow patterns that eventually generate the strongest scalar gradients or the finest local dissipation...