Issue 
A&A
Volume 624, April 2019



Article Number  A65  
Number of page(s)  11  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201834455  
Published online  11 April 2019 
Scaling laws of quasiperiodic pulsations in solar flares
^{1}
Department of Physics, University of Warwick, Coventry CV4 7AL, UK
email: am.broomhall@warwick.ac.uk
^{2}
St. Petersburg Branch, Special Astrophysical Observatory, Russian Academy of Sciences, 196140 St. Petersburg, Russia
Received:
18
October
2018
Accepted:
25
February
2019
Context. Quasiperiodic pulsations (QPPs) are a common feature of solar flares, but there has previously been a lack of observational evidence to support any of the theoretical models that might explain the origin of these QPPs.
Aims. We aimed to determine if there are any relationships between the QPP period and other properties of the flaring region, using a previously assembled sample of flares with QPPs. If any relationships exist, then these can be compared with scaling laws for the theoretical QPP mechanisms.
Methods. To obtain the flaring region properties, we made use of the Atmospheric Imaging Assembly (AIA) 1600 Å and Helioseismic and Magnetic Imager (HMI) data. The flare ribbons are visible in AIA 1600 Å images, and the positive and negative magnetic polarity ribbons can be distinguished and the magnetic properties determined in the HMI magnetograms. The ribbon properties calculated in this study were the ribbon separation distance, area, total unsigned magnetic flux, and average magnetic field strength. Only the flares that occurred within ±60° of the solar disc centre were included, which meant a sample of 20 flares with 22 QPP signals.
Results. Positive correlations were found between the QPP period and the ribbon properties. The strongest correlations were with the separation distance and magnetic flux. Because these ribbon properties also correlate with the flare duration and because the relationship between the QPP period and flare duration may be influenced by observational bias, we also made use of simulated data to determine whether artificial correlations were introduced. These simulations show that although QPPs cannot be detected for certain combinations of QPP period and flare duration, this does not introduce an apparent correlation.
Conclusions. There is evidence of relationships between the QPP period and flare ribbon properties, and in the future, the derived scaling laws between these properties can be compared to equivalent scaling laws for theoretical QPP mechanisms.
Key words: Sun: activity / Sun: flares / Sun: oscillations / Sun: radio radiation / Sun: UV radiation / Sun: Xrays, gamma rays
© ESO 2019
1. Introduction
Quasiperiodic pulsations (QPPs) are frequently observed in solar and stellar flare light curves (recent publications include Brosius et al. 2016; Hayes et al. 2017; Li et al. 2017; Kolotkov et al. 2018; Doyle et al. 2018; Jackman et al. 2019) and remain to be fully understood. Several different mechanisms for the generation of QPP signals have been proposed (for recent reviews, see Nakariakov et al. 2016; Van Doorsselaere et al. 2016; McLaughlin et al. 2018), but there has previously been a lack of observational evidence to favour one or more of these mechanisms above the others. Limitations of the spatial and temporal resolution of solar flare data, along with saturation effects that many instruments suffer from during flares, mean that although a QPP signal may be clearly visible in Sunasastar data, plasma motions that might be associated with the QPPs are not usually resolved in the imaging observations. Despite the lack of spatial information on the QPPs themselves, the highquality data from instruments that have been observing the Sun over the past few years have allowed systematic statistical studies of QPPs in solar flares to be made (Simões et al. 2015; Inglis et al. 2016; Pugh et al. 2017b; Dominique et al. 2018). By examining large populations, the general properties of QPPs can be determined, and these can be compared with the properties of QPPs that are produced by the various theoretical mechanisms.
An example of observational QPP characteristics that can be compared with theoretical mechanisms are scaling laws between the QPP and flaring region properties. For the QPP mechanisms based on standing magnetohydrodynamic (MHD) oscillations of the flaring coronal loop, the QPP period is expected to scale linearly, with the loop length based on the theory of MHD modes of a straight cylinder (Nakariakov & Melnikov 2009). In some limits (such as long wavelengths), this relationship might break down, however (Nakariakov et al. 2012). For other mechanisms, no linear relationship between the QPP period and loop length may exist. For example, simulations performed by Ruan et al. (2018), based on the mechanism proposed by Fang et al. (2016), showed oscillations with a period of 25 s in the soft Xray light curve of the flare. Based on this period, the authors suggested that the oscillations were most likely the result of a fast sausage mode excited by motions associated with a KelvinHelmholtz instability. They found that if the sausage mode were present in the whole flaring loop, then the phase speed of the mode would be greater than the external Alfvén speed, which is not consistent with MHD wave theory. Therefore they proposed that the waves may be reflected before reaching the loop footpoints as a result of the change in loop cross section, and hence the associated characteristic spatial scale may be smaller than the loop length.
A linear relationship between the QPP period and an associated characteristic spatial scale is expected for some mechanisms that are not based on standing MHD oscillations (e.g. Takasao & Shibata 2016; Takahashi et al. 2017), but these spatial scales may not have a linear relationship with observable spatial scales such as the loop length. Takasao & Shibata (2016) found that their “magnetic tuning fork” mechanism produced a period that was related to the magnetic field strength as P ∝ B^{−2.1}, and compared this to the expected scaling of P ∝ B^{−0.43} for QPPs produced by slow magnetoacoustic oscillations of the flaring loop. Although it is not currently possible to directly measure the magnetic field in the corona, the magnetic field strength inside a coronal loop might be expected to scale with the magnetic field measured at the loop footpoints in the photosphere.
Other mechanisms exist that would result in a QPP period independent of the flaring region properties (e.g. Nakariakov et al. 2006), or that do not yet make testable predications of relationships between flaring region properties and the QPP period (e.g. Thurgood et al. 2017).
This study continues the work of Pugh et al. (2017b), making use of the same sample of flares with QPPs. While Pugh et al. (2017b) made the first attempt at linking QPP periods with spatial scales, the present work improves on this by focussing on the spatial scales of the flaring regions, rather than on the spatial scales of the active region as a whole. In addition, we use simulated data to explore the potential for observational bias to influence correlations found between the QPP and flare properties. Details of the data used are given in Sect. 2, while the analysis methods are described in Sect. 3. Section 4 gives the results of the search for relationships between the QPP and flare properties, along with the determination for how observational bias can affect the QPP periods that are detected for different flare durations. Finally, conclusions are given in Sect. 5.
2. Observations
The sample used for this study is made up of flares from the longlived active region NOAA 12172/12192/12209, which was present on the Sun between September and November 2014. During this time interval, the Sun was well observed, therefore data could be taken from multiple instruments in order to give the best opportunity of detecting QPP signals. Xray timeseries data were taken from the Geostationary Operational Environmental Satellite (GOES) Xray sensor (XRS), the Extreme ultraviolet SpectroPhotometer (ESP) channel of the Extreme ultraviolet Variability Experiment (EVE) on board the Solar Dynamics Observatory (SDO), the Fermi Gammaray Burst Monitor (GBM), and the Detector of the Roentgen and Gammaray Emissions (DRGE) instrument on board Vernov. In addition, microwave timeseries data from Nobeyama Radioheliograph (NoRH) were used. Characteristics of the data from these instruments, along with details of how data uncertainties were estimated for those instruments whose data did not include uncertainties, are described in Pugh et al. (2017b).
Spatial properties of the flare ribbons were determined using data from the Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI), both on board SDO. From AIA (Lemen et al. 2012), the 1600 Å images were used to identify the flare ribbons, while the HMI (Scherrer et al. 2012) lineofsight magnetogram images were used to distinguish the positive and negative magnetic polarity flare ribbons, and estimate magnetic fluxes and field strengths in the region. AIA data cubes were constructed with a cadence of 24 s, and the time interval used covered the impulsive phase of the flare (i.e., from the GOES flare start time to the GOES peak time). The AIA images can be subject to blooming when the CCD is saturated, which often occurs during flares, therefore frames with visible blooming were removed. For the HMI images, a single frame was used for each flare because there was little change in the magnetograms over the course of the impulsive phases of the flare, and these frames were taken at about the times of the flare peaks. All images were used with the full resolution of 4096 × 4096 pixels.
3. Data analysis
The sample of flares with QPPs used in this study was taken from Pugh et al. (2017b), who found that 37 out of 181 flares had a QPP signal with a peak in the power spectrum above the 95% global confidence level. The method used to detect the QPPs did not require any detrending of the data, and involved calculating power spectrum confidence levels that accounted for data uncertainties and the presence of red noise. Further details of the method can be found in Pugh et al. (2017a), and an example of one of the flares with detected QPPs is shown in Fig. 1.
Fig. 1. Left: time derivatives of a section of flare 092 (see Table 1) observed by GOES/XRS. The top panel shows the 1–8 Å emission and the bottom panel the 0.5–4 Å emission. Right: corresponding power spectra, where the red solid lines are broken powerlaw fits to the spectra, the red dotted lines represent the 95% confidence levels, and the red dashed lines the 99% levels. There is a significant peak in both power spectra at a period of . 
Because lineofsight effects must be accounted for when measuring spatial properties of the flares, only flares that occurred within around ±60° of the solar disc could be included in this study. This gave a sample of 20 flares with 22 QPP signals, and the properties of this sample are given in Table 1.
3.1. Flare ribbon properties
The flare ribbon properties were determined using a similar approach as in Toriumi et al. (2017), using the AIA 1600 Å and HMI lineofsight magnetogram data. For all AIA and HMI frames, the lineofsight projection effect of features closer to the solar limb that appeared to be smaller than when closer to the disc centre was first accounted for. This was done by differentially rotating the frames to the time when the active region was approximately at the central meridian, using the IDL Solar Software routine drot_map. These times were taken to be 2014 September 26 18:30:16 UT for the active region’s first crossing of the solar disc as NOAA 12172, 2014 October 23 14:00:40 UT for the second crossing as NOAA 12192, and 2014 November 19 15:00:40 for the third crossing as NOAA 12209. Next, the images were cropped so that only the region containing the active region and the full extent of all flare ribbons was included. The coordinates used for the cropped regions, in the form [x_{0}, x_{1}, y_{0}, y_{1}] and in the units of arcseconds for the active region’s three crossings of the solar disc, are [−250, 250, −432, −147], [−334, 334, −568, −120], and [−551, 551, −541, 16], respectively.
To identify the flare ribbons, a threshold brightness for the AIA 1600 Å data was first determined. A period of time when no active regions and very few magnetic features were present on the solar disc was chosen, between 2018 March 28 21:00:00 and 21:30:00 UT. After averaging the AIA 1600 Å frames during this period and cropping the frames to the region [−550, 550, −550, 16], similar to the regions defined above, the mean and standard deviation were calculated. Based on these values, the brightness threshold was defined as the mean plus 40 times the standard deviation. For each AIA 1600 Å frame of each flare, flare ribbon pixels were defined as those with values greater than this threshold. Following the method of Toriumi et al. (2017), a ribbon composite was then formed by selecting the pixels that exceeded the brightness threshold in any of the AIA 1600 Å frames. This reduces the possibility that background noise is included as part of the ribbons because a high brightness threshold is used, while ensuring that the full extent of the flare ribbons is represented. Although there is risk of overestimating the ribbon area using this approach if the ribbon separation grows with time, the flares from the active region used in this study were mostly confined and therefore do not exhibit much growth of the ribbon separation distance (Thalmann et al. 2015, Sun et al. 2015).
To separate the positive and negative magnetic polarity flare ribbons, for each flare the ribbon pixels were overlaid on the corresponding HMI magnetogram. The ribbons were then separated based on whether they were above a positive or negativepolarity region in the HMI data. Figure 2 shows an example of how the flare ribbons appear in the AIA 1600 Å data for one of the flares in Table 1, along with the corresponding HMI magnetogram and the composite flare ribbon contours. HMI magnetograms and flare ribbon contours for all flares in the sample are shown in Fig. A.1.
Fig. 2. Left: AIA 1600 Å images averaged over the impulsive phase of flare 092, showing the full extent of the flare ribbons. Right: HMI magnetogram showing the host active region at the time of the flare peak in the GOES 1–8 Å waveband. The red and blue contours show the positions of the composite flare ribbons with positive and negative magnetic polarity, respectively. The green line joins the geometric centroids of the positive and negativepolarity ribbon components. 
Flares with QPPs used in this study.
To calculate the flare ribbon areas, S_{ribbon}, for each of the ribbon composite pixels the area on the solar surface that the pixel would correspond to if it were located at the disc centre was multiplied by a cosine correction factor, to account for the fact that different pixels correspond to different surface areas as a result of the spherical nature of the Sun. These corrected pixel areas were then summed to obtain the total area of the flare ribbons. The ribbon separation distance, d_{ribbon}, was estimated from the great circle distance between the geometric centroids of the positive and negativepolarity components of the flare ribbons. The total unsigned magnetic flux, Φ_{ribbon}, below the ribbons was defined as Φ_{ribbon} = ∫_{Sribbon}B ds, and the average magnetic field strength as B_{ribbon} = Φ_{ribbon}/S_{ribbon}. The calculated flare ribbon properties are given in Table B.1.
3.2. Simulated flare light curves
In order to explore the potential for observational bias to result in apparent relationships between the QPP periods and flare properties in the observations, a set of simulated flare light curves was generated. For these simulated flares, 16 flare durations, QPP periods, and QPP durations were used, which were all uniformly distributed in log space. This meant that a total of 4096 simulated flares were created, each with a different combination of the 16 QPP periods, QPP durations, and flare durations. The flares were given a time cadence of 1 s, with log durations evenly spaced between 1.300 and 3.505 (or between around 20 and 3200 s in linear space), log QPP periods between 0.7 and 1.9 (or 5 and 79 s in linear space), and log QPP durations between 1.0 and 2.5 (or 10 and 316 s in linear space). These parameter spaces were chosen to be similar to the observed flare and QPP parameters reported by Pugh et al. (2017b). For the basic flare time profile, we made use of the following expression from Gryciuk et al. (2017):
where
I(t) is the flare intensity as a function of time, and A, B, C, and D are arbitrary parameters. For this set of simulated flares, these arbitrary parameters were set to A = 1, B = t_{flare}/15, C = t_{flare}/10, and D = 3/t_{flare}.
The QPP signal was modelled as a Gaussian modulated sinusoid, as described by the following expression:
where I_{QPP}(t) is the QPP signal intensity as a function of time, A_{QPP} is the signal amplitude (set to 0.1 for all flares), P is the period, t_{0} is set to be the flare peak time, and τ_{G} is the Gaussian decay time.
After normalising Eq. (1) by the peak intensity of the flare, so that the resulting flare peak intensity was 1, it was summed with Eq. (3) and normally distributed white noise with a standard deviation of 0.02. One of the simulated flares is shown in Fig. 3 as an example.
Fig. 3. One of the simulated flare light curves, with a flare duration of 213 s, a QPP period of 12 s, and a QPP decay time of 100 s. 
After we generated the set of simulated flares, we attempted to detect the QPPs using a simplified version of the method used in Pugh et al. (2017b) because there were no uncertainties on the simulated data. For each flare, a power spectrum was first generated, and a broken powerlaw model with the following form was fitted to the spectrum:
where is the spectral power of the fit as a function of frequency, f_{break} is the frequency at which the powerlaw break occurs, α and β are powerlaw indices, and c is a constant.
Next the confidence levels were calculated. Noise in the power spectrum follows a chisquared distribution with two degrees of freedom (d.o.f.), and the corresponding probability density is
Therefore the probability of having a value Z in the power spectrum that is greater than some threshold z is
In order to calculate the “global” confidence level, the number of data points in the spectrum must be accounted for. For a power spectrum sampled at N independent frequencies, the probability is equivalent to
where ϵ_{N} is the falsealarm probability, and the approximation holds when is ϵ_{N} small (Chaplin et al. 2002). Equating Eqs. (6) and 7 and rearranging gives
Finally, the fact that the power spectrum is not normalised must be accounted for because the above expression will only give the confidence level for a power spectrum dominated by white noise and with a mean value of 1. For the flare power spectra, the noise is distributed around the broken power law, so the confidence level can be found from , where J_{j} is the observed spectral power at frequency f_{j}. The 95% global confidence level (corresponding to ϵ_{N} = 0.05) was used as the detection threshold for the QPP signals. In addition, only significant spectral peaks with periods greater than four times the time cadence and less than a quarter of the duration of the time series were included because outside of this range, it is unclear whether a periodic signal can be reliably detected, and the same constraint was applied in Pugh et al. (2017b).
4. Results and discussion
4.1. Comparing QPP and flare properties
Figures 4–7 show scatter plots of the detected QPP period, P, with the flare ribbon separation distance d_{ribbon}, total unsigned magnetic flux Φ_{ribbon}, average magnetic field strength B_{ribbon}, and area S_{ribbon}, respectively. All plots show evidence of a positive correlation, and powerlaw fits give the following relationships:
Fig. 4. QPP periods plotted against the separation of the geometric centroids of the two flare ribbons. The differently coloured points correspond to the different instruments used to detect the QPP signals, and the grey line is a linear fit in log space. The Pearson correlation is 0.53 with a pvalue of 0.01, and the Spearman correlation is 0.64 with a pvalue of 0.001. 
Fig. 5. QPP periods plotted against the total unsigned magnetic flux of the photosphere below the flare ribbons. The Pearson correlation is 0.59 with a pvalue of 0.004, and the Spearman correlation is 0.68 with a pvalue of 0.0005. 
Fig. 6. QPP periods plotted against the average magnetic field strength of the photosphere below the flare ribbons. The Pearson correlation is 0.47 with a pvalue of 0.03, and the Spearman correlation is 0.50 with a pvalue of 0.02. 
Fig. 7. QPP periods plotted against the area of the flare ribbon. The Pearson correlation is 0.58 with a pvalue of 0.005, and the Spearman correlation is 0.59 with a pvalue of 0.004. 
These relationships are all nonlinear, therefore the Spearman rank correlation coefficient is a more suitable measure of the correlation than the Pearson correlation coefficient if the properties are considered in linear space. Based on this, the strongest correlations are with the total unsigned magnetic flux (with a Spearman correlation of 0.68) and the ribbon separation distance (with a Spearman correlation of 0.64). If the Pearson correlation coefficients are calculated based on the logarithm of the properties, then values of 0.68, 0.73, 0.57, and 0.65 are obtained for Figs. 4–7, respectively.
Pugh et al. (2017b) also found a correlation between the QPP period and the total duration of the flares as measured in the GOES 1–8 Å waveband. Because it is sometimes difficult to determine exactly when a flare ends in the GOES/XRS data, a more reliable estimate of the flare duration may be the time between when the GOES 1–8 Å flux starts to increase and when the 1–8 Å flux reaches its peak, which roughly corresponds to the impulsive phase of the flare. Therefore plots of the QPP period against the impulsive phase duration of the flare, t_{flare}, are shown in Fig. 8, where the left plot is for the full sample of flares with QPPs identified by Pugh et al. (2017b) and the right plot is for the subset of flares used in this study. These still show positive correlations, and powerlaw fits to the two versions of the QPP period versus flare duration plots gave the following relationships:
Fig. 8. Two variants of the QPP period plotted against the flare impulsive phase duration. Left: full set of 37 flares with QPPs found by Pugh et al. (2017b). The Pearson correlation is 0.53 with a pvalue of 0.0002, and the Spearman correlation is 0.50 with a pvalue of 0.0006. Right: plot includes only those 20 flares that occurred within ±60° of the solar disc centre, and hence is the same sample of flares used in Figs. 4–7. The Pearson correlation coefficient is 0.56 with a pvalue of 0.006, and the Spearman correlation coefficient is 0.65 with a pvalue of 0.001. 
There is potential for observational bias to limit which QPP signals can be detected, as mentioned by Pugh et al. (2017b), because for longperiod QPPs in a shortduration flare there would not be a sufficient number of flare intensity pulses present for the signal to have a significant peak in the power spectrum. Toriumi et al. (2017) showed that the flare ribbon properties correlate strongly with the flare duration, hence if there is potential for observational bias to affect the relationship between the QPP period and flare duration, then there is also potential for the relationships between the QPP period and ribbon properties to be affected. Figure 9 shows plots similar to those presented by Toriumi et al. (2017), but for the flare sample used in this study, and powerlaw fits to the plots give the following relationships:
Fig. 9. Flare impulsive phase durations plotted against the flare ribbon separation distance (left), total magnetic flux (middle), and area (right). The grey lines show linear fits in log space. Left: the Pearson correlation is 0.80 with a pvalue of 9 × 10^{−6}, and the Spearman correlation is 0.84 with a pvalue of 1 × 10^{−6}. Middle: the Pearson correlation is 0.93 with a pvalue of 5 × 10^{−10}, and the Spearman correlation is 0.92 with pvalue of 1 × 10^{−9}. Right: the Pearson correlation is 0.83 with a pvalue of 1 × 10^{−6}, and the Spearman correlation is 0.88 with a pvalue of 9 × 10^{−8}. 
Although a different estimate was used for the flare duration, these relationships and the correlation coefficients given in Fig. 9 are consistent with those reported by Toriumi et al. (2017). The results of Toriumi et al. (2017) are also consistent with other studies, which show that the flare duration is correlated with the GOES peak Xray flux (Veronig et al. 2002), and the peak Xray flux is correlated with the flare ribbon area (Kazachenko et al. 2017).
4.2. Checking for observational biases using simulated data
Figure 10 shows results from the first set of simulated data, described in Sect. 3.2. The parameters input into the model flares are represented by the grey points, while the black points represent the parameters of the flares with QPP signals that were detected by the method described in Sect. 3.2. The majority of detections match the uniform in logspace sampling of periods and durations that were input into the simulated flares (i.e., the black points mostly lie above the grey points), but a few inaccurate detections are visible, where the black points are between the grey points. This is because we used the 95% confidence level as the detection threshold. The important aspect of these plots is the distribution of the black points. As expected, observational bias has resulted in triangular shapes to the plots in Fig. 10, with an absence of points in the longperiod, shortduration regions. There is no absence of points in the shortperiod, longduration regions, however, therefore the correlation coefficients are small.
Fig. 10. Simulated flare QPP period plotted against the flare duration (left) and QPP signal duration (right). The grey points represent the input parameters of the simulated flares with QPPs, while the black points represent the simulated flares with QPP signals that were detected using the method described in Sect. 3.2. Left: the Pearson correlation is 0.23 with a pvalue of 0.0, and the Spearman correlation is 0.30 with a pvalue of 4 × 10^{−27}. Right: the Pearson correlation is 0.31 with a pvalue of 0.0, and the Spearman correlation is 0.30 with a pvalue of 3 × 10^{−27}. 
Pugh et al. (2017b) showed a positive correlation between the QPP period and the duration of the QPP signal, defining the duration as the time interval that gave the most significant peak in the power spectrum. Fitting a linear model gave the relationship
where P is the period and τ is the QPP signal duration time. To determine whether this relationship could result in an apparent dependence between the period and flare duration, another set of simulated flares was generated assuming the relationship to be real. For this set of flares, the same ranges of periods and flare durations were used as before, but the QPP duration was instead set to depend on the period (rearranging Eq. (18) to give τ = 1.30P^{1.61}). Because a Gaussian decay time was used for the QPPs in the simulated flares, whereas for the set of real flares the QPP duration was taken to be the time interval that gave the most prominent peak in the power spectrum, here we assumed τ ≈ 2τ_{G}. When we attempt to detect the QPPs in the simulated flares and then produce a period versus flare duration plot for the flares with detected QPPs, we obtain the plot shown in Fig. 11. The correlation coefficients for this plot are fairly small, therefore this is evidence against the possibility that observational bias and a relationship between the QPP period and QPP duration could result in an apparent relationship between the QPP period and flare duration, with no physical basis behind the relationship. Hence the results of analysing these simulated data support the idea that the observed QPP period versus flare duration relationship is, at least in part, a true relationship.
Fig. 11. Simulated flare QPP period plotted against the flare duration for the case where the QPP signal duration is set to depend on the flare duration. The grey points represent the input parameters of the simulated flares with QPPs, while the black points represent the simulated flares with QPP signals that were detected using the method described in Sect. 3.2. The Pearson correlation is 0.18 with a pvalue of 0.1, and the Spearman correlation is 0.21 with a pvalue of 0.07. 
5. Conclusions
We performed a comparative analysis of QPPs and properties of flare ribbons in the longlived superactive region NOAA 12172/12192/12209. An important feature of this active region is that the hosted flares lack a substantial increase in ribbon separation distance, thus making the process of determining ribbon properties more straightforward.
This study has revealed correlations between the QPP period and all four flare ribbon properties, determined using the AIA 1600 Å and HMI data, as well as the duration of the flare impulsive phase. These ribbon properties are the area, separation distance, total unsigned magnetic flux, and average magnetic field strength. Because there is potential for observational bias to affect the relationship observed between the QPP period and flare duration and the flare duration strongly correlates with the flare ribbon properties, tests with simulated flare light curves were performed. These suggest that the correlations cannot be explained by observational bias alone, and therefore support the idea that they are real. In addition, the strongest correlation is between the QPP period and total unsigned magnetic flux in the flare ribbons, not the flare duration.
The results of this study open up the possibility for revealing the mechanisms behind QPPs. The obtained empirical scaling laws of the QPP period with the flare parameters do not seem to match any of the current proposed theoretical mechanisms that have comparable scaling laws, however. For example, on one hand, the positive correlation of the period with the ribbon separation distance is consistent with a similar property of standing kink oscillations (Goddard et al. 2016). On the other hand, the observed increase in period with the magnetic field strength is not typical of kink oscillations; the increase in field strength should lead to an increase in Alfvén and kink speeds, and hence a decrease in period. Likewise, the period of sausage oscillations decreases as the fast magnetoacoustic speed (and hence the magnetic field) increases, while the period is only weakly dependent on the oscillating loop length (which would relate to the ribbon separation) (e.g. Nakariakov et al. 2012). Kink and sausage oscillations should not be disregarded as possible mechanisms for QPPs, however, because both kink and sausage oscillation periods are also determined by the mass of the oscillating loop, which may affect the scaling. In addition, for some promising mechanisms, the relationships between the QPP period and the flaring region parameters have not yet been addressed. This is especially applicable to the mechanisms based on repetitive spontaneous magnetic reconnection (see Van Doorsselaere et al. 2016; McLaughlin et al. 2018 for recent reviews). Thus, the empirically determined scaling presented in this paper demonstrates the need for dedicated theoretical modelling of this scaling for various QPP mechanisms, to allow them to be validated.
In addition, while this study is the first successful attempt at observing these scaling laws, the obtained scaling laws are not yet definitive and some possible shortcomings should be considered. Firstly, a larger sample size should be used to improve the accuracy and precision of the scaling laws. Secondly, because the QPP detection method used is based on the periodogram, only stationary and weakly nonstationary oscillatory signals in the timeseries data could be detected. The effects of strong modulation of the QPP amplitude or period were not accounted for, therefore the detection method is conservative (see Pugh et al. 2017b for a discussion). Another possible shortcoming is the lack of any classification of the QPPs, the need for which has been stressed in Nakariakov et al. (2019). QPPs of different classes may be caused by different mechanisms that correspond to different scaling laws between the period and flare parameters. The identification of different classes of QPPs would allow us to search for scaling laws typical for the different QPP classes, and verify the hypotheses of their mechanisms. Nevertheless, we consider the obtained scaling laws as encouraging findings that should stimulate development of the theoretical models and further observational studies.
Acknowledgments
The authors would like to acknowledge the ISSI International Team “Quasiperiodic Pulsations in Stellar Flares: a Tool for Studying the SolarStellar Connection” led by AMB for many productive discussions, and are grateful to I. N. Myagkova and A. V. Bogomolov for providing the data from Vernov, along with the GOES/XRS, SDO/EVE, SDO/HMI, SDO/AIA, Fermi/GBM, and Nobeyama teams for providing the rest of the data used. AMB acknowledges support from the Royal Society International Exchanges grant IEC\R2\170056, and VMN acknowledges support from the STFC consolidated grant ST/P000320/1 and the Russian Foundation for Basic Research grant No. 175280064 BRICSA. Flare information was provided courtesy of SolarMonitor.org.
References
 Brosius, J. W., Daw, A. N., & Inglis, A. R. 2016, ApJ, 830, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., Elsworth, Y., Isaak, G. R., et al. 2002, MNRAS, 336, 979 [NASA ADS] [CrossRef] [Google Scholar]
 Dominique, M., Zhukov, A. N., Dolla, L., Inglis, A., & Lapenta, G. 2018, Sol. Phys., 293, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Doyle, J. G., Shetye, J., Antonova, A. E., et al. 2018, MNRAS, 475, 2842 [NASA ADS] [CrossRef] [Google Scholar]
 Fang, X., Yuan, D., Xia, C., Van Doorsselaere, T., & Keppens, R. 2016, ApJ, 833, 36 [NASA ADS] [CrossRef] [Google Scholar]
 Goddard, C. R., Nisticò, G., Nakariakov, V. M., & Zimovets, I. V. 2016, A&A, 585, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gryciuk, M., Siarkowski, M., Sylwester, J., et al. 2017, Sol. Phys., 292, 77 [NASA ADS] [CrossRef] [Google Scholar]
 Hayes, L. A., Gallagher, P. T., McCauley, J., et al. 2017, J. Geophys. Res. (Space Phys.), 122, 9841 [NASA ADS] [CrossRef] [Google Scholar]
 Inglis, A. R., Ireland, J., Dennis, B. R., Hayes, L., & Gallagher, P. 2016, ApJ, 833, 284 [NASA ADS] [CrossRef] [Google Scholar]
 Jackman, J. A. G., Wheatley, P. J., Pugh, C. E., et al. 2019, MNRAS, 482, 5553 [NASA ADS] [CrossRef] [Google Scholar]
 Kazachenko, M. D., Lynch, B. J., Welsch, B. T., & Sun, X. 2017, ApJ, 845, 49 [NASA ADS] [CrossRef] [Google Scholar]
 Kolotkov, D. Y., Pugh, C. E., Broomhall, A.M., & Nakariakov, V. M. 2018, ApJ, 858, L3 [NASA ADS] [CrossRef] [Google Scholar]
 Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Li, D., Zhang, Q. M., Huang, Y., Ning, Z. J., & Su, Y. N. 2017, A&A, 597, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 McLaughlin, J. A., Nakariakov, V. M., Dominique, M., Jelínek, P., & Takasao, S. 2018, Space Sci. Rev., 214, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., & Melnikov, V. F. 2009, Space Sci. Rev., 149, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., Foullon, C., Verwichte, E., & Young, N. P. 2006, A&A, 452, 343 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nakariakov, V. M., Hornsey, C., & Melnikov, V. F. 2012, ApJ, 761, 134 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., Pilipenko, V., Heilig, B., et al. 2016, Space Sci. Rev., 200, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., Kolotkov, D. Y., Kupriyanova, E. G., et al. 2019, Plasma Phys. Controlled Fusion, 61, 014024 [NASA ADS] [CrossRef] [Google Scholar]
 Pugh, C. E., Broomhall, A.M., & Nakariakov, V. M. 2017a, A&A, 602, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pugh, C. E., Nakariakov, V. M., Broomhall, A.M., Bogomolov, A. V., & Myagkova, I. N. 2017b, A&A, 608, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ruan, W., Xia, C., & Keppens, R. 2018, A&A, 618, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207 [Google Scholar]
 Simões, P. J. A., Hudson, H. S., & Fletcher, L. 2015, Sol. Phys., 290, 3625 [NASA ADS] [CrossRef] [Google Scholar]
 Sun, X., Bobra, M. G., Hoeksema, J. T., et al. 2015, ApJ, 804, L28 [Google Scholar]
 Takahashi, T., Qiu, J., & Shibata, K. 2017, ApJ, 848, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Takasao, S., & Shibata, K. 2016, ApJ, 823, 150 [NASA ADS] [CrossRef] [Google Scholar]
 Thalmann, J. K., Su, Y., Temmer, M., & Veronig, A. M. 2015, ApJ, 801, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Thurgood, J. O., Pontin, D. I., & McLaughlin, J. A. 2017, ApJ, 844, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Toriumi, S., Schrijver, C. J., Harra, L. K., Hudson, H., & Nagashima, K. 2017, ApJ, 834, 56 [NASA ADS] [CrossRef] [Google Scholar]
 Van Doorsselaere, T., Kupriyanova, E. G., & Yuan, D. 2016, Sol. Phys., 291, 3143 [NASA ADS] [CrossRef] [Google Scholar]
 Veronig, A., Temmer, M., Hanslmeier, A., Otruba, W., & Messerotti, M. 2002, A&A, 382, 1070 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Additional figure
Fig. A.1. HMI magnetograms for all flares listed in Table 1, showing the host active region at the time of the flare peak in the GOES 1–8 Å waveband. The red and blue contours show the positions of the composite flare ribbons with positive and negative magnetic polarity, respectively. The green lines join the geometric centroids of the positive and negativepolarity ribbon components. 
Fig. A.1. continued. 
Appendix B: Additional table
Summary of the flare properties.
All Tables
All Figures
Fig. 1. Left: time derivatives of a section of flare 092 (see Table 1) observed by GOES/XRS. The top panel shows the 1–8 Å emission and the bottom panel the 0.5–4 Å emission. Right: corresponding power spectra, where the red solid lines are broken powerlaw fits to the spectra, the red dotted lines represent the 95% confidence levels, and the red dashed lines the 99% levels. There is a significant peak in both power spectra at a period of . 

In the text 
Fig. 2. Left: AIA 1600 Å images averaged over the impulsive phase of flare 092, showing the full extent of the flare ribbons. Right: HMI magnetogram showing the host active region at the time of the flare peak in the GOES 1–8 Å waveband. The red and blue contours show the positions of the composite flare ribbons with positive and negative magnetic polarity, respectively. The green line joins the geometric centroids of the positive and negativepolarity ribbon components. 

In the text 
Fig. 3. One of the simulated flare light curves, with a flare duration of 213 s, a QPP period of 12 s, and a QPP decay time of 100 s. 

In the text 
Fig. 4. QPP periods plotted against the separation of the geometric centroids of the two flare ribbons. The differently coloured points correspond to the different instruments used to detect the QPP signals, and the grey line is a linear fit in log space. The Pearson correlation is 0.53 with a pvalue of 0.01, and the Spearman correlation is 0.64 with a pvalue of 0.001. 

In the text 
Fig. 5. QPP periods plotted against the total unsigned magnetic flux of the photosphere below the flare ribbons. The Pearson correlation is 0.59 with a pvalue of 0.004, and the Spearman correlation is 0.68 with a pvalue of 0.0005. 

In the text 
Fig. 6. QPP periods plotted against the average magnetic field strength of the photosphere below the flare ribbons. The Pearson correlation is 0.47 with a pvalue of 0.03, and the Spearman correlation is 0.50 with a pvalue of 0.02. 

In the text 
Fig. 7. QPP periods plotted against the area of the flare ribbon. The Pearson correlation is 0.58 with a pvalue of 0.005, and the Spearman correlation is 0.59 with a pvalue of 0.004. 

In the text 
Fig. 8. Two variants of the QPP period plotted against the flare impulsive phase duration. Left: full set of 37 flares with QPPs found by Pugh et al. (2017b). The Pearson correlation is 0.53 with a pvalue of 0.0002, and the Spearman correlation is 0.50 with a pvalue of 0.0006. Right: plot includes only those 20 flares that occurred within ±60° of the solar disc centre, and hence is the same sample of flares used in Figs. 4–7. The Pearson correlation coefficient is 0.56 with a pvalue of 0.006, and the Spearman correlation coefficient is 0.65 with a pvalue of 0.001. 

In the text 
Fig. 9. Flare impulsive phase durations plotted against the flare ribbon separation distance (left), total magnetic flux (middle), and area (right). The grey lines show linear fits in log space. Left: the Pearson correlation is 0.80 with a pvalue of 9 × 10^{−6}, and the Spearman correlation is 0.84 with a pvalue of 1 × 10^{−6}. Middle: the Pearson correlation is 0.93 with a pvalue of 5 × 10^{−10}, and the Spearman correlation is 0.92 with pvalue of 1 × 10^{−9}. Right: the Pearson correlation is 0.83 with a pvalue of 1 × 10^{−6}, and the Spearman correlation is 0.88 with a pvalue of 9 × 10^{−8}. 

In the text 
Fig. 10. Simulated flare QPP period plotted against the flare duration (left) and QPP signal duration (right). The grey points represent the input parameters of the simulated flares with QPPs, while the black points represent the simulated flares with QPP signals that were detected using the method described in Sect. 3.2. Left: the Pearson correlation is 0.23 with a pvalue of 0.0, and the Spearman correlation is 0.30 with a pvalue of 4 × 10^{−27}. Right: the Pearson correlation is 0.31 with a pvalue of 0.0, and the Spearman correlation is 0.30 with a pvalue of 3 × 10^{−27}. 

In the text 
Fig. 11. Simulated flare QPP period plotted against the flare duration for the case where the QPP signal duration is set to depend on the flare duration. The grey points represent the input parameters of the simulated flares with QPPs, while the black points represent the simulated flares with QPP signals that were detected using the method described in Sect. 3.2. The Pearson correlation is 0.18 with a pvalue of 0.1, and the Spearman correlation is 0.21 with a pvalue of 0.07. 

In the text 
Fig. A.1. HMI magnetograms for all flares listed in Table 1, showing the host active region at the time of the flare peak in the GOES 1–8 Å waveband. The red and blue contours show the positions of the composite flare ribbons with positive and negative magnetic polarity, respectively. The green lines join the geometric centroids of the positive and negativepolarity ribbon components. 

In the text 
Fig. A.1. continued. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.