Recent fully nonlinear, kinetic three-dimensional simulations of magnetic reconnection [1] evolve structures and exhibit dynamics on multiple scales, in a manner reminiscent of turbulence. These simulations of reconnection are among the first to be performed at sufficient spatio-temporal resolution to allow formal quantitative analysis of statistical scaling which we present here. We find that the magnetic field fluctuations generated by reconnection are anisotropic, have non-trivial spatial correlation and exhibit the hallmarks of finite range fluid turbulence; they have non-Gaussian distributions, exhibit Extended Self-Similarity in their scaling and are spatially multifractal. Furthermore, we find that the field J · E is also multifractal, so that magnetic energy is converted to plasma kinetic energy in a manner that is spatially intermittent. This suggests that dissipation in this sense in collisionless reconnection on kinetic scales has an analogue in fluid-like turbulent phenomenology, in that it proceeds via multifractal structures generated by an intermittent cascade. Magnetic reconnection is a fundamental process that converts magnetic energy into various forms of plasma kinetic energy. It is thought to occur in a variety of space, astrophysical and laboratory applications, with parameter regimes spanning from collisional to highly collisionless plasmas (e.g. see Ref.[2] and references therein). While many studies have focused on laminar initial conditions, it is now widely recognized that the influence of turbulence remains a major uncertainty. Depending on the application, the turbulence may arise from a spectrum of instabilities within the reconnection layer or from pre-existing magnetic fluctuations in the ambient plasma. Within the magnetohydrodynamic (MHD) model, there has been progress on both ideas -either by starting with a laminar current sheet to explore instabilities [3][4][5][6] or by directly driving turbulence [7][8][9][10][11].Moving beyond the MHD model into kinetic regimes, most research has focused on initially laminar current sheets within a variety of descriptions [12] including two fluid, hybrid and fully kinetic simulations, which allow a complete description of the electron physics responsible for breaking the frozen-flux condition in collisionless parameter regimes [13,14]. As larger kinetic simulations became possible, one surprising result was that the nonlinear evolution of reconnection produced extended electron-scale current sheets, with half-thickness on the order of the electron inertial length and lengths that can extend beyond the ion inertial scale [15][16][17][18][19]. These predictions have since been confirmed in spacecraft observations [20]. While the precise details depend on the strength of the guide field (i.e. magnetic shear angle), it has been demonstrated that electron pressure anisotropy plays a key role in setting up and driving these layers [21]. Thus the existence of these structures is now well accepted and a variety of two-dimensional (2D) kinetic simul...