海南地区电离层Es特性研究

利用中国科学院空间科学与应用研究中心海南电离层观测站DPS-4电离层测高仪2002年3月至2005年2月的观测数据,对太阳活动下降期间海南地区上空发生的偶发E层(Es)的逐年变化和季节变化进行了统计分析.主要结果表明,在这3年中海南地区Es事件出现次数逐年增加;Esr,Esq出现最频繁,其次是Esr,而Esc和Esf出现次数较少,出现次数最少的是Esh;Es事件在夏季出现最为频繁,秋季和冬季次之,春季最为稀少;各季节不同类型Es事件的出现率也不同,其中Esl和Esq一般在各季节占主导地位,Esr其次,Esc,Esf和Esh出现次数较少,尤其是Esf在秋季几乎没有出现.这些结果对于进一步研究Es产生的物理机制提供了探测基础,同时对于电离层空间天气预报模型的建立也具有重要的意义.      

武汉地区电离层TEC和NmF2及板厚的季节变化

通过利用武汉电离层观测站(114.4°E,30.6°N)1980-1990年对E8T-Ⅱ卫星信标的法拉第旋转测量的TEC(电子浓度总含量)数据,以及由测高仪测量的1980-1990年间的f0F2(F2层临界频率)数据,研究了武汉地区TEC,NmF2(最大电子浓度)和板厚的季节变化,同时比较了IRI和武汉单站模式在预测NmF2季节性方面的有效性.武汉单站模式在预测NmF2季节性变化方面优于IRI模式.      

Estimating sea level rise around Australia using a new approach to account for low frequency climate signals

Regional sea level studies help to identify the vulnerable areas to the sea level rise. This study investigates the impact of climate modes on sea level variations and trends around Australia using altimetry data, climate indices, and sea level records from tide gauge stations. Here, we show that the sea level variations are negatively correlated with climate indices to the north and west of Australia. The spectral analyses of the climate indices and tide gauge data suggest that a low frequency signal with a period of 11 years emerges during the mid 1980s. Since the 25-year length of the satellite altimetry record is yet too short to detect low frequency signals, their effect on the estimation of regional sea level trend is unknown. Therefore, we estimate the sea level trend with consideration of this signal and using a two-step method. All signals with periods shorter than 7.5 years are first removed from sea level time series and then the trend is estimated using the parametric model that includes the 11-year signal. The skill of the parametric model in explaining the variations in sea level anomaly validates the presence of the 11-year signal detected in the spectrograms of the tide gauge data and climate indices. The average sea level trend for the study area is estimated as 3.85 ± 0.15 mm/year.     

Tide variation monitoring based improved GNSS-MR by empirical mode decomposition

Global Navigation Satellite System multipath reflectometry (GNSS-MR) technology has great potential for monitoring tide level changes. GNSS-MR techniques usually extract signal-to-noise ratio (SNR) residual sequences using quadratic polynomials; however, such algorithms are affected considerably by satellite elevation angles. To improve the stability and accuracy of an SNR residual sequence, this study proposed an SNR signal decomposition method based on empirical mode decomposition (EMD). First, the SNR signal is decomposed by EMD, following which the SNR residual sequence is obtained by combining the corresponding intrinsic mode function with the frequency range of the coherent signal. Second, the Lomb–Scargle spectrum is analyzed to obtain the frequency of the SNR residual sequence. Finally, the SNR frequency is converted into the tide height. To verify the validity of the SNR residual sequence obtained by EMD, the algorithm performance was assessed using multigroup satellite elevation angle intervals with measured data from two station, SC02 in the United States and RSBY in Australia. Experimental results demonstrated that the accuracy of the improved algorithm was improved in the low-elevation range. The improved algorithm demonstrated high monitoring accuracy, and the effective number was not less than 80% of the total in SC02, which means it could effectively describe the trend of the tide with accuracy of approximately 10 cm, meanwhile, the RMS error of RSBY could be reduced by 30 cm, to the maximum extent. The EMD method effectively expands the range of available GNSS-MR elevations, avoids the loss of effective information, enhances considerably the utilization rate of GNSS data, and improves the accuracy of GNSS-MR tide level monitoring.     

Effects of sector structure of the interplanetary magnetic field on the upper mesosphere–lower thermosphere dynamics

In this paper we study the influence of the interplanetary magnetic field (IMF) polarity changes caused by the Earth passing through the IMF sector boundary on the dynamic processes taking place in neutral atmosphere within the altitude interval of the upper mesosphere–lower thermosphere (83–101 km). The analysis has revealed the influence of the IMF sector structure on dynamics of the upper mesosphere–lower thermosphere. There has been a significant seasonal variation of the wind reaction to the IMF polarity changes observed. The influence of the IMF polarity changes on neutral atmosphere dynamics within the altitude range of 83–101 km is most pronounced in the zonal component of neural wind when the IMF polarity changes from negative to positive in all the seasons except for spring and when IMF polarity changes from positive to negative – in spring only.     

Dependence of the 27-day variation of cosmic rays on the global magnetic field of the Sun

We show that the higher range of the heliolongitudinal asymmetry of the solar wind speed in the positive polarity period (A > 0) than in the negative polarity period (A < 0) is one of the important reasons of the larger amplitudes of the 27-day variation of the galactic cosmic ray (GCR) intensity in the period of 1995–1997 (A > 0) than in 1985–1987 (A < 0). Subsequently, different ranges of the heliolongitudinal asymmetry of the solar wind speed jointly with equally important corresponding drift effect are general causes of the polarity dependence of the amplitudes of the 27-day variation of the GCR intensity. At the same time, we show that the polarity dependence is feeble for the last unusual minimum epoch of solar activity 2007–2009 (A < 0); the amplitude of the 27-day variation of the GCR intensity shows only a tendency of the polarity dependence. We present a three dimensional (3-D) model of the 27-day variation of GCR based on the Parker’s transport equation. In the 3-D model is implemented a longitudinal variation of the solar wind speed reproducing in situ measurements and corresponding divergence-free interplanetary magnetic field (IMF) derived from the Maxwell’s equations. We show that results of the proposed 3-D modeling of the 27-day variation of GCR intensity for different polarities of the solar magnetic cycle are in good agreement with the neutron monitors experimental data. To reach a compatibility of the theoretical modeling with observations for the last minimum epoch of solar activity 2007–2009 (A < 0) a parallel diffusion coefficient was increased by ∼40%.     

Study of the 27-day variations of the galactic cosmic ray intensity and anisotropy

We demonstrate that the general features of the radial and azimuthal components of the anisotropy of galactic cosmic rays can be studied by the harmonic analysis method using data from an individual neutron monitor with cut off rigidity <5 GV. In particular, we study the characteristics of the 27-day (solar rotation period) variations of the galactic cosmic ray intensity and anisotropy, solar wind velocity, interplanetary magnetic field strength and sunspot number. The amplitudes of the 27-day variations of the galactic cosmic ray anisotropy are greater, and the phases more clearly established, in A > 0 polarity periods than in A < 0 polarity periods at times of minimum solar activity. The phases of the 27-day variations of the galactic cosmic rays intensity and anisotropy are opposite with respect to the similar changes of the solar wind velocity in A > 0 polarity periods. No significant dependence of the amplitude of the 27-day variation of the galactic cosmic ray anisotropy on the tilt angle of the heliospheric neutral sheet is found. Daily epicyclegrams obtained by Chree’s method show that the 27-day variations of the galactic cosmic ray anisotropy during A > 0 polarity periods follow elliptical paths with the major axes oriented approximately along the interplanetary magnetic field. The paths are more irregular during A < 0 polarity periods.     

Characteristics of the seasonal variation of the global tropopause revealed by COSMIC/GPS data

Using the Global Navigation Satellite System (GNSS) radio occultation observations from Formosa Satellite mission-3/Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) from 2007 to 2012, the climatological characteristics of the global tropopause was studied, with the following features identified. The overall results generally agree with previous studies. The tropopause has an obvious zonal structure, with more zonal characteristics in the Southern Hemisphere than the Northern Hemisphere. The vertical shape of the tropopause is sharp in the tropics and broad in the sub-tropical latitudes, with the sharpest latitudinal gradient in the mid-latitudes of both hemispheres. The global tropopause exists in a large range between 8 km and 17 km (or between 100 hPa and 340 hPa). The highest tropopause is over the South Asian monsoon regions for the entire year. The spatial structure of the tropopause in the polar region is of concentric structure, with an altitude between 7.5 km and 10 km. It is more symmetric in the Antarctic than the Arctic. Differing from other places, the height of the tropopause in the Antarctic is higher in winter as opposed to summer. The tropopause has distinct seasonal variability, especially in polar regions.     

Variations of geomagnetic cutoff rigidity in the southern hemisphere close to 70°W (South-Atlantic Anomaly and Antarctic zones) in the period 1975–2010

We report the existence of rapid variations in (effective) geomagnetic cutoff rigidity (Rc) between the equatorial and Antarctic zones adjacent to the Andes Mountains, revealed by the variation rate of geomagnetic cutoff rigidity (VRc) in the period 1975–2010. Our analysis is based on empirical records and theoretical models of the variations in cosmic rays and on the structure of geomagnetic fields. These have given us a different view of variations in Rc in time and space along the 70°W meridian, where secular variations in the geomagnetic field are strongly influenced by the proximity of the South Atlantic Magnetic Anomaly (SAMA), one of the most important characteristics of the terrestrial magnetic field that affects our planet, close from the equator to the 50°S parallel and from South America to South Africa. The VRc presents rapid changes in mid-latitudes where SAMA exerts its influence despite the existence of smooth changes in the geomagnetic field. This shows that these changes occur mainly in the spatial configuration, rather than in the temporal evolution of Rc. The analysis was performed using measurements from the Chilean Network of Cosmic Rays and Geomagnetism Observatories, equipped with BF-3 and latest generation He-3 neutron monitors, Fluxgate magnetometers, geomagnetic reference field (IGRF) and Tsyganenko 2001 model (just for completeness).     

基于流场/声爆耦合伴随方程的超声速公务机声爆优化

基于自主研发的大规模并行结构化网格CFD求解器PMB3D以及并行化伴随方程求解器PADJ3D，开展了流场/声爆伴随方程的求解研究。首先采用标准算例，对内部CFD代码PMB3D软件和声爆预测代码进行了声爆计算可信度验证，以及声爆强度对近场声压梯度的校核。针对并行环境下多块对接网格的近场声压提取操作的复杂性，提出了"包围盒"的方法实现并行环境下近场声压装配单元编号、网格块编号以及对应的进程编号确定，基于声爆计算坐标将并行传递的数据进行一维排序，为声爆预测、伴随方程以及梯度求解提供输入条件。通过线性插值雅克比矩阵实现均匀坐标系梯度信息向非均匀坐标转换，并进一步根据结构化网格特征提出了插值原则，简化了近场声压转换雅克比矩阵的变分。通过装配单元记录，实现声爆强度对流场守恒变量的变分结果向各个进程装配，将装配结果作为流场伴随方程的右端项实现流场声爆耦合伴随方程的求解。此外，对小型超声速公务机开展了声爆优化，对比分析了设计前后的声压及其频谱特性。