Impact of local and global factors and meteorological parameters in temporal variation of atmospheric potential gradient

In this research work, we have performed comparative diurnal variations of atmospheric Potential Gradient (PG) of fair-weather days by using the data of three stations installed in Northern, Pakistan for the year 2018. We investigated the impact of both local and global factors and meteorological parameters in the diurnal variation of atmospheric Potential Gradient on the annual and seasonal time scale. We observed two peaks, primary and secondary. This is because of the land-based measurements of annual and seasonal variations. The annual average curve of Potential Gradient of all three stations: Islamabad (CES), Muzaffarabad (MZF), and Balakot (BKT) demonstrated a notable deviation from the standard oceanic Carnegie curve. The atmospheric Potential Gradient variations are due to numerous meteorological factors e.g., air pollution, humidity, aerosol particles, fog, and temperature. Among three stations, the MZF station is located in highland (mountainous) and it demonstrated a higher atmospheric Potential Gradient. We further differentiate the results of our three stations with global results for authenticity and observed coherence between them. In addition, a positive correlation of fair-weather Potential Gradient is observed with temperature and a notable correlation between relative humidity and atmospheric Potential Gradient for all the three observatories.     

晴天大气电场时空变化特征分析

利用子午工程5个大气电场观测站点近10年的近地面晴天大气电场观测数据,在年变化、季节变化、日变化三种不同时间尺度上进行了对比分析。结果表明:不同纬度站点的日平均晴天大气电场峰谷类型不同,且部分站点的波峰出现了逐年左移或右移的趋势;在电场幅值变化方面,位于中低纬度地区的站点呈现出逐年减小的变化特征,而位于中高纬度的站点呈现出逐年增加的变化趋势,且这种年变化均是线性的;多元回归分析表明,最大波峰出现时间与地理经度呈负相关,而与地理纬度呈正相关;各站点在近10年中均未出现明显的纬度效应;冬季的晴天大气电场日平均值水平较高,夏季较低,且各季节的日平均晴天大气电场最小值及最大值均近似呈正态分布;最小值及最大值在年度和季节两个时间尺度上的变化规律是基本一致的。这些研究结果揭示了晴天大气电场在不同时间尺度上的变化特征。      

西藏阿里地区的晴天大气电场特征

晴天大气电场是一个区域大气电场变化的基准场,是大气电场特性的研究基础。对于具有不同地质条件、地理环境和地形特点的区域,晴天大气电场特征具有明显区别。利用在西藏阿里地震台（80.12°E,32.51°N）安装的一台场磨式大气电场仪于2021年10月10日至2021年11月10日观测到的数据,通过对大气电场数据建模与分析,结合气象条件,得到西藏阿里地区23天晴天条件下的平均卡耐基曲线。经过平滑处理去除仪器噪声的波动,得到西藏阿里地区的晴天条件下的标准卡耐基曲线,进而与北京市昌平区十三陵观测台站（116.23°E,40.25°N）同期的晴天大气电场特征进行对比,并对其电场特性差异和原因进行了讨论。研究结果对分析高原地区晴天大气电场特征具有重要的参考价值和科学意义。      

TEC derived from some GPS stations in Nigeria and comparison with the IRI and NeQuick models

Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.     

Periodic variations of ionization rate in the low and middle atmosphere

The cosmic ray ionization source functions which were obtained using a simplified extensive air shower model are used to calculate the eleven year cycle, seasonal and diurnal variations of ionization rate in the low and middle atmosphere. The ionization source function, as a function of the penetrating depth and the energy of cosmic ray particles, is the ionization rate per unit depth for a unit flux of incoming cosmic ray particles with certain energy.The calculation of the eleven year cycle variation of ionization rate in the low and middle atmosphere due to the modulation of galactic cosmic ray intensity by solar activity shows that the amplitude is larger at a higher magnetic latitude and is generally larger at higher altitudes. The relative amplitude of fluctuation of the ionization peak value (at altitudes near 15 km) is up to 45% in the magnetic polar region. The ionization rate, due to the seasonal variation of the atmospheric density, varies from several per cent below the ionization peak to several tens per cent above the peak. This seasonal variation of ionization rate reaches 35% at 70 km. The diurnal variation of atmospheric densities caused by atmospheric tidal oscillation can produce a diurnal variation of the ionization rate to an amplitude of several per cent at altitudes above 40 km. The diurnal oscillation is less than 1% below 35 km.     

Validation of IRI-2012 TEC model over Ethiopia during solar minimum (2009) and solar maximum (2013) phases

This paper investigates the capacity of the latest version of the International Reference Ionosphere (IRI-2012) model in predicting the vertical Total Electron Content (vTEC) over Ethiopian regions during solar minimum (2009) and solar maximum (2013) phases. This has been carried out by comparing the IRI-2012 modeled and experimental vTEC inferred from eight ground based dual frequency GPS (Global Positioning System) receivers installed recently at different regions of the country. In this work, the diurnal, monthly and seasonal variation in the measured vTEC have been analyzed and compared with the IRI-2012 modeled vTEC. During the solar minimum phase, the lowest and highest diurnal peak of the experimental vTEC are observed in July and October, respectively. In general, the diurnal variability of vTEC has shown minimum values around 0300 UT (0600 LT) and maximum values between around 1000 and 1300 UT (1300 and 1600 LT) during both solar activity phases. Moreover, the maximum and minimum monthly and seasonal mean hourly vTEC values are observed in October and July and in the March equinox and June solstice, respectively. It is also shown that the IRI-2012-model better predicts the diurnal vTEC in the time interval of about 0000–0300 UT (0300–0600 LT) during the solar minimum phase. However, the model generally overestimates the diurnal vTEC except in the time interval of about 0900–1500 UT (1200–1800 LT) during the solar maximum phase. The overall result of this work shows that the diurnal vTEC prediction performance of the model is generally better during the solar minimum phase than during solar maximum phase. Regarding the monthly and seasonal prediction capacity of the model, there is a good agreement between the modeled and measured monthly and seasonal mean hourly vTEC values in January and December solstice, respectively. Another result of the work depicts that unlike the GPS–TEC the IRI-2012 TEC does not respond to the effect resulted from geomagnetic storms.     

Behavior of the equivalent slab thickness over three European stations

The total electron content (TEC) derived from the global positioning system (GPS) and the F2-layer peak electron density obtained from Digisonde data have been used to study the diurnal, seasonal and solar activity variations of the ionospheric equivalent slab thickness (τ) over three European stations located at Pruhonice (50.0°N, 15.0°E), Ebro (40.8°N, 0.5°E) and El Arenosillo (37.1°N, 353.3°E). The diurnal variation of the τ is characterized by daytime values lower than nighttime ones for all seasons at low solar activity while daytime values larger than nighttime characterizes the diurnal variation for summer at high solar activity. A double peak is noticeable at dusk and at dawn, better expressed for winter at low solar activity. The seasonal variations of τ depend on local time and solar activity, the daytime values of τ increases from winter to summer whereas nighttime values of τ show the opposite. The effect of the solar activity on τ depends on local time and season, there being very sensitive for winter nighttime values of τ. The results of this study are compared with those presented by other authors.     

海南地区电离层闪烁观测与GISM模式预测的比较分析

为了获得全球电离层闪烁模式GISM在中国低纬地区预测的精度和可靠性,利用海南三亚GPS电离层闪烁监测系统一年的观测数据与模式预测结果进行对比分析.结果表明,在太阳活动低年,GISM能较好地反映海南地区电离层闪烁的季节变化、日变化和空间分布特性;对于季节变化,模式与观测结果在中等强度闪烁条件下较为相符,而在强闪烁和弱闪烁条件下有不同程度的偏差;模式预测的闪烁日变化与实际观测基本一致,但在闪烁发生率出现最大值的时间上模式预测要滞后约1 h左右;在电离层闪烁发生率的空间分布上,模式预测与实际观测较为相符,即海南地区南面电离层闪烁发生率高于北面.      

Variations of Thermosphere Density Based on Joint Analysis of In-Situ Measurement Data From Shenzhou Spacecrafts

Based on the measurements made by Atmospheric Density Detectors (ADDs) onboard Chinese spacecraft Shenzhou 2-4, the variations of thermosphere density are revealed. During the quiet period, the density at spacecraft altitude of 330～410km exhibited a dominant diurnal variation, with high value on dayside and low value on nightside. The ratio of the diurnal maximum density to the minimum ranged from 1.7 to 2.0. The ratio shows a positive correlation with the level of solar activity and a negative correlation with the level of geomagnetic activity. When a geomagnetic disturbance comes, the atmospheric density at the altitude of 330～410km displayed a global enhancement. For a strong geomagnetic disturbance, the atmospheric density increased by about 56%, and reached its maximum about 6～7 hours after the geomagnetic disturbance peak. The density asymmetry was also observed both in the southern and northern hemisphere during the geomagnetic disturbance peak.      

北京地区电离层电子含量的某些特性

从对1980年至今的北京电离层电子含量的常规观测资料分析发现:电离层电子含量(IEC)的周日变化曲线除了绝大多数呈单峰变化外, 还有非扰动状态下的周日双峰变化、夜增长现象以及日出后不久的次极大正午之后的主极大等.本文着意于报道北京地区的IEC双峰特性和夜增现象及其随太阳活动的周期变化, 并与其它报道进行比较.      

Global structure,seasonal and interannual variability of the eastward propagating tides seen in the SABER/TIMED temperatures (2002–2007)

The present paper is focused on the global spatial (altitude and latitude) structure, seasonal and interannual variability of the most stable in amplitude and phase eastward propagating diurnal and semidiurnal tides with zonal wavenumbers 2 and 3 derived from the SABER/TIMED temperatures for full 6 years (January 2002–December 2007). The tidal results are obtained by an analysis method where the tides (migrating and nonmigrating) and the planetary waves (zonally travelling, zonally symmetric and stationary) are simultaneously extracted from the satellite data. It has been found that the structures of the eastward propagating diurnal tides with zonal wavenumbers 3 and 2 change from antisymmetric with respect to the equator below ∼85 km height, to more symmetric above ∼95 km. The seasonal behavior of the DE3 is dominated by annual variation with maximum in August–September reaching average (2002–2007) amplitude of ∼15 K, while that of the DE2 by semiannual variation with solstice maxima and with average amplitude of ∼8 K. These tides revealed some interannual variability with a period of quasi-2 years. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 2 in the southern hemisphere (SH) is dominated by annual variation with maximum in the austral summer (November–January) while that in the northern hemisphere (NH) by semiannual variation with equinoctial maxima. The SE2 maximizes near 115 km height and at latitude of ∼30° reaching an average amplitude of ∼6 K. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 3 in both hemispheres indicates a main maximum during June solstice and a secondary one during December solstice. The tide maximizes near 110–115 km height and at a latitude of ∼30° reaching an average amplitude of ∼4.8 K in the SH and ∼4 K in the NH. The tidal structures of the two eastward propagating semidiurnal tides are predominantly antisymmetric about the equator.     

Variability study of the crest-to-trough TEC ratio of the equatorial ionization anomaly around 120°E longitude

In this paper, latitudinal profiles of the vertical total electron content (TEC) deduced from the dual-frequency GPS measurements obtained at ground stations around 120°E longitude were used to study the variability of the equatorial ionization anomaly (EIA). The present study mainly focuses on the analysis of the crest-to-trough TEC ratio (TEC-CTR) which is an important parameter representing the strength of EIA. Data used for the present study covered the time period from 01 January, 1998 to 31 December, 2004. An empirical orthogonal function analysis method is used to obtain the main features of the TEC-CTR’s diurnal and seasonal variations as well as its solar activity level dependency. Our results showed that: (1) The diurnal variation pattern of the TEC-CTR at 120°E longitude is characterized by two remarkable peaks, one occurring in the post-noon hours around 13–14 LT, and the other occurring in the post-sunset hours around 20–21 LT, and the post-sunset peak has a much higher value than the post-noon one. (2) Both for the north and south crests, the TEC-CTR at 120°E longitude showed a semi-annual variation with maximum peak values occurring in the equinoctial months. (3) TEC-CTR for the north crest has lower values in summer than in winter, whereas TEC-CTR for the south crest does not show this ‘winter anomaly’ effect. In other words, TEC-CTR for both the north and south crests has higher values in the northern hemispheric winter than in the northern hemispheric summer. (4) TEC-CTR in the daytime post-noon hours (12–14 LT) does not vary much with the solar activity, however, TEC-CTR in the post-sunset hours (19–21 LT) shows a clear dependence on the solar activity, its values increasing with solar activity.     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

Variation of ionospheric total electron content in Indian low latitude region of the equatorial anomaly during May 2007–April 2008

The ionospheric total electron content (TEC), derived by analyzing dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Varanasi (geomagnetic latitude 14°, 55′N, geomagnetic longitude 154°E) is studied. Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the solar minimum period from May 2007 to April 2008. It is found that the daily maximum TEC near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semiannual variation is seen with two maxima occurring in both spring and autumn. Statistical studies indicate that the variation of EIA crest in TEC is poorly correlated with Dst-index (r = −0.03) but correlated well with Kp-index (r = 0.82). The EIA crest in TEC is found to be more developed around 12:30 LT.     

Seasonal and diurnal behavior of NO2 column density in Antarctica as observed by Mark IV,Brewer Ozone Spectrophotometer Number 153

The Meteorology Department of India has been measuring vertical column density of NO₂ at Maitri (70.7°S, 11.7°E), Antarctica since July 1999 using a Mark IV, Brewer Ozone Spectrophotometer. Maitri is situated at the south of the Antarctic circle. An analysis of 6 years of data shows that NO₂ column has seasonal variation with a maximum value during summer. It is also found that during the period when sun does not set, the NO₂ column exhibits a diurnal variation, with a peak around noon and lower values in the morning and afternoon hours. Using a simple steady-state chemical reaction scheme, an attempt has been made to explain these features.     

Annual,seasonal and diurnal variations of integrated water vapor using GPS observations over Hyderabad,a tropical station

This paper presents annual, seasonal and diurnal variations of integrated water vapor (IWV) derived from Global Positioning System (GPS) measurements for a tropical site, Hyderabad (17.4° N, 78.46° E). The zenith wet delay (ZWD) due to the troposphere has been computed using GPS observations and collocated meteorological data. ZWD is converted to IWV with very little added uncertainty. Mean monthly IWV values show maximum in July (~50 kg m⁻²) and minimum in December (~15 kg m⁻²). Fast Fourier Transform (FFT) and Harmonic analyses methods have been adopted to extract amplitudes and phases of diurnal (24 h), semi-diurnal (12 h) and ter-diurnal (8 h) oscillations which yielded comparable results. Amplitude of the 24 h component is observed to be maximum in spring whereas 12 h and 8 h components maximize in summer. A cross-correlation study between available daily IWV values and corresponding surface temperatures over one year produced a good correlation coefficient (0.44). The correlation obtained for different seasons got reduced to 0.25, 0.02, −0.39 and 0.21 for winter, spring, summer and autumn seasons respectively. The correlation between IWV and rainfall is poor. The coefficients obtained for the whole year is 0.05 and −0.13 for the rainy season.     

Assessment of improvement of the IRI model over Ethiopia for the modeling of the variability of TEC during the period 2013–2016

This paper discusses the monthly and seasonal variation of the total electron content (TEC) and the improvement of performance of the IRI model in estimating TEC over Ethiopia during the solar maximum (2013–2016) phase employing as reference the GPS derived TEC data inferred from four GPS receivers installed in different regions of Ethiopia; Assosa (geog 10.05°N, 34.55°E, Geom. 7.01°N), Ambo (8.97°N, 37.86°E, Geom. 5.42°N), Nazret (8.57°N, 39.29°E, Geom. 4.81°N) and Arba Minch (6.06°N, 37.56°E, Geom. 2.62°N). The results reveal that, in the years 2013–2016, the highest peak GPS-derived diurnal VTEC is observed in the March equinox in 2015 over Arba Minch station. Moreover, both the arithmetic mean GPS-derived and modelled VTEC values, generally, show maximum and minimum values in the equinoctial and June solstice months, respectively in 2014–2015. However, in 2013, the minimum and maximum arithmetic mean GPS-derived values are observed in the March equinox and December solstice, respectively. The results also show that, even though overestimation of the modelled VTEC has been observed on most of the hours, all versions of the model are generally good to estimate both the monthly and seasonal diurnal hourly VTEC values, especially in the early morning hours (00:00–03:00?UT or 03:00–06:00?LT). However, it has also been shown that the IRI 2007 and IRI 2012 versions generally perform best in matching the diurnal GPS derived TEC values as compared to that of the IRI 2016 version. In addition, the IRI 2012 version with IRI2001 option for the topside electron density shows the highest overestimation of the VTEC as compared to the other options. None of the versions of the IRI model are proved to be able to capture the effects of geomagnetic storms.     

1979—1980年中国地区地磁日变场幅度的空间分布及其逐日变化

对静日变化资料进行滤波和付氏分析,取其合成日变幅,得到的季节平均规律及统计分布更真实地反映了S_q日变化的空间分布特征.资料的计算表明,S_q垂直分量日变幅空间分布的极值位置与水平分量的结果不一致,特别是在冬季.此外,静日变化的周日、半周日谐波的振幅、相位也存在着明显的季节变化,而且夏季周日与半周日的初位相总是偏大.      

Ionospheric behavior of the F2 peak parameters foF2 and hmF2 at Hainan and comparisons with IRI model predictions

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with hourly time interval resolution for the diurnal variation, obtained with DPS-4 digisonde observations at Hainan (19.4°N, 109.0°E) are used to study the low latitude ionospheric variation behavior. The observational results are compared with the International Reference Ionospheric Model (IRI) predictions. The time period coverage of the data used for the present study is from March 2002 to February 2005. Our present study showed that: (1) In general, IRI predictions using CCIR and URSI coefficients follow well the diurnal and seasonal variation patterns of the experimental values of foF2. However, CCIR foF2 and URSI foF2 IRI predictions systematically underestimate the observed results during most time period of the day, with the percentage difference ΔfoF2 (%) values changing between about −5% and −25%, whereas for a few hours around pre-sunrise, the IRI predictions generally overestimate the observational ones with ΔfoF2 (%) sometimes reaching as large as ∼30%. The agreement between the IRI results and the observational ones is better for the year 2002 than for the other years. The best agreement between the IRI results and the observational ones is obtained in summer when using URSI coefficients, with the seasonal average values of ΔfoF2 (%) being within the limits of ±10%. (2) In general, the IRI predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the observational results. However, when using the measured M(3000)F2 as input, the diurnal variation pattern of hmF2 given by IRI2001 has a much better agreement with the observational one with the detailed fine structures including the pre-sunrise and post-sunset peaks reproduced reasonably well. The agreement between the IRI predicted hmF2 values using CCIR M(30,000)F2 option and the observational ones is worst for the afternoon to post-midnight hours for the high solar activity year 2002. During daytime hours the agreement between the hmF2 values obtained with CCIR M(30,000)F2 option and the observational ones is best for summer season. The discrepancy between the observational hmF2 and that obtained with CCIR M(30,000)F2 option stem from the CCIR M(3000)F2 model, which does not produce the small scale structures observed in the measured M(3000)F2.     

Stability of the seasonal variations in diurnal and semidiurnal components of mid-latitude F2 layer parameters

We report work utilizing 15-min resolution ionospheric data obtained with DPS-4 digisonde in 2003–2011 to study the seasonal variations in amplitudes and phases of the most powerful spectral components of the F2 layer critical frequency (foF2) and peak height (hmF2) fluctuations over Irkutsk (52.5°N, 104.0°E). We show that fluctuations of both parameters contain quasi-harmonic components with periods of T_n = 24/n h (n = 1–7). The number of distinct spectral peaks varies from 3 in summer to 7 in winter. Amplitude and phase characteristics of the diurnal (n = 1) and semidiurnal (n = 2) components is studied using the data sets extracted from the original data sets with band-pass filter. It has been found that the amplitudes of diurnal/semidiurnal foF2 and diurnal hmF2 components are maximum in winter and minimum in summer. Amplitudes of the diurnal components vary gradually; those of the foF2 semidiurnal one, abruptly, thus forming a narrow winter maximum in November–January. The phase (local time of maximum) of the diurnal foF2 component increases gradually by 4–6 h from winter to summer. The phase of the semidiurnal foF2 component is nearly stable in winter/summer and sharply decreases (increases) by 2–3 h near the spring (autumn) equinox. The phase of the diurnal component of hmF2 (local time of minimum) varies slightly between 1130 and 1300 LT; that of the semidiurnal one decreases (increases) by 4–6 h from January to March (from September to November). The results obtained show that the main features of seasonal variations in the diurnal and semidiurnal components of the mid-latitude F2 layer parameters recur consistently during the solar activity growth and decline phases.