Comparison of observed ionospheric vertical TEC over the sea in Indian region with IRI-2016 model

The vertical ionospheric TEC values obtained from GAGAN grid based ionospheric delay correction values over the sea in the Indian equatorial region have been compared with the corresponding values derived from the International Reference Ionosphere model, IRI-2016. The objective of this work is to study the deviation of the vertical TEC derived from the IRI model from ground truths over the sea for different conditions. This will serve the basic intention of assessing the candidature of the IRI model as an alternative ionospheric correction model in navigation receivers in terms of accuracy. We have chosen different solar activity periods, seasons, geomagnetic conditions, locations etc. for our comparison and analysis. The TEC values by the IRI-2016 were compared with the actual measured values for the given conditions and errors were obtained. The measured vertical TEC values at the ionospheric grid points were derived from the GAGAN broadcast ionospheric delay data and used as reference. The IRI model with standard internal functions was used in estimating the TEC at the same ionospheric grid points. The errors in the model derived values are statistically analysed. Broadly, the results show that, for the Indian sector over the sea, the IRI model performs better on quiet days in off equatorial regions, particularly in the northern region. The overall performance degrades for other conditions with the model generally underestimating the true TEC values and most severely in the equatorial region. The performance is worst in this region for the disturbed days of the equinoctial period. The comparison study is also done with the TEC data measured directly by dual frequency GPS receivers. The results were found to be in general agreement with those obtained by comparing the model with GAGAN broadcast data as reference. This study will be useful in considering the IRI-2016 model for real time estimates of TEC as an alternative to the current parametric model in a satellite navigation receiver in absence of other options.     

Ionospheric total electron content response to annular solar eclipse on June 21, 2020

The effects of physical events on the ionosphere structure is an important field of study, especially for navigation and radio communication. The paper presents the spatio-temporal ionospheric TEC response to the recent annular solar eclipse on June 21, 2020, which spans across two continents, Africa and Asia, and 14 countries. This eclipse took place on the same day as the June Solstice. The Global Navigation Satellite System (GNSS) based TEC data of the Global Ionosphere Maps (GIMs), 9 International GNSS Service (IGS) stations and FORMOSAT-7/COSMIC-2 (F7/C2) were utilized to analyze TEC response during the eclipse. The phases of the TEC time series were determined by taking the difference of the observed TEC values on eclipse day from the previous 5-day median TEC values. The results showed clear depletions in the TEC time series on June 21. These decreases were between 1 and 9 TECU (15–60%) depending on the location of IGS stations. The depletions are relatively higher at the stations close to the path of annular eclipse than those farther away. Furthermore, a reduction of about ?10 TECU in the form of an equatorial plasma bubble (EPB) was observed in GIMs at ~20° away from the equator towards northpole, between 08:00–11:00 UT where its maximum phase is located in southeast Japan. Additionally, an overall depletion of ~10% was observed in F7/C2 derived TEC at an altitude of 240 km (hmF2) in all regions affected by the solar eclipse, whereas, significant TEC fluctuations between the altitudes of 100 km ? 140 km were analyzed using the Savitzky-Golay smoothing filter. To prove TEC depletions are not caused by space weather, the variation of the sunspot number (SSN), solar wind (V_SW), disturbance storm-time (Dst), and Kp indices were investigated from 16th to 22nd June. The quiet space weather before and during the solar eclipse proved that the observed depletions in the TEC time series and profiles were caused by the annular solar eclipse.     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

Scaling analysis and model estimation of solar corona index

A monthly average solar green coronal index time series for the period from January 1939 to December 2008 collected from NOAA (The National Oceanic and Atmospheric Administration) has been analysed in this paper in perspective of scaling analysis and modelling. Smoothed and de-noising have been done using suitable mother wavelet as a pre-requisite. The Finite Variance Scaling Method (FVSM), Higuchi method, rescaled range (R/S) and a generalized method have been applied to calculate the scaling exponents and fractal dimensions of the time series. Autocorrelation function (ACF) is used to find autoregressive (AR) process and Partial autocorrelation function (PACF) has been used to get the order of AR model. Finally a best fit model has been proposed using Yule-Walker Method with supporting results of goodness of fit and wavelet spectrum. The results reveal an anti-persistent, Short Range Dependent (SRD), self-similar property with signatures of non-causality, non-stationarity and nonlinearity in the data series. The model shows the best fit to the data under observation.     

Application of a modified generalized sail model on a solar sail with wrinkles

In this paper, we present an analysis of effect of wrinkles on the solar sail performance. We describe different analytical, semi-analytical and numerical approaches to the calculation of general large-scale curvature of a solar sail as well as parameters of so-called wrinkled domains, and introduce the impact of such wrinkles on the thrust and torque of the solar sail. Finally, we present a model of an optically-orthotropic surface for such non-ideal sail, providing a connection with the Generalized Sail Model, and other solar sail thrust models.     

Non-stationarity and cross-correlation effects in the MHD solar activity

The analysis of turbulent processes in sunspots and pores which are self-organizing long-lived magnetic structures is a complicated and not yet solved problem. The present work focuses on studying such magneto-hydrodynamic (MHD) formations on the basis of flicker-noise spectroscopy using a new method of multi-parametric analysis. The non-stationarity and cross-correlation effects taking place in solar activity dynamics are considered. The calculated maximum values of non-stationarity factor may become precursors of significant restructuring in solar magnetic activity. The introduced cross-correlation functions enable us to judge synchronization effects between the signals of various solar activity indicators registered simultaneously.     

Investigating the Coronal Mass Ejections associated with DH type-II radio bursts and solar flares during the ascending phase of the solar cycle 24

We studied a set of 74 CMEs, with shedding the light on the halo-CMEs (HCMEs), that are associated with decametric – hectometric (DH) type-II radio bursts (1–16?MHz) and solar flares during the period 2008–2014. The events were classified into 3 groups (disk, intermediate, and limb events) based on their longitudinal distribution.We found that the events are mostly distributed around 15.32° and 15.97° at the northern and southern solar hemispheres, respectively. We found that there is a clear dependence between the longitude and the CME’s width, speed, acceleration, mass, and kinetic energy. For the CMEs’ widths, most of the events were HCMEs (～62%), while the partial HCMEs comprised ～35% and the rest of events were CMEs with widths less than 120°. For the CMEs’ speeds, masses, and kinetic energies, the mean values showed a direct proportionality with the longitude, in which the limb events had the highest speeds, the largest masses, and the highest kinetic energies. The mean peak flux of the solar flares for different longitudes was comparable, but the disk flares were more energetic. The intermediate flares were considered as gradual flares since they tended to last longer, while the limb flares were considered as impulsive flares since they tended to last shorter.A weak correlation (R?=?0.32) between the kinetic energy of the CMEs and the duration of the associated flares has been noticed, while there was a good correlation (R?=?0.76) between the kinetic energy of the CMEs and the peak flux of the associated flares. We found a fair correlation (R?=?0.58) between the kinetic energy of the CMEs and the duration of the associated DH type-II radio bursts.     

Latitudinal distribution of solar coronal active regions

We studied the cyclic evolution of the latitudinal distribution of solar coronal active regions based on daily images from SOHO EIT for the period 1995–2017. Fully automated software was used, which included the following steps: initial preparation of images in the data series, normalization of histograms and correction of limb brightening, segmentation of images using threshold intensity values obtained from their histograms, scanning of segmented images in heliographic coordinates and obtaining profiles of latitudinal distribution of coronal active regions for each image of the data series. From the output data, we obtained a temporary change in the latitudinal distribution profiles and the migration of activity centers on the solar disk. From the period of minimum activity to the next minimum in both hemispheres, activity centers begin to migrate from high latitudes towards the equator. At the same time, the general center of activity repeatedly changes the direction of migration. The latitudinal distribution of the so-called presence factor of coronal active regions closely resembles the magnetic butterfly diagram, which proves their direct causal relationships. Variations in the presence factor of coronal active regions are correlated with cyclic variations in the sunspot daily numbers.     

Plasma polarization electric field derived from radio sounding of solar wind acceleration region with spacecraft signals

Presented is the analytical approximation of averaged solar wind velocity radial dependence in the solar wind acceleration region at heliolatitudes below 60° under low and moderate solar activity. This empirical approximation is based on the data of radio sounding of the solar corona with radio signals from various spacecraft. Deduced is an equation connecting the solar wind velocity radial dependence and the radial dependence of solar wind plasma polarization electric field intensity. This allows constructing a semi-empirical radial dependence of plasma polarization electric field corresponding to the empirical radial dependence of solar wind velocity. Main properties of the semi-empirical dependence, which is based on radio sounding data, are described.     

Solar sailing in realistic mode,based on a locally-optimal control laws

Ballistic design of solar sailing missions in the solar system is composed of defining the design parameters, the control programs, and the trajectories that provide performance goals of a flight. The use of a solar sail spacecraft imposes specific restrictions on mission parameters that include the degradation limit on the flight duration, the maximum temperature of solar sail's surface, the minimum distance from the Sun, the maximum angular velocity of the spacecraft's rotation and others.Many authors considered the impact of these restrictions on the design of the mission separately, but they used a sophisticated method of finding the exact optimal motion control or applied the most straightforward laws of motion control. This paper uses local-optimal control laws at the complete mathematical models of motion and functioning of solar sail spacecraft to describe a technique of designing interplanetary missions. The described method avoids the need to obtain an accurate optimal solution to the control problem and does not cause significant computational difficulties.