Ionospheric precursor of a destructive earthquake that occurred on April 6, 2009 at L’Aquila (Italy)

On the basis of the 15-min data from a series of ground-based vertical ionospheric sounding stations, a study of variations of the foF2 critical frequency before the strong earthquake (M = 6.3) that occurred on April 6, 2009 at L’Aquila (Italy) was carried out. The earthquake epicenter was located 85 km north-eastward from Rome. Approximately 20 h prior to the earthquake, a well-pronounced statistically significant effect of foF2 increase relative to the average background for magnetically quiet days was observed for almost 1.5 h at the Rome ionospheric station. In this case, at control stations distanced from the earthquake epicenter, no statistically significant deviations of foF2 from the background values were detected during the same observations period. This fact provides grounds for consideration of the foF2 increase observed at Rome station as a possible ionospheric precursor of this earthquake.     

Modelling F2-layer seasonal trends and day-to-day variability driven by coupling with the lower atmosphere

This paper presents results from the TIME-GCM-CCM3 thermosphere–ionosphere–lower atmosphere flux-coupled model, and investigates how well the model simulates known F2-layer day/night and seasonal behaviour and patterns of day-to-day variability at seven ionosonde stations. Of the many possible contributors to F2-layer variability, the present work includes only the influence of ‘meteorological’ disturbances transmitted from lower levels in the atmosphere, solar and geomagnetic conditions being held at constant levels throughout a model year.In comparison to ionosonde data, TIME-GCM-CCM3 models the peak electron density (NmF2) quite well, except for overemphasizing the daytime summer/winter anomaly in both hemispheres and seriously underestimating night NmF2 in summer. The peak height hmF2 is satisfactorily modelled by day, except that the model does not reproduce its observed semiannual variation. Nighttime values of hmF2 are much too low, thus causing low model values of night NmF2. Comparison of the variations of NmF2 and the neutral [O/N₂] ratio supports the idea that both annual and semiannual variations of F2-layer electron density are largely caused by changes of neutral composition, which in turn are driven by the global thermospheric circulation.Finally, the paper describes and discusses the characteristics of the F2-layer response to the imposed ‘meteorological’ disturbances. The ionospheric response is evaluated as the standard deviations of five ionospheric parameters for each station within 11-day blocks of data. At any one station, the patterns of variability show some coherence between different parameters, such as peak electron density and the neutral atomic/molecular ratio. Coherence between stations is found only between the closest pairs, some 2500 km apart, which is presumably related to the scale size of the ‘meteorological’ disturbances. The F2-layer day-to-day variability appears to be related more to variations in winds than to variations of thermospheric composition.     

中低纬地区电离层对CIR和CME响应的统计分析

本文利用中低纬日本地区(131°E,35°N)GPS-TEC格点化数据,分析了2001—2009年间109个共转相互作用区(CIR)事件、45个日冕物质抛射(CME)事件引起的地磁扰动期间电离层的响应.结果表明,电离层暴的类型随太阳活动的变化而有不同的变化,CIR事件引发的电离层正相暴、正负双相暴多发生在太阳活动下降年,负相暴多发生在高年,负正双相暴多发生在低年;CME事件引发的电离层正相暴和负相暴多发生在高年.CIR和CME引发的不同类型的电离层暴的季节性差异不大,在夏季多发生正负双相暴.电离层暴发生时间相对地磁暴的时延大部分在-6~6h之间,但CIR引发的电离层暴时延范围更广,在-12~24h之间,而CME引发的电离层暴时延主要在-6~6h之间.中低纬的电离层暴多发生在主相阶段,其中CIR引发的双相暴也会发生在初相阶段.电离层负暴多发生在AE最大值为800~1200nT之间.CIR引起的电离层扰动持续时间较长,一般在1~6天左右,而CME引起的电离层扰动持续时间一般在1~4天左右.     

赤道地区电离层foF2在第23/24太阳活动周极低年期间创造了极低纪录?

太阳辐射是电离层的电离源,强烈地调制电离层的变化.探索不同太阳辐射水平下的电离层状态,有助于认识电离层演变及其内在的基本物理过程.太阳活动在2008—2009年处于有记录以来的极低水平,研究电离层在此期间的变化及与其它太阳活动低年的差异是一个有益的课题.本文利用位于美洲扇区磁赤道地区Jicamarca台站(12.0°S, 283.2°E; dip 0.28°)测高仪观测的电离层F2层临界频率foF₂数据探讨赤道地区foF₂的行为.分别对第22/23太阳活动周低年(1996—1997)和第23/24活动周低年(2008—2009)的月中值、季节中值和滑动年均值进行分析,确认相比上一个太阳低年而言,在2008—2009年foF₂滑动年均值和不同季节中值在各个地方时均降低,而月中值存在降低和升高.对foF₂的时间尺度特性的分析发现,在本太阳周低年foF₂长时间尺度分量下降,而短时间尺度分量呈现不一致的变化.我们认为,现有文献报道给出2008—2009年与以往太阳活动低年对比结果不一致有可能归因于所用分析方法关注的时间尺度不相同.     

Correlation analysis of vertical total electron content

It is necessary to model and analyze the ionospheric effects due to a direct relationship between Global Positioning System (GPS) applications and changes in the ionosphere. In order to monitor these changes, the ionosphere can be represented by the vertical total electron content (VTEC) which can be used to analyze ionospheric conditions from a variety of stations. In this study, 21 stations were used to carry out analysis and estimation of VTEC. Three days during a geomagnetic storm, namely, 7, 8, and 9 January 2005, are chosen for investigation. In addition, the de-correlation time of the VTEC was estimated to define ionospheric variations in time using autocorrelation analysis. The de-correlation time of the ionosphere is based on correlation times estimated by using autocorrelation functions. From the high-latitude stations, the mean of the correlation times decreased from 8 to 6 epochs during a storm. In this time period, it was found from the station results that the ionosphere was more affected at the high-latitude than at the mid-latitude region.     

Derivation and test of ionospheric activity indices from real-time ionosonde observations in the European region

New ionospheric activity indices are derived from automatically scaled online data from several European ionosonde stations. These indices are used to distinguish between normal ionospheric conditions expected from prevailing solar activity and ionospheric disturbances caused by specific solar and atmospheric events (flares, coronal mass ejections, atmospheric waves, etc.). The most reliable indices are derived from the maximum electron density of the ionospheric F 2-layer expressed by the maximum critical frequency foF 2. Similar indices derived from ionospheric M(3000)F 2 values show a markedly lower variability indicating that the changes of the altitude of the F 2-layer maximum are proportionally smaller than those estimated from the maximum electron density in the F 2-layer. By using the ionospheric activity indices for several stations the ionospheric disturbance level over a substantial part of Europe (34°N–60°N; 5°W–40°E) can now be displayed online.     

Hour-to-hour variability of the ionosphere: An application of the classical multidimensional scaling method

The classical multidimensional scaling(MDS) method is introduced and applied in the study of the hour-to-hour ionospheric variability based on the ionospheric fo F2 observed at three ionosonde stations in East-Asia in 2002 and 2007. Results from the matrix eigen decompositions indicate that the annual part of the ionospheric variation is large in middle latitude and solar maximum period(2002) while low in the low latitude and solar minimum period(2007). The connectivity maps of the hour-to-hour ionospheric variability based on MDS method show some common diurnal features. The ionospheric connectivity between adjacent hours near noon hours and near midnight hours is high. The ionospheric connectivity between adjacent hours near sunrise hours and near sunset hours is poor, especially for the sunrise hours. Also there are latitudinal and solar activity dependences in this kind of connectivity. These results revealed from the ionospheric connectivity maps are useful physically and in practice for the ionospheric forecasting on the hour-to-hour scale.     

Spatial correlation of foF2 and vTEC under quiet and disturbed ionospheric conditions: A case study

The spatial variations of the ionospheric F2-layer vertical incidence critical frequency (foF2) and GPS-derived vertical total electron content (vTEC) under geomagnetically quiet and disturbed days are examined using measurements from the latitudinal and longitudinal chains of ionospheric stations and GPS receivers over the European area. Plots produced for January 2005 are used to discuss temporal structures in terms of the prevailing solar-terrestrial conditions. Then the line trends procedure has been applied to simultaneous data collected from a limited number of measuring stations during quiet monthly median ionospheric conditions as well as during the storm period of 16–23 January 2005. The procedure is explained involving an application of the least squares method to define latitudinal and longitudinal dependence of foF2 and vTEC at different locations. Examples of coefficients of determination thereby produced show that the linear regression equations are very helpful in predicting longitudinal and latitudinal vTEC and foF2 variations during the quiet as well as disturbed ionospheric conditions.     

基于陆态网络GPS数据的电离层空间天气监测与研究

中国大陆构造环境监测网络(简称陆态网络)是以全球卫星导航定位系统(GNSS)为主,辅以多种空间观测技术,实时动态监测大陆构造环境变化,探求其对资源、环境和灾害的影响的地球科学综合观测网络.基于陆态网络约200个基准站的GPS观测数据,本文探讨了其在电离层空间天气监测与研究方面的应用.包括磁暴期间电离层暴扰动形态,大尺度电离层行进式扰动,太阳耀斑引起的电离层骚扰和低纬电离层不规则体结构等.研究结果表明:陆态网络布局合理,观测数据质量良好,完全可用于中国及周边地区电离层空间天气监测与研究,为进一步开展我国电离层空间天气预警和预报奠定了观测基础.     

Application of predictive regional ionosphere model to medium range RTK positioning

Real Time Kinematic (RTK) GPS positioning over longer distances requires a support of atmospheric (ionospheric and tropospheric) corrections, since the atmospheric errors decorrelate with the growing distances and cannot be completely eliminated by double differencing of the satellite observations. Currently, the most commonly used approach is to derive the atmospheric corrections at the reference station network and provide them in real time to the roving receiver. Another solution, proposed here, is to use predictive atmospheric models in order to derive the atmospheric corrections. This paper presents the test results of the performance assessment of the predictive ionosphere model (UWM-IPM) application to medium-range RTK positioning. The rover data collected within 25 to 67 km from the closest reference station were processed in the kinematic mode with the support of the ionospheric corrections derived from the UWM-IPM model. The RTK solution was derived in both single-and multi-baseline modes, and compared to the two reference solutions obtained without the ionospheric corrections. All numerical tests were carried out using the MPGPS software developed in cooperation with The Ohio State University; a recent extension to the software, developed at the University of Warmia and Mazury in Olsztyn, introduces the predictive ionosphere model to the RTK solution. The test results are very promising, and indicate that predicted ionosphere corrections can effectively support medium-range RTK positioning, and allow for fast ambiguity resolution over distances of several tens of kilometers under moderate ionospheric conditions.     

Ionospheric behavior over Europe during the solar eclipse of 3 October 2005

An annular eclipse occurred over Europe in the morning hours of 3 October 2005. The well-defined obscuration function of the solar radiation during the eclipse provided a good opportunity to study the ionospheric/thermospheric response to solar radiation changes. Since the peak electron density behavior of the ionospheric F2 layer follows the local balance of plasma production, loss and transport, the ionospheric plasma redistribution processes significantly affect the shape of the electron density profile. These processes are discussed here based on a comparison of vertical incidence sounding (VS) and vertical total electron content (TEC) data above-selected ionosonde stations in Europe. The equivalent slab thickness, derived with a time resolution of 10 min, provides relatively good information on the variation of the electron density profile during the eclipse. The computations reveal an increased width of the ionosphere around the maximum phase. As indicated by the available measurements over Spain, the photo production is significantly reduced during the event leading to a slower increase of the total ionization in comparison with the neighboring days. The supersonic motion of the Moon's cool shadow through the atmosphere may generate atmospheric gravity waves that propagate upward and are detectable as traveling ionospheric disturbances at ionospheric heights. High-frequency (HF) Doppler shift spectrograms were recorded during the eclipse showing a distinct disturbance along the eclipse path. Whereas the ionosonde measurements at the Ebro station/Spain in the vicinity of the eclipse path reveal the origin of the wave activity in the lower thermosphere below about 180 km altitude, the similar observations at Pruhonice/Czech Republic provide arguments to localize the origin of the abnormal waves in the middle atmosphere well below the ionospheric heights. Although ionosonde and HF Doppler measurements show enhanced wave activity, the TEC data analysis does not, which is an indication that the wave amplitudes are too small for detecting them via this interpolation method. The total ionization reduces up to about 30% during the event. A comparison with similar observations from the solar eclipse of 11 August 1999 revealed a quite different ionospheric behavior at different latitudes, a fact that needs further investigation.     

欧洲地区电离层峰值电子密度逐日变化的相关距离研究

本文采用欧洲22个台站的电离层F^₂层峰值电子密度N_mF^₂,分析了其逐日变化分量的相关距离S,着重研究了S的周日变化、季节变化及其随太阳活动和地磁活动的变化.首先用指数型函数模式来拟合任意两站间电离层逐日变化的相关系数R随间距d的变化,由此估算出逐日变化的相关距离S.详细研究了S在不同的季节（春季、夏季、秋季和冬季）,不同的地磁活动（平静和扰动）及不同的太阳活动（低、中和高）随世界时的变化（周日变化）.结果表明：（１）S的范围一般为400～1600 km;（２）S值在白天比夜间大;（３）S值具有季节变化,夏季最大,冬季最小,春秋季差异不大;（４）S值在地磁扰动时比平静时大;（５）当太阳活动低时,S值在日落到正午间要比太阳活动中或高时明显偏小,而在正午到日落间则与太阳活动中或高时差异不大.根据以上结果,我们认为：（１） 太阳辐射对电离层逐日变化的影响是大尺度的,并在白天和太阳活动高时大于晚上和太阳活动低时;（２） 地磁活动的影响也是大尺度的;（３）气象活动的影响是相对小尺度的,且逐日变化具有季节性.本文从相关尺度分析的角度,证实了电离层逐日变化来源于太阳辐射、地磁活动和气象活动因素的论断.     

Studies of the ionospheric D-region using partial reflections in spring 2004 at middle and high latitudes

Results of the observations of the ionospheric effects of two solar flares in April 2004 performed using partial reflections are presented. The studies were carried out at the measuring facilities located in different latitudinal regions: at Vasil’sursk station in the Nizhni Novgorod region and at Tumannyi station in the Murmansk region. The quantitative estimates of the electron density in the polar and midlatitude D region under quiet conditions and during solar flares were obtained. The correlation between rapid variations in electron concentration at heights of about 80 km at these stations was found and it was shown that during solar flares the electron density at heights of 60–70 km corresponds to the intensity of the X-ray flux in the range of 0.5–3 Å, which points to the action of the linear law of recombination in the ionospheric D region.     

GPS phase scintillation observed over a high-latitude Antarctic station during solar minimum

Radio wave scintillation due to the presence of ionospheric irregularities can cause fading and phase variation of L-band navigation signals such as those used by the Global Positioning System (GPS). The high-latitude regions experience increased scintillation events under geomagnetically disturbed conditions, particularly during solar maximum. This paper presents phase scintillation measurements observed at the South African Antarctic polar research station during solar minimum conditions. In addition, a multi-instrument approach is shown in which the scintillation events are associated with auroral electron precipitation. Furthermore, it is shown that external energy injection into the ionosphere is an important factor in the development of irregularities producing scintillation.     

Comparison of variations in the <Emphasis Type="Italic">F</Emphasis>2 layer critical frequency before the strong Italian earthquakes and during the magnetic storm: Data of the stations for ground-based vertical sounding of the ionosphere

The possible seismic ionospheric effects before two strong earthquakes in Italy, which took place near Rome station (January 7, 1962, M = 6.0; April 5, 1998, M = 5.3), have been studied using the data of several European stations for ground-based vertical sounding of the ionosphere. An analysis of the behavior of the foF2 critical frequencies for the earthquake that occurred on January 7, 1962, made it possible to simultaneously consider the effects, originating during earthquake preparation under quiet geomagnetic conditions and caused by the magnetic strom that began on January 10, 1962, using a unified technique. It has been indicated that the distinguished precursor effects of two considered earthquakes with magnitudes not more than 6.0 are localized in the region with an epicentral distance of about 400 km, whereas the disturbances caused by the magnetic storm are of a planetary character. The amplitude of foF2 seismic ionospheric disturbances is much smaller than the maximal variations caused by the magnetic storm; however, the absolute value of this amplitude is a factor of 1.5 as large as the standard deivation from the monthly median.     

基于北斗、GLONASS和GPS系统的中低纬电离层特性联合探测

基于GNSS(Global Navigation Satellite Systems)的发展,我们利用具有北斗、GLONASS和GPS三系统信号接收功能的接收机观测的数据,结合电离层总电子含量(Total Electron Content, TEC)的反演算法,提取出GNSS三系统观测的电离层TEC;同时,将GNSS三系统获取的TEC应用到电离层TEC地图、行进式扰动、不规则体结构和电离层的太阳耀斑响应等方面的研究中,这也是首次使用三种GNSS系统数据对电离层进行联合探测研究.研究结果表明,增加了北斗系统的GNSS三系统在研究中国地区电离层TEC地图、周日变化、逐日变化,行进式扰动以及电离层的实时监测等方面较单系统的GPS具有明显的优势.     

Variability of parameters of the <Emphasis Type="Italic">F</Emphasis>2-layer maximum in the quiet midlatitude ionosphere under low solar activity: 1. Statistical properties

Results of statistical analysis of the properties of variability of F2-layer maximum parameters (critical frequency foF2 and the height hmF2) in quiet midlatitude ionosphere under low solar activity in the daytime (1000–1500 LT) and nighttime (2200–0300 LT) hours are presented on the basis of Irkutsk station data for 2007–2008. It is found that the distribution density of δfoF2 could be presented as consisting of two distinctly different normal laws of this distribution, one of which corresponds to weak (|δfoF2| < 10%) fluctuations in foF2 and the other corresponds to strong (30% > |δfoF2| > 10%) fluctuations. Weak fluctuations in foF2 to a substantial degree are related to ionospheric variability at times less of than 1–3 h and determine the δfoF2 variability in the daytime hours. Strong fluctuations in foF2 are mainly related to day-to-day variability of the ionosphere at a fixed local time, the variability increasing by approximately a factor of 3 during the transition from day to night and determining the δfoF2 variability in the nighttime hours. The distribution density of ΔhmF2 is close to the normal distribution law. An interpretation of the different character of the distribution densities of δfoF2 and ΔhmF2 is given.     

Variability of the ionosphere

Hourly foF2 data from over 100 ionosonde stations during 1967–89 are examined to quantify F-region ionospheric variability, and to assess to what degree the observed variability may be attributed to various sources, i.e., solar ionizing flux, meteorological influences, and changing solar wind conditions. Our findings are as follows. Under quiet geomagnetic conditions (K_p<1), the 1-σ (σ is the standard deviation) variability of N_max about the mean is approx. ±25–35% at ‘high frequencies’ (periods of a few hours to 1–2 days) and approx. ±15–20% at ‘low frequencies’ (periods approx. 2–30 days), at all latitudes. These values provide a reasonable average estimate of ionospheric variability mainly due to “meteorological influences” at these frequencies. Changes in N_max due to variations in solar photon flux, are, on the average, small in comparison at these frequencies. Under quiet conditions for high-frequency oscillations, N_max is most variable at anomaly peak latitudes. This may reflect the sensitivity of anomaly peak densities to day-to-day variations in F-region winds and electric fields driven by the E-region wind dynamo. Ionospheric variability increases with magnetic activity at all latitudes and for both low and high frequency ranges, and the slopes of all curves increase with latitude. Thus, the responsiveness of the ionosphere to increased magnetic activity increases as one progresses from lower to higher latitudes. For the 25% most disturbed conditions (K_p>4), the average 1-σ variability of N_max about the mean ranges from approx. ±35% (equator) to approx. ±45% (anomaly peak) to approx. ±55% (high-latitudes) for high frequencies, and from approx. ±25% (equator) to approx. ±45% (high-latitudes) at low frequencies. Some estimates are also provided on N_max variability connected with annual, semiannual and 11-year solar cycle variations.     

Temporal variability of the physical and chemical water characteristics at a coastal monitoring observatory: Station ENSENADA

The temporal variability of the physical and chemical conditions of coastal waters off Ensenada, Baja California (Mexico) was characterized. A historical analysis was made based on 11 years (1998–2008) of temperature and salinity data records measured quarterly by IMECOCAL, along a transect perpendicular to the coast (CalCOFI line 100). Moreover, the physical and chemical conditions at a coastal monitoring observatory called station ENSENADA were described using a 2-year data series (October 2006–November 2008) obtained with improved temporal resolution. The historical analysis of line 100 showed marked seasonal variability in the thermohaline conditions associated with fluctuations in the flow of the equatorward California Current and the poleward California Undercurrent, as well as with coastal upwelling events whose magnitude and frequency increase towards spring–summer. Interannual variability was also observed, related to warm and/or cold ENSO phases that modify the characteristics of the water column in this coastal region. The most striking characteristics of the interannual variability at station ENSENADA were La Niña conditions recorded from summer 2007 to mid 2008. During this cold ENSO phase, temperature, salinity, dissolved oxygen, density, and dissolved inorganic carbon data revealed the anomalous presence of subsurface water at the surface layers in spring 2008. Results suggest that the coastal observatory is sensitive to the temporal variability of hydrographic conditions on shelf coastal waters (<50 km) off Ensenada in the northern BC region. Consequently, station ENSENADA would be a good location to high-frequency monitors the oceanographic conditions of the transitional region between tropical/subtropical and subarctic systems of the California Current System.     

DIAS effective sunspot number as an indicator of the ionospheric activity level over Europe

DIAS (European Digital Upper Atmosphere Server) effective sunspot number — R12eff was recently introduced as a proxy of the ionospheric conditions over Europe for regional ionospheric mapping purposes. Although a pre-processing step for the real-time update of the Simplified Ionospheric Regional Model (SIRM) to real-time conditions, R12eff is available in real time by DIAS system (http://dias.space.noa.gr) for independent operational use. In this paper we discuss the efficiency of R12eff to specify ionospheric conditions over Europe. For this purpose, the diurnal R12eff’s reference pattern was determined on monthly basis and for different solar cycle phases. The deviation of the real-time R12eff estimates from the reference values, ΔR12eff was found to be highly correlated with the foF2 storm-time disturbances, especially during large scale effects indicating that DIAS-R12eff can provide a reliable estimator of the ionospheric activity level over a substantial part of Europe and a powerful tool for ionospheric specification applications.