Temporal evolution analysis of storm-enhanced density during an intense magnetic storm on March 2015

An intense Storm Enhancement Density (SED) event with the magnetic storm occurred on 17–24 March 2015 has been investigated. The morphological character of the SED during different phase of the magnetic storm is examined and compared with the non-storm time. Three intensity indexes, i.e., “general” SED index, “heavy” SED index and “severe” SED index, are defined to represent the intensity of SED respectively represented by the numbers of the ionospheric total electron content (TEC) grids with TEC > 60 TECu, TEC > 80 TECu and TEC > 100 TECu. The temporal evolution of the SED intensity indexes during a time span covering the non-storm time and the magnetic storm time have also been investigated. The SED exhibits a shape with two parallel slender troughs in the middle and low latitudes during the non-storm time and then gradually develops into an ellipse structure as the development of magnetic storm. The intensity of SED and the fluctuation of the TEC evolution are generally corresponding to the fluctuation of Dst index. The analyzing results enrich our understanding of the temporal and spatial evolution of the ionospheric SED.     

Preliminary analysis of the quality and positioning performance of BDS-3 global interoperable signal B1C&B2a

BeiDou-3 Navigation Satellite System (BDS-3) satellites are equipped with the new generation GNSS signals B1C and B2a, which support the interoperability with GPS and Galileo systems. In this study, the pseudo-range multipath error and carrier phase observation noise of the BDS-3 B1C and B2a signals were evaluated based on zero baseline measurements from the day of year (DOY) 113 to 116, 2020. Further, the precision and performance of the single point positioning (SPP) and precise point positioning (PPP) are assessed at 9 stations. This assessment manifests that the standard deviations (STDs) of the pseudo-range multipath error are about 0.09 ~ 0.22 m, while STDs of the carrier phase observation noise are about 0.075 mm. For the single-frequency SPP, its positioning precision is about 2.03 ~ 4.85 m and 3.29 ~ 10.73 m at the 99.99% confidence level in horizontal and vertical directions, respectively, while the dual-frequency SPP precision is about 1.92 ~ 8.02 m and 4.81 ~ 12.77 m in horizontal and vertical directions, respectively. For the daily static PPP, the convergence time is about 20 ~ 30 min, while the daily positioning precision can reach 1.38 ~ 4.42 cm and -1.31 ~ 4.34 cm in horizontal and vertical directions, respectively. In general, the quality and the SPP and PPP performance of the BDS-3 B1C&B2a signals are comparable to the GPS and Galileo.     

Observation of aphelion cloud belt over Martian tropics,its evolution,and associated dust distribution from MCS data

The present study uses five Martian years of observations from Mars Climate Sounder onboard Mars Reconnaissance Orbiter for investigating the Aphelion Cloud Belt (ACB) over the tropics. Analysis of zonal mean water ice column opacity suggests that the spatial extension of the ACB is mainly confined over the tropics and mid-latitudes (-20 – 40°N) during L_S ~ 45 – 135° (L_S = 0° signifies northern spring equinox). The ACB is seen primarily in the nighttime only due to the truncation of the daytime profile observations at significantly higher altitudes (at ~30 km). Zonal mean ice extinction profiles show ACB’s altitudinal range within ~10 – 40 km, and the existence of a thin cloud band in the absence of a thick ACB during aphelion season. Three phases of the ACB could be identified as the formation phase during L_S = 45 – 75° (phase 1), the peak phase during L_S = 76 – 105° (phase 2), and the decaying phase during L_S = 106 – 135° (phase 3). Observation of the cloud latitude belt shows a northward movement starting from phase 2, prominent over regions nearby Lunae Planum and Xanthe Terra. During this phase, the top level of thick clouds within the ACB decreases to ~20 km in the southern hemisphere, while it increases a little over the northern hemisphere (NH). The decreasing tendency continues in phase 3 over the entire region ?10 – 10°N, and the thick cloud base moves higher over the NH, though the vertical depth of it becomes narrower than phase 2. Temperature profiles do not show any noticeable influence on the northward evolution of the ACB. However, the study at a regional level indicates a possible association of upper tropospheric dustiness with the ACB’s evolution. The mechanism is evident in the correlation analysis mostly at an altitude range of ~18 – 35 km. The migrating semidiurnal tide (SMD) as a proxy of dust or water ice forcing, and the calculated upper tropospheric dust radiative heating, shows an apparent northward movement of their peak amplitude within the three phases of the ACB. This match between the spatiotemporal variations of the SMD and the water ice was not observed previously. However, the correlating behavior seems to be prominent in the areas nearby Lunae Planum and Xanthe Terra and the upper-tropospheric region of the atmosphere.     

Installation and preliminary data analysis of Penang magnetic data Acquisition system (MAGDAS) in Malaysia

MAGDAS PEN was established on 19th September 2019 as one of the MAGDAS observatory arrays located at Universiti Sains Malaysia (USM) (5.15°, 100.50°). The main objective of the MAGDAS project is to monitor global electromagnetic and the ambient plasma density in the geospace environment. This installation has contributed to a better understanding of the Sun-Earth coupling system. This paper presents the installation process of one of the MAGDAS magnetometers named YU-8 T magnetic sensor and provides a preliminary analysis of geomagnetic HDZ components amplitude-time that was observed at PEN station. A one-month HDZ-geomagnetic field data was processed from 1st November to 30th November 2019. The daily variations with a consistent pattern in delta H geomagnetic field components are observed throughout the day with eastward electric field effects that are observed during solar peak hours. The delta H-component gradually increases around 0700LT and reaches the maximum reading at 1300LT with a range of value ~ 40-70nT. The value slowly decreases that started from 1400LT until the night time. The reading during the night time shows a constant variation with magnitude varies in between ?10nT to + 10nT. The average H-component value of the night time is used as the baseline for the observation system. Overall, the observed trends portray a good sign of solar quiet field and Sq with no solar-terrestrial disturbances.     

Prediction of geosynchronous electron fluxes using an artificial neural network driven by solar wind parameters

There are hundreds of satellites operating at the geosynchronous (GEO) orbit where relativistic electrons can cause severe damage. Thus, predicting relativistic electron fluxes is significant for spacecraft safety. In this study, using GOES satellite data during 2011–2020, we propose two neural network models with two hidden layers to predict geosynchronous relativistic electron fluxes at two energy channels (>0.8 MeV and > 2 MeV). The number of input neurons of the two channels (>0.8 MeV and > 2 MeV) are determined to be 36 and 44, respectively. The > 0.8 MeV model has 22 and 9 neurons in the hidden layers, while the > 2 MeV model has 25 and 15 neurons in the hidden layers. The input parameters include the north–south component of the interplanetary magnetic field, solar wind speed, solar wind dynamic pressure and solar wind proton density. Through the analysis of different time delays, we determine that the optimal time delays of two energy channels (>0.8 MeV and > 2 MeV) are 8 days and 10 days, respectively. The training set and validation set (Jan 2011-Dec 2018) are divided by the 10-fold cross-validation method, and the remaining data (Jan 2019-Feb 2020) is used to analyze the model performance as a test set. The prediction results of both energy channels show good agreement with satellite observations indicated by low RMSE (～0.3 cm^-2sr^-1s^?1), high PE (～0.8) and CC (～0.9). These results suggest that only using solar wind parameters is capable of obtaining reasonable predictions of geosynchronous relativistic electron fluxes.     

An engineering model of Titan surface winds for Dragonfly landed operations

Winds near the ground on Titan for the Dragonfly landing site (near Selk crater, 10°N) for the mid-2030s (Titan late southern summer, L_s ~ 310°) are estimated for mission design purposes. Prevailing winds due to the global circulation are typically 0.5 m/s, and do not exceed 1 m/s. Local terrain-induced flows such as slope winds appear to be similarly capped at 1 m/s. At various landing sites and times, these two contributions will vectorially combine to yield steady winds (for part of a Titan day, Tsol) of up to 2.0 m/s, but typically less – the slope wind component will be small in the mid-morning. In early afternoon, as on Earth and Mars, solar-driven convection in the planetary boundary layer will cause wind fluctuations of the order of 0.1 m/s, varying with a typical timescale of ~1000 s. Occasionally this convection organizes into coherent ‘dust devil’ vortices: detectable vortices with speeds of 1 m/s are predicted about once per Titan day. We have introduced the convective velocity scale combined with the advection time of PBL cells as a metric to derive the frequency of occurrence of gusts associated with convective vortices (‘dust devils’). Maximum possible vortex winds on Titan of 2.8 m/s may be expected only once per 40 Tsols, and define the maximum wind (4.8 m/s at 10 m height) that Dragonfly must tolerate without damage. The applicability of different wind combinations, scaled to the height of relevant Dragonfly components above the ground (e.g. the maximum corresponds to 3.9 m/s at 1.3 m height) by a logarithmic wind profile, to Dragonfly design and operations are discussed.     

The ultraviolet sky: An overview from the GALEX surveys

The Galaxy Evolution Explorer (GALEX) has performed the first surveys of the sky in the ultraviolet (UV). Its legacy is an unprecedented database with more than 200 million source measurements in far-UV (FUV) and near-UV (NUV), as well as wide-field imaging of extended objects, filling an important gap in our view of the sky across the electromagnetic spectrum. The UV surveys offer unique sensitivity for identifying and studying selected classes of astrophysical objects, both stellar and extra-galactic. We examine the overall content and distribution of UV sources over the sky, and with magnitude and color. For this purpose, we have constructed final catalogs of UV sources with homogeneous quality, eliminating duplicate measurements of the same source. Such catalogs can facilitate a variety of investigations on UV-selected samples, as well as planning of observations with future missions. We describe the criteria used to build the catalogs, their coverage and completeness. We included observations in which both the far-UV and near-UV detectors were exposed; 28,707 fields from the All-Sky Imaging survey (AIS) cover a unique area of 22,080 square degrees (after we restrict the catalogues to the central 1° diameter of the field), with a typical depth of ∼20/21mag (FUV/NUV, in the AB mag system), and 3008 fields from the Medium-depth Imaging Survey (MIS) cover a total of 2251 square degrees at a depth of ∼22.7mag. The catalogs contain ∼71 and ∼16.6 million sources, respectively. The density of hot stars reflects the Galactic structure, and the number counts of both Galactic and extra-galactic sources are modulated by the Milky Way dust extinction, to which the UV data are very sensitive.     

A Gemini observation of the anomalous X-ray pulsar 1RXS J170849-400910

The anomalous X-ray pulsars (AXPs) represent a growing class of neutron stars discovered at X-ray energies. While the nature of their multi-wavelength emission mechanism is still under debate, evidence has been recently accumulating in favor of their magnetar nature. Their study in the optical and infrared (IR) wavelengths has recently opened a new window to constrain the proposed models. We here present a brief overview of AXPs and our Gemini-South observation of 1RXS J170849-400910, which is a relatively bright AXP discovered with ROSAT and later found to be an 11 s X-ray pulsar by ASCA. The observation was taken with the near-IR imager Flamingos in J (1.25 μm), H (1.65 μm), and K_s (2.15 μm). We confirm the recent detection by (ApJ, 589, L93–L96) of a source coincident with the CHANDRA source (candidate ‘A’). Our derived magnitudes of J = 20.6 (0.2), H = 18.6 (0.2), and K_s = 17.1 (0.2) are consistent with those derived by (ApJ, 589, L93–L96), and indicate that if this source is indeed the IR counterpart to 1RXS J170849-400910, then there is no evidence of variability from this AXP. However, given the lack of IR variability and the relatively high IR to X-ray flux of this source when compared to the other AXPs, we conclude that this source is unlikely the counterpart of the AXP, and that the other source (candidate ‘B’) within the CHANDRA error circle should not be ruled out as the counterpart. Further monitoring of these sources and a deep observation of this complex field are needed to confirm the nature of these sources and their association with the AXP.     

垂直激波条件下能量粒子加速机制的模拟研究

在具有湍动的磁场和垂直激波条件下对大量测试粒子的轨迹进行了数值计算,研究了激波强度和粒子初始能量对于粒子穿越激波的平均能量变化的影响,分析了漂移加速(SDA)在不同条件下对粒子加速的贡献,并给出了一个与数值结果相符合的漂移加速理论公式△E=amvivup(1-1/s).结果表明,加入磁场湍流后,垂直激波条件下粒子仍主要受到漂移加速作用,而基于粒子引导中心的耗散漂移加速理论在此条件下失效.      

Two years of INTEGRAL result on X-ray binaries

INTEGRAL is the ESA lead International Gamma-Ray Astrophysical Laboratory, successfully launched the 17th October 2002 from Baikonur with a Proton vehicle. In view of the high sensitivity of the two γ-ray instruments IBIS and SPI and their capability to provide at the same time image, spectra and time profiles of all the sources in their wide field of view, a key project was approved as “Core Programme” to obtain deep observations of the Galactic Centre (GCDE) and to exploit regular scan of the whole Galaxy Plane since the beginning of the mission. This paper will briefly review the main astrophysical results obtained in the field of high energy Galactic sources with the INTEGRAL/IBIS γ-ray Imager onboard INTEGRAL, and make a comparison with the previous scenario depicted by the BeppoSAX and RXTE results.     

Accuracy analysis of GNSS-IR snow depth inversion algorithms

In recent years, with the continuous development of Global Navigation Satellite System (GNSS), it has been applied not only to navigation and positioning, but also to Earth surface environment monitoring. At present, when performing GNSS-IR (GNSS Interferometric Reflectometry) snow depth inversion, Lomb-Scargle Periodogram (LSP) spectrum analysis is mainly used to calculate the vertical height from the antenna phase center to the reflection surface. However, it has the problem of low identification of power spectrum analysis, which may lead to frequency leakage. Therefore, Fast Fourier Transform (FFT) spectrum analysis and Nonlinear Least Square Fitting (NLSF) are introduced to calculate the vertical height in this paper. The GNSS-IR snow depth inversion experiment is carried out by using the observation data of P351 station in PBO (Plate Boundary Observatory) network of the United States from 2013 to 2016. Three algorithms are used to invert the snow depth and compared with the actual snow depth provided by the station 490 in the SNOTEL network. The observations data of L1 and L2 bands are respectively used to find the optimal combination between different algorithms further to improve the accuracy of GNSS-IR snow depth inversion. For L1 band, different snow depths correspond to different optimal algorithms. When the snow depth is less than 0.8 m, the inversion accuracy of NLSF algorithm is the highest. When the snow depth is greater than 0.8 m, the inversion accuracy of FFT algorithm is higher. Therefore, according to the different snow depth, a combined algorithm of NLSF + FFT is proposed for GNSS-IR snow depth inversion. Compared with the traditional LSP algorithm, the inversion accuracy of the combined algorithm is improved by 10%. For L2 band data, the results show that the accuracy of snow depth inversion of various algorithms do not change with the variations of snow depth. Among the three single algorithms, the inversion accuracy of FFT algorithm is better than that of LSP and NLSF algorithms.     

Signatures of magnetic reconnection in solar radio observations?

Magnetic reconnection occurs during eruptive processes (flares, CMEs) in the solar corona. This leads to a change of magnetic connectivity. Nonthermal electrons propagate along the coronal magnetic field thereby exciting dm- and m-wave radio burst emission after acceleration during reconnection or other energy release processes in heights of some Mm to ⩾700 Mm. We summarize the results of some case studies which can be interpreted as radio evidence of magnetic reconnection: under certain conditions, simple spectral structures (pulsation pulses, reverse drift bursts) are formed by simultaneously acting but widely spaced radio sources. Narrowband spikes are emitted as a side-effect during large-scale coronal loop collisions. In dynamic radio spectra, the lower fast mode shock formed in the reconnection outflow appears as type II burst-like but nondrifting emission lane. It has been several times observed at the harmonic mode of the local plasma frequency between 250 and 500 MHz and at heights of ≈200 Mm.     

Interplanetary proton cumulated fluence model update

Solar particle events leading to important increase of particle fluxes at energies of order of magnitude ranging from MeV to GeV constitute an important hazard for space missions. They may lead to effects seen in microelectronics or damage to solar cells and constitute a potential hazard for manned missions. Cumulative damage is commonly expressed as a function of fluence which is defined as the integral of the flux over time. A priori deterministic estimates of the expected fluence cannot be made because over the time scale of a space mission, the fluence can be dominated by the contribution of a few rare and unpredictable high intensity events. Therefore, statistical approaches are required in order to estimate fluences likely to be encountered by a space mission in advance. This paper extends work done by Rosenqvist et al. [Rosenqvist, L., Hilgers, A., Evans, H., Daly, E., Hapgood, M., Stamper, R., Zwickl, R., Bourdarie, S., Boscher, D. Toolkit for updating interplanetary proton-cumulated fluence models. J. Spacecraft Rockets, 42(6), 1077–1090, 2005] to describe an updated predictive engineering model for the proton interplanetary fluence with energies >30 MeV. This model is derived from a complete list of solar proton fluences based on data from a number of calibrated sources covering almost three solar cycles.     

Comparative evaluation of five global gravity models over a part of the Bay of Bengal

Several global gravity models (GGMs) are freely available in the public domain, which can be utilised to study the earth's gravity field in almost every part of the globe. The present study compared the free-air gravity anomalies calculated from the five GGMs EGM2008, EIGEN6C4, GECO, XGM2019e_2159, and SGG-UGM-2 archived by the International Centre for Global Earth Models (ICGEM) with respect to shipborne gravity in the Bay of Bengal. The average correlation coefficient and covariance are ～ 96 % and ～ 450mGal². The mean difference between the shipborne and the modelled gravity is ? 5 mGal. Relatively higher amplitude gravity differences observed at the continental-oceanic transition, the 85°E and Ninetyeast ridges, and the western basin are possibly due to high gradient, dominant density contrasts, and rugged topography. The average standard deviation and root-mean-square-error (RMSE) of the differences are ～ 6.5 mGal and ～ 7.5 mGal. A significantly lower standard deviation and RMSE found for the models generated at higher degree/order compared to lower degree/order is due to diminishing omission error of the GGMs with increasing degrees of truncation. The spectral analysis and coherence estimation of the modelled gravity demonstrate excellent correspondence for anomalies wider than ～ 25 km. The agreement between anomaly amplitudes and shapes and calculated statistics indicates that the accuracy and resolution of the modelled gravity data are certainly good enough for regional-scale studies aiming to map Moho topography and mantle structures. However, the delineation of shorter wavelength features originating from the shallow-level basement/sedimentary might be uncertain and requires further validations. The present study confirms that all five models show promising results in terms of their accuracy and resolution. Moreover, the SGG-UGM-2 and XGM2019e_2159 models compare favourably with respect to the GECO, EIGEN6C4 and EGM2008 models in the Bay of Bengal.     

Observation bio-effect of SW-global solar radiation in Ilorin in the tropics

The annual variation of global (total) solar radiation measured over four (4) years (1995–1998) at Ilorin (8°32′N, 4°34′E) was studied and the bio-effects of the variation is herein reported. The weekly averages of the solar radiation flux were obtained and plotted. The profile of the weekly averages of the insolation shows two notable “wells” corresponding to the harmattan dry season and the rainy season. It was deduced that the duration of the seasons was directly proportional to the sizes of the “wells”, while the depths of the “wells” correspond to the severity of the causes of the “wells”. The fall and rise rates of the depth of the edges of the wells are inversely proportional to the ecological effects of the variation of the radiation. The rates of decrease and increase of the radiation flux affect the depths of the “wells” and were found to be rapid and shallow for the “Harmattan Well”, and gentle and deep for the “Rain Well”. The width (duration), rate of change and depth of the “Harmattan Well” bring about stress for both animals and plants.     

The evolution of a complex solar radio burst corresponding to special configuration of microwave sources

A complex radio burst associated with periodic (∼1 and 6 min) pulsations and several kinds fine structures, e.g., normal- and reverse-drifting type III bursts, zebra patterns, and slowly drifting structure was observed with the radio spectrometers (1.0–2.0, 2.6–3.8, 5.2–7.6, and 0.65–1.5 GHz) at the National Astronomical Observatories of China (NAOC) in Beijing and Yunnan on 19 October 2001. In combination with the images of 17 and 34 GHz from NoRH and the magnetograms from MDI we reveal the existence and evolution of preexisting and new emerging sources, and find the horseshoe-shaped structure of microwave sources intensity during the late phase of the burst. Through the detailed comparison of the evolution of each source with the time profiles of radio bursts corresponding to these sources we indicate that the intimate correlation between the microwave sources evolution and the generation of the radio burst associated fine structures. Some fine structures can be considered as the MHD turbulence and plasma emission mechanism, based on the anisotropic beam instability and hybrid waves generations. From the characteristics of observations we may presume that the coronal magnetic structures should contain an extended coronal loop system and multiple discrete electrons acceleration/injection sites. The mechanisms of this complex radio burst are deal with the incoherent gyrosynchrotron emission from the trapped electrons and the coherent plasma emission from the non trapped electrons.     

Study of ionospheric D region changes during solar flares using MF radar measurements

During solar flares, the X-ray radiation suddenly increases, resulting in an increase in the electron density of the atmospheric D region and a strong absorption of short-wave radio waves. Based on Langfang medium frequency (MF) radar, this paper analyzed the variation characteristics of D region in the lower ionosphere from 62 km to 82 km. The analysis focused on multiple C-level and M-level solar flare events before and after the large-scale flare event at 11:53 (UT) on September 6, 2017. The results show that it is difficult to detect the electron density over 70 km in Langfang during solar flares, but the electron density value can be obtained as low as 62 km, and the stronger the flare intensity, the lower the detectable electron density height. Besides, the equal electron density height, the received power of X and O waves will also be significantly reduced during the flares, and the reduction of equal electron density height has a weak linear relationship with flare intensity.     

Influence of different data fusion methods on the accuracy of three-dimensional displacements field

Combining multi-source data can improve the accuracy and the spatial resolution of the three-dimensional (3-D) displacements field. How to effectively integrate multi-source data to obtain high-precision and high spatial resolution 3-D displacements field is worthy of further study. The stochastic model and fusion model of integrating multi-source data affect the accuracy of data fusion. In this paper, based on the least squares method, the effects of different stochastic models and data fusion models on the 3-D displacements field’s accuracy are studied. The optimal method for estimating large-scale 3-D displacements field from integrated InSAR, leveling and GPS measurements is obtained. Then we realize the integrating InSAR, leveling and GPS measurements to obtain the high-precision 3-D displacements velocity field in Tianjin (China) from 2016 to 2018. The results are validated with GPS measurements at 6 independent stations, with the root mean squares (RMS) residuals of the discrepancies being 2.39 mm/yr, 2.54 mm/yr and 2.83 mm/yr in eastern, northern and vertical directions, respectively. By comparing different stochastic models, the 3-D displacements field obtained from multi-source data is optimized by the variance component estimation-least squares method, which is better than weighted least squares (WLS) method. By comparing different data fusion models, the accuracy of the horizontal displacements velocity is better than that of interpolated GPS results. The horizontal displacements component has a great influence on the vertical displacements velocity accuracy in the process of acquiring the 3-D displacements velocity by integrating InSAR, GPS and leveling measurements. This study provides a reference method for integrating multi-source data to obtain 3-D displacements field. This method effectively utilizes the advantages of GPS, InSAR and leveling measurements, and extends the limitations of single technical in describing surface-time scale applications. The 3-D displacements information with a large spatial scale and high spatial resolution provide a reliable data basis for studying the crustal movement and its dynamic mechanism.     

Collaborative VLA,SOHO and RHESSI observations of evolving sources of energy release in the corona above active regions

Very Large Array (VLA) observations at 20 and 91 cm wavelength are compared with data from the SOHO (EIT and MDI) and RHESSI solar missions to investigate the evolution of decimetric Type I noise storms and Type III bursts and related magnetic activity in the photosphere and corona. The combined data sets provide clues about the mechanisms that initiate and sustain the decimetric bursts and about interactions between thermal and nonthermal plasmas at different locations in the solar atmosphere. On one day, frequent, low-level hard X-ray flaring observed by RHESSI appears to have had no clear affect on the evolution of two closely-spaced Type I noise storm sources lying above the target active region. EIT images however, indicate nearly continuous restructuring of the underlying EUV loops which, through accompanying low-level magnetic reconnection, might give rise to nonthermal particles and plasma turbulence that sustain the long-lasting Type I burst emission. On another day, the onset of an impulsive hard X-ray burst and subsequent decimetric burst emission followed the gradual displacement and coalescence of a small patch of magnetic magnetic polarity with a pre-existing area of mixed magnetic polarity. The time delay of the impulsive 20 and 91 cm bursts by up to 20 min suggests that these events were unlikely to represent the main sites of flare electron acceleration, but instead are related to the rearrangement of the coronal magnetic field after the main flare at lower altitude. Although the X-ray flare is associated with the decimetric burst, the brightness and structure of a long-lasting Type I noise storm from the same region was not affected by the flare. This suggests that the reconfiguration of the coronal magnetic fields and the subsequent energy release that gave rise to the impulsive burst emission did not significantly perturb that part of the corona where the noise storm emission was located.     

Improving regional ionospheric TEC mapping based on RBF interpolation

Radial basis function (RBF) interpolation with multi-quadric is developed to perform ionospheric total electron content (TEC) mapping for the Chinese region between 15°N ~ 40°N and 100°E ~ 125°E. TEC measurements from the Centre for Orbit Determination in Europe (CODE) covering the solar maximum year 2011 are used to investigate the performance of the proposed RBF interpolation method. The differences between the RBF interpolated TEC and the CODE TEC are within 0.5 TECU and the root mean square error (RMSE) is very small when 49 data points are used. The maximum difference is ~5 TECU and the error is less than 1 TECU with 25 samples. Our study suggests that a random distribution of measurement points gives smaller RMSEs than a homogenous distribution when the number of sample points is low. The study indicates that RBF interpolation offers a powerful and reliable tool for ionospheric TEC mapping.