One candidate-model for heating the solar corona is magnetic reconnection that embodies Ohmic dissipation of current sheets. When numerous small-scaled magnetic reconnection occur, then it is possible to heat the corona. Due to the limitations of current instrumentation, nanoflares cannot be resolved. But their importance is evaluated via statistics by finding the power-law index of the energy distribution. This method is however biased due to technical and physical reasons. We aim to overcome limitations imposed by observations and statistical analysis. This way, we will identify, and study the small scale impulsive events. We employ a three-dimensional magnetohydrodynamic (3D-MHD) simulation using the \bifrost code. We also employ a new technique to identify the evolution of 3D Joule heating events in the corona. Then, we derive parameters describing the heating events in these locations. We report on the identification of heating events. We obtain the distribution of duration, released energy, and volume. We also find weak power-law correlation between these parameters. In addition, we extract information about geometrical parameters of 2D slices of 3D events, and about the evolution of resolved Joule heating compared to the total Joule heating and the magnetic energy in the corona. Even though the energy power index is less than 2, when classifying the energy release into three categories regarding with respect to the energy release (pico-, nano-, and micro-events), we find that nano-events release $82 \ \%$ of the resolved energy. This fraction corresponds to an energy flux larger than the one needed to heat the corona. It seems that the most popular population is the one containing short-lived, with small spatial extend events.

The solar coronal heating problem is an open question since 1939. One proposed model for the transport and release of mechanical energy generated in the sub-phorospheric layers and photosphere is the nanoflare model that incorporates Ohmic heating which releases a part of the energy stored in the magnetic field via magnetic reconnection. The problem with the verification of this model is that we cannot resolve observationally small scale events. Histograms of observable characteristics of flares, show powerlaw behavior, for both energy release rate, size and total energy. Depending on the powerlaw index of the energy release, nanoflares might be an important candidate for coronal heating; we seek to find that index. In this paper, we employ a numerical 3D-MHD simulation produced by the numerical code Bifrost, and a new technique to identify the 3D heating events at a specific instant. The quantity we explore is the Joule heating, which is explicitly correlated with the magnetic reconnection because depends on the curl of the magnetic field. We are able to identify 4136 events in a volume $24 \times 24 \times 9.5 \ \textrm{Mm}^3$ (i.e. $768 \times 786 \times 331$ grid cells) of a specific snapshot. We find a powerlaw slope of the released energy per second, and two powerlaw slopes of the identified volume. The identified energy events do not represent all the released energy, but of the identified events, the total energy of the largest events dominate the energy release. Most of the energy release happens in the lower corona, while heating drops with height. We find that with a specific identification method that large events can be resolved into smaller ones, but at the expense of the total identified energy releases. The energy release which cannot be identified as an event favours a low energy release mechanism.