A method for monitoring the precipitation of magnetospheric particles comprising the steps of: - detecting charged magnetospheric particles by at least one particles detector installed on board at least one satellite vehicle in orbit associating the detected particles with respective detection data, said detection data including a time and a geomagnetic longitude at which the detection of the particle occurred; - processing the detection data, comprising estimating or measuring a kinetic energy of the detected magnetospheric particles and associating the respective estimate or measurement of kinetic energy with the detected magnetospheric particles; - obtaining a first count value N H correlated to the number of charged particles, detected in a period of time, which is associated with a relatively higher estimate or measurement of kinetic energy included in a first energy range; - obtaining a second count value N L correlated to the number of charged particles, detected in said period of time, which is associated with a relatively lower estimate or measurement of kinetic energy included in a second energy range; - calculating a value of a relative variation of the second count value N L or of the first count value N H; - determining that in the aforementioned period of time an impulsive event of precipitation of charged particles in the magnetosphere occurred - MPP event - comparing the aforementioned value of the variation with a threshold value; - assigning to the MPP event a geomagnetic longitude and a time in which the MPP event occurred; - defining one or more groups of MPP events, each group comprising MPP events which occurred in a time range at the same geomagnetic longitude or at relatively close geomagnetic longitudes; - identifying one or more of the defined groups of MPP events as indicative of an activity of terrestrial origin, such as for example a pre-seismic o seismic activity, based on the number of MPP events included in the group and/or based on the values of the variations calculated in the calculating step; the method being characterized in that: said step of calculating a value of a relative variation comprises calculating a value of a relative variation of the second count value N L with respect to the first count value N H. - detecting charged magnetospheric particles by at least one particles detector installed on board at least one satellite vehicle in orbit associating the detected particles with respective detection data, said detection data including a time and a geomagnetic longitude at which the detection of the particle occurred; - processing the detection data, comprising estimating or measuring a kinetic energy of the detected magnetospheric particles and associating the respective estimate or measurement of kinetic energy with the detected magnetospheric particles; - obtaining a first count value N H correlated to the number of charged particles, detected in a period of time, which is associated with a relatively higher estimate or measurement of kinetic energy included in a first energy range; - obtaining a second count value N L correlated to the number of charged particles, detected in said period of time, which is associated with a relatively lower estimate or measurement of kinetic energy included in a second energy range; - calculating a value of a relative variation of the second count value N L or of the first count value N H; - determining that in the aforementioned period of time an impulsive event of precipitation of charged particles in the magnetosphere occurred - MPP event - comparing the aforementioned value of the variation with a threshold value; - assigning to the MPP event a geomagnetic longitude and a time in which the MPP event occurred; - defining one or more groups of MPP events, each group comprising MPP events which occurred in a time range at the same geomagnetic longitude or at relatively close geomagnetic longitudes; - identifying one or more of the defined groups of MPP events as indicative of an activity of terrestrial origin, such as for example a pre-seismic o seismic activity, based on the number of MPP events included in the group and/or based on the values of the variations calculated in the calculating step;