GPS-TEC variations during low solar activity period (2007–2009) at Indian low latitude stations

The paper is based on the ionospheric variations in terms of vertical total electron content (VTEC) for the low solar activity period from May 2007 to April 2009 based on the analysis of dual frequency signals from the Global Positioning System (GPS) satellites recorded at ground stations Varanasi (Geographic latitude 25°16′ N, Longitude 82°59′ E), situated near the equatorial ionization anomaly crest and other two International GNSS Service (IGS) stations Hyderabad (Geographic latitude 17°20′ N, longitude 78°30′ E) and Bangalore (Geographic latitude 12°58′ N, longitude 77°33′ E) in India. We describe the diurnal and seasonal variations of total electron content (TEC), and the effects of a space weather related event i.e. a geomagnetic storm on TEC. The mean diurnal variation during different seasons is brought out. It is found that TEC at all the three stations is maximum during equinoctial months (March, April, September and October), and minimum during the winter months (November, December, January and February), while obtaining intermediate values during summer months (May, June, July and August). TEC shows a semi-annual variation. TEC variation during geomagnetic quiet as well as disturbed days of each month and hence for each season from May 2007 to April 2008 at Varanasi is examined and is found to be more during disturbed period compared to that in the quiet period. Monthly, seasonal and annual variability of GPS-TEC has been compared with those derived from International Reference Ionosphere (IRI)-2007 with three different options of topside electron density, NeQuick, IRI01-corr and IRI 2001. A good agreement is found between the GPS-TEC and IRI model TEC at all the three stations.     

Effect of magnetic activity on ionospheric time delay at low latitude

The purpose of this work is to investigate the effect of magnetic activity on ionospheric time delay at low latitude Station Bhopal (geom. lat. 23.2°N, geom. long. 77.6°E) using dual frequency (1575.42 and 1227.60 MHz) GPS measurements. Data from GSV4004A GPS Ionospheric Scintillation and TEC monitor (GISTM) have been chosen to study these effects. This paper presents the results of ionospheric time delay during quiet and disturbed days for the year 2005. Results show that maximum delay is observed during quiet days in equinoxial month while the delays of disturbed period are observed during the months of winter. We also study the ionospheric time delay during magnetic storm conditions for the same period. Results do not show any clear relationship either with the magnitude of the geomagnetic storm or with the main phase onset (MPO) of the storm. But most of the maximum ionospheric time delay variations are observed before the main phase onset (MPO) or sudden storm commencement (SSC) as compared to storm days.     

Longitude dependent response of the GPS derived ionospheric ROTI to geomagnetic storms

The local time dependent effects of geomagnetic storm on the ionospheric TEC and Rate of change of TEC Index (ROTI) are studied here using the GPS data for four different low latitude stations: Ogaswara, Japan (24.29?°N, 153.91?°E; Geomagnetic: 17.21?°N, 136.16?°W); Surat, India (21.16?°N, 72.78?°E; Geomagnetic: 12.88?°N, 146.91?°E); Bogota, Colombia (4.64?°N, ?74.09?°E; Geomagnetic: 14.42?°N, 1.67?°W); and Kokee park Waimea, Hawaii, US (22.12?°N, ?159.67?°E; Geomagnetic: 22.13?°N, 91.19?°W). The solar wind velocity and geomagnetic indices: Dst, Kp and IMF Bz are utilized to validate the geomagnetic storms registered during the years 2011 and 2012. Using the GPS based TEC data and computed values of ROTI, the storm induced ionospheric irregularities generation and inhibition has been studied for all stations. The present study suggests that, the F-region irregularities of a scale length of few kilometers over the magnetic equator are locally affected by geomagnetic storms. This study also shows a good agreement (70–84 %) with the Aaron’s criteria (Aarons, Radio Sci., 26:1131–1149, 1991; Biktash, Ann. Geophys., 19:731–739, 2004) as significant absence and enhancement of ROTI was found to be influenced by the local time of the negative peak of Dst index association.     

Study of total electron content variations over equatorial and low latitude ionosphere during extreme solar minimum

In recent years with the advancement in satellite based navigational applications, study of Total Electron Content (TEC) has gained significant importance. It is well known that due to dynamical behaviour of equatorial and low latitude ionosphere, the levels of ionization is relatively high herein. The sustained decrease in solar extreme ultraviolet radiations during the current minimum is greater than any in recent history. This gives us the opportunity to study the observations of global positioning system total electron content (GPS-TEC) dual frequency signals from the GPS satellites continuously recorded at Trivandrum (an equatorial station) and Delhi (a low latitude station) during the extremely low solar activity period from January 2007 to June 2009. This study illustrates the diurnal, seasonal and annual variations of TEC during the extended solar minimum period. This study also investigates the behaviour of daytime ionosphere around spring and autumn equinoxes at low solar activity period. The results clearly reveal the presence of equinoctial asymmetry which is more pronounced at equatorial station Trivandrum. The diurnal variation of TEC shows a short-lived day minimum which occurs between 0500 to 0600 LT at both the stations. Delhi TEC values show its steep increase and reach at its peak value between 1200 and 1400 LT, while at the equator the peak is broad and occurs around 1600 LT. Further, the daily maximum TEC ranges from about 5 to 40 TEC units at Trivandrum and about 10 to 40 TEC units at Delhi, which correspond to range delay variations of about 1 to 8 m at the GPS L1 frequency of 1.575 GHz. The Maximum values of TEC were observed during spring equinox rather than autumn equinox, showing presence of semi annual variation at both the locations. The minimum values of TEC were observed during the summer solstice at Trivandrum indicating the presence of winter anomaly at equatorial region while Delhi TEC values were minimum during winter solstice showing absence of winter anomaly. Also the TEC values at both the locations have been decreasing since 2007 onwards exhibit good positive correlation with solar activity.     

On accounting for temperature and pressure in determining the Rayleigh scattering in the earth’s atmosphere

We calculated the variations of Rayleigh optical depth with changes of pressure and temperature for three observation sites: Simferopol (φ = 44°57′N, λ = 34°8′E, h = 265 m above sea level), Nauchny (φ = 44°43′N, λ = 34°3′E, h = 583 m), and Ai-Petry meteorological station (φ = 44°24′N, λ = 34°6′E, h = 1180 m).     

The phase ratios of the color index: Mapping of two regions of the near side of the Moon

From the ground-based colorimetry performed for two surface regions of the near side of the Moon, images of the phase ratio of the color index C(600 nm/470 nm) have been built for the phase angles between 2° and 95°. It has been found that for phase angles smaller than α ∼ 40°–50°, the color index of the highlands grows with the phase quicker than that of the mare regions. For larger phase angles, α > 50°, a reverse situation is observed. The laboratory data on the spectrophotometry of the lunar samples confirm the peculiarities found in the phase dependence of color. The influence of multiple scattering on the phase dependence of the color of the mare and highland regions of the Moon are discussed.     

Latitude Dependence of the Variations of Sunspot Group Numbers (SGN) and Coronal Mass Ejections (CMEs) in Cycle 23

The 12-month running means of the conventional sunspot number Rz, the sunspot group numbers (SGN) and the frequency of occurrence of Coronal Mass Ejections (CMEs) were examined for cycle 23 (1996 – 2006). For the whole disc, the SGN and Rz plots were almost identical. Hence, SGN could be used as a proxy for Rz, for which latitude data are not available. SGN values were used for 5° latitude belts 0° – 5°, 5° – 10°, 10° – 15°, 15° – 20°, 20° – 25°, 25° – 30° and > 30°, separately in each hemisphere north and south. Roughly, from latitudes 25° – 30° N to 20° – 25° N, the peaks seem to have occurred later for lower latitudes, from latitudes 20° – 25° N to 15° – 20° N, the peaks are stagnant or occur slightly earlier, and then from latitudes 15° – 20° N to 0° – 5° N, the peaks seem to have occurred again later for lower latitudes. Thus, some latitudinal migration is suggested, clearly in the northern hemisphere, not very clearly in the southern hemisphere, first to the equator in 1998, stagnant or slightly poleward in 1999, and then to the equator again from 2000 onwards, the latter reminiscent of the Maunder butterfly diagrams. Similar plots for CME occurrence frequency also showed multiple peaks (two or three) in almost all latitude belts, but the peaks were almost simultaneous at all latitudes, indicating no latitudinal migration. For similar latitude belts, SGN and CME plots were dissimilar in almost all latitude belts except 10° – 20° S. The CME plots had in general more peaks than the SGN plots, and the peaks of SGN often did not match with those of CME. In the CME data, it was noticed that whereas the values declined from 2002 to 2003, there was no further decline during 2003 – 2006 as one would have expected to occur during the declining phase of sunspots, where 2007 is almost a year of sunspot minimum. An inquiry at GSFC-NASA revealed that the person who creates the preliminary list was changed in 2004 and the new person picks out more weak CMEs. Thus a subjectivity (overestimates after 2002) seems to be involved and hence, values obtained before and during 2002 are not directly comparable to values recorded after 2002, except for CMEs with widths exceeding 60°.     

A Study on the Reliability of the Ionospheric VTEC and Satellite DCB Derived from a Regional GPS Network

According to the requirement of high-precision satellite navigation, we have introduced the method for the quasi-realtime monitoring of variations of the regional ionospheric total electron content (TEC) and GPS satellite differential code bias (DCB), based on the dual-frequency carrier-phase smoothed pseudorange data obtained from a regional GPS network in China. Especially, we have studied the feasibility of retrieving DCB independently from the regional GPS networks with different sizes. For this purpose, 3 regional networks based on the countrywide GPS stations are investigated. The comparisons of the computed DCB and VTEC(vertical total electron content) with those of CODE(the Center for Orbit Determination in Europe) indicate that in order to realize a reliable quasi-realtime measurement of DCB by a regional network, there is a certain requirement on the size of the regional network, and that the relative accuracies of the quasi-realtime VTEC and DCB measured by using a Chinese GPS regional network can reach 2.0TECu and 0.25 ns, respectively.     

A synthesis map of the sky at 34.5 MHz

This paper describes a wide-field survey made at 34.5 MHz using GEETEE,1 the low frequency telescope at Gauribidanur (latitude 13°36′12′′N). This telescope was used in the transit mode and by per forming 1-D synthesis along the north-south direction the entire observable sky was mapped in a single day. This minimized the problems that hinder wide-field low-frequency mapping. This survey covers the declination range of-50° to + 70° (- 33° to +61° without aliasing) and the complete 24 hours of right ascension. The synthesized beam has a resolution of 26′ x 42′ sec (δ- 14°. 1). The sensitivity of the survey is 5 Jy/beam (1σ). Special care has been taken to ensure that the antenna responds to all angular scale structures and is suitable for studies of both point sources and extended objects This telescope is jointly operated by the Indian Institute of Astrophysics, Bangalore and the Roman Research Institute, Bangalore.     

Response of low latitude ionosphere to the geomagnetic storm of May 30

Response of low latitude ionosphere to the geomagnetic storm of May 30, 2005 in the Indian longitude sector has been investigated by using the GPS data recorded at three stations namely, Udaipur, Hyderabad and Bengaluru. The event is noteworthy due to the fact that the Z component of interplanetary magnetic field (IMF Bz) remained southward for about 10 hours, coincident with the local day time for the Indian longitude sector, along with significantly higher values of AE and ASY-H indices. However, we neither found any evidence for the presence of long lasting storm time electric fields nor could we infer episodes of eastward-westward penetration of electric fields under steady southward IMF Bz and unsteady ring current conditions. On the storm day, the maximum enhancement in the total electron content has been found to be about 60%. The ionosonde observations also showed increased critical frequency (foF₂) and the height (h_PF₂) of the F layer. The foF₂ was enhanced by ∼60% which is consistent with the enhancement in total electron content. The slow rise and long duration enhancement of h_PF₂ and foF₂ have been attributed to the upwelling by the meridional neutral winds, caused by continuous energy inputs at higher latitudes. The poleward expansion of the equatorial ionization anomaly has also been observed on May 30. On May 31, the following day of the storm, significantly suppressed anomaly with near absence of its northern crest in the Indian longitude sector, revealed the effect of storm induced disturbance dynamo electric fields.     

NGC 4650 A: A nearly edge-on ring galaxy?

The peculiar galaxy NGC 4650 A (α=12^h 42^m. 1; = δ—40° 26′; 1950·0) has been studied by means of direct and spectral observations with the ESO 3·6-m telescope. It is interpreted as a prolate, elliptical galaxy surrounded by a warped ring of H II regions, dust and stars. The distance is 47 Mpc (H ₀=55 km s⁻¹ Mpc⁻¹). The ring is seen nearly edge-on (inclination 85°) and it rotates. It has a diameter of about 21 kpc and is bluer than the elliptical galaxy for which the (M/L _v) ratio is ∼12 in solar units. The observed configuration may be the result of interaction with the nearby galaxy, NGC 4650.     

Mathematical modelling of ionospheric TEC from Turkish permanent GNSS Network (TPGN) observables during 2009–2017 and predictability of NeQuick and Kriging models

The present study examines the ionospheric Total Electron Content (TEC) variations in the lower mid-latitude Turkish region from the Turkish permanent GNSS network (TPGN) and International GNSS Services (IGS) observations during the years 2009 to 2017. The corresponding vertical TEC (VTEC) predicted by Kriging and NeQuick-2 models are evaluated to realize their efficacy over the country. We studied the diurnal, seasonal and spatial pattern of VTEC variation and tried to estimate by a new mathematical model using the long term of 9 years VTEC data. The diurnal variation of VTEC demonstrates a normal trend with its gradual enhancement from dawn to attain a peak around 09:00–14.00 UT and reaching the minimum level after 22.00 UT. The seasonal behavior of VTEC indicates a strong semi-annual variation of VTEC with maxima in September equinox followed by March equinox and minima in June solstice followed by December solstice. Also, the spatial variation in VTEC depicts a meaningful longitudinal/latitudinal pattern altering with seasons. It decreases longitudinally from the west to the east during March equinox and June solstice increases with latitude. The comparative analysis among the GNSS-VTEC, Kriging, NeQuick and the proposed mathematical model are evaluated with the help one way ANOVA test. The analysis shows that the null hypothesis of the models during storm and quiet days are accepted and suggesting that all models are statistically significantly equivalent from each other. We believe the outcomes from this study would complement towards a relatively better understanding of the lower mid-latitude VTEC variation over the Turkish region and analogous latitudes over the globe.     

Ionospheric total electron content seismo-perturbation after Japan’s March 11, 2011, M=9.0 Tohoku earthquake under a geomagnetic storm; a nonlinear principal component analysis

Nonlinear principal component analysis (NLPCA) is implemented to analyze the spatial pattern of total electron content (TEC) anomalies 3 hours after Japan’s Tohoku earthquake that occurred at 05:46:23 on 11 March, 2011 (UTC) (M _w =9). A geomagnetic storm was in progress at the time of the earthquake. NLPCA and TEC data processing were conducted on the global ionospheric map (GIM) for the time between 08:30 to 09:30 UTC, about 3 hours after this devastating earthquake and ensuing tsunami. Analysis results show stark earthquake-associated TEC anomalies that are widespread, and appear to have been induced by two acoustic gravity waves due to strong shaking (vertical acoustic wave) and the generation of the tsunami (horizontal Rayleigh mode gravity wave). The TEC anomalies roughly fit the initial mainshock and movement of the tsunami. Observation of the earthquake-associated TEC anomalies does not appear to be affected by a contemporaneous geomagnetic storm.     

Slow and Fast Rotating Coronal Holes from Geomagnetic Indices

The evolution of the 27-day recurrence in the series of two solar indices (Wolf number WN and 10.7 cm radio flux F) and two geomagnetic indices (Dst and ζ, variance of the geomagnetic field recorded at a magnetic observatory) have been studied over the 1957 – 2007 time span. Spectral energies contained in two period domains (25 – 27.3 and 27.3 – 31 days), designated as E ₁ and E ₂, have been computed. Whereas the evolution of E ₁ is the same for the four indices, that of E ₂ is essentially different for WN and F on the one hand, Dst and ζ on the other hand. Some general conclusions on the dynamics of the solar outer layers are inferred from these results. First the solar activity, as measured by WN, and when averaged over a few years, evolves in the same way whatever the latitude. Second, two families of coronal holes (CHs) are identified; the rapidly and the slowly rotating CHs evolve quite differently.     

A general catalogue of ephemerides and apparent orbits for visual binaries

Summary Informations on 736 pairs of visual binaries are given in the form of the Tables. The General Catalogue gives ephemerides (t,θ°,ρ″) in 20 years (1984–2003) for each pair with the figures of its apparent elliptical orbit where the positions of secondary component relative to the primary one at different epochs are indicated.The General Catalogue contains four parts: Part one — Source and Grading of Orbit; Part two — Ephemérides and Atlas of Apparent Orbits; Part three — Classifications of Visual Binary Stars; Part four — Indexes of Visual Binary Stars. The principle of calculation and the statistical data are presented in this paper. There are seven statistical tables, giving the elements distribution of true and apparent orbits, grading distribution of orbits, number distribution with different physical property of component of binary star. The number of binary stars in anyone constellation, the number of binary stars brighter than 6^m.5. The number of binary stars nearer than 25 parsec.     

Automated search for star clusters in large multiband surveys: II. Discovery and investigation of open clusters in the galactic plane

Automated search for star clusters in J, H, K _s data from 2MASS catalog has been performed using the method developed by Koposov et al. (2008). We have found and verified 153 new clusters in the interval of the galactic latitude −24° < b < 24°. Color excesses E(B − V), distance moduli and ages were determined for 130 new and 14 yet-unstudied known clusters. In this paper, we publish a catalog of coordinates, diameters, and main parameters of all the clusters under study. A special web-site available at has been developed to facilitate dissemination and scientific usage of the results.     

Accuracy and Limitations of Fitting and Stereoscopic Methods to Determine the Direction of Coronal Mass Ejections from Heliospheric Imagers Observations

Using data from the Heliospheric Imagers (HIs) onboard STEREO, it is possible to derive the direction of propagation of coronal mass ejections (CMEs) in addition to their speed with a variety of methods. For CMEs observed by both STEREO spacecraft, it is possible to derive their direction using simultaneous observations from the twin spacecraft and also, using observations from only one spacecraft with fitting methods. This makes it possible to test and compare different analysis techniques. In this article, we propose a new fitting method based on observations from one spacecraft, which we compare to the commonly used fitting method of Sheeley et al. (J. Geophys. Res. 104, 24739, 1999). We also compare the results from these two fitting methods with those from two stereoscopic methods, focusing on 12 CMEs observed simultaneously by the two STEREO spacecraft in 2008 and 2009. We find evidence that the fitting method of Sheeley et al. (J. Geophys. Res. 104, 24739, 1999) may result in significant errors in the determination of the CME direction when the CME propagates outside of 60°±20° from the Sun – spacecraft line. We expect our new fitting method to be better adapted to the analysis of halo or limb CMEs with respect to the observing spacecraft. We also find some evidence that direct triangulation in the HI fields-of-view should only be applied to CMEs propagating approximatively toward Earth (± 20° from the Sun – Earth line). Last, we address one of the possible sources of errors of fitting methods: the assumption of radial propagation. Using stereoscopic methods, we find that at least seven of the 12 studied CMEs had a heliospheric deflection of less than 20° as they propagated in the HI fields-of-view, which, we believe, validates this approximation.     

Predicting the Amplitude of a Solar Cycle Using the North – South Asymmetry in the Previous Cycle: II. An Improved Prediction for Solar Cycle 24

Recently, using Greenwich and Solar Optical Observing Network sunspot group data during the period 1874 – 2006, Javaraiah (Mon. Not. Roy. Astron. Soc. 377, L34, 2007: Paper I), has found that: (1) the sum of the areas of the sunspot groups in 0° – 10° latitude interval of the Sun’s northern hemisphere and in the time-interval of −1.35 year to +2.15 year from the time of the preceding minimum of a solar cycle n correlates well (corr. coeff. r=0.947) with the amplitude (maximum of the smoothed monthly sunspot number) of the next cycle n+1. (2) The sum of the areas of the spot groups in 0° – 10° latitude interval of the southern hemisphere and in the time-interval of 1.0 year to 1.75 year just after the time of the maximum of the cycle n correlates very well (r=0.966) with the amplitude of cycle n+1. Using these relations, (1) and (2), the values 112±13 and 74±10, respectively, were predicted in Paper I for the amplitude of the upcoming cycle 24. Here we found that the north – south asymmetries in the aforementioned area sums have a strong ≈44-year periodicity and from this we can infer that the upcoming cycle 24 will be weaker than cycle 23. In case of (1), the north – south asymmetry in the area sum of a cycle n also has a relationship, say (3), with the amplitude of cycle n+1, which is similar to (1) but more statistically significant (r=0.968) like (2). By using (3) it is possible to predict the amplitude of a cycle with a better accuracy by about 13 years in advance, and we get 103±10 for the amplitude of the upcoming cycle 24. However, we found a similar but a more statistically significant (r=0.983) relationship, say (4), by using the sum of the area sum used in (2) and the north – south difference used in (3). By using (4) it is possible to predict the amplitude of a cycle by about 9 years in advance with a high accuracy and we get 87±7 for the amplitude of cycle 24, which is about 28% less than the amplitude of cycle 23. Our results also indicate that cycle 25 will be stronger than cycle 24. The variations in the mean meridional motions of the spot groups during odd and even numbered cycles suggest that the solar meridional flows may transport magnetic flux across the solar equator and potentially responsible for all the above relationships. The author did a major part of this work at the Department of Physics and Astronomy, UCLA, 430 Portola Plaza, Los Angeles, CA 90095-1547, USA.     

Resonances of low orders in the planetary system of HD37124

The full set of published radial velocity data (52 measurements from Keck + 58 ones from ELODIE + 17 ones from CORALIE) for the star HD37124 is analysed. Two families of dynamically stable high-eccentricity orbital solutions for the planetary system are found. In the first one, the outer planets c and d are trapped in the 2/1 mean-motion resonance. The second family of solutions corresponds to the 5/2 mean-motion resonance between these planets. In both families, the planets are locked in (or close to) an apsidal corotation resonance. In the case of the 2/1 MMR, it is an asymmetric apsidal corotation (with the difference between the longitudes of periastra Δω ~ 60°), whereas in the case of the 5/2 MMR it is a symmetric antialigned one (Δω = 180°). It remains also possible that the two outer planets are not trapped in an orbital resonance. Then their orbital eccentricities should be relatively small (less than, say, 0.15) and the ratio of their orbital periods is unlikely to exceed 2.3 − 2.5.     

Large-scale mapping of the IRAS 0042 + 5530 region in the 12CO (J = 1−0) and 13CO (J = 1−0) molecular lines

We present new radio observations of molecular lines in the region of high mass star formation, namely G122.0-7.1. A large-scale map of the emission observed in the ¹²CO (J = 1−0) and ¹³CO (J = 1−0) lines covers the area of 15′ × 9′, revealing two dense regions. The molecular bipolar outflows have been resolved in ASO1 region. It is associated with the known candidate YSO nearby IRAS 0042 + 5530. Also, a new dense region has been discovered in the North-Western part of the G122.0-7.1 at a distance of 5′ from IRAS 0042 + 5530. Its position is close to the peak of 4850 MHz emission. The text was submitted by the authors in English.