Perseids 1901-2100

Introduction
Perseids are one of the most famous meteor showers. Perhaps the only other shower can rival them in popularity - the Leonids in November known with their phenomenal storms. The reason is in regular and very high level of activity as well as summer time in northern hemisphere (where Perseids are best seen) when the shower is active. Comfortable weather allows very large number of observers watch it.
The parent comet of the shower is 109P Swift-Tuttle having orbital period about 135 years. The comet passed its last perihelion in December 1992, which caused considerable increase in Perseid activity in the 1st half of 1990s, up to 600 meteors per hour. Primary maximum, produced by the fresh comet material was being observed till the end of 1990s, gradually decreasing to 100-120 meteors per hour. It had ceased to appear since 2000. In general the future shower activity is to decline as the parent comet is now moving into the outer layers of the Solar system. It doesn't deny however the possibility of bursts in separate years.
Principally when a comet came close to the Sun (and to the Earth so far), our planet meets more dust ejected by it and we see more meteors. But this picture is a very big simplification of the real situation as during each perihelion passage the comet ejects material which soon forms itself into a dust trail. These trails continue moving in space and can keep their structure during a long time. So far the stream consists of a number of trails with slightly differing orbits. This means that the shower activity can vary greatly from year to year.
Can we compute the movement of these trails and make a prediction of shower activity? We can if we have enough volume of information. Such computations are most frequntly done for the Leonids as this shower is studied much better than the Perseids, for which much less predictions were done. For example, Esko Lyytinen issued a prediction for the Perseids 2004, namely for the encounter with 1 rev. trail at 20:50 UT 11 August. Around this time Perseid activity increased up to about ZHR=170 (ZHR - see below), which was considerably higher than usual levels of maximum (ZHR=100-120).
As in case with other showers, evolutuon of Perseid stream is traced with the modelling very good. However, the Perseids have the following distinctive feature: the main part of trails passes too far outside of the Earth orbit and so far observed Perseid activity is almost totally caused by bachground material. Now we have only one more or less clearly fixed case of trail encounter - already mentioned Perseids 2004.
Particles ejected by the comet form lengthy trails. One of the reasons is radiation pressure force, which acts parallel with gravitational force. The latter is dependent on a particle mass, i.e. it is proportional to the third power of particle radius. The outcrying radiation pressure is defined by the second power of particle radius. So far the influence of radiation pressure is the more the less size of a particle is. Its influence is equivalent to the diminishing of gravitational constant G. So it increases the orbital period of particles, and the tinier a particle is, the more it is continuously retarded from larger particles after their ejection from the comet. This process therefore leads to the formation of lengthy comet trails.
Meteor modelling is done through computation of orbital evolution of particles ejected by the comet with different velocities in directions tangential to the comet trajectory at the moment of perihelion. In the reality, of course, particles are ejected not only at the point of perihelion, but also during several months around it. However, comets are close to the perihelion during quite a little time comparing to their overall orbital period and main perturbations happen around their aphelions, so when comets are closer to the Sun newly ejected particles are moving very close to them in a compact dust cloud. This is the reason we can take that cloud as completely ejected at the point of perihelion.
Speaking of directions in which particles are ejected we can say that, again, in the reality they are ejected far not only in tangential directions, but in all possible ones. However, ejection velocities (from 0 to 100 m/s, and the overwhelming majority of real ejections - from 0 to 20 m/s) are negligibly small comparing to the own comet velocity (from 30 to 40 km/s) near the Earth's orbit), ejected particles have only slightly changed orbits and don't "fly away in all directions". Radial part of ejection velocity defines only thickness of a trail, which usually reaches several hundred of thousand kilometers. The shape of the trail is defined by tangential part of ejection velocity.
And the last. Non-gravitational forces are often not taken into consideration in meteor calculations, as is in our case. However, some of them, say, radiation pressure, can be considered indirectly. As far as this kind of force works as diminishing of gravitational constant G, this is equivalent to increase of ejection velocity which could be easily accounted in the model. So this non-gravitational force, as many others doesn't change the configuration of trails, but leads to shifting of particles with different masses along them.
As spoken previously, Perseid trails modelling allowed to prepare very good predictions of shower activity in the previous years. More serious problem is prediction of outburst intensity - how strong the maximum could be. For such predictions special empirical models were elaborated (the single possible way in this case) but as before for their improvement new observations are very necessary.
The paper presents the description of assumed past and future Perseid activity for the period 1901-2100. Computations were done for each year in this period, and, as the results are presented in 2006, they are truly predictive for 2007-2100 years, for the rest of them "postpredictions" were compiled. Also, as the models, used in computations are based after all on observations of real activity in the past, we will not compare each postpredition with real Perseid activity in respective years.

Computation characteristics
This paper presents the results of the Perseid meteor stream simulation aimed to prediction of its meteor activity in 1901-2100. The simulation was made for trails of latest 7 revolutions, i.e, from the 1992 trail. The Author used the program "Comet's Dust 2.0" by S. Shanov and S. Dubrovsky to calculate orbital elements of ejected meteor particles. To estimate expected ZHRs for different encounters the model described in [4] was used with some Author's alterations made in order to adopt the model for ejection velocity (v) instead of da0 (difference in a-semimajor axis) and to turn the model from the Leonid stream (for which it was originally created) to the Perseids. The computation considered only gravitational forces, however, the results are on the whole in good accordance with these of other researchers. The prediction includes all encounters found within interval +/-0.007 a.u. The following parts of trails were computed: the first 5 rev. trails for ejection velocities [-50;100] m/s, 6-7 rev. trails - [-30;50] m/s.
Predictions for different years of Perseid activity are divided into decades, for each year with meaningful encounters a table describing them is given. A typical sample of such table is shown in Table A1.
Table А1
2000
encounters with trail
trailyearrD-rEVejfM(fMD)sol.long.Max. timeZHRexReliability
rev.-а.е.m/s-°UTmeteors-
...........................
where 2000 - described year of Perseid activity; trail is the number of revolutions of given trail; year - the year of trail formation; rD-rE - the distance between the Earth's orbit and trail particles orbits (positive value means the descending node of trail is outside the Earth's orbit, negative - the node is inside it); Vej - velocity of particles ejection on the given trail part (positive values means the particles were ejected against the comet movement, negavive - particles were ejected along the comet movement); fM(fMD) - the characteristic of longitudinal density of the trail, it is derived from the time interval between passages of minimal distance to the Earth's orbit by particles with different ejection velocities; sol. long. - solar longitude corresponding to the maximum; Max. time - assumed time of maximum; ZHRex - ZHR expected to be produced by the trail (it should be noted, that this ZHR is "pure", i.e. it is given without consideration of background meteor activity. So far, this ZHRex value should be added to expected background activity at the time of given trail maximum to obtain the total expected level of activity; Reliabitity - an estimation of reliability of the given maximum.
We have to give some details on reliability. For each encounter with a trail the tables give the Author's estimations of its reliability (on 5-degree scale) and for traditional maximum - the Author's estimation of its intensity (also on 5-degree scale). The values of these estimations have the following interpretation:
А. Reliability of encounters with trails:
4 - very high reliability. The enhancement will almost exactly take place, its actual maximum time should differ from predicted one up to no more than several minutes and its intensity should be quite close to the predicted value.
3 - high reliability. The enhancement will almost exactly take place, but its actual maximum time can differ from predicted one up to 15-20 minutes and its intensity can substantially, in 2-3 times, differ from the predicted value.
2 - moderate reliability. The enhancement will very likely take place, but there is small possibility of its absense. Its actual maximum time can differ from the predicted one up to many tens of minutes and its intensity can differ from the predicted value in many times.
1 - low reliability. There is substantial possibility of enhancement appearance, but its absence is also quite likely. Its actual maximum time can differ from the predicted one up to many tens of minutes and its intensity can differ from the predicted value in many times.
0 - very low reliability. The enhancement appearance is possible, but its absence is more likely. Its actual maximum time can differ from the predicted one up to many hours and its intensity can differ from the predicted value in many times.
Each decade will be accompanied with general comments on what variations of background Perseids activity are to be expected in different years of the described decade. In case of presence meaningful encounters with trails they also would be commented.
As it is impossible to display in tables all nuances and detales for each prediction, information given in comments has priority-driven character comparing to the data in tables. It doesn't mean that it will be in contradiction with tables (at least, not in the Author's point of view), but formal figures far not always can tell the reader how to understand prediction correcly.

Orbit of the comet 109P in 1901-2100
The Author used initial orbital elements of 109P, starting from perihelion of 950 and up to 1992, presented by Yau, Yeomans and Weissman [2]. Orbital elements of 109P in the period 1901-2100, as well as values of minimal distances to the Earth orbit for these elements and relative solar longitudes are given in Table A2.
Table А2
time of perihelionqeAOPNodeimin. dist.sol. long.
-а.е.-°°°AU°
1992.12.12.323940.95821750.9635892153.00138139.44442113.42658-0.0009423139.43815
Orbital elements are given for the Epoch J2000. The following of them are denoted with symbols: q - perihelion distance; e - eccentricity; AOP - argument of perihelion; Node - longitude of ascending node; i - inclination. Positive value of minimal distance means that point of such minimum lies outside the Earth orbit, and negative value means that this point is inside the Earth orbit.