The Polar Cap (PC) indices: Relations to solar wind parameters and global magnetic activity

The polar geomagnetic activity resulting from solar wind–magnetosphere interactions can be characterized the Polar Cap (PC) indices, PCN and PCS. PC index values are derived from polar magnetic variations calibrated on a statistical basis such that the index approximate values in units of mV/m of the interplanetary “geo-effective” (or “merging”) electric field (E_M) conveyed by the solar wind. The timing and amplitude relations of the PC index to solar wind plasma and magnetic field parameters are reported. The solar wind effects are parameterized in terms of the geo-effective electric field (E_M) and the dynamical pressure (P_DYN). The PC index has a delayed and damped response to E_M variations and display saturation-like effects for E_M values exceeding 10 mV/m. Steady or slowly varying levels of solar wind dynamical pressure have little or no impact on the PC index above the effects related to E_M for which the solar wind velocity is also a factor. Sharp increases in the dynamical pressure generate impulsive variations in the PC index comprising a initial negative impulse of 5–10 min duration followed by a positive impulse lasting 10–20 min. Typical amplitudes of both the negative and the positive impulses are 0.2–0.5 units. A sharp decrease in the pressure produces the inverse sequence of pulses in the PC index. Auroral substorm activity represented by the AL index level has a marked influence on the average PC/E_M level at the transition from very quiet (AL0 nT) to disturbed conditions while more or less disturbed conditions (AL<−100 nT) have no systematic effect on the average PC/E_M values. At distinct substorm events the PC/E_M ratio has a minimum (0.8) in the pre-onset phase at around 20 min before substorm onset. The average ratio gradually increases in the expansion phase to reach a maximum value (1.1) at around 40 min after substorm onset (or 20 min after the largest (negative) peak in AL). At substorm recovery during the next 2 h the PC/E_M ratio decreases. Finally, we report on the application of polar magnetic variations to model the disturbance storm time (Dst) index development during magnetic storms by using the PC index as a source function to quantify the energy input to the ring current representing accumulated storm energy and characterized by the Dst index.     

Some evidence of ground power enhancements at frequencies of global magnetospheric modes at low latitude

A statistical analysis of the power spectra of the geomagnetic field components H and D for periods ranging between 3 min and 1 h was conducted at a lowlatitude observatory (Aquila, L = 1.6) at the minimum and maximum of the solar cycle. For both components, during daytime intervals, we found evidence of power enhancements at frequencies predicted for global modes of the Earths magnetosphere and occasionally observed at auroral latitudes in the F-region drift velocities (approximately at 1.3, 1.9, 2.6, and 3.4 mHz). Nighttime observations reveal a relative low frequency H enhancement associated with the bay occurrence together with a peak in the H/D power ratio which sharply emerges at 1.2 mHz in the premidnight sector. The strong similarity between solar minimum and maximum suggests that these modes can be considered permanent magnetospheric features. A separate analysis on a two-month interval shows that the observed spectral characteristics are amplified by conditions of high-velocity solar wind.     

The height of the maximum ionospheric electron density over the MU radar

Ionospheric F₂-layer peak height h_mF₂ variations, as measured over 1986–1995 by the MU radar (34.85°N, 136.1°E) and as calculated with a theoretical model, are discussed. The diurnal variations of the measured peak height for different seasons and levels of solar activity are compared with those estimated from ionosonde M3000F₂ and IRI predictions. Also given are the measured ion drift velocities and meridional neutral winds needed to understand the dynamic behavior of the F₂-layer. It is found that: (1) h_mF₂ is generally higher during periods of the solar maximum than during periods of the solar minimum, and higher in summer than in winter; (2) for the solar maximum, h_mF₂ drops markedly in the morning and in the afternoon, while, for the solar minimum, the h_mF₂ minimum occurs in the morning during summer and usually in the afternoon during winter. In general, the measured h_mF₂ is well reproduced by our model when we use the observed drift velocities and plasma temperatures as inputs. Our modeling study shows that the neutral wind contributes strongly to the diurnal variation of h_mF₂ in winter by lowering the ionization layer by day, particularly for the solar maximum; it also helps to enlarge the day–night difference of h_mF₂ in summer. The northward electromagnetic drifts that usually cancel the neutral wind effect have only a minor effect for the location of the MU radar. Other features of the observed h_mF₂ variations, e.g., the solar maximum–minimum difference, the summer–winter difference, and the morning and afternoon drops, are explained by the basic processes of O⁺ production, loss and diffusion, as influenced by the atomic oxygen concentration and neutral and plasma temperatures.     

Equatorial ionospheric response to the 10.7 cm radio flux over two sunspot cycles (1969-1991)

It is evident that fluctuations in a standard ionospheric parameter, the minimum (virtual) height (h^′F) of the equatorial F-region in the African (Ouagadougou), Asian (Manila) and American (Huancayo) longitudinal sectors, closely resemble changes in solar activity as deduced from the 10.7 cm solar flux index (S), over two solar cycles (1969-91). The monthly median hourly value of h^′F, particularly in the post-sunset period (18-20 LT), are positively correlated with the monthly average S. The value of h^′F can be deduced from an empirical formula: h^′F=0.68S+218.3, with the correlation coefficient (r) between h^′F and S being 0.78. The diurnal distribution of r during daytime (06-14 LT) was radically different for the African and Asian longitudinal sectors during 1980-1991, with the most pronounced difference in the post-noon period (12-14 LT) when the correlation coefficients r for the Asian and African sectors are 0.8 and 0.2, respectively. Thus, the daytime F-region in the African sector responded far less to changes in solar activity than the Asian F-region during this cycle. This longitudinal anomaly was however absent in the preceding cycle (1969-1979) when the African and Asian sectors were both characterised by low daytime and pronounced post-sunset correlation coefficient r. The American sector appears to have a high correlation coefficient r in daytime increasing to a small maximum in the post-sunset interval. The post-sunset enhancement in r is a characteristic feature for equatorial stations only (corrected geomagnetic latitude <10^°).     

Solar EUV Flux Variations in 11–Year Solar Activity Cycle According to Measurements On-Board PROGNOZ Satellite Series, INTERBALL-1 and ELECTRO Satellites during 1978–1997 Period

The paper presents data on flux variations of ionizing EUV radiation for wavelengths shorter than 130 nm and in L-alpha hydrogen line on the basis of measurements performed on-board a Russian Earth-orbiting and interplanetary spacecrafts in the period since 1978 to 1997. These data were received by using common heritage instrumentation installed on four PROGNOZ satellites, INTERBALL-1, PHOBOS, and the geosynchronous meteorological satellite ELECTRO. Though the time series of observations is not continuous, it gives a possibility to estimate flux variations in certain periods of time and during 11–year cycle. During quiet periods of solar activity (except for flares) a level of diurnal variations is rather low, it being equal to about 3%. Radiation changes in L-alpha line from maximum to minimum in an 11–year cycle of solar activity achieves 200% and increases in the range of short waves. Data measured onboard the Russian satellites are compared with data obtained by Solar Mesospheric Explorer.     

Day-time variations of the equivalent electron concentration in the non-sounded lower ionosphere

Summary Based on data on the lowest reflected frequencyf _min and on information on the lower and upper boundaries of the non-sounded lower ionosphere, an equivalent electron concentration for all concentrations below the correspondingf _min was determined. Day-time variations of the equivalent concentration are investigated, confirming that there is a cosine relation to the solar zenith angle. The power index of that relation has an outlined seasonal course with a maximum in April and October, while the absolute seasonal minimum is during the winter (the summer minimum is slightly outlined). The mean yearly values of the index are almost constant:n _N 0.5 for solaractivity,I ₁₅₀₀ to 115.10^–22 W Hz^–1 m². During higher activityn _N changes correspondingly toI ₁₅₀₀ according to relation (12). The variations ofn _N during high solar activity show that the altitude gradient and temperature gradient in the low ionosphere are becoming proportional toI ₁₅₀₀ when the solar x-ray radiation exceeds a certain level. The results obtained confirm the reliability of the method developed for employingf _min in aeronomic investigations.     

Recurrence of flare energy releases in solar active regions (Cycle 23)

The following results have been achieved in this work. The distribution of the recurrence times of solar flare events is generally lognormal. The typical flare recurrence times at the cycle 23 minimum and maximum are different: the average times (100–200 min) are most typical of the maximum; at the same time, the minimum is simultaneously characterized by short (several tens of seconds) and long (from several hundreds to a thousand of minutes). The minimal flare recurrence time tends to decrease in an active region with increasing sunspot group area in this region. The average flare recurrence times in an active region have typical values of 120^m, 210^m, 300^m, 400^m, and 530^m, which is close to the typical periods of long-period sunspot oscillations. The total number of flares in an active region depends on the sunspot area in this region and the flare energy release rate.     

Beginning of the new solar cycle (cycle 24) in the large-scale open solar magnetic field

It is proposed to determined minimums of the 11-year solar cycles based on a minimal flux of the large-scale open solar magnetic field. The minimal fluxes before the finished cycle 23 (Carrington rotation CR 1904) and the started cycle 24 (CR 2054, April 2007) were equal to 1.8 × 10²² and 1.2 × 10²² μs, respectively. The long-term tendency toward an approach to a deep minimum of solar activity is confirmed. On the assumption that magnetic flux variations from minimums to maximums are proportional to each other, the anticipated value of the maximal Wolf number during cycle 24 is estimated as W _max = 80.     

Solar eclipse of August 1, 2008, above Kharkov: 2. Observation results of wave disturbances in the ionosphere

Quasi-periodic variations in the power of incoherent scattered signals, caused by wave disturbances of the electron concentration in the ionosphere, are analyzed for the day of a partial solar eclipse and for a background day. The windowed and adaptive Fourier transforms and the wavelet transform are used for spectral analysis. The spectral parameters of the wave disturbances at altitudes of 100–500 km in the 10–120 min period range differed significantly on the day of the solar eclipse and on the background day. Variations in the spectrum began near the onset of the phase of maximum disk occultation and continued no less than 2 h. The amplitude of time variations N was 2 × 10⁹–4 × 10¹⁰ m^?3, and the relative amplitude was 0.10–0.15. Wave disturbances have been compared for five solar eclipses; the comparison shows a noticeable variation in the spectrum of the wave disturbances during these events.     

Reaction of the critical frequency of the <Emphasis Type="Italic">F</Emphasis>2 layer to a sharp depletion in atmospheric pressure

Changes in the critical frequencies of the F2 layer at several midlatitude stations of ionospheric vertical sounding during a sharp depletion in atmospheric pressure under quiet solar and geomagnetic conditions are analyzed. It is shown that in such periods, the observed foF2 values differ from the mean values by approximately 10–15% and the deviations from the mean could be both negative (in the daytime hours) and positive (at night). Such variations in foF2 could be referred to the known class of ionospheric disturbances observed under a quiet geomagnetic situation, that is, to the so-called “Q-disturbances.” Analysis of wavelet spectra of foF2 variations shows the presence in the F region of oscillations of various periods (from 0.5 to 10 days). The decrease in the amplitude of daily variations during pressure depletion is found. Presumably, the observed effect is caused by the dynamic impact of waves formed in the lower atmosphere on the ionospheric F2 layer.     

Solar wind velocity distribution on the heliospheric current sheet during Carrington rotations 1787-1795

The solar wind velocity distribution in the heliosphere is best represented using a v-map, where velocity contours are plotted in heliographic latitude-longitude coordinates. It has already been established that low-speed regions of the solar wind on the source surface correspond to the maximum bright regions of the K-corona and the neutral line of the coronal magnetic field. In this analysis, v-maps on the source surface for Carrington rotations (CRs) 1787-1795, during 1987, have been prepared using the interplanetary scintillation measurements at Research Institute of Atmospherics (RIA), Nagoya Univ., Japan. These v-maps were then used to study the time evolution of the low-speed (\leq450 km s⁻¹) belt of the solar wind and to deduce the distribution of solar wind velocity on the heliospheric current sheet. The low-speed belt of the solar wind on the source surface was found to change from one CR to the next, implying a time evolution. Instead of a slow and systematic evolution, the pattern of distribution of solar wind changed dramatically at one particular solar rotation (CR 1792) and the distributions for the succeeding rotations were similar to this pattern. The low-speed region, in most cases, was found to be close to the solar equator and almost parallel to it. However, during some solar rotations, they were found to be organised in certain longitudes, leaving regions with longitudinal width greater than 30^○ free of low-speed solar wind, i.e. these regions were occupied by solar wind with velocities greater than 450 km s⁻¹. It is also noted from this study that the low-speed belt, in general, followed the neutral line of the coronal magnetic field, except in certain cases. The solar wind velocity on the heliospheric current sheet (HCS) varied in the range 300–585 km s⁻¹ during the period of study, and the pattern of velocity distribution varied from rotation to rotation.     

Carriers of magnetism and vertical inhomogeneities in ultramafic-mafic rocks of the Kondyor Massif (Aldan shield)

Dunite samples from a borehole drilled in the platiniferous concentrically-zoned Kondyor Massif are studied by electron spin resonance (ESR). The spectrum profiles, relative intensities I, and volume magnetic susceptibilities κ are analyzed. These values experience strong irregular variations, sometimes by an order of magnitude, in the upper and medium parts of the column, at depths from 100 to 400 m; and the variations decay at greater depths. The magnetic properties of the samples are determined by iron (II) ions in the olivine lattice and by iron (III) ions in the magnetite and pyrrhotite microphases and in the products of breakdown of the solid solution: chromiferous magnetite, chromoferrite, etc. The I and κ values are directly related: κ_max = 27.8 × 10⁻³ SI units, κ_min=2.63 × 10⁻³ SI units, and κ_mean = 12.7 × 10⁻³ SI units. The maximum κ values are found in the zones with elevated contents of magnetite and pyrrhotite particles, and the minimum ones, in zones with few medium and small clusters with Fe³⁺ ions. The uneven distributions of solid solutions and magnetic phases over depths are suggested to be related to the disturbances in the conditions of crystallization.     

Comparison of model electron densities and temperatures with Millstone Hill observations during undisturbed periods and the geomagnetic storms of 16–23 March and 6–12 April 1990

Measurements of F-region electron density and temperature at Millstone Hill are compared with results from the IZMIRAN time-dependent mathematical model of the Earths ionosphere and plasmasphere during the periods 16–23 March and 6–12 April 1990. Each of these two periods included geomagnetically quiet intervals followed by major storms. Satisfactory agreement between the model and the data is obtained during the quiet intervals, provided that the recombination rate of O⁺(⁴S) ions was decreased by a factor of 1.5 at all altitudes during the nighttime periods 17–18 March, 19–20 March, 6–8 April and 8–9 April in order to increase the NmF2 at night better to match observations. Good model/data agreement is also obtained during the storm periods when vibrationally excited N₂ brings about factor-of-2-4 reductions in daytime NmF2. Model calculations are carried out using different expressions for the O⁺ – O collision frequency for momentum transfer, and the best agreement between the electron-density measurements and the model results is obtained when the CEDAR interim standard formula for the O⁺ – O collision frequency is used. Deviations from the Boltzmann distribution for the first five vibrational levels of NI were calculated. The calculated distribution is highly non-Boltzmann at vibrational levels j > 2, and the Boltzmann distribution assumption results in the increase of 10–30% in calculated NmF2 during the storm-time periods. During the March storm at solar maximum the model results obtained using the EUVAC solar flux model agree a little better with the observations in comparison with the EUV94 solar flux model. For the April storm period of moderate solar activity the EUV94X model results agree better with the observations in comparison to the EUVAC model.     

Compressional wave events in the dawn plasma sheet observed by Interball-1

Compressional waves with periods greater than 2 min (about 10–30 min) at low geomagnetic latitudes, namely compressional Pc5 waves, are studied. The data set obtained with magnetometer MIF-M and plasma analyzer instrument CORALL on board the Interball-1 are analyzed. Measurements performed in October 1995 and October 1996 in the dawn plasma sheet at −30 R_E ≤ X_GSM and |Z_GSM| ≤ 10 R_E are considered. Anti-phase variations of magnetic field and ion plasma pressures are analyzed by searching for morphological similarities in the two time series. It is found that longitudinal and transverse magnetic field variations with respect to the background magnetic field are of the same order of magnitude. Plasma velocities are processed for each time period of the local dissimilarity in the pressure time series. VeloCity disturbances occur mainly transversely to the local field line. The data reveal the rotation of the veloCity vector. Because of the field line curvature, there is no fixed position of the rotational plane in the space. These vortices are localized in the regions of anti-phase variations of the magnetic field and plasma pressures, and the vortical flows are associated with the compressional Pc5 wave process. A theoretical model is proposed to explain the main features of the nonlinear wave processes. Our main goal is to study coupling of drift Alfven wave and magnetosonic wave in a warm inhomogeneous plasma. A vortex is the partial solution of the set of the equations when the compression is neglected. A compression effect gives rise to a nonlinear soliton-like solution.     

Ring current and auroral electrojets in connection with interplanetary medium parameters during magnetic storm

The relationship between the auroral electrojet indices (AE) and the ring current magnetic field (DR) was investigated by observations obtained during the magnetic storm on 1–3 April 1973. During the storm main phase the DR development is accompanied by a shift of the auroral electrojets toward the equator. As a result, the standard AE indices calculated on the basis of data from auroral observatories was substantially lower than the real values (AE’). To determine AE’ during the course of a storm main phase data from subauroral magnetic observatories should be used. It is shown that the intensity of the indices (AE’) which take into account the shift of the electrojets is increased substantially relative to the standard indices during the storm main phase. AE’ values are closely correlated with geoeffective solar wind parameters. A high correlation was obtained between AE’ and the energy flux into the ring current during the storm main phase. Analysis of magnetic field variations during intervals with intense southward IMF components demonstrates a decrease of the saturation effect of auroral electrojet currents if subauroral stations magnetic field variations are taken into account. This applies both to case studies and statistical data. The dynamics of the electrojets in connection with the development of the ring current and of magnetospheric substorms can be described by the presence (absence) of saturation for minimum (maximum) AE index values during a 1-h interval. The ring current magnetic field asymmetry (ASY) was calculated as the difference between the maximum and minimum field values along a parallel of latitude at low latitudes. The ASY value is closely correlated with geoeffective solar wind parameters and simultaneously is a more sensitive indicator of IMF B_z variations than the symmetric ring current. ASY increases (decreases) faster during the main phase (the recovery phase) than DR. The magnetic field decay at low latitudes in the recovery phase occurs faster in the afternoon sector than at dusk.     

Effective barometric admittance and gravity residuals

In the analysis of surface gravity signals that may originate from the Earth's core, the step of correcting for the atmospheric pressure fluctuations is one that must be done carefully. We apply two techniques for determining the local, or effective, barometric admittance function between simultaneous observations of surface gravity and pressure. The first is a frequency domain fit that computes the admittance on a band-by-band basis. Using data from both the Canadian and French superconducting gravimeters we determined that the magnitude of the local, or background, admittance increases smoothly and monotonically from about 0.2 μgal mbar⁻¹ at long periods (> 10 days) to about 0.35 μgal mbar⁻¹ at frequencies greater than 3 cycles per day (c.p.d.); the phase lag is within a few degrees of 180°. By comparison, the effective admittances of the large-scale harmonics of the solar heating tide (S₁---S₇) are much smaller, between 0.1 and 0.3 μgal mbar⁻¹, for most of the harmonics of a day. In the second approach we fit a symmetrical time domain admittance function having lengths between 1 and 19 h using both a standard least-squares fit to a white noise residual and a new, and clearly superior, fit assuming a brown noise residual. Both time and frequency domain approaches give comparable results and contribute to a significant lowering of the residual level in non-tidal bands.     

Solar activity imprints in tree ring width from Chile (1610–1991)

We have investigated the solar activity signal in tree ring data from two locations in Chile. The tree ring time series extended over a period of ∼400 yr. Spectral and wavelet analysis techniques were employed. We have found evidence for the presence of the solar activity Schwabe (∼11 yr), Hale (∼22 yr), fourth-harmonic of the 208-yr Suess cycle (∼52 yr) and Gleissberg (∼80 yr) cycles. The Gleissberg cycle of tree ring data is in anti-phase with solar activity. Wavelet and cross-wavelet techniques revealed that the periods found are intermittent, possibly because solar activity signals observed in tree rings are mostly due to solar influence on local climate (rainfall, temperature, and cloud cover) where trees grow up. Further, cross-wavelet analysis between sunspot and tree ring time series showed that the cross power around the 11 yr solar cycle is more significant during periods of high solar activity (grand maximum) than during periods of low solar activity (grand minimum). As Glaciar Pio XI is practically at the Pacific Ocean level, the tree-ring response may be stronger due to the heating of the Pacific Ocean water following an increase of the solar radiation incidence rather than at the higher altitudes of Osorno region.     

Observations of the response time of high-latitude ionospheric convection to variations in the interplanetary magnetic field using EISCAT and IMP-8 data

We have combined ∼300 h of tristatic measurements of the field-perpendicular F region ionospheric flow measured overhead at Tromsø by the EISCAT UHF radar, with simultaneous IMP-8 measurements of the solar wind and interplanetary magnetic field (IMF) upstream of the Earth’s magnetosphere, in order to examine the response time of the ionospheric flow to changes in the north-south component of the IMF (B_z). In calculating the flow response delay, the time taken by field changes observed by the spacecraft to first effect the ionosphere has been carefully estimated and subtracted from the response time. Two analysis methods have been employed. In the first, the flow data were divided into 2 h-intervals of magnetic local time (MLT) and cross-correlated with the “half-wave rectifier” function V²B_s, where V is the solar wind speed, and B_s is equal to IMF B_z if the latter is negative, and is zero otherwise. Response delays, determined from the time lag of the peak value of the cross-correlation coefficient, were computed versus MLT for both the east-west and north-south components of flow. The combined data set suggests minimum delays at ∼1400 MLT, with increased response times on the nightside. For the 12-h sector centred on 1400 MLT, the weighted average response delay was found to be 1.3 ± 0.8 min, while for the 12-h sector centred on 0200 MLT the weighted average delay was found to increase to 8.8 ± 1.7 min. In the second method we first inspected the IMF data for sharp and enduring (at least ∼5 min) changes in polarity of the north-south component, and then examined concurrent EISCAT flow data to determine the onset time of the corresponding enhancement or decay of the flow. For the case in which the flow response was timed from whichever of the flow components responded first, minimum response delays were again found at ∼1400 MLT, with average delays of 4.8 ± 0.5 min for the 12-h sector centred on 1400 MLT, increasing to 9.2 ± 0.8 min on the nightside. The response delay is thus found to be reasonably small at all local times, but typically ∼6 min longer on the nightside compared with the dayside. In order to make an estimate of the ionospheric information propagation speed implied by these results, we have fitted a simple theoretical curve to the delay data which assumes that information concerning the excitation and decay of flow propagates with constant speed away from some point on the equatorward edge of the dayside open-closed field line boundary, taken to lie at 77° magnetic latitude. For the combined cross-correlation results the best-fit epicentre of information propagation was found to be at 1400 MLT, with an information propagation phase speed of 9.0 km s⁻¹. For the combined event analysis, the best-fit epicentre was also found to be located at 1400 MLT, with a phase speed of 6.8 km s⁻¹.     

Spectral energy contributions of quasi-periodic oscillations (2–35 days) to the variability of the f oF2

The relative contributions of quasi-periodic oscillations from 2 to 35 days to the variability of f_oF2 at middle northern latitudes between 42°N and 60°N are investigated. The f_oF2 hourly data for the whole solar cycle 21 (1976–1986) for four European ionospheric stations Rome (41.9°N, 12.5°E), Poitiers (46.5°N, 0.3°E), Kaliningrad (54.7°N, 20.6°E) and Uppsala (59.8°N, 17.6°E) are used for analysis. The relative contributions of different periodic bands due to planetary wave activity and solar flux variations are evaluated by integrated percent contributions of spectral energy for these bands. The observations suggest that a clearly expressed seasonal variation of percent contributions exists with maximum at summer solstice and minimum at winter solstice for all periodic bands. The contributions for summer increase when the latitude increases. The contributions are modulated by the solar cycle and simultaneously influenced by the long-term geomagnetic activity variations. The greater percentage of spectral energy between 2 to 35 days is contributed by the periodic bands related to the middle atmosphere planetary wave activity.     

Statistical analysis of earthquake activity at Etna volcano (March 1981 eruption)

Seismic activity that preceded, accompanied, and followed the 17–23 March 1981 Etnean eruption has been statistically analyzed.On the grounds of both time evolution of seismicity and catalogue completeness, three time intervals have been defined (12 February–2 March, 12–17 March, 19–31 March) and for each of these periods both the b coefficient of theGutenberg-Richter's (1956) relationship and the E parameter of the cluster size (Shlien andToksoz, 1970) have been calculated.No significant variations were observed between the first and third periods, while lower values of bothb andE coefficients were found in the second one. These findings might indicate that changes in the seismicity features occur just before the eruption start.Small but fast variations in the stress field acting on the volcano might originate this type of seismic activity, while the importance of the tectonic control on volcanic phenomena seems to be confirmed.