On mid-term periodicities of high-latitude solar activity

The mid-term periodicities of polar faculae are studied separately for the total disk, northern and southern hemispheres of the Sun for a time interval from 1951 August to 1998 December. Apart from the 11-year Schwabe cycle which is the fundamental period and is found in all of the three time series, the following prominent results are found: (1) the rotational periodicity of solar activity at high latitudes is approximately from 28 to 32 days; (2) a large number of quasi-periods appearing in low-latitude solar activity (annual variation, 1.3–1.7 years, quasi-biennial oscillation, and 4–5 years) also exist in polar faculae; (3) the periodicities on both hemispheres are not identical.     

A morphological study of CMEs using wavelet analysis

The principal objective of this study is to analyze structures of Coronal Mass Ejections (CMEs) using a wavelet technique. We use data measured with the SOHO/LASCO coronographs C2–C3, EIT and STEREO COR1A–B, COR2A–B. We have found that different structures show up in a CME at different spatial scales of wavelets. We also study the orientation of the most intense flux and find the orientation of the structures in the CME.     

Comparison of GPS TEC measurements with IRI-2007 TEC prediction over the Kenyan region during the descending phase of solar cycle 23

This paper presents an analysis of the Total Electron Content (TEC) derived from the International GNSS Service receiver (formerly IGS) at Malindi (2.9°S, 40.1°E), Kenya for the periods 2004–2006 during the declining phase of solar cycle 23. The diurnal, monthly and seasonal variations of the TEC are compared with TEC from the latest International Reference Ionosphere model (IRI-2007). The GPS–TEC exhibits features such as an equatorial noon time dip, semi-annual variations, Equatorial Ionization Anomaly and day-to-day variability. The lowest GPS–TEC values are observed near the June solstice and September equinox whereas largest values are observed near the March equinox and December solstice. The mean GPS–TEC values show a minimum at 03:00 UT and a peak value at about 10:00 UT. These results are compared with the TEC derived from IRI-2007 using the NeQuick option for the topside electron density (IRI–TEC). Seasonal mean hourly averages show that IRI-2007 model TEC values are too high for all the seasons. The high prediction primarily occur during daytime hours till around midnight hours local time for all the seasons, with the highest percentage deviation in TEC of more 90% seen in September equinox and lowest percentage deviation in TEC of less than 20% seen in March equinox. Unlike the GPS–TEC, the IRI–TEC does not respond to geomagnetic storms and does overestimate TEC during the recovery phase of the storm. While the modeled and observed data do correlate so well, we note that IRI-2007 model is strongly overestimating the equatorial ion fountain effect during the descending phase of solar cycle, and this could be the reason for the very high TEC estimations.     

Attempt to identify a source mechanism of Mercury’s sodium exosphere by a spectrometer using Fabry–Perot etalon

The Mercury’s Sodium Atmosphere Spectral Imager (MSASI) on BepiColombo (BC) will address a range of fundamental scientific questions pertaining to Mercury’s exosphere. The measurements will provide new information on regolith–exosphere–magnetosphere coupling as well as new understanding of the dynamics governing the exosphere bounded by the planetary surface, the solar wind and interplanetary space. MSASI is a high-dispersion visible spectrometer working in the spectral range around sodium D2 emission (589 nm). A tandem Fabry–Perot etalon is used to achieve a compact design. We presents a design of the spectral analyzer using Fabry–Perot interferometer. We conclude that: (1) The MSASI optical design is practical and can be implemented without new or critical technology developments; (2) The thermally-tuned etalon design is based on concepts, designs and materials that have good space heritage.     

Modeling and observations of the north–south ionospheric asymmetry at low latitudes at long deep solar minimum

Using the physics based model SUPIM and FORMOSAT-3/COSMIC electron density data measured at the long deep solar minimum (2008–2010) we investigate the longitude variations of the north–south asymmetry of the ionosphere at low latitudes (±30° magnetic). The data at around diurnal maximum (12:30–13:30 LT) for magnetically quiet (Ap ? 15) equinoctial conditions (March–April and September–October) are presented for three longitude sectors (a) 60°E–120°E, (b) 60°W–120°W and (c) 15°W–75°W. The sectors (a) and (b) have large displacements of the geomagnetic equator from geographic equator but in opposite hemispheres with small magnetic declination angles; and sector (c) has large declination angle with small displacement of the equators; vertical E × B drift velocities also have differences in the three longitude sectors. SUPIM investigates the importance of the displacement of the equators, magnetic declination angle, and E × B drift on the north–south asymmetry. The data and model qualitatively agree; and indicate that depending on longitudes both the displacement of the equators and declination angle are important in producing the north–south asymmetry though the displacement of the equators seems most effective. This seems to be because it is the displacement of the equators more than the declination angle that produces large north–south difference in the effective magnetic meridional neutral wind velocity, which is the main cause of the ionospheric asymmetry. For the strong control of the neutral wind, east–west electric field has only a small effect on the longitude variation of the ionospheric asymmetry. Though the study is for the long deep solar minimum the conclusions seem valid for all levels of solar activity since the displacement of the equators and declination angle are independent of solar activity.     

Planetary wave coupling of the atmosphere–ionosphere system during the Northern winter of 2008/2009

This paper presents the global spatial (latitude and altitude) structure and temporal variability of the ∼23-day ionospheric zonally symmetric (s = 0) planetary wave (PW) seen in the Northern winter of 2008/2009 (October 2008–March 2009). It is shown that these ∼23-day ionospheric oscillations are forced from PWs propagating from below. The COSMIC ionospheric parameters f_oF2 and h_mF2 and electron density at fixed altitudes and the SABER temperatures were utilized in order to define the waves which are present simultaneously in the atmosphere and ionosphere. The long-period PWs from the two data sets have been extracted through the same data analysis method. The similarity between the lower thermospheric ∼23-day (s = 0) temperature PW and its ionospheric electron density response provides valuable and strong experimental evidence for confirming the paradigm of atmosphere–ionosphere coupling.     

Effects of sector structure of the interplanetary magnetic field on the upper mesosphere–lower thermosphere dynamics

In this paper we study the influence of the interplanetary magnetic field (IMF) polarity changes caused by the Earth passing through the IMF sector boundary on the dynamic processes taking place in neutral atmosphere within the altitude interval of the upper mesosphere–lower thermosphere (83–101 km). The analysis has revealed the influence of the IMF sector structure on dynamics of the upper mesosphere–lower thermosphere. There has been a significant seasonal variation of the wind reaction to the IMF polarity changes observed. The influence of the IMF polarity changes on neutral atmosphere dynamics within the altitude range of 83–101 km is most pronounced in the zonal component of neural wind when the IMF polarity changes from negative to positive in all the seasons except for spring and when IMF polarity changes from positive to negative – in spring only.     

Morphology of the G condition occurrence over Irkutsk

Using Irkutsk digisonde data obtained in 2003–2011, a morphological analysis of the G condition occurrence has been made. The G condition was found to occur during daylight hours in summer; in winter, it is extremely rare, and its appearance is associated with intense magnetic storms. In the years of moderate solar activity, the G condition is most frequently registered at Kp ? 4, in the forenoon. During low solar activity, it can be observed under quiet geomagnetic conditions; in most cases, local time of its appearance shifts to afternoon hours. The highest percentage of the G condition occurrence (7.7–6.4%) was recorded in June and July 2008 when the levels of solar and geomagnetic activity were abnormally low.     

Generation of secondary waves due to intensive large-scale AGW traveling

By using data from GPS receivers we detected huge-amplitude solitary large-scale traveling acoustic-gravity waves (LS AGW) which manifested themselves as perturbations of total electron content (TEC) of duration of about 40 min. Originated in the auroral area after significant alterations of geomagnetic field intensity during geomagnetic storms on 29–30 October 2003, LS disturbances propagated with a velocity about 1000–1200 m/s and caused generation of secondary small-scale (SS) waves with time period of 2–10 min. Such SS structure followed the solitary intensive AGW at a distance more than 4000 km. However, we observed such phenomenon only within the territory with high values of “vertical” TEC and steep gradients of TEC. Apparently, these conditions are necessary for generation of SS waves due to propagation of LS AGW.     

Interseismic slip rate of the Garze–Yushu fault belt in the Tibetan Plateau from C-band InSAR observations between 2003 and 2010

InSAR time series techniques can provide high-spatial resolution deformation fields across an active fault belt, even for zones with heavy vegetation coverage. An interseismic deformation map across the Garze–Yushu fault belt in the Tibetan Plateau, ∼300 km by ∼100 km, is derived from C-band Envisat/ASAR imagery collected between 2003 and 2010. Unlike previous research, we obtain a lookup figure which relates the slip rate with the fault locking depth, the dip angle and the rake angle. The estimated slip rate changes significantly with the locking depth and the rake angle, but relatively little with the dip angle. When considering the focal mechanism solutions of historical earthquake along the Garze–Yushu fault, the interseismic slip rate of the Garze–Yushu fault is close to a value of 6.4 mm/yr, which is between the highest (18.2 mm/yr) and the lowest (3.1 mm/yr) slip rate from GPS estimations, but slightly less than the minimum value (∼ 7 mm/yr) from the geological estimations. The earthquake recurrence interval on the Yushu part of Garze–Yushu fault equals 272 yr, and the April 14, 2010 Mw 6.9 earthquake has not completely released the accumulated strain energy between 1738 and 2010.     

Observation of Mercury’s sodium tail using Fabry–Perot Interferometer

We observed sodium emission from Mercury’s atmosphere using a Fabry–Perot Interferometer at Haleakala Observatory on June 14, 2006. The Fabry–Perot Interferometer was used as a wavelength-tunable filter. The spectra of the surface reflection were subtracted from the observed spectra because sodium emission is contaminated by the surface reflection of Mercury. The image obtained in our observation shows the sodium exosphere extended to the anti-solar direction. The lifetime of a sodium atom was estimated to be 1.6 × 10⁴ to 1.9 × 10⁵ s with an error by a factor of 3–4.     

Prediction of horizontal gas–solid flows under different gravitational fields

In this paper the performance of horizontal pneumatic conveying under different gravity environments is evaluated. An Euler–Lagrange approach validated versus ground experiments is employed to predict the relevant particle variables such as particle mass flux, mean conveying and fluctuating velocities in terrestrial, lunar and micro-gravity conditions. Gravity reduced computations predict a reduction in the global particle–wall collision frequency. Also, in the case of low wall roughness and small particle mass loading, reduction of gravity acceleration implies an increase of particle–wall collision frequency with the upper wall of the channel affecting greatly the particle mass flux profile. In the case of high wall roughness and/or high particle-to-fluid mass loading (i.e., around 1.0) particle conveying characteristics are similar in the three gravity conditions evaluated. This is due to the fact that both, wall roughness and inter-particle collisions reduce gravitational settling. However, the influence of gravity on the additional pressure loss along the channel due to the conveying of the particles is much reduced.     

Ion-acoustic waves in unmagnetized collisionless weakly relativistic plasma of warm-ion and isothermal-electron using time-fractional KdV equation

Collisionless unmagnetized plasma consisting of a mixture of warm ion-fluid and isothermal-electron is considered, assuming that the ion flow velocity has a weak relativistic effect. The reductive perturbation method has been employed to derive the Korteweg–de Vries (KdV) equation for small – but finite-amplitude electrostatic ion-acoustic waves in this plasma. The semi-inverse method and Agrawal’s method lead to the Euler–Lagrange equation that leads to the time fractional KdV equation. The variational-iteration method given by He is used to solve the derived time fractional KdV equation. The calculations show that the fractional order may play the same rule of higher order dissipation in KdV equation to modulate the soliton wave amplitude in the plasma system. The results of the present investigation may be applicable to some plasma environments, such as space-plasmas, laser-plasma interaction, plasma sheet boundary layer of the earth’s magnetosphere, solar atmosphere and interplanetary space.     

The Earth’s magnetosphere response to interplanetary medium conditions on January 21–22, 2005 and on December 14–15, 2006

The Earth’s magnetosphere response to interplanetary medium conditions on January 21–22, 2005 and on December 14–15, 2006 has been studied. The analysis of solar wind parameters measured by ACE spacecraft, of geomagnetic indices variations, of geomagnetic field measured by GOES 11, 12 satellites, and of energetic particle fluxes measured by POES 15, 16, 17 satellites was performed together with magnetospheric modeling based in terms of A2000 paraboloid model. We found the similar dynamics of three particle populations (trapped, quasi-trapped, and precipitating) during storms of different intensities developed under different external conditions: the maximal values of particle fluxes and the latitudinal positions of the isotropic boundaries were approximately the same. The main sources caused RC build-up have been determined for both magnetic storms. Global magnetospheric convection controlled by IMF and substorm activity driven magnetic storm on December 14–15, 2006. Extreme solar wind pressure pulse was mainly responsible for RC particle injection and unusual January 21, 2005 magnetic storm development under northward IMF during the main phase.     

Reminiscences of cosmic ray research in Mexico

Cosmic ray research in Mexico dates from the early 1930s with the work of the pioneering physicist, Manuel Sandoval Vallarta and his students from Mexico. Several experiments of international significance were carried out during that period in Mexico: they dealt with the geomagnetic latitude effect, the north–south and west–east asymmetry of cosmic ray intensity, and the sign of the charge of cosmic rays. The international cosmic ray community has met twice in Mexico for the International Cosmic Ray Conferences (ICRC): the fourth was held in Guanajuato in 1955, and the 30th took place in Mérida, in 2007. In addition, an international meeting on the Pierre Auger Collaboration was held in Morelia in 1999, and the International Workshop on Observing UHE Cosmic Rays took place in Metepec in 2000. A wide range of research topics has been developed, from low-energy Solar Energetic Particles (SEP) to the UHE. Instrumentation has evolved since the early 1950s, from a Simpson type neutron monitor installed in Mexico City (2300 m asl) to a solar neutron telescope and an EAS Cherenkov array, (within the framework of the Auger International Collaboration), both at present operating on Mt. Sierra La Negra in the state of Puebla (4580 m asl). Research collaboration has been undertaken with many countries; in particular, the long-term collaboration with Russian scientists has been very fruitful.     

Effect of solar activity on the repetitiveness of some meteorological phenomena

In this paper we research the relationship between solar activity and the weather on Earth. This research is based on the assumption that every ejection of magnetic field energy and particles from the Sun (also known as Solar wind) has direct effects on the Earth’s weather. The impact of coronal holes and active regions on cold air advection (cold fronts, precipitation, and temperature decrease on the surface and higher layers) in the Belgrade region (Serbia) was analyzed. Some active regions and coronal holes appear to be in a geo-effective position nearly every 27 days, which is the duration of a solar rotation. A similar period of repetitiveness (27–29 days) of the passage of the cold front, and maximum and minimum temperatures measured at surface and at levels of 850 and 500 hPa were detected. We found that 10–12 days after Solar wind velocity starts significantly increasing, we could expect the passage of a cold front. After eight days, the maximum temperatures in the Belgrade region are measured, and it was found that their minimum values appear after 12–16 days. The maximum amount of precipitation occurs 14 days after Solar wind is observed. A recurring period of nearly 27 days of different phases of development for hurricanes Katrina, Rita and Wilma was found. This analysis confirmed that the intervals of time between two occurrences of some particular meteorological parameter correlate well with Solar wind and A index.     

Regular features of the polar ionosphere characteristics from Digisonde measurements over Norilsk

Regular features of the polar ionosphere have been studied using its local empirical model of the electron density distribution in the bottomside ionosphere. The local empirical model was derived from the hand-scaled ionogram data recorded by DPS-4 Digisonde at Norilsk, Russia (69.4N, 88.1E; 60N GLAT, 166E GLON) for a 6-year period from December, 2002 to December, 2008. The paper describes the technique used to build the local empirical model and discusses its diurnal, seasonal, and solar activity specifications in comparison with the standard IRI-2007 climatological model for the same period of time, long-term observations from the European Incoherent Scatter UHF radar (1988–1999), and the high-latitude ionosondes data. Primary focus of the paper is behavior of the three F2 layer parameters: the F2 peak density (NmF2), the peak height (hmF2) and the bottomside thickness (B0). Special emphasis of the paper is the analysis of the winter anomaly manifestation at Norilsk and the peculiar diurnal–seasonal behavior of hmF2 under low solar activity, named as a “polar day effect”.     

Generation of vertical–horizontal and horizontal–horizontal gravity gradients using stochastically modified integral estimators

The Earth’s gravity field modelling is an ill-posed problem having a sensitive solution to the error of data. Satellite gravity gradiometry (SGG) is a space technique to measure the second-order derivatives of geopotential for modelling this field, but the measurements should be validated prior to use. The existing terrestrial gravity anomalies and Earth gravity models can be used for this purpose. In this paper, the second-order vertical–horizontal (VH) and horizontal–horizontal (HH) derivatives of the extended Stokes formula in the local north-oriented frame are modified using biased, unbiased and optimum types of least-squares modification. These modified integral estimators are used to generate the VH and HH gradients at 250 km level for validation purpose of the SGG data. It is shown that, unlike the integral estimator for generating the second-order radial derivative of geopotential, the system of equations from which the modification parameters are obtained is unstable for all types of modification, with large cap size and high degree, and regularization is strongly required for solving the system. Numerical studies in Fennoscandia show that the SGG data can be estimated with an accuracy of 1 mE using an integral estimator modified by a biased type least-squares modification. In this case an integration cap size of 2.5° and a degree of modification of 100 for integrating 30′ × 30′ gravity anomalies are required.     

Two-manoeuvres transfers between LEOs and Lissajous orbits in the Earth–Moon system

The purpose of this work is to compute transfer trajectories from a given Low Earth Orbit (LEO) to a nominal Lissajous quasi-periodic orbit either around the point L₁ or the point L₂ in the Earth–Moon system. This is achieved by adopting the Circular Restricted Three-Body Problem (CR3BP) as force model and applying the tools of Dynamical Systems Theory.     

Lunar apex–antapex cratering asymmetry as an impactor recorder in the Earth–Moon system

The degree of apex–antapex cratering asymmetry of a synchronously rotating satellite primarily depends on the mean encounter velocity of impactors with respect to the planetary system and the orbital velocity of the satellite. This means that we can estimate the mean encounter velocity of impactors by observing the apex–antapex cratering asymmetry, if the relationship between these is known. To apply this technique to the Moon, we attempt to derive the relationship between the mean encounter velocity of impactors and the degree of the lunar cratering asymmetry as a function of time, considering the temporal variation in the lunar orbital velocity during the last 4.0 Gyr. We used the cratering asymmetry of Zahnle et al. [Zahnle, K., Schenk, P., Sobieszczyk, S. et al. Differential cratering of synchronously rotating satellites by ecliptic comets. Icarus 153, 111–129, 2001] to obtain the relationship. Applying this relationship enables us to estimate the impactor’s velocity of the Earth–Moon system from an investigation of the spatial distribution of lunar craters. Furthermore, we re-evaluate the cratering asymmetry’s influence on lunar cratering chronology.