Characteristics of temporal variability of long-duration bursts of high-energy radiation associated with thunderclouds on the Tibetan plateau

From 1998 to 2017, neutron monitors located at an altitude of 4300 m on the Tibetan plateau detected 127 long-duration bursts of high-energy radiation in association with thunderclouds. These bursts typically lasted for 10 to 40 minutes, and 89\% of them occurred between 10:00 and 24:00 local time. They were also found to be more likely to occur at night, especially during 18:00$-$06:00 local time period. The observed diurnal and seasonal variations in burst frequency were consistent with the frequencies of lightning and precipitation on the Tibetan plateau. Based on 19 years of data, the present study suggests that an annual variation in burst frequency has a periodicity of $\sim$16 years and a lag of $\sim$3 years relative to solar activity.


Introduction
In the 1920s, C. T .R. Wilson proposed that electrons in thunderclouds are accelerated to high energy in strong electric fields and produce X-rays by bremsstrahlung (Wilson 1924(Wilson , 1925)).Since then, many observations have shown that lightning and thunderclouds generate such X-ray signals, confirming that they are indeed powerful particle accelerators (Dwyer et al. 2012, and references therein).Beyond the Wilson prediction, bremsstrahlung gamma rays with a few tens of MeV have been detected from lightning and thunderclouds.The detection of the high-energy photons suggests that lightning and thunderclouds can accelerate electrons beyond a few tens of MeV inside them.Theoretical works (e.g., Gurevich et al. 1992;Dwyer 2003Dwyer , 2007;;Babich et al. 2004Babich et al. , 2010) ) have pointed out that the production of those high-energy radiations is due to relativistic runaway electrons.Seed electrons produced by, for example, cosmic rays are thought to be Page 2 of 14 Tsuchiya et al. Progress in Earth and Planetary Science (2024) 11:26 accelerated to relativistic energies when electric fields exceed a critical threshold and become such runaway electrons.Solar activity influences the heliospheric magnetic field, resulting in variations in the intensity of galactic cosmic rays reaching the Earth in an approximately 11-year cycle.Therefore, it is thought that solar activity, by affecting the flux of galactic cosmic rays reaching the Earth, also impacts the production of runaway electrons.However, it remains unclear whether cosmic rays or solar activity plays a role in producing the runaway electrons responsible for high-energy bremsstrahlung gamma rays.
Bursts of high-energy radiation during thunderstorms can be categorized into two types when classified simply by their duration.One type is characterized by intense short-duration bursts lasting for microseconds to milliseconds.The other type is a long-duration burst, which lasts from a few seconds to a few tens of minutes and exhibits relatively fainter characteristics compared to the short-duration bursts.Numerous measurements of these short-duration bursts have been taken by space-and ground-based detectors and revealed several characteristics including the source altitude and their relation to lightning.In contrast to short-duration bursts, longduration bursts such as gamma-ray glows or terrestrial ground enhancements (TGEs) have been observed less frequently until recently.However, observations made in the coastal area of the Japan Sea (Torii et al. 2002;Tsuchiya et al. 2007aTsuchiya et al. , 2011;;Wada et al. 2021a), at high mountains (Torii et al. 2009;Tsuchiya et al. 2009;Chilingarian et al. 2010Chilingarian et al. , 2011;;Tsuchiya et al. 2012) and at high aircraft altitude (Østgaard et al. 2019;Kochkin et al. 2021) have led to increased detection of these long-duration bursts and have drawn attention to their characteristics, including their duration and their correlation with meteorological phenomena.
The increasing amount of data on long-duration bursts have allowed us to compile a comprehensive catalog of such phenomena and to gain insight into their characteristics at various observation sites, including high mountain regions (Chilingarian et al. 2019) and coastal areas of the Japan Sea (Wada et al. 2021b).However, more data are needed to explore whether there is a correlation between long-duration bursts and cosmic rays and to understand the meteorological and geographical factors that contribute to the generation of these bursts.By studying these conditions, we can identify commonalities among long-duration bursts and develop a deeper understanding of their production mechanisms.To this end, we will present observational results obtained by neutron monitors (NM) at Yangbajing on the Tibetan plateau (TP) during 1998−2017 and then discuss the temporal variability of long-duration bursts.

Overview
The Yangbajing NM has been in operation at the Yangbajing cosmic-ray observatory for more than two decades.The observatory is located at an altitude of 4300 m (corresponding to an atmospheric depth of ∼ 600 g cm −2 ) in the TP.In addition, electric field mills (EFMs) have been deployed at the observatory in 2010 to study possible correlations between variations in count rate measured by a detector such as the Yangbajing NM and electric fields (EFs) (Amenomori et al. 2013).To clarify the points of the present work, we briefly explain their characteristics below.

Neutron monitor
The main purpose of the Yangbajing NM is to investigate the ion acceleration at the solar surface and to monitor galactic cosmic-ray intensity associated with the solar activity (Miyasaka et al. 2005;Tsuchiya et al. 2007b).In the present study, the NM is repurposed to examine the production mechanism of long-duration bursts of highenergy radiation associated with thunderclouds.Hereafter, the term "long-duration burst" is employed to refer to the high-energy radiation detected by the Yangbajing NM in association with thunderclouds.
As shown in Fig. 1, the Yangbajing NM consists of 28 NM-type BF 3 counters surrounded by lead blocks and polyethylene in the standard NM64 configuration.
Each NM unit has one BF 3 counter with dimensions of 190 cm in length and 15 cm in diameter.The total area of the Yangbajing NM is 32 m 2 .The outermost polyeth- ylene serves the purpose of moderating incoming fast neutrons, while the embedded lead blocks breed the moderated neutrons via (n, xn) reactions to enhance the detection efficiency for neutrons (x represents a number).When entering the BF 3 counters, the moderated neutrons undergo the (n, α ) reaction: n + 10 B → α + 7 Li .The cross sections of this reaction become larger as the Fig. 1 Arrangement of the Yangbajing NM in a room and a cross-sectional view of one unit, measured in centimeters neutron energy approaches the thermal energy, thus allowing efficient detection of neutrons in thermal to epithermal energy ranges.To further improve neutron detection efficiency, the isotopic composition of 10 B in each BF 3 counter is enriched to 98% from a natural iso- topic ratio of ∼ 20%.Charged particles produced by the above reactions lose energy due to the ionization loss in gas, resulting in anode signals in the counter.A data acquisition system collects the output signals from each counter and registers the number of signals every one second.
In our previous study (Tsuchiya et al. 2012), we determined the proportions of individual particles that contribute to the total count rate of the Yangbajing NM.Secondary cosmic-ray neutrons and protons dominate the total count rate with fractions of 83 and 12%, respectively.In contrast, leptons, including gamma rays and electrons, contribute to it with only 5% [see also Clem and Dorman (2000)].The previous study also indicated that the same type of NMs as the Yangbajing NM have a small but non-negligible detection efficiency for gamma rays through photonuclear reactions with lead blocks ( γ 206,207,208), with a threshold gamma-ray energy of 7 MeV (IAEA 2000).Given these insights, we demonstrated that bremsstrahlung photons were responsible for the observed count increases during thunderstorms in the Yangbajing NM, rather than neutrons produced by photonuclear reactions in the atmosphere [e.g., Enoto et al. (2017)].Therefore, for the present study as well, we assume that count increases of the Yangbajing NM can be attributed primarily to bremsstrahlung photons with energy higher than 7 MeV.

Electric field mill
Electric field data play a crucial role in understanding the association of observed increases with thunderclouds and lightning.For this purpose, two BOLTECK EFM-100 commercial electric field mills have been operational at the observatory since 2010, positioned 25 m apart.One is situated at ground level, while the other is mounted atop an experimental building.The vertical separation between them is approximately 3.4 m.The signals from both EFMs are transmitted to computers within the building via optical cables.These signals are recorded at an interval of every 0.1 s, providing field strength below 100 kV m −1 , with a measurement resolution of 20 V m −1 .

Analysis
To identify long-duration bursts, we performed several steps of analysis of data collected by the Yangbajing NM from 1998 to 2017.First, we made a 5-minute count history corrected for atmospheric pressure variations for each day.Pressure fluctuations typically induce count changes within ±2 %.Second, we constructed an occur- rence histogram of NM counts in every 5-minute count history, N 5 min , and fitted a Gaussian function to this his- togram to determine the mean value N BG and the stand- ard deviation σ 1day of the Gaussian function every day.The resulting σ 1day served as an estimate of the uncer- tainty in the 5-minute counts for each day.Third, we calculated the statistical significance, σ s , for each N 5min using the formula; This allowed us to assess whether the daily NM data contained count bins with σ s > 4 .The probability of encoun- tering a bin with a significance of 4 or higher by chance is 5 × 10 −5 .Therefore, out of the 288 bins in each 5-minute count history, the number of bins with a significance of 4 σ or greater is only 0.014.Among the N 5min counts, we defined the most prominent one as the peak count, N p , in each 5-minute count history.This data set, which we refer to as the "potential burst set", contains 135 events that represent possible long-duration bursts associated with thunderclouds.
From the potential data set, we extracted a 5-minute count history over a time span of ± 2 h centered on the peak time " t p " when N p was recorded.This time interval of ± 2 h is sufficient to cover the entire duration of each burst and to estimate the background level of each burst.To define the burst duration " t d ", we identified a con- tinuous time range in which the statistical significance σ s remained above 1 before and after the peak time t p .By excluding all count bins in this t d , two-time regions are left in every extracted count history to estimate the background.The background level was determined by fitting count bins in the two time regions with a quadratic function.The background level in t d was then estimated by interpolating the derived quadratic function.
After completing the above analysis, an eye scan was performed to detect any suspicious events that may have been caused by irregular factors such as power failures.As a result of this scan, eight events were identified as suspicious events in the potential burst set.Subsequently, a total of 127 long-duration bursts were selected from the NM data collected between 1998 and 2017.This data set is referred to as the final burst set.None of the events in the final burst set was found to be associated with solar energetic particle events, including ground-level enhancements (Vainio et al. 2017;Cohen and Mewaldt 2018), suggesting that such a solar particle event does not contribute to the production of long-duration bursts. (1) The current analysis is limited in the detection of faint bursts, because it only considers those with > 4 σ .This selection criterion is one of the reasons for the detection of only 127 events during the long observation period (1998−2017).This frequency is also likely influenced by the higher NM threshold of 7 MeV and the limited spatial illumination of a long-duration burst.Gamma-ray detectors using NaI and BGO scintillators can readily capture long-duration bursts with energies of a few MeV or less (Chilingarian et al. 2019;Wada et al. 2021b), while these low-energy gamma rays are not detectable by the Yangbajing NM.In addition, it is known that a long-duration burst illuminates only a limited area (0.5-1 km) on the ground (Tsuchiya et al. 2011;Wada et al. 2019;Tsurumi et al. 2023).At the altitude of 4300 m, the mean free path of gamma rays with energy of 10−50 MeV is approximately 700−900 m.Therefore, it would be more difficult to clearly detect long-duration bursts originating significantly farther away from NMs than the mean free path.

Example of time profiles of long-duration bursts
Figure 2 depicts examples of the 5-minute count histories in the final burst set.It is evident that the individual count histories exhibit significant count increases.Table 1 presents several key characteristics for each of 127 longduration bursts, including its net count increase denoted as " N tot " and the related statistical significance.The N tot is calculated by subtracting the interpolated background from the total observed count in t d .The observed count increases are superimposed on gradual count changes caused primarily by solar diurnal variation, which typically ranges within ±0.3 % in amplitude.This variation arises from the anisotropy of galactic cosmic rays below several tens of GeV in interplanetary space.
Unlike the gradual changes induced by the solar diurnal variation, EF variations result in faster changes in the NM count rate.Figure 3 shows examples of the relationship between EF changes and observed count increases.The EF intensity remains at ∼ 0.1 kV m −1 under fair weather conditions, with its sign following atmospheric physics Fig. 2 Examples of the NM count rate histories per 5 min in the final burst set.Each panel is labeled with the date of the event acquisition.The horizontal axis depicts universal time (UT).Each error shown corresponds to a statistical 1 σ .Red solid lines represent estimated background levels, while two black dashed lines indicate a range of ± σ 1day from the background level conventions.A positive value denotes an EF directed downward from the sky to the ground and vice versa.The EF is found to exhibit significant variations lasting for 1−1.5h.While the count increases in the events (Fig. 3) seem to coincide with periods of relatively large positive EFs, the peak times of these count increases do not necessarily align with the peak times of the positive EFs.A possible explanation could be that the alterations in the EF within a thundercloud, responsible for electron acceleration, may not precisely synchronize with the observed variations in the EF at the ground surface.Intense, sharp peaks are also found in EF variations (e.g., 20150827 and 20160617) and indicate the presence of lightning.These lightning events causing significant peaks in EF variations do not appear to impact the observed long-duration bursts.This is probably because the location of the lightning observed by the EFM is far from NMs.Therefore, when investigating the relationship between EF and longduration bursts, it would be necessary to specify not only the EF intensity but also the distance to a thundercloud or lightning.This will be our future task.
All other events detected after the EFM installations are also found to have occurred under similar EF conditions (Additional file 1), suggesting a common underlying mechanism to produce a long-duration burst.Typically, a thundercloud is considered to have a tripole  charge structure, comprising a lower positive charge region (LPCR), a middle negative charge region, and an upper positive charge region (Cooray 2015).Electrons, located between the LPCR and the middle negative charge region, undergo the electric-field acceleration toward the LPCR and emit high-energy photons through bremsstrahlung directed to the ground.Based on measurements taken in the TP at an altitude of 4500 m above sea level (a.s.l.), it was found that the cloud base altitude is generally 1 km above the ground surface (Qie et al. 2005).In our previous study (Tsuchiya et al. 2012), we examined a long-duration burst that occurred on 22 July 2010 (listed as 20100722 in Table 1), and found that about 90% of the NM count enhancement was caused by bremsstrahlung photons emitted by electrons accelerated in thunderclouds.This estimation assumed a cloud base altitude of typically about km above the ground surface.Actually, the constitution ratio of these photons varies according to the distance between the cloud base and the ground, ranging from 85% at a distance of 5 km to 95% at a distance of 0.3 km (Tsuchiya et al. 2012).Therefore, the present bursts exhibit constitution ratios of 85−95% if the cloud-base altitude is in between 0.3 and 5 km.

Duration
Figure 4 presents the duration distribution for the final burst set.Most of the observed duration ranged between 10 and 40 min, occasionally extending to ∼ 60 min.This trend agrees with the typical life cycle of thunderclouds, 30−60 min, during the presummer and summer seasons.Their mature stage, which is likely to produce lightning, generally persists for 15-30 min (Ferrier and Houze 1989).
The ARGO-YBJ experiment, operating ∼ 200 m away from the NM location, observed variations of the airshower rate during 2012 (Axikegu et al. 2022).These variations lasted 14 to 55 min, and are in good agreement with our results.Count enhancements with similar duration have been identified on Mt.Aragats (3250 m a.s.l.) and Mt.Lomnický Štít (2634 m a.s.l.), using multiple detectors such as NMs, a NaI detector, and a plastic scintillation detector.The count increases obtained by NMs (Chilingarian et al. 2010(Chilingarian et al. , 2011) ) lasted ∼ 10 min, and most enhancements of the NaI count in energy > 3 MeV lasted for a few minutes to 50 min, with rare occasions extending to 120 min (Chilingarian et al. 2019).Similarly, those obtained at Mt. Lomnický Štít lasted for 1-15 min (Chum et al. 2020).The typical durations within 50 min are in agreement with the present result.However, the > 120 min duration has never been observed at the Yangbajing observatory.This discrepancy in burst duration may be attributed to two factors.One is the different meteorological conditions between Mt.Aragats (3250 m a.s.l) and Yangbajing (4300 m a.s.l.): for example, successive generations of cells in thunderclouds might be more pronounced on Mt.Aragats.This difference may result in different life cycles for thundercloud-charged layers, where electron acceleration takes place.The other is a difference in the energy threshold for photons of 7 MeV for the Yangbajing NM and 3 MeV for their NaI detector.This difference in energy thresholds could also lead to different duration behaviors.
In contrast to high-altitude observations made during the pre-summer or summer seasons, observations conducted at the Japan Sea in winter seasons have revealed an alternative perspective on long-duration bursts, which generally last only a few minutes (Torii et al. 2002;Tsuchiya et al. 2007aTsuchiya et al. , 2011;;Wada et al. 2021b).According to Kitagawa and Michimoto (1994), tripole charge structures become apparent during the mature stages of winter thunderclouds.These tripole structures persist for less than 10 min during early or late winter and for even shorter duration, less than a few minutes during midwinter.Thus, it is reasonable to postulate that differences in the life cycles of thunderclouds are responsible for the noticeable difference in burst duration among long-duration bursts during the summer and winter seasons.

Diurnal variation
The distribution of the onset times for the final burst set is shown in Fig. 5.This distribution is clearly divided into two sections: one covering the period from 00:00 to 10:00 local time (LT) and the other covering 10:00 to 24:00 LT.LT is calculated as 6 h + universal time (UT), consider- ing the geographical longitude of 90 • .In the total number of observed bursts, 11% of them began at 00:00-10:00 LT and the rest occurred at 10:00-24:00 LT.The latter time region contains the period of increased convective activity in convective clouds over the TP as reported by Uyeda et al. (2001) which showed that strong updrafts generated by land surfaces heated due to solar radiation in TP result in the rapid development of convective clouds during 09:00−13:00 LT in the premonsoon seasons and 12:00−16:00 LT in the monsoon seasons.In addition, Qie et al. (2003) investigated lightning activity between 1998 and 2002 in TP using the lightning image sensor onboard the Tropical Rainfall Measuring Mission.Their investigation indicates that 6% and 78% of lightning flashes in the diurnal lightning activity occur in 00:00−10:00 LT and 12:00−19:00 LT, respectively.Therefore, a clear correlation is evident between the distribution of the observed burst start times and the temporal pattern of lightning flashes.Wada et al. (2021b) examined the onset time distribution of 70 long-duration bursts detected in the coastal area of the Japan Sea over 4 winter seasons from 2015 to 2018.They divided the data into two parts: the daytime spanning 06:00−18:00 LT and the nighttime encompassing the remaining time period.They showed that the number of events in daytime is 15 and significantly smaller than 55, the number of events in nighttime.The present study also shows a comparable trend, indicating that the number of events during the day, 43, is less than the number of events during the night, 84, with a

Seasonal variation
Figure 6 shows the seasonal variation in the occurrences number of the final burst set.Among them, 122 events (97% of the total 127 events) were observed between May and September, corresponding to the premonsoon and monsoon seasons in TP.The seasonal trend obtained is similar to that of precipitation and lightning activity within the TP region (Qie et al. 2003).This similarity suggests that long-duration bursts tend to occur during the rainy season when lightning is more frequent, as expected.
According to the EF measurements from March 2010 to February 2012 (Amenomori et al. 2013), the total number of thundercloud occurrences at Yangbajing was 158 days.During this period, the Yangbajing NM detected only seven long-duration bursts (Table 1).Interestingly, Amenomori et al. (2013) showed that thunderclouds were also observed in 54 days during winter seasons from December to March.However, the Yangbajing NM did not detect any long-duration bursts during this winter period, except only one recorded between 1998 and 2017 (20040203 in Table 1).This infrequent detection during winter suggests that the EF intensity due to thunderclouds may not be the sole factor governing the production of long-duration bursts.
The ARGO-YBJ experiment, conducted near the NM location, detected twenty count enhancements during thunderstorms in 2012 between April and October (Axikegu et al. 2022).This period coincides with the present result.Notably, four (20120714, 20120804, 20120911, and 20121010) of the six long-duration bursts detected in 2012 (Table 1) align with the dates recorded by the ARGO-YBJ experiment.The effective area of the Yangbajing NM (32 m 2 ) is smaller than that of the ARGO-YBJ detector ( ∼ 110 × 100 m 2 ).A larger effective area is likely to detect more long-duration bursts with narrow illumination areas.Hence, the difference in effective area is considered one of the main factors contributing to the variation in the number of events detected in 2012 (6 for the Yangbajing NM and 20 for the ARGO-YBJ experiment).Muraki et al. (2004) studied a seasonal variation using a large-area meson monitor installed at Mt. Norikura (2770 m a.s.l.) in Japan.The meson monitor comprises scintillation detectors, and primarily detects charged particles but also has some sensitivity to neutral particles such as neutrons and gamma rays.Analyzing data obtained from October 1990 to January 2002, they found that the derived seasonal variation matches the pattern of thunderstorm occurrence.The seasonal variation observed with the Yangbajing NM is similar to their result.This is probably because the two areas are influenced by the Asian monsoon.Chilingarian et al. (2019) also showed a seasonal variation, with a notable peak in May, which may probably be due to different meteorological conditions compared to those observed in the TP and at Mt. Norikura.et al. 1992;Babich et al. 2007;Carlson et al. 2008;Gurevich and Karashtin 2013).As is well known, variations in sunspot number reflect solar activity that has the 11-year cycle, and is inversely correlated with variations of galactic cosmic-ray intensity at Earth.Thus, it is natural to conjecture that solar activity may influence the occurrence of long-duration bursts.
To quantify the annual variation of the occurrence rate in the final burst set, we employed the following functions for fitting: (1) a constant function, (2) a linear function, and (3) a periodic function.Functions (1) and ( 2) are denoted f c (t) = a 0 and f l (t) = a 0 + a 1 t , respec- tively.Here, a 0 and a 1 represent a constant and a slope, respectively.The variable t shows the elapsed year since 1998.The function (1) assumes that solar activity or cosmic rays exert time-independent effects on the observed number of long-duration bursts.Function (2) assumes a where A, T and φ denote a normalization constant, a period, and a phase, respectively.As listed in Table 2, the resultant χ 2 /ν values ( ν is the degree of freedom) indicate that f p (t) provides the best fit to reproduce the derived annual variation of the final burst set.This result implies the potential presence of periodicity in the annual occurrence rate of the observed long-duration bursts.
For comparison, we evaluated periods of annual variations in the sunspot number and cosmic-ray flux as well by fitting the following function to the variations This function considers a clear decrease (increase) in the amplitude of the sunspot number variation (cosmic-ray flux) over the years.The resulting curves of f sun (t) and f cr (t) are drawn in Fig. 7, and the parameters determined are listed in Table 2.The evaluated periods (T) for the sunspot number and cosmic-ray flux deviate from the period estimated from the final burst set by 2.2σ and 2.7σ (Table 2).In addition, by using the derived periods (T) and phases ( φ ) in Table 2, we can calculate the peak years of f p and f sun , and the year of the minimum of f cr , each measured from 1998, as 6.9 ± 1.3 years, 3.58 ± 0.13 years, and 4.756 ± 0.003 years, respectively.From these values, possible slight lags of f p with respect to f sun and f cr were estimated as 3.3 ± 1.3 years and 2.1 ± 1.3 years, respectively.While the present analysis suggests a potential connection between long-duration bursts and solar activity or cosmic rays, we need more data over a longer (2) f p (t) = a 0 + A sin 2π T t + φ , (3) f sun,cr (t) = a 0 + a 1 t + A sin 2π T t + φ .
time period before definitively confirming the periodicity and delay.

Conclusions
Long-duration bursts were detected by the Yangbajing NM located at 4300 m a.s.l.during the period of 1998−2017.The detected bursts have a typical duration of 10 to 40 min.This duration distribution suggests that electrons undergo continuous acceleration to 10 MeV or higher within a few tens of minutes.In the coastal area of the Japan Sea, long-duration bursts have also been observed in winter seasons and typically last for a few minutes.The difference in duration would result from the distinct life cycles of winter and summer thunderclouds.
The seasonal and diurnal variations in the occurrence number of long-duration bursts agreed well with those variations of lightning and precipitation on the TP.This agreement justifies the selection and subsequent analysis of the Yangbajing NM data.The present work, based on 19-year data, suggests the possibility that the annual variation of long-duration bursts has a periodicity of ∼ 16 years and a lag of ∼ 3 years relative to solar activity.To gain further insight into the observed annual variability, it would be necessary to conduct longer-term high-energy radiation observations coupled with comprehensive meteorological measurements.
The present work highlights the importance of NM data for understanding the temporal variability and underlying production mechanisms of long-duration bursts of high-energy radiation associated with thunderclouds.Since the 1950s, numerous NMs have been operational at various observatories around the world (Bütikofer 2018).These NM data would allow us to unravel the real connections between long-duration bursts, cosmic rays, and solar activity.

Fig. 3
Fig. 3 Comparison between the BG-subtracted NM count rate and electric field (solid line).The abscissa depicts universal time (UT).The errors quoted from the NM counts are statistical 1 σ

Fig. 4
Fig. 4 Duration distribution of long-duration bursts in the final burst set

Fig. 6 Fig. 5
Fig. 6 Distribution of months in which events occurred in the final burst set (filled histogram).Filled squares and circles indicate lightning activity and precipitation, respectively, during the years 1998−2002 in the central part of the TP(Qie et al. 2003)

Fig. 7
Fig. 7 Top: annual occurrence variation of the final data set (red histogram).The black histogram shows an annual variation in sunspot number expressed in %.Solid red and black curves denote a periodic function of f p (t) and f sun (t) , respectively, while dashed lines represent f c (t) and f l (t).The error bar represents statistical 1 σ.Bottom: the annual variation of cosmic-ray flux measured by neutron monitors located at Jungfraujoch.The NM data were obtained from the following https:// cidas.isee.nagoya-u.ac.jp/ WDCCR/ index.html AbbreviationsTPTibetan plateau NM Neutron monitor EFM Electric field mill

Table 1
The parameters determined for individual events in the final burst set

Table 1
(continued) a Start time b End time c Significance of N tot

Table 2
The fitting parameters determined for each function a ν is degree of freedom.