Negative VTEC with spherical harmonic expansion model under different solar activity conditions

The Earth’s ionosphere can be described by a spherical harmonic (SH) expansion up to a specific degree. However, there exist negative vertical total electron content (VTEC) values in the global ionosphere map (GIM) with the SH expansion model. In this contribution, we specifically investigated the negative VTEC values that are induced by the SH expansion model and validated the performance of the inequality-constrained least squares (ICLS) method in eliminating the negative VTEC values. The GPS data from 2004 to 2017 was selected to cover one solar cycle and the experiments under different solar activity conditions were analyzed. The results in our work show that the occurrence of the negative VTEC values is attributed to the deficiency of the SH expansion model when the VTEC itself is small instead of the unevenly distribution of the GNSS stations. The negative VTEC values appear periodically in the temporal domain, showing apparently one year and half year periods. During one year, two peaks in June and December can be observed in the time series of the negative VTEC values. The number of negative VTEC values in June is obvious larger than that in December. During one solar cycle, the number of negative VTEC values under quiet solar activity condition is obvious larger than that under strong solar activity condition. In the spatial domain, the appearance of the negative VTEC values is strongly related with the movement of the subsolar point. In the latitude of the subsolar point has the largest magnitude, the negative values will appear on the opposite hemisphere and the further from the subsolar point the more negative values. The maximum number of the negative VTEC values in the southern hemisphere appears in June, while the peak value in the northern hemisphere appears in December. The maximum number of negative VTEC values in the southern hemisphere is generally larger than that in the northern hemisphere. In addition, the negative VTEC values are distributed both at middle latitude and high latitude in the southern hemisphere, while they are mainly distributed at high latitude in the northern hemisphere. When the ICLS method is used, the negative VTEC values can be eliminated efficiently and it has nearly no influence on the positive VTEC values. The ICLS method can also improve the receiver’s differential code bias (DCB) and significantly decrease the unreasonable negative slant TEC (STEC) values along the lines of sight. Using the final GIM product of the Jet Propulsion Laboratory (JPLG) as a reference, the root mean square (RMS) of the ICLS solution shows maximum 25%, 20% and 45% improvement relative to the least squares (LS) solution at northern high latitude, southern middle latitude and southern high latitude, respectively.     

SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

Aspects related to variability of radiative cooling by NO in lower thermosphere,TEC and O/N2 correlation,and diffusion of NO into mesosphere during the Halloween storms

Nitric Oxide is a very important trace species which plays a significant role acting as a natural thermostat in Earth’s thermosphere during strong geomagnetic activity. In this paper, we present various aspects related to the variation in the NO Infrared radiative flux (IRF) exiting the thermosphere by utilizing the TIMED/SABER (Thermosphere Ionosphere Mesosphere Energetics and Dynamics/ Sounding of the Atmosphere using Broadband Emission Radiometry) observational data during the Halloween storm which occurred in late October 2003. The Halloween storm comprised of three intense-geomagnetic storms. The variability of NO infrared flux during these storm events and its connection to the strength of the geomagnetic storms were found to be different in contrast to similar super storms. The connection between the quantum of energy outflux from the upper atmosphere into space in terms of NO IRF and the duration of storms is established. The NO radiative cooling, and the closely correlated depletion in O/N₂ ratio are controlled by the Joule heating intensity (proxied by AE-index). The collisional excitation rate of NO, calculated using the modelled datasets of WACCM-X (Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension), correlates well with the observed pattern of radiative emission by NO. Observational datasets from TIMED/GUVI (Global Ultra-Violet Imager) and MIT Haystack observatory madrigal GNSS (Global navigation satellite system) total electron content (TEC) database shows that the TEC and O/N₂ enhancement in low-mid northern hemispheric latitudes are mainly controlled by the z-component of Interplanetary magnetic field (IMF-B_z). The penetration of eastward electric field during the storm events is found to be responsible for the overall enhancement of TEC. The contribution of enhanced day-side TEC in observed variation of O/N₂ ratio by GUVI is also reported. It is also seen that during substorms related events the night-time polar region experiences more cooling due to NO than the daytime polar region. The connections between the mid- and low-latitude enhancement in NO IRF with the propagation of LSTIDs (Large-scale traveling ionospheric disturbances) in combination with the O/N₂ variability, and the altitudinal variation in NO flux with the progression of the storm is also investigated. This study presents the evidence on the role of diffusion processes in the large scale enhancement of NO in the mesospheric altitudes.     

考虑界面损伤的立式贮存药柱结构完整性研究

为了进一步结合实际分析固体火箭发动机药柱在立式贮存条件下的结构完整性，考虑推进剂/衬层界面损伤模式在复杂应力条件下具有多样性。以某型固体火箭发动机为例，与常规将衬层设置为粘接单元相比，模型在推进剂与绝热层之间设置粘接接触。对固体火箭发动机在立式贮存环境时经历固化降温、充气内压和重力载荷联合作用下有无界面损伤时的发动机进行仿真分析。结果表明:界面损伤的存在导致推进剂/绝热层界面这个薄弱环节更危险；该型固体火箭发动机药柱在充气内压增大过程中在人工脱粘层根部部位应力呈先增大后减小趋势；在充气内压达到0.085MPa之前，推进剂与绝热层之间考虑界面损伤时，推进剂在垂直于轴向的靠近人工脱粘层根部部分更容易损伤，之后则推进剂垂直于轴向的初始点更容易损伤。该结论可以为固体火箭发动机结构完整性精确仿真提供一定的指导。     

Ionospheric storm due to solar Coronal mass ejection in September 2017 over the Brazilian and African longitudes

Coronal mass ejection (CME) occurs when there is an abrupt release of a large amount of solar plasma, and this cloud of plasma released by the Sun has an intrinsic magnetic field. In addition, CMEs often follow solar flares (SF). The CME cloud travels outward from the Sun to the interplanetary medium and eventually hits the Earth’s system. One of the most significant aspects of space weather is the ionospheric response due to SF or CME. The direction of the interplanetary magnetic field, solar wind speed, and the number of particles are relevant parameters of the CME when it hits the Earth’s system. A geomagnetic storm is most geo-efficient when the plasma cloud has an interplanetary magnetic field southward and it is accompanied by an increase in the solar wind speed and particle number density. We investigated the ionospheric response (F-region) in the Brazilian and African sectors during a geomagnetic storm event on September 07–10, 2017, using magnetometer and GPS-TEC networks data. Positive ionospheric disturbances are observed in the VTEC during the disturbed period (September 07–08, 2017) over the Brazilian and African sectors. Also, two latitudinal chains of GPS-TEC stations from the equatorial region to low latitudes in the East and West Brazilian sectors and another chain in the East African sector are used to investigate the storm time behavior of the equatorial ionization anomaly (EIA). We noted that the EIA was disturbed in the American and African sectors during the main phase of the geomagnetic storm. Also, the Brazilian sector was more disturbed than the African sector.