Introducing thresholds to analyze time series of emission from the Sun enables a new and simple definition of solar flare events, and their interoccurrence times. Rescaling time by the rate of events, the waiting and quiet time distributions both conform to scaling functions that are independent of the intensity threshold over a wide range. The scaling functions are well described by a two parameter function, with parameters that depend on the phase of the solar cycle. For flares identified according to the current, standard definition, similar behavior is found.PACS numbers: 96.60. Rd, 05.45.Tp, 05.65.+b The solar corona is a very high Reynolds number turbulent plasma producing intermittent bursts of radiation. Plasma forces twist the coronal magnetic fields until stresses are suddenly released, an avalanching process governed by magnetic reconnection [1]. The released magnetic energy induces radiative emission that can be detected as a flare. Flares exhibit scale invariant statistics. For instance, the probability distribution of flare energies is a power law spanning more than eight orders of magnitude [2,3], similar to the Gutenberg-Richter law for earthquakes. The distribution of magnetic concentration sizes on the photosphere is also scale invariant, and the coronal magnetic network embodies a scale-free network [4,5]. In fact, a model of self-organized criticality (SOC) with avalanches of reconnecting flux tubes reproduces the observed scale-free network structure [4,6].As part of the debate on the characterization of magnetohydrodynamic turbulence in this regime [1,4,5,6,7,8,9], interest has focused on comparing interoccurrence times between flares with those in models of SOC. Analyses of flare catalogs have indicated scale invariance of the waiting times, but the behavior was found to vary with the phase of the solar cycle [10] and with the methods used to analyze the catalogs. (See e.g. Ref. [10,11].) The prior belief that avalanches occur with Poissonian waiting times in the well-known BTW sandpile model [12] (giving an exponential distribution of waiting times) argued against the SOC hypothesis [8]. However, including a finite detection threshold leads to a power law distribution of quiet times even for the BTW model [13], when durations and quiet times are measured with the same clock. Since the turnover time scale for flux to be regenerated in the corona is of the order of ten hours [14], while the correlated waiting time intervals between flares can extend up to years, the physical mechanism(s) responsible for these correlations resides in the turbulent convective region beneath the photosphere that generates magnetic flux and drives it into the corona. Systematic studies of the temporal pattern of flares can give insight into the dynamics of magnetic flux in the convective region, or the solar dynamo, which is difficult to observe directly.Here we show that the interoccurrence times between flares has a hierarchical scaling structure when flares are defined as intervals during which the emissi...