Antarctic meteor observations using the Davis MST and meteor radars

This paper presents the meteor observations obtained using two radars installed at Davis (68.6°S, 78.0°E), Antarctica. The Davis MST radar was installed primarily for observation of polar mesosphere summer echoes, with additional transmit and receive antennas installed to allow all-sky interferometric meteor radar observations. The Davis meteor radar performs dedicated all-sky interferometric meteor radar observations. The annual count rate variation for both radars peaks in mid-summer and minimizes in early Spring. The height distribution shows significant annual variation, with minimum (maximum) peak heights and maximum (minimum) height widths in early Spring (mid-summer). Although the meteor radar count rate and height distribution variations are consistent with a similar frequency meteor radar operating at Andenes (69.3°N), the peak heights show a much larger variation than at Andenes, while the count rate maximum-to-minimum ratios show a much smaller variation. Investigation of the effects of the temporal sampling parameters suggests that these differences are consistent with the different temporal sampling strategies used by the Davis and Andenes meteor radars. The new radiant mapping procedure of [Jones, J., Jones, W., Meteor radiant activity mapping using single-station radar observations, Mon. Not. R. Astron. Soc., 367(3), 1050–1056, doi: 10.1111/j.1365-2966.2006.10025.x, 2006] is investigated. The technique is used to detect the Southern delta-Aquarid meteor shower, and a previously unknown weak shower. Meteoroid speeds obtained using the Fresnel transform are presented. The diurnal, annual, and height variation of meteoroid speeds are presented, with the results found to be consistent with those obtained using specular meteor radars. Meteoroid speed estimates for echoes identified as Southern delta-Aquarid and Sextantid meteor candidates show good agreement with the theoretical pre-atmospheric speeds of these showers (41 km s⁻¹ and 32 km s⁻¹, respectively). The meteoroid speeds estimated for these showers show decreasing speed with decreasing height, consistent with the effects of meteoroid deceleration. Finally, we illustrate how the new radiant mapping and meteoroid speed techniques can be combined for unambiguous meteor shower detection, and use these techniques to detect a previously unknown weak shower.     

Possible physical properties of the flow in the vicinity of the heliopause

An interface between the fully ionized hydrogen plasma of the solar wind (SW) and the partially ionized hydrogen gas flow of the local interstellar medium (LISM) is formed as a region where there is a strong interaction between these two flows. The interface is bounded by the solar wind termination shock (TS) and the LISM bow shock (BS) and is separated on two regions by the heliopause (HP) separating the solar wind and charged component of the LISM (plasma component below). The BS is formed due to the deceleration of the supersonic LISM flow relative to the solar system. Regions of the interface between the TS and HP and between the HP and BS were in literature named as the inner and outer heliosheaths, respectively. An investigation of the structure and physical properties of the heliosheath is at present especially interested due to the fact that Voyager-1 and Voyager-2 have crossed the TS in December 2004 (Burlaga, L.F., Ness, N.F., Acuna, M.Y., et al. Crossing the termination shock into the the heliosheath. Magnetic fields. Science 309, 2027–2029, 2005; Fisk, L.A. Journey into the unknown beyond. Science 309, 2016–2017, 2005; Decker, R.B., Krimigis, S.M., Roelof, E.C., et al. Voyager 1 in the foreshock, termination shock and heliosheath. Science 309, 2020–2024, 2005; Stone, E.C., Cummings, A.C., McDonald, F.B., et al. Voyager 1 explores the termination shock region and the heliosheath beyond. Science 309, 2017–2020, 2005) and in September 2007 (Jokipii, J.R. A shock for Voyager 2. Nature 454, 38–39, 2008; Gurnett, D.A., Kurth, W.S. Intense plasma waves at and near the solar wind termination shock. Nature 454, 78–80, 2008. doi: 10.1038/nature07023; Wang, L., Lin, R.P., Larson, D.E., Luhmann, J.G. Domination of heliosheath pressure by shock-accelerated pickup ions from observations of neutral atoms. Nature 454, 81–83, 2008. doi: 10.1038/nature07068.14; Burlaga, L.F., Ness, N.F., Acuna, M.H., et al. Magnetic fields at the solar wind termination shock. Nature 454, 75–77, 2008. doi: 10.1038/nature07029; Richardson, J.D., Kasper, J.C., Wang, C., et al. Cool heliosheath plasma and deceleration of the upstream solar wind at the termination shock. Nature 454, 63–66, 2008. doi: 10.1038/nature07024; Stone, E.C., Cummings, A.C., McDonald, F.B., et al. An asymmetric solar wind termination shock. Nature 454, 71–74, 2008. doi: 10.1038/nature07022; Decker, R.B., Krimigis, S.M., Roelof, E.C., et al. Mediation of the solar wind termination shock by non-thermal ions. Nature 454, 67–70, 2008. doi: 10.1038/nature 07030), respectively, and entered to the inner heliosheath.     

Two-spacecraft observations of reconnection at the magnetopause: Model results and data comparison

We revisit an example of “quasi-steady” magnetic reconnection at the dayside magnetopause on February 11, 1998, observed by Equator-S and Geotail at the dawnside magnetopause. Phan et al. [Phan, T.D. et al., 2000. Extended magnetic reconnection at the Earth’s magnetopause from detection of bi-directional jets. Nature 404, 848–850.] reported oppositely directed jets at these spacecrafts and inferred a length of the reconnection line of about 38R_E. Pinnock et al. [Pinnock, M., Chisham, G., Coleman, I.J., Freeman, M.P., Hairston, M., Villain, J.-P., 2003. The location and rate of dayside reconnection during an interval of southward interplanetary magnetic field. Ann. Geophys. 21, 1467–1482.] used measurements from SuperDARN radars to show that the reconnection electric field was variable. Here we complement this work by obtaining snapshots of the reconnection electric field from the in situ observations. To do this, we apply a reconstruction method based on a model of compressible Petschek-type magnetic reconnection. This independent method uses magnetic field observations as input data to calculate the reconnection electric field. We obtain average values of E_rec in the range of 0.4–2.4 mV/m. Further we infer a distance perpendicular to the reconnection line of 0.4–0.6R_E. The model results are compared with the two studies mentioned above. It thus appears that while the transfer of momentum for this event is indeed large-scale, the actual rate depends on the time it is measured.     

Investigation of precursory phenomena in the ionosphere,atmosphere and groundwater before large earthquakes of M > 6.5

Precursory phenomena in the ionosphere, atmosphere and groundwater before large earthquakes (M > 6.5) are extensively investigated toward the earthquake prediction. Upward tornado type seismic clouds occurred near the epicenter associated with strong LF-VLF radio noises from lightning discharges in the evening of January 9, 1995 [Yamada, T., Oike, K. On the increase of electromagnetic noises before and after the 1995 Hyogo-Ken Nanbu earthquake. In: Hayakawa M. (Ed.), Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes. TERRAPUB, Tokyo, pp. 417–427, 1999] and anomalous foEs increases up to 10 MHz were detected at Shigaraki, 90 km of the epicenter and at Kokubunji, 500 km east of the epicenter [Ondoh, T. Anomalous sporadic-E layers observed before M7.2 Hyogo-ken Nanbu earthquake; Terrestrial gas emanation model. Adv. Polar Upper Atmos. Res. 17, 96–108, 2003; Ondoh, T. Anomalous sporadic-E ionization before a great earthquake, Adv. Space Research 34, 1830–1835, 2004] associated with strong ELF noises from lightning discharges in the daytime on January 15, 1995 [Hata, M., Fujii, T., Takumi, I. EM precursor of large-scale earthquakes in Japan, in: Abstracts of International Workshop on Seismo Electromagnetics (IWSE 2005), Univ. Electro-Communications, Chofu, Tokyo, Japan, March 15–17, pp. 182–186, 2005] before the M7.2 Hyogoken–Nanbu earthquake of January 17, 1995. The anomalous foEs increases occurred at epicentral distances within 500 km that are the same as those of the terrestrial gas emanations along active faults before large earthquakes [King, C.-Y. Gas geochemistry applied to earthquake prediction: An overview. J. Geophys. Res. 91(B12), 12269–12281, 1986]. The anomalous foEs increases seem to be a seismic precursor because geomagnetic and solar conditions were very quiet all day on January 15,1995 and the normal foEs in Japanese winter is below 6 MHz. No significant pre-seismic geomagnetic field variation was detected at epicentral distance of 100 km before this earthquake [Ondoh, T., Hayakawa, M. Anomalous occurrence of sporadic-E layers before the Hyogoken–Nanbu earthquake, M7.2 of January 17, 1995. In: Hayakawa, M. (Ed.), Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes, TERRAPUB, Tokyo, pp. 629–639, 1999; Ondoh, T., Hayakawa, M. Seismo discharge model of anomalous sporadic E ionization before great earthquakes. In: Hayakawa, M., O.A. Molchanov, (Eds.), Seismo Electromagnetics: Lithosphere–Atmosphere–Ionosphere Couplings, TERRAPUB, Tokyo, pp. 385–390, 2002; Ondoh. T., Hayakawa, M. Synthetic study of precursory phenomena of the M7.2 Hyogo-ken Nanbu earthquake. Phys. Chem. Earth 31, 378–388, 2006]. The foF2 decrease and h’F increase occurred before the M7.8 Hokkaido Nansei-Oki earthquake of July 12,1993 in a geomagnetic quiet period [Ondoh, T. Ionospheric disturbances associated with great earthquake of Hokkaido southwest coast, Japan of July 12, 1993. Phys. Earth Planet. Interiors. 105, 261–269, 1998; Ondoh, T. Seismo ionospheric phenomena. Adv. Space Res. 26, 8, 1267–1272, 2000]. Characteristic phase changes at terminator times of Omega 10.2 kHz waves passing 70 km of the epicenter extended toward darker local times by 1 h for 3 days before this earthquake due to lowering of the wave reflection height or ion density increases in the D region [Hayakawa, M., Molchanov, O. A., Ondoh, T., Kawai, E. The precursory signature effect of the Kobe earthquake on VLF subionospheric signals. J. Commun. Res. La., 43, 00. 169–180, 1996]. The radon concentration in the atmosphere over Ashiya fault, Kobe [Yasuoka, Y., Shinogi, M. Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe, Japan, earthquake. Health Phys. 72(5), 759–761, 1997] and in the groundwater at 17 m well in Nishinomiya, Japan [Igarashi, G., Saeki, S., Takahata, N., Sumikawa, K., Tasaki, S., Sasaki, Y., Takahashi, M., Sano, Y. Ground-water radon anomaly before the Kobe earthquake in Japan. Science 269, 60–61, 1995] had gradually increased since 2 months before the M7.2 earthquake, increased suddenly in December 1994, and rapidly returned to the normal low level of October, 1994 [Yasuoka, Y., Shinogi, M. 1997. Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe. Japan, earthquake. Health Phys. 72(5), 759–761.]. Radon concentration changes in the groundwater before the M 7.0 Izu-Oshima-kinkai earthquake, Japan on January 14, 1978 [Wakita, H., Nakamura, Y., Notsu, K., Noguchi, M., Asada, T. 1980. Radon anomaly: a possible precursor of the 1978 Izu-Oshima-kinkai earthquake. Science 207, 882–883] and the M6.8 Chengkung earthquake, Taiwan on December 10, 2003 [Kuo, T., Fan, K., Chen, W., Kuochen, H., Han, Y., Wang, C., Chang, T., Lee, Y. Radon anomaly at the Antung Hot Spring before the Taiwan M6.8 Chengkung earthquake. Proceedings, Thirty-First Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, January 30–February 1, 2006, SGP-TR-179, 2006] are also investigated to find common features of the groundwater radon concentration changes before large earthquakes (M > 6.5) in comparison with those before the M7.2 Hyogoken–Nanbu earthquake. Groundwater radon concentrations before the 3 large earthquakes had shown common characteristic changes of gradually initial ones from the normal level since about 2 months before the earthquake onsets, rapid decreases down to the minimum, and quick increases up to the maximum at 7–20 days before the earthquake onsets, respectively. These are very useful characteristics of pre-seismic radon anomaly for the earthquake prediction or warning. Promising observations toward the earthquake prediction are also discussed.     

Plasma and electromagnetic wave simulations of meteors

Every day billions of meteoroids impact and disintegrate in the Earth’s atmosphere. Current estimates for this global meteor flux vary from 2000 to 200,000 tons per year, and estimates for the average velocity range between 10 km/s and 70 km/s. The basic properties of this global meteor flux, such as the average mass, velocity, and chemical composition remain poorly constrained. We believe much of the mystery surrounding the basic parameters of the interplanetary meteor flux exists for the following reason, the unknown sampling characteristics of different radar meteor observation techniques, which are used to derive or constrain most models. We believe this arises due to poorly understood radio scattering characteristics of the meteor plasma, especially in light of recent work showing that plasma turbulence and instability greatly influences meteor trail properties at every stage of evolution. We present our results on meteor plasmas simulations of head echoes using particle in cell (PIC) ions, which show that electric fields strongly influence early stage meteor plasma evolution, by accelerating ions away from the meteoroid body. We also present the results of finite difference time domain electromagnetic simulations (FDTD), which can calculate the radar cross section of the simulated meteor plasmas. These simulations have shown that the radar cross section depends in a complex manner on a number of parameters. These include the angle between radar and meteor entry, a large dependence on radar frequency, which shows that for a given meteor plasma size and density, the reflectivity as a function of probing radar frequency varies, but typically peaks below 100 MHz.     

Effective non-vertical and apparent cutoff rigidities for a cosmic ray latitude survey from Antarctica to Italy in minimum of solar activity

In this paper we will report the results of the computation of cutoff rigidities of vertical and non-vertical incident cosmic ray particles. Non-vertical effective cutoff rigidities have been computed by tracing particle trajectories through the “real” geomagnetic magnetic field comprising the International Geomagnetic Reference Field model (IGRF95, IAGA Division 5 Working Group 8, 1996: Sabaka, T.J., Langel, R.A., Baldwin, R.T., Conrad, J.A. The geomagnetic field, 1900–1995, including the large scale fields from magnetospheric sources and NASA candidate models for the 1995 IGRF revision. J. Geomag. Geoelect. 49, 157–206, 1997.) and the Tsyganenko [Tsyganenko, N.A. A magnetospheric magnetic field model with a warped tail current sheet. Planet. Space Sci. 37, 5–20, 1989.] magnetosphere model. The computation have been done for the backward route (from Antarctica to Italy) of the Italian Antarctic ship survey 1996–1997, for geographic points corresponding to the daily average coordinates of the ship; for zenith angles 15°, 30°, 45° and 60°, and azimuth angles from 0° to 360° in steps of 45°. By means of the obtained non-vertical cutoffs the apparent cutoff rigidities have been calculated. The information on integral multiplicities of secondary neutrons detected by the neutron monitor in dependence of the zenith angle of incoming primary cosmic ray particles have also been used. This information is based on the theoretical calculations of meson-nuclear cascades of primary protons with different rigidities arriving to the Earth’s atmosphere at the zenith angles of 0°, 15°, 30°, 45°, 60° and 75°. The difference between the computed apparent and vertical cutoff rigidities reaches ∼1 GV at rigidities >7–8 GV. At rigidities of 10–16 GV, the difference between the apparent and vertical cutoff rigidities is larger than that obtained earlier by Clem et al. [Clem, J.M., Bieber, J.W., Duldig, M., Evenson, P., Hall, D., Humble, J.E. Contribution of obliquely incident particles to neutron monitor counting rate. J. Geophys. Res. 102, 26919–26926, 1997.] and Dorman et al. [Dorman, L.I., Villoresi, G., Iucci, N., Parisi, M., Tyasto, M.I., Danilova, O.A., Ptitsyna, N.G. Cosmic ray survey to Antarctica and coupling functions for neutron component near solar minimum (1996–1997), 3. Geomagnetic effects and coupling functions. J. Geophys. Res. 105, 21047–21056, 2000.].     

Cosmic ray modulation with a Fisk-type heliospheric magnetic field and a latitude-dependent solar wind speed

The effect of a latitude-dependent solar wind speed on a Fisk heliospheric magnetic field [Fisk, L. A. Motion of the footpoints of heliospheric magnetic field lines at the Sun: implications for recurrent energetic particle events at high heliographic latitudes. J. Geophys. Res. 101, 15547–15553, 1996] was first discussed by Schwadron and Schwadron and McComas [Schwadron, N.A. An explanation for strongly underwound magnetic field in co-rotating rarefaction regions and its relationship to footpoint motion on the the sun. Geophys. Res. Lett. 29, 1–8, 2002. and Schwadron, N.A., McComas, D.J. Heliospheric “FALTS”: favored acceleration locations at the termination shock. Geophys. Res. Lett. 30, 41–1, 2003]. Burger and Sello [Burger, R.A., Sello, P.C. The effect on cosmic ray modulation of a Parker field modified by a latitudinal-dependent solar wind speed. Adv. Space Res. 35, 643–646, 2005] found a significant effect for a simplified 2D version of a latitude-dependent Fisk-type field while Miyake and Yanagita [Miyake, S., Yanagita, S. The effect of a modified Parker field on the modulation of the galactic cosmic rays. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, 445–448, 2007] found a smaller effect. The current report improves on a previous attempt Hitge and Burger [Hitge, M., Burger, R.A. The effect of a latitude-dependent solar wind speed on cosmic-ray modulation in a Fisk-type heliospheric magnetic field. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, pp. 449–450, 2007] where the global change in the solar wind speed and not the local speed gradient was emphasized. The sheared Fisk field of Schwadron and McComas [Schwadron, N.A., McComas, D.J. Heliospheric “FALTS”: Favored acceleration locations at the termination shock. Geophys. Res. Lett. 30, 41–1, 2003.) is similar to the current Schwadron–Parker hybrid field. Little difference is found between the effects of a Parker field and a Schwadron–Parker hybrid field on cosmic-ray modulation, in contrast to the results of Burger and Sello and Miyake and Yanagita [Burger, R.A., Sello, P.C. The effect on cosmic ray modulation of a Parker field modified by a latitudinal-dependent solar wind speed. Adv. Space Res. 35, 643–646, 2005 and Miyake, S., Yanagita, S. The effect of a modified Parker field on the modulation of the galactic cosmic rays. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, pp. 445–448, 2007]. The two-dimensional approximation used by these authors is therefore inadequate to model the complexities of the actual three-dimensional field. We also show that a Fisk-type field with a latitude-dependent solar wind speed (Schwadron–Parker hybrid field) decreases both the relative amplitude of recurrent cosmic ray intensity variations and latitude gradients and yields similar constants of proportionality for these quantities as for the constant solar wind speed case.     

Midday bottomside electron density profiles during moderate solar activity and comparison with IRI-2001

Bottomside electron density (Ne-h) profiles during midday (10–14 h) are analyzed using modern digital ionosonde observations at a low-middle latitude station, New Delhi (28.6N, 77.2E, dip 42.4N), for the period from January 2003 to December 2003, pertaining to moderate solar activity (MSA). Each individual profile is normalized with respect to the peak height and density (hmF2, NmF2) of the F2-region. These profiles are compared with those obtained from the International Reference Ionosphere (IRI-2001) model. Bilitza [Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 36 (2), 261–275, 2001] using both the options namely: Gulyaeva’s model [Gulyaeva, T.L. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7 (6) 39–48, 1987] and B0 Tab. option [Bilitza, D., Radicella, S.M., Reinisch, B.W., Adeniyi, J.O., Mosert Gonzalez, M.E., Zhang, S.R., Obrou, O. New B0 and B1 models for IRI. Adv. Space Res. 25 (1), 89–95, 2000]. The study reveals that during summer and equinox, the IRI model with B0 Tab. option in general, produces better agreement with the observed median profiles, while the IRI predictions using Gulyaeva’s option, overestimate the electron density distribution at all the heights below the F2-peak. However, during winter, in general, the IRI model, using both the options, reveals shows fairly good agreement with the observations.     

Energetic neutral atom (ENA) and charged particle periodicities in Saturn’s magnetosphere

The Magnetospheric Imaging Instrument (MIMI) on the Cassini spacecraft has observed energetic neutral atoms (ENA) and charged particles at Saturn from mid-2004 to the present. The particles often but not always reveal striking periodic behavior that seems to depend on the type of particle and spacecraft location. When subjected to a Lomb periodogram analysis, energetic electrons (>150 keV) exhibited strong frequency peaks near 10.80 h (the nominal or “base” period of Saturn kilometric radiation) during 2006–2008, but essentially no periodicity during 2005. The electron periodograms also show pronounced “double” frequency peaks in 2007 and 2008. Energetic protons (3–26 keV) show strong peaks near the same period for 2005–2007, but none for 2008. Oxygen ions at the same energies display strong peaks for 2005 and 2006, but not for 2007 and 2008. By projecting the ENA images onto Saturn’s equatorial plane or onto a plane perpendicular to the equatorial plane and then summing the data in the appropriate dimension, “strip” images can be constructed from which a time history can be derived. These time histories of ENA emissions are also subjected to a Lomb periodogram analyses. The energetic hydrogen neutrals (20–50 keV) exhibited periodic behavior only during 2007, while energetic oxygen neutrals (64–144 keV) displayed a strong SKR-like period in 2005 and 2006 but not for 2007 or 2008. Some of this behavior may be due to changing spacecraft aspect relative to the ENA emissions, and some of it may be real. This periodic behavior may be consistent with a rotating anomaly that “flashes” brightly in the midnight-to-dawn sector once per 10.8 h, with the flash parameters depending on particle species and energy.     

Gravity waves and wind-shear in the MLT at 23°S

We have used the technique suggested by Hocking [Hocking, W. A new approach to momentum flux determinations using SKiYMET meteor radars. Ann. Geophys. 23, 2005.] to derive short period wind variances in the 80–100 km region from meteor radar data. We find that these fluctuating winds, assumed to correspond to gravity waves and turbulence, are closely correlated with the vertical shear of the horizontal tidal winds. This close correlation suggests that in situ wind shear may be a major source of gravity waves and turbulence in the MLT. If this is the case, gravity waves generated in the troposphere and propagating up to the MLT region, generally assumed to constitute an important influence on the climatology of the region, may be a less important source of energy and momentum in the 80–100 km region than has been hitherto believed.     

Variations of the solar acoustic high-degree mode frequencies over solar cycle 23

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present variations of the solar acoustic mode frequencies caused by the solar activity cycle. High-degree (100 < ? < 900) solar acoustic modes were analyzed using global helioseismology analysis techniques over most of solar cycle 23. We followed the methodology described in details in [Korzennik, S.G., Rabello-Soares, M.C., Schou, J. On the determination of Michelson Doppler Imager high-degree mode frequencies. ApJ 602, 481–515, 2004] to infer unbiased estimates of high-degree mode parameters ([see also Rabello-Soares, M.C., Korzennik, S.G., Schou, J. High-degree mode frequencies: changes with solar cycle. ESA SP-624, 2006]). We have removed most of the known instrumental and observational effects that affect specifically high-degree modes. We show that the high-degree changes are in good agreement with the medium-degree results, except for years when the instrument was highly defocused. We analyzed and discuss the effect of defocusing on high-degree estimation. Our results for high-degree modes confirm that the frequency shift scaled by the relative mode inertia is a function of frequency and it is independent of degree.     

Electric field and ion density anomalies in the mid latitude ionosphere: Possible connection with earthquakes?

The problem of earthquake prediction has stimulated the search for a correlation between seismic activity and ionospherical anomalies. We found observational evidence of possible earthquake effects in the near-equatorial and low latitude ionosphere; these ionospheric anomalies have been proposed by Gousheva et al. [Gousheva, M., Glavcheva, R., Danov, D., Angelov P., Hristov, P., Influence of earthquakes on the electric field disturbances in the ionosphere on board of the Intercosmos-Bulgaria-1300 satellite. Compt. Rend. Acad. Bulg. Sci. 58 (8) 911–916, 2005a; Gousheva, M., Glavcheva, R., Danov, D., Angelov, P., Hristov, P., Kirov, B., Georgieva, K., Observation from the Intercosmos-Bulgaria-1300 satellite of anomalies associated with seismic activity. In: Poster Proceeding of 2nd International Conference on Recent Advances in Space Technologies: Space in the Service of Society, RAST ‘2005, June 9–11, Istanbul, Turkey, pp. 119–123, 2005b; Gousheva, M., Glavcheva, R., Danov, D., Angelov, P., Hristov, P., Kirov, B., Georgieva, K., Satellite monitoring of anomalous effects in the ionosphere probably related to strong earthquakes. Adv. Space Res. 37 (4), 660–665, 2006]. This paper presents new results from observations of the quasi-static electric field and ion density on board INTERCOSMOS-BULGARIA-1300 satellite in the mid latitude ionosphere above sources of moderate earthquakes. Data from INTERCOSMOS-BULGARIA-1300 satellite and seismic data (World Data Center, Denver, Colorado, USA) for magnetically quiet and medium quiet days are juxtaposed in time-space domain. For satellite’s orbits in the time period 15.09–01.10.1981 an increase in the horizontal and vertical components of the quasi-static electric field and fluctuations of the ion density are observed over zones of forthcoming seismic events. Some similar post effects are observed too. The emphasis of this paper is put on the anomalies which specify the mid latitude ionosphere. The obtained results contain important information because they confirm our previous results for near-equatorial and low latitude regions.     

长时间流星不均匀体母体特征个例研究

流星体坠入地球大气烧蚀电离产生流星等离子体尾迹,在等离子体不稳定性过程作用下产生流星不均匀体.利用光学视频和无线电雷达在低纬三亚开展流星体烧蚀和流星不均匀体综合探测结果,发展了一种获取流星不均匀体母体（流星体）特征参数的方法,并对2015年12月双子座流星雨期间观测的一次长持续时间流星不均匀体事例进行了分析,得到了其母体速度、质量和轨道等特征,结果显示产生这次流星不均匀体的流星体速度和轨道等具有双子座流星特点.该方法可应用于流星不均匀体及其母体特征研究.      

廊坊中频雷达流星观测及初步结果

中频雷达用来开展夜间100km高度以上的流星观测,获得流星随时间、高度、方位的分布情况及流星体速度、流星辐射点、流星余迹径向速度等参数,其探测数据可用于流星天文学、中层大气动力学等领域的研究.利用2017年11月16日12:00UT-22:00UT期间廊坊观测站（39.4°N,116.7°E）的中频雷达数据,首次开展了中国中纬度地区夜间流星观测实验,共检测到94个流星回波信号,集中分布在97～115km高度范围内,平均高度为106.5km,计算得到了流星回波的双极扩散系数、方位分布等相关参数,并与国外中频雷达流星探测结果进行了初步比较.      

Global dynamics of high area-to-mass ratios GEO space debris by means of the MEGNO indicator

In this paper we provide an extensive analysis of the global dynamics of high-area-to-mass ratios geosynchronous (GEO) space debris, applying a recent technique developed by Cincotta and Simó [Cincotta, P.M., Simó, C.Simple tools to study global dynamics in non-axisymmetric galactic potentials–I. Astron. Astrophys. (147), 205–228, 2000.], Mean Exponential Growth factor of Nearby Orbits (MEGNO), which provides an efficient tool to investigate both regular and chaotic components of the phase space.     

Asteroid and comet hazard: Identification problem of observed space objects with the parental bodies

This article focuses on the genetic identification of observed small cosmic bodies with alleged parental bodies; namely, comets, asteroids and meteoroid swarms. There is a problem of the upper D-value limit as a measure of proximity between the orbits of the bodies in the five-dimensional phase space (Southworth and Hawkins, 1963). In the study of genetic relationships of the comet and meteor complexes, the D value is usually taken as equal to 0.2 for all meteor showers. However, the upper D limit should be investigated for each meteoroid complex. For example, such investigation was performed for the Taurid meteor complex (Porub?an et al., 2006). In this paper, the upper D-criterion limit value was investigated for the Perseid meteor shower. The 1862 III Swift–Tuttle comet is its parental comet.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

Amateur observations of solar eclipses and derivation of scientific data

This work presents the educational approach of using total solar eclipse occurrences as a scientific process learning aid. The work reviews the basic scientific aims and experiments included in the observational programs “Total solar eclipse 1999 and 2006” (Stoev, A., Kiskinova, N., Muglova, P. et al. Complex observational programme of the Yuri Gagarin Public Astronomical Observatory and STIL, BAS, Stara Zagora Department for the August 11, 1999 total solar eclipse, in: Total Solar Eclipse 1999 – Observational Programmes and Coordination, Proceedings, Recol, Haskovo, pp. 133–137, 1999a (in Bulgarian); Stoeva, P.V., Stoev, A.D., Kostadinov, I.N. et al. Solar Corona and Atmospheric Effects during the March 29, 2006 Total Solar Eclipse, in: 11th International Science Conference SOLAR–Terrestrial Influences, Sofia, November 24–25, pp. 69–72, 2005).     

The variability of the proton cosmic ray flux on the Sun’s way around the galactic center

The discovery of direct evidences for the acceleration of high energetic particles at the shell supernova remnant RXJ1713.7-3946 underlined the need to calculate the cosmic ray (CR) distribution in the Galaxy on a spatial grid fine enough to resolve the changes in the CR density due to these kind of objects. It was shown before by Büsching et al. [Büsching, I., Kopp, A., Pohl, M., Schlickeiser, R., Perrot, C., Grenier, I. Cosmic-ray propagation properties for an origin in supernova remnants. ApJ 619, 314–326, 2005] that the discrete nature (both in space and time) of super novae (SN) as sources of Galactic CR leads to CR spectra changing in space and time, resulting in a range of possible CR spectra at a given location in the Galaxy. As the most frequent SN types Ib and II are found within spiral arms, one can expect a significant difference of the range of possible spectra in and outside spiral arms. We investigate the variability of the local interstellar CR proton spectrum during the motion of the Sun in and out of spiral arms in its journey around the Galactic center. Using the code described by Büsching et al. [Büsching, I., Kopp, A., Pohl, M., Schlickeiser, R., Perrot, C., Grenier, I. Cosmic-ray propagation properties for an origin in supernova remnants. ApJ 619, 314–326, 2005], the proton CR density in the Galaxy is calculated with high spatial and temporal resolution (75 pc in galactocentric radius r and azimuth φ at the position of the Sun and 20 pc in z (perpendicular to the galactic plane), with a time step of 1 kyr), assuming stochastically distributed point sources with a probability distribution that resembles the spiral structure of our Galaxy. We find the averaged CR flux outside the spiral arm to be about 50% of that inside a spiral arm. We further find spatial and temporal variations of the CR flux inside spiral arms.     

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”.