On the GPS-based ionospheric perturbation after the Tohoku earthquake of March 11, 2011

Based on the data from the GPS receiving networks in Japan and America which have a high time resolution (2 min), two-dimensional (2D) distributions of the variations in the ionospheric total electron content (TEC) are constructed both close to and far from of the epicenter of the submarine earthquake of March 11, 2011 in Japan. Above the epicenter, a diverging multi-period disturbance appears after the main shock due to the acoustic gravity waves. Far from the epicenter, the wave trains associated with the tsunamigenic atmospheric internal gravity waves are revealed. These atmospheric waves significantly advance the arrival of the tsunami signal initially on the Hawaiian islands and then on the western coast of North America. The presence of the tsunami precursor in the form of atmospheric gravity waves is supported by the numerical calculations and by the analysis of the dispersion relation for the waves in the atmosphere. The detected ionospheric responses close and far from the epicenter can be used in the early tsunami warning systems.     

Performance of different TEC models to provide GPS ionospheric corrections

The existence of a worldwide international GPS service (IGS) permanent network of dual-frequency receivers makes the computation of global ionospheric maps (GIMs) of total electron content (TEC) feasible. The GIMs computed by the IGS Associate Analysis Centers on a daily basis and by other kinds of forecast GIMs, which can be computed from, for instance, the international reference ionosphere (IRI) model, and the GPS broadcast models in the navigation message, can be applied to a broad diversity of fields, for instance as, navigation and time transfer.In this context, the performance of different kinds of models are presented in order to determine the accuracy of the different GIM. This is carried out by comparison with the TOPEX data that provides an independent and precise (at the level of few TECU) vertical TEC determination over the oceans and seas. Thus, the obtained accuracies, in terms of global relative error, ranging from 54% corresponding to the GPS broadcast model, to about 41% corresponding to IRI climatological model, and to less than 30% corresponding to GPS data driven models.     

Laurentia crustal motion observed using TOPEX/POSEIDON radar altimetry over land

A new method to estimate the vertical crustal motion from satellite altimetry over land was developed. The method was tested around Hudson Bay, where the observed vertical motion is largely caused by the incomplete glacial isostatic adjustment (GIA) as a result of the Laurentide ice sheet deglaciation since the last glacial maximum (LGM). Decadal (1992–2003) TOPEX/POSEIDON radar altimetry data over land surfaces were used. The results presented here are improved compared to a previous study (Lee, H., Shum, C.K., Kuo, C.Y., Yi, Y., Braun, A., 2008. Application of TOPEX altimetry for solid Earth deformation studies. Terr. Atmos. Ocean. Sci. 19, 37–46. doi:10.3319/TAO.2008.19.1-2.37(SA).) which estimated vertical motion only over relatively flat land surfaces (standard deviation of the height variation <40 cm). In this study, we extended the concept of traditional 1-Hz (one-per-frame) radar altimeter ocean stackfiles to build 10-Hz (10-per-frame) land stackfiles over Hudson Bay land regions, and succeeded in obtaining vertical motion estimates over much rougher surfaces (standard deviation of the height variation <2 m). 90-m C-band Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) is used as a reference surface to select an optimal waveform retracker, to correct surface gradient errors, and to calculate land surface anomalies. Here, we developed an alternative retracker, called the modified threshold retracker, resulting in decadal vertical motion time series over a 1500 km by 1000 km region covering northern Ontario, northeastern Manitoba, and the Great Lakes region which is at the margin of the former Laurentide ice sheet. The average of the estimated uncertainties for the vertical motion is 2.9 mm/year which is comparable to 2.1 mm/year of recent GPS solutions. The estimated vertical motion is compared with other geodetic observations from GPS, tide gauge/altimetry, GRACE, and several GIA models. The data agree best with the laterally varying 3D GIA model, RF3S20 (β = 0.4) whereas the combination of land altimetry solution with other measurements match best with the models RF3S20 (β = 0.0) or RF3S20 (β = 0.2) in terms of mean and standard deviation of the differences. It is anticipated that this innovative technique could potentially be used to provide additional constraints for GIA model improvement, and be applied to other geodynamics studies.     

Laurentia crustal motion observed using TOPEX/POSEIDON radar altimetry over land

A new method to estimate the vertical crustal motion from satellite altimetry over land was developed. The method was tested around Hudson Bay, where the observed vertical motion is largely caused by the incomplete glacial isostatic adjustment (GIA) as a result of the Laurentide ice sheet deglaciation since the last glacial maximum (LGM). Decadal (1992–2003) TOPEX/POSEIDON radar altimetry data over land surfaces were used. The results presented here are improved compared to a previous study (Lee, H., Shum, C.K., Kuo, C.Y., Yi, Y., Braun, A., 2008. Application of TOPEX altimetry for solid Earth deformation studies. Terr. Atmos. Ocean. Sci. 19, 37–46. doi:10.3319/TAO.2008.19.1-2.37(SA).) which estimated vertical motion only over relatively flat land surfaces (standard deviation of the height variation <40 cm). In this study, we extended the concept of traditional 1-Hz (one-per-frame) radar altimeter ocean stackfiles to build 10-Hz (10-per-frame) land stackfiles over Hudson Bay land regions, and succeeded in obtaining vertical motion estimates over much rougher surfaces (standard deviation of the height variation <2 m). 90-m C-band Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) is used as a reference surface to select an optimal waveform retracker, to correct surface gradient errors, and to calculate land surface anomalies. Here, we developed an alternative retracker, called the modified threshold retracker, resulting in decadal vertical motion time series over a 1500 km by 1000 km region covering northern Ontario, northeastern Manitoba, and the Great Lakes region which is at the margin of the former Laurentide ice sheet. The average of the estimated uncertainties for the vertical motion is 2.9 mm/year which is comparable to 2.1 mm/year of recent GPS solutions. The estimated vertical motion is compared with other geodetic observations from GPS, tide gauge/altimetry, GRACE, and several GIA models. The data agree best with the laterally varying 3D GIA model, RF3S20 (β = 0.4) whereas the combination of land altimetry solution with other measurements match best with the models RF3S20 (β = 0.0) or RF3S20 (β = 0.2) in terms of mean and standard deviation of the differences. It is anticipated that this innovative technique could potentially be used to provide additional constraints for GIA model improvement, and be applied to other geodynamics studies.     

First complementary observations by ionospheric tomography, the EISCAT Svalbard radar and the CUTLASS HF radar

Experimental results are presented from ionospheric tomography, the EISCAT Svalbard radar and the CUTLASS HF radar. Tomographic measurements on 10 October 1996, showing a narrow, field-aligned enhancement in electron density in the post-noon sector of the dayside auroral zone, are related to a temporal increase in the plasma concentration observed by the incoherent scatter radar in the region where the HF radar indicated a low velocity sunwards convection. The results demonstrate the complementary nature of these three instruments for polar-cap ionospheric studies.     

The height of the maximum ionospheric electron density over the MU radar

Ionospheric F₂-layer peak height h_mF₂ variations, as measured over 1986–1995 by the MU radar (34.85°N, 136.1°E) and as calculated with a theoretical model, are discussed. The diurnal variations of the measured peak height for different seasons and levels of solar activity are compared with those estimated from ionosonde M3000F₂ and IRI predictions. Also given are the measured ion drift velocities and meridional neutral winds needed to understand the dynamic behavior of the F₂-layer. It is found that: (1) h_mF₂ is generally higher during periods of the solar maximum than during periods of the solar minimum, and higher in summer than in winter; (2) for the solar maximum, h_mF₂ drops markedly in the morning and in the afternoon, while, for the solar minimum, the h_mF₂ minimum occurs in the morning during summer and usually in the afternoon during winter. In general, the measured h_mF₂ is well reproduced by our model when we use the observed drift velocities and plasma temperatures as inputs. Our modeling study shows that the neutral wind contributes strongly to the diurnal variation of h_mF₂ in winter by lowering the ionization layer by day, particularly for the solar maximum; it also helps to enlarge the day–night difference of h_mF₂ in summer. The northward electromagnetic drifts that usually cancel the neutral wind effect have only a minor effect for the location of the MU radar. Other features of the observed h_mF₂ variations, e.g., the solar maximum–minimum difference, the summer–winter difference, and the morning and afternoon drops, are explained by the basic processes of O⁺ production, loss and diffusion, as influenced by the atomic oxygen concentration and neutral and plasma temperatures.     

First incoherent scatter radar observations of ionospheric heating on the second electron gyro-harmonic

We report first results from a unique experiment performed at the HIPAS ionospheric modification facility in conjunction with the Poker Flat incoherent scatter radar in Alaska. High-power radio waves at 2.85 MHz, which corresponds to the second electron gyro-harmonic at ~245 km altitude, were transmitted into the nighttime ionosphere. Clear evidence of F-region ionospheric electron temperature enhancements were found, for the first time at this pump frequency, maximizing when the pump frequency is close to the second gyro-harmonic and double resonance. This is consistent with previous pump-enhanced artificial optical observations. We estimate the plasma heating efficiency to be approximately double that for higher pump frequencies.     

Usage of regularization techniques for effective electron collision frequency determination from oblique ionospheric sounding

Computer modelling is used to investigate the possibility of determining ionospheric parameters from slightly oblique ionospheric soundings, using absorption data for decametric radio waves of different polarization. It is shown that with mean square measurement errors of 0.5 dB, and using regularization algorithms to solve the inverse problems, electron collision frequency profiles can be obtained for the night F-region with errors of less than 30%. Both temperatures of electrons and neutrals are also determined to within 10%.     

Investigation by backscatter radar of artificial irregularities produced in ionospheric plasma heating experiments

The artificial ionospheric turbulence occurs in the ionosphere illuminated by high power HF radio waves. There are a lot of irregularities stretched along the geomagnetic field in this region. The investigation of the artificially disturbed ionospheric region is based on the reception of back scattered signals (BSS) which permits the basic parameters of this region to be estimated and its inhomogeneous structure to be described.Experiments were carried out using ‘Sura’ heating facility in the frequency range of ? = 4.7–9 MHz (ordinary mode) with the effective radiated power P_eff = 50–70 MW beamed vertically upwards. The most important dynamic and statistical BSS characteristics (the built-up time, the relaxation and autocorrelation times, the BSS amplitude spatial correlation function and power spectra) were measured using probe waves in the frequency range of that made it possible to obtain the basic parameters of the artificial irregularities. The model representation of a disturbed region in a form of a periodic structure gives a possibity to evaluate the scale of the structure, the whole size of the disturbance and its power and to calculate the main BSS characteristics.     

重力场恢复中的基于星载GPS的低轨卫星简化动力学定轨方法研究

低轨重力卫星轨道的精确确定是获得精密地球重力场模型的前提, 而精密重力场模型又是获得高精度定轨结果的保证.本文简述了利用卫星重力方法恢复地球重力场及简化动力学方法确定低轨卫星轨道的数学模型,并简单分析和比较现有的几种重力场模型.用CHAMP实测数据,结合现有的重力场模型,系统分析、研究了不同阶次、不同重力场模型对低轨卫星定轨精度的影响;研究了不同间隔的随机速度脉冲在简化动力学方法中对模型误差的吸收、调节作用.计算结果表明,在定轨中,选择合理阶数的、较精确的重力场模型及合理间隔的随机脉冲参数,不但可以提高计算效率,更能提高定轨精度.     

Calibration of hydrological models in ungauged basins based on satellite radar altimetry observations of river water level

This paper proposes a new orientation to address the problem of hydrological model calibration in ungauged basin. Satellite radar altimetric observations of river water level at basin outlet are used to calibrate the model, as a surrogate of streamflow data. To shift the calibration objective, the hydrological model is coupled with a hydraulic model describing the relation between streamflow and water stage. The methodology is illustrated by a case study in the Upper Mississippi Basin using TOPEX/Poseidon (T/P) satellite data. The generalized likelihood uncertainty estimation (GLUE) is employed for model calibration and uncertainty analysis. We found that even without any streamflow information for regulating model behavior, the calibrated hydrological model can make fairly reasonable streamflow estimation. In order to illustrate the degree of additional uncertainty associated with shifting calibration objective and identifying its sources, the posterior distributions of hydrological parameters derived from calibration based on T/P data, streamflow data and T/P data with fixed hydraulic parameters are compared. The results show that the main source is the model parameter uncertainty. And the contribution of remote sensing data uncertainty is minor. Furthermore, the influence of removing high error satellite observations on streamflow estimation is also examined. Under the precondition of sufficient temporal coverage of calibration data, such data screening can eliminate some unrealistic parameter sets from the behavioral group. The study contributes to improve streamflow estimation in ungauged basin and evaluate the value of remote sensing in hydrological modeling. Copyright © 2011 John Wiley & Sons, Ltd.     

基于修正K_p指数的中纬区单站电离层暴预报方法研究

对电离层暴进行有效预报具有很重要的实际应用价值,但目前这一工作还存在着许多方面有待改进.本文假设电离层扰动的季节性规律占主导地位,借鉴Kutiev and Muhtarov构建修正K_p指数表示地磁活动的方法,建立了一个适用于中纬地区Yakutsk单台站电离层暴预报的线性经验模型.检验表明,该模型对Yakutsk站夏季和春秋季电离层扰动的预测效果明显优于IRI中的STORM模型,模型改善度平均可达35%左右;能很好地反映磁暴连续发生的情况下δfoF₂分阶段下降的形态;能反映出一定的电离层正相扰动及其持续时间,但是正扰动的预报精度有待进一步提高.     

The COMPASS-2 satellite and the ground-based LOIS vector sensing radar facility as novel tools for ionospheric plasma diagnostics

To give a more detailed and complete understanding of physical plasma processes that govern the solar–terrestrial space, and to develop qualitative and quantitative models of the magnetosphere–ionosphere–thermosphere coupling, it is necessary to design and build the next generation of instruments for space diagnostics and monitoring. Novel ground-based wide-area sensor networks, such as the LOFAR Outrigger In Scandinavia (LOIS, LOFAR: Low Frequency Array) radar facility, comprising wide band, and vector-sensing radio receivers with full three-dimensional polarization coverage, and multi-spacecraft plasma diagnostics should help solve outstanding problems of space physics and describe long-term environmental changes. The new digital radio frequency analyzer (RFA) on board the low-orbiting COMPASS-2 satellite was designed to monitor and investigate the ionospheric plasma properties. This two-point ground-based and topside ionosphere-located space plasma diagnostic can be a useful new tool for monitoring and diagnosing turbulent plasma properties. The RFA on board the COMPASS-2 satellite is the first in a series of experiments which is planned to be launched into the near-Earth environment.The main purpose of this presentation is to describe new advanced diagnostic techniques of the near-Earth space plasma and point out the scientific challenges of the COMPASS-2 and LOIS experiments.     

A statistical survey of dayside pulsed ionospheric flows as seen by the CUTLASS Finland HF radar

Nearly two years of 2-min resolution data and 7- to 21-s resolution data from the CUTLASS Finland HF radar have undergone Fourier analysis in order to study statistically the occurrence rates and repetition frequencies of pulsed ionospheric flows in the noon-sector high-latitude ionosphere. Pulsed ionospheric flow bursts are believed to be the ionospheric footprint of newly reconnected geomagnetic field lines, which occur during episodes of magnetic flux transfer to the terrestrial magnetosphere - flux transfer events or FTEs. The distribution of pulsed ionospheric flows were found to be well grouped in the radar field of view, and to be in the vicinity of the radar signature of the cusp footprint. Two thirds of the pulsed ionospheric flow intervals included in the statistical study occurred when the interplanetary magnetic field had a southward component, supporting the hypothesis that pulsed ionospheric flows are a reconnection-related phenomenon. The occurrence rate of the pulsed ionospheric flow fluctuation period was independent of the radar scan mode. The statistical results obtained from the radar data are compared to occurrence rates and repetition frequencies of FTEs derived from spacecraft data near the magnetopause reconnection region, and to ground-based optical measurements of poleward moving auroral forms. The distributions obtained by the various instruments in different regions of the magnetosphere were remarkably similar. The radar, therefore, appears to give an unbiased sample of magnetopause activity in its routine observations of the cusp footprint.     

Spatial observations by the CUTLASS coherent scatter radar of ionospheric modification by high power radio waves

Results are presented from an experimental campaign in April 1996, in which the new CUTLASS (Co-operative UK twin-located Auroral Sounding System) coherent scatter radar was employed to observe artificial field aligned irregularities (FAI) generated by the EISCAT (European Incoherent SCATter) heating facility at Tromso, Norway. The distribution of back-scatter intensity from within the heated region has been investigated both in azimuth and range with the Finland component of CUTLASS, and the first observations of artificial irregularities by the Iceland radar are also presented. The heated region has been measured to extend over a horizontal distance of 170 ± 50 km, which by comparison with a model of the heater beam pattern corresponds to a threshold electric field for FAI of between 0.1 and O.OlV/m. Differences between field-aligned and vertical propagation heating are also presented.     

The synthesis of travelling ionospheric disturbance (TID) signatures in HF radar observations using ray tracing

Characteristic signatures are often observed in HF radar range-time-intensity plots when travelling ionospheric disturbances (TIDs) are present. These signatures, in particular the variation of the F-region skip distance, have been synthesised using a ray tracing model. The magnitude of the skip variation is found to be a function of the peak electron density perturbation associated with the TID and radar frequency. Examination of experimental observations leads to an estimate of the peak electron density perturbation amplitude of around 25% for those TIDs observed by the CUTLASS radar system. The advantage of using the skip variation over the radar return amplitude as an indicator of density perturbation is also discussed. An example of a dual radar frequency experiment has been given. The investigation of the effect of radar frequency on the observations will aid the optimisation of future experiments.     

极区电离层对IMF Bz 4次快速转向的响应——EISCAT/ESR雷达观测

利用EISCAT VHF和EISCAT Svalbard（ESR）雷达观测数据,对2003年2月12日IMF Bz分量4次快速方向转换期间,极区电离层,尤其是极尖/极隙区的响应特征进行了分析研究.随着IMF Bz方向的多次快速变化,地面雷达观测到极尖/极隙区所在位置随着开放-闭合磁力线边界在纬度方向上来回移动.在此期间,极区电离层等离子体水平对流多次反向,表现出与IMF Bz分量强的负相关性.进一步分析表明：极区磁层-电离层系统在日侧对IMF极性变化的平均响应时间约为3 min.     

First observation of quasi-2-day oscillations in ionospheric plasma frequency at fixed heights

The existence and development of the quasi-2-day oscillations in the plasma frequency variations of the F region at northern middle latitudes are investigated. A new approach to study the quasi-2-day oscillations is presented, using a methodology that allows us to do such a study at fixed heights. The hourly values of plasma frequency at fixed heights, from 170 km to 220 km at 10 km step, obtained at the Observatori de lEbre station (40.8°N, 0.5°E) during 1995 are used for analysis. It is found that quasi-2-day oscillations exist and persisted in the ionospheric plasma frequency variations over the entire year 1995 for all altitudes investigated. The dominant period of oscillation ranges from 42 to 56 h. The amplitude of oscillation is from 0.1 MHz to 1 MHz. The activity of the quasi-2-day oscillation is better expressed during the summer half year when several enhancements, about 15–30 days in duration, were observed. The largest enhancements of the oscillation occurred during early June, July and early August; i. e., near and after the summer solstice when the 2-day wave in the middle neutral atmosphere typically displays its largest activity in the Northern Hemisphere. The results obtained may help us understand better the possible influencing mechanisms between the 2-day wave in the middle neutral atmosphere and the ionospheric quasi-2-day oscillations.     

CUTLASS Finland radar observations of the ionospheric signatures of flux transfer events and the resulting plasma flows

The CUTLASS Finland radar has been run in a two-beam special scan mode, which offered excellent temporal and spatial information on the flows in the high-latitude ionosphere. A detailed study of one day of this data revealed a convection reversal boundary (CRB) in the CUTLASS field of view (f.o.v) on the dayside, the direction of plasma flow either side of the boundary being typical of a dawn-cell convection pattern. Poleward of the CRB a number of pulsed transients are observed, seemingly moving away from the radar. These transients are identified here as the ionospheric signature of flux transfer events (FTEs). Equatorward of the CRB continuous backscatter was observed, believed to be due to the return flow on closed field lines. The two-beam scan offered a new and innovative opportunity to determine the size and velocity of the ionospheric signatures associated with flux transfer events and the related plasma flow pattern. The transient signature was found to have an azimuthal extent of 1900 ± 900 km and an poleward extent of 250 km. The motion of the transient features was in a predominantly westward azimuthal direction, at a velocity of 7.5 ± 3 km.