Incoherent scatter radar observations of AGW/TID events generated by the moving solar terminator

Observations of traveling ionospheric disturbances (TIDs) associated with atmospheric gravity waves (AGWs) generated by the moving solar terminator have been made with the Millstone Hill incoherent scatter radar. Three experiments near 1995 fall equinox measured the AGW/TID velocity and direction of motion. Spectral and cross-correlation analysis of the ionospheric density observations indicates that ST-generated AGWs/TIDs were observed during each experiment, with the more-pronounced effect occurring at sunrise. The strongest oscillations in the ionospheric parameters have periods of 1.5 to 2 hours. The group and phase velocities have been determined and show that the disturbances propagate in the horizontal plane perpendicular to the terminator with the group velocity of 300–400 m s^–1 that corresponds to the ST speed at ionospheric heights. The high horizontal group velocity seems to contradict the accepted theory of AGW/TID propagation and indicates a need for additional investigation.     

Seasonal variation of space–time parameters of internal gravity waves at Kharkiv (49°30′N, 36°51′E)

Internal gravity waves (IGWs) over Kharkiv (49°30′N, 36°51′E) have been studied by using an automatic goniometer of a meteor radar (AG MR), with its antenna directed to the East. The AG MR carries out Doppler measurements of the radial drift velocity and the position of the reflecting area of meteor trails. In order to obtain information about IGW parameters, an algorithm was used that processes the AG MR wind data by dividing the measured volume into several sub-volumes, and by carrying out wavelet analysis of wind variations in these areas. Results of 1 year (1987) of AG MR data show a good qualitative correspondence with literature results. There is an increase of IGW activity during the change of the zonal background wind in spring and autumn. It was found that the mean IGW horizontal phase velocity and predominating horizontal propagation direction change with season. However, the mean period (1.5 h), the dominant amplitude (30 m/s) and the vertical phase velocity do not strongly vary in the course of the year. The mean vertical and horizontal wavelengths were found to be 40 and 250 km, respectively.     

Determining parameters of large-scale traveling ionospheric disturbances of auroral origin using GPS-arrays

The intention in this paper is to investigate the form and dynamics of large-scale traveling ionospheric disturbances (LS TIDs) of auroral origin. We have devised a technique for determining LS TID parameters using GPS-arrays whose elements can be selected from a large set of GPS stations forming part of the International GPS Service network. The method was used to determine LS TID parameters during a strong magnetic storm of September 25, 1998. The North-American sector where many GPS stations are available, and also the time interval 00:00–06:00 UT characterized by a maximum value of the derivative Dst were used in the analysis. The study revealed that this period of time was concurrent with the formation of the main ionospheric trough with a conspicuous southward wall in the range of geographic latitudes 50–60° and the front width of no less than 7500 km. The auroral disturbance-induced large-scale solitary wave with a duration of about 1 h and the front width of at least 3700 km propagated in the equatorward direction to a distance of no less than 2000–3000 km with the mean velocity of about 300 m/s. The wave front behaved as if it ‘curled’ to the west in longitude where the local time was around afternoon. Going toward the local nighttime, the propagation direction progressively approximated an equatorward direction.     

大尺度电离层行扰的GPS观测

利用日本境内的高空间分辨率的双频GPS台站资料,观测研究了发生于2000年7月中旬太阳强活动期间的一次大尺度电离层行扰. 结果表明:在7月15日11:00UT-1:00UT期间观测区域的电离层中出现了大尺度电离层行扰. 在15:00UT之前,扰动周期为2h左右,在15:00UT以后,扰动周期为1h左右;总电子含量扰动幅度的变化范围约为1-2TECU;通过对15:00-17:00UT之间总电子含量扰动曲线同相位点的分析,发现这期间的电离层行扰的扰动速度约为600-700m/s,扰动波长在2200km左右,扰动传播的方向几乎沿着经线从高纬向低纬传播. 该行扰与此次强太阳活动有直接的关系,因其发生在7月15日的磁暴急始之前数小时,因此与磁暴本身没有因果关系,应与磁暴之前先期到达地球空间的高能质子流有关.     

耗散大气中水平不均匀风场对内重力波传播的影响

采用包括耗散的射线跟踪方法，计算了在水平不均匀风场作用下，不同尺度重力波从对流层直至220km观测高度的传播，结果表明，垂直于重力波传播方向的风以及风剪切能够引起波射线的折射，从而导致重力波明显偏离初始传播方向.在强顺风场作用下，由于风场引起的捕获，大量重力波不能传播到观测高度.由于风场引起的多普勒频移，小周期的重力波在弱顺风条件下能够传播到观测高度.由于反射作用，强逆风场不支持周期低于约18min的较高频重力波的传播.而在弱逆风作用下，大部分中尺度范围重力波都能够传播到观测高度.本文统计了武汉电离层观象台的TID观测数据随热层风场的分布，统计结果与模拟结果符合较好.     

Estimating the Characteristics of Traveling Ionospheric Disturbances from Vertical Incidence Ionograms within a Compound Parabolic Layer Model

The characteristic U-shaped traces (cusps) in ionograms are associated with the passage of traveling ionospheric disturbances (TIDs), which lead to horizontal electron density gradients in the ionosphere and, therefore, to off-angle reflections in radio sounding. A new form of representation is considered for daily ionospheric sounding data. A compound parabolic layer model is proposed, which allows analytical calculation of ray paths to speed up the “homing-in” of the rays. Changes in the shape of the trace in the ionogram due to varying the TID characteristics are examined. A discussion is given of the possibilities for estimating TID characteristics from digitized vertical sounding ionograms.     

Gravity wave signatures at mesopause heights

The ALOMAR SOUSY Radar operated at 53.5 MHz has been used in a five-beam configuration to study dynamical processes at gravity wave periods in the summer polar mesosphere. A case study of a gravity wave with a period of about 9 min is presented and analysed in some detail. The three-dimensional wave number vector is determined from the phase information of the 9-min velocity oscillations obtained in all beam directions and all range gates. The horizontal wavelengths in the north and east direction are estimated to be about 60 and 50 km, respectively. The echo power variations, simultaneously observed at different beam pointing positions, are investigated using cross-correlation analysis. The results show that these wave-associated variations lead to horizontal and vertical deformation of the echoing layers. The horizontal wavelength in the north-south direction, derived from the echo power modification, is in good agreement with the one obtained from the velocity analysis.     

与日夜交替线移动相关的中尺度电离层扰动GPS网观测

本文选取由国际GPS服务中心(IGS)提供的北美中纬度地区GPS网TEC观测数据,通过多通道最大熵频谱分析方法研究了2005年（太阳活动低年）地磁活动平静期间日出和日落时由于日夜交替线的移动而激发的中尺度电离层扰动(MSTID),并统计分析其季节变化特性.结果表明：(1) 日出或日落期间,在中纬度地区经常观测到由日夜交替线移动激发的中尺度电离层扰动.扰动主要沿日夜交替线运动方向传播,平均持续时间约2～3 h;振幅在0.2～0.8 TECU之间,水平波长,水平相速度和周期分别为300±150 km,150±80 m/s和25±15 min;(2) 由日夜交替线移动激发的中尺度电离层扰动在春秋分出现率较少;在夏季,扰动在黄昏时出现率达最大值,在日出后少量出现;而冬季则日出后的扰动效应更为明显.分析表明,在中纬地区,这种扰动出现率随季节的变化与不同季节的日出日落时刻太阳EUV辐射通量变化过程的快慢,以及电离层中离子损失过程快慢有关.     

An experimental study of wave packets in travelling ionospheric disturbances

The complex demodulation method was employed to analyse observations of TIDs carried out by means of multifrequency HF Doppler sounding at vertical incidence over Almaty (76°55′E, 43°15′N). Two types of behaviour of the amplitude and phase in the band-pass filtered Doppler shift oscillations were revealed. These were packet-like wave structures and quasi-stochastic TIDs. A study of the phase behaviour inside wave packets showed that most of them had a quasi-monochromatic structure, some of them had oscillations with a period which increased with time, and none had oscillations with a period which decreased with time. Phase jumps between nearby packets were explained as a result of interference between two adjacent packets having similar oscillation periods. A method of calculating the correlation coefficient for the wave packets propagating on a background of quasi-stochastic TIDs was developed. This showed that wave packets (having periods between 7 and 50 min) travelled in the vertical plane practically without any loss of coherence, while quasi-stochastic TIDs had a short coherence length; the velocity of the quasi-stochastic TIDs was noticeably larger than that of the wave packets.     

Analysis of gravity waves in the tropical middle atmosphere over La Reunion Island (21°S, 55°E) with lidar using wavelet techniques

The capabilities of the continuous wavelet transform (CWT) and the multiresolution analysis (MRA) are presented in this work to measure vertical gravity wave characteristics. Wave properties are extracted from the first data set of Rayleigh lidar obtained between heights of 30 km and 60 km over La Reunion Island (21°S, 55°E) during the Austral winter in 1994 under subtropical conditions. The altitude-wavelength representations deduced from these methods provide information on the time and spatial evolution of the wave parameters of the observed dominant modes in vertical profiles such as the vertical wavelengths, the vertical phase speeds, the amplitudes of temperature perturbations and the distribution of wave energy. The spectra derived from measurements show the presence of localized quasi-monochromatic structures with vertical wavelengths <10 km. Three methods based on the wavelet techniques show evidence of a downward phase progression. A first climatology of the dominant modes observed during the Austral winter period reveals a dominant night activity of 2 or 3 quasi-monochromatic structures with vertical wavelengths between 1/2 km from the stratopause, 3/4 km and 6/10 km observed between heights of 30 km and 60 km. In addition, it reveals a dominant activity of modes with a vertical phase speed of –0.3 m/s and observed periods peaking at 3/4 h and 9 h. The characteristics of averaged vertical wavelengths appear to be similar to those observed during winter in the southern equatorial region and in the Northern Hemisphere at mid-latitudes.     

Observations of short-period waves with an all-sky camera in the infrared airglow above Yakutsk

The parameters of internal gravity waves detected based on the variations in the hydroxyl molecule emission are statistically analyzed. The wave structures were registered with an all-sky infrared camera at Maimaga optical station (? = 63° N, λ = 129.5° E). The data obtained in the winter period of 1998–2002 are analyzed. In total, 162 waves, the majority of which propagated westward, were recorded. The wavelengths vary from 15.4 to 100 km (the average value is ～31 km); the observed horizontal phase velocities change from 19 to 166 m/s (the average value is ～60 m/s), and the estimated periods are 9–90 min (the average value is ～11 min). The statistical characteristics of the waves do not differ from those of similar waves at middle and low latitudes. The azimuthal dependence of the wave propagation direction is consistent with the theory of wave filtration by a background wind in the middle atmosphere. Probable sources of the waves are mountain ranges located at a distance of 200 km east of the observation site. Somewhat greater values of the mean wavelength and wave propagation velocities than those recorded at lower latitudes may be due to the lower loss of energy and velocity of the waves during their propagation from the source to the mesosphere, although other causes are not ruled out. Ripple-type waves have the same direction of propagation as band-type waves.     

Statistics of gravity waves in the lower stratosphere over Beijing based on high vertical resolution radiosonde

Atmospheric gravity waves, with small to medium scales, prevail in the atmosphere and have global ef- fects. Many researches show that gravity waves are the main source that causes the variation of wind and temperature field in the stratosphere, and that the break-up of upward propagating gravity waves is the dominant sources of small scale turbulent and mixing processes in the middle atmosphere. Theories and ob- servations indicate that the redistribution of momen- tum, caused by the generati…     

Three-dimensional structure of mesoscale eddies in the western tropical Pacific as revealed by a high-resolution ocean simulation

The three-dimensional structure of mesoscale eddies in the western tropical Pacific(6°S–20°N, 120°E–150°E)is investigated using a high-resolution ocean model simulation. Eddy detection and eddy tracking algorithms are applied to simulated horizontal velocity vectors, and the anticyclonic and cyclonic eddies identified are composited to obtain their three-dimensional structures. The mean lifetime of all long-lived eddies is about 52 days, and their mean diameter is 147 km. Two typical characteristics of mesoscale eddies are revealed and possible dynamic explanations are analyzed. One typical characteristic is that surface eddies are generally separated from subthermocline eddies along the bifurcation latitude(~13°N) of the North Equatorial Current in the western tropical Pacific, which may be associated with different eddy energy sources and vertical eddy energy fluxes in subtropical and tropical gyres. Surface eddies have maximum swirl velocities of 8–9 cm s~(-1) and can extend to about 1500 m depth. Subthermocline eddies occur below 200 m, with their cores at about 400–600 m depth, and their maximum swirl velocities can reach 10 cm s~(-1). The other typical characteristic is that the meridional velocity component of the eddy is much larger than the zonal component. This characteristic might be due to more zonal eddy pairs(two eddies at the same latitude),which is also supported by the zonal wavelength(about 200 km) in the high-frequency meridional velocity component of the horizontal velocity.     

Storm observations by remote sensing and influences of gustiness on ocean waves and on generation of rogue waves

The impact of the gustiness on surface waves under storm conditions is investigated with focus on the appearance of wave groups with extreme high amplitude and wavelength in the North Sea. During many storms characterized by extremely high individual waves measured near the German coast, especially in cold air outbreaks, the moving atmospheric open cells are observed by optical and radar satellites. According to measurements, the footprint of the cell produces a local increase in the wind field at sea surface, moving as a consistent system with a propagation speed near to swell wave-traveling speed. The optical and microwave satellite data are used to connect mesoscale atmospheric turbulences and the extreme waves measured. The parameters of open cells observed are used for numerical spectral wave modeling. The North Sea with horizontal resolution of 2.5?km and with focus on the German Bight was simulated. The wind field “storm in storm,” including moving organized mesoscale eddies with increased wind speed, was generated. To take into account the rapid moving gust structure, the input wind field was updated each 5?min. The test cases idealized with one, two, and four open individual cells and, respectively, with groups of open cells, with and without preexisting sea state, as well the real storm conditions, are simulated. The model results confirm that an individual-moving open cell can cause the local significant wave height increase in order of meters within the cell area and especially in a narrow area of 1–2?km at the footprint center of a cell (the cell's diameter is 40–90?km). In a case of a traveling individual open cell with 15?m·s^?1 over a sea surface with a preexisting wind sea of and swell, a local significant wave height increase of 3.5?m is produced. A group of cells for a real storm condition produces a local increase of significant wave height of more than 6?m during a short time window of 10–20?min (cell passing). The sea surface simulation from modeled wave spectra points out the appearance of wave groups including extreme individual waves with a period of about 25?s and a wavelength of more than 350?m under the cell's footprint. This corresponds well with measurement of a rogue wave group with length of about 400?m and a period of near 25?s. This has been registered at FiNO-1 research platform in the North Sea during Britta storm on November 1, 2006 at 04:00 UTC. The results can explain the appearance of rogue waves in the German Bight and can be used for ship safety and coastal protection. Presently, the considered mesoscale gustiness cannot be incorporated in present operational wave forecasting systems, since it needs an update of the wind field at spatial and temporal scales, which is still not available for such applications. However, the scenario simulations for cell structures with appropriate travel speed, observed by optical and radar satellites, can be done and applied for warning messages.     

A study of atmospheric gravity waves and travelling ionospheric disturbances at equatorial latitudes

A global coupled thermosphere-ionosphere-plasmasphere model is used to simulate a family of large-scale imperfectly ducted atmospheric gravity waves (AGWs) and associated travelling ionospheric disturbances (TIDs) originating at conjugate magnetic latitudes in the north and south auroral zones and subsequently propagating meridionally to equatorial latitudes. A fast dominant mode and two slower modes are identified. We find that, at the magnetic equator, all the clearly identified modes of AGW interfere constructively and pass through to the opposite hemisphere with unchanged velocity. At F-region altitudes the fast AGW has the largest amplitude, and when northward propagating and southward propagating modes interfere at the equator, the TID (as parameterised by the fractional change in the electron density at the F2 peak) increases in magnitude at the equator. The amplitude of the TID at the magnetic equator is increased compared to mid-latitudes in both upper and lower F-regions with a larger increase in the upper F-region. The ionospheric disturbance at the equator persists in the upper F-region for about 1 hour and in the lower F-region for 2.5 hours after the AGWs first interfere, and it is suggested that this is due to enhancements of the TID by slower AGW modes arriving later at the magnetic equator. The complex effects of the interplays of the TIDs generated in the equatorial plasmasphere are analysed by examining neutral and ion winds predicted by the model, and are demonstrated to be consequences of the forcing of the plasmasphere along the magnetic field lines by the neutral air pressure wave.     

北极地区低平流层惯性重力波的观测研究

南极地区重力波活动有大量报道,相对而言,北极地区重力波的研究还很少.本文利用极区Ny-Alesund站点（78.9°N,11.9°E）无线电探空仪从2012年4月1日到2017年3月31日共5年的观测数据,统计分析了北极地区低平流层惯性重力波的特征.观测显示,月平均纬向风在20 km以下盛行东向风,再随着高度增加,逐渐呈现出半年振荡现象.对流层顶高度在5~13 km范围内变化,其月平均高度显示出年循环,最高出现在夏季,约为10 km,最低出现在冬季,约为8.5 km.对流层和低平流层月平均温度都显示出明显的年周期变化,这与中低纬度观测结果有所不同.结合Lomb-Scargle谱分析和矢端曲线方法,估算了准单色惯性重力波参数.个例研究表明,低平流层惯性重力波呈现出远离源区的自由传播性质.统计结果显示,惯性重力波的水平和垂直波长分别集中在50~450 km和1~4 km范围内,本征频率集中在1~2.5倍惯性频率间,这些值都比中低纬度观测值稍小.垂直方向本征相速度主要集中在-0.3~0 m·s^-1,而纬向和经向本征相速度集中在-40~40 m·s^-1之间.在5年的观测中,大约91.5%的惯性重力波向上传播.在冬季和早春,由于极地平流层极涡活动,激发出向下传播的惯性重力波,因此,向下传播的比例上升到相应月份的20%左右.由于低层大气盛行的东向风的滤波效应,低平流层大部分惯性重力波向西传播.波能量呈现出明显的年周期变化,最大值在冬季、最小值在夏季,与北半球中低纬度观测结果一致,表明北半球重力波活动普遍冬季强、夏季弱.     

Vertical plasma drift velocities in the polarization jet observation by ground Doppler measurements and driftmeters on DMSP satellites

Vertical and horizontal plasma drifts are investigated during the polarization jet (PJ) detection in the F2 ionospheric layer based on the Doppler measurements at the Yakutsk meridian chain of subauroral ionospheric stations. It is shown that the velocities of vertical and horizontal drifts are significantly higher than the background motion during PJ observation periods. The ionospheric plasma motion direction changes from upward to downward on the polar edge of the main ionospheric trough. Doppler measurements on the DPS-4 ionosondes are compared with the simultaneous measurements of the plasma drift on the DMSP satellites during their passage near the Yakutsk meridian. The two kinds of measurements are in good agreement with each other. During the magnetic storm of June 23, 2005, by measurements of the DMSP satellites, the velocities of upward plasma flows were 1.0–1.4 km/s at a satellite altitude of 850 km. In the ionospheric F region, this speed corresponds to 150 m/s. According to satellite measurements, the westward drift velocity reached 2.5 km/s. The development of the polarization jet in the ionosphere was accompanied by a tenfold decrease in the electron density in 15–30 min.     

Modeling the dynamic process of tsunami earthquake by liquid-solid coupling model

Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock     

Wind tunnel investigation of horizontal and vertical sand fluxes of ascending and descending sand particles in aeolian sand transport

The horizontal and vertical sand mass fluxes in aeolian sand transport are investigated in a wind tunnel by PTV (particle tracking velocimetry). According to the particle velocity and volume fraction of each individual particle from PTV images, the total horizontal sand mass flux, the horizontal mass fluxes of ascending and descending sand particles, and upward and downward vertical sand fluxes are analyzed. The results show that the horizontal mass fluxes of ascending and descending sand particles generally decrease with the increase of height and can be described by an exponential function above about 0.03 m height. At the same friction velocity, the decay heights of the total horizontal sand mass flux and the horizontal mass fluxes of ascending and descending sand particles are very similar. The proportion of horizontal mass flux of ascending sand particles is generally about 0.3–0.42, this means the horizontal mass flux of descending sand particles makes an important contribution to the total horizontal sand mass flux. Both the upward and downward vertical sand mass fluxes generally decrease with height and they are approximately equal at the same height and friction velocity. The relation between upward (or downward) vertical sand mass flux and horizontal sand mass flux can be described by a power function. The present study is used to help understand the transport of ascending and descending sand particles. Copyright © 2016 John Wiley & Sons, Ltd.     

Quantitative models of the fallout and dispersal of tephra from volcanic eruption columns

A theoretical model of clast fallout from convective eruption columns has been developed which quantifies how the maximum clast size dispersal is determined by column height and wind strength. An eruption column consists of a buoyant convecting region which rises to a heightH _B where the column density equals that of the atmosphere. AboveH _B the column rises further to a heightH _T due to excess momentum. BetweenH _T andH _B the column is forced laterally into the atmosphere to form an upper umbrella region. Within the eruption column, the vertical and horizontal velocity fields can be calculated from exprimental and theoretical studies and consideration of mass continuity. The centreline vertical velocity falls as a nearly linear function over most of the column's height and the velocity decreases as a gaussian function radially away from the centreline. Both column height and vertical velocity are strong functions of magma discharge rate. From calculations of the velocity field and the terminal fall velocity of clasts, a series of particle support envelopes has been constructed which represents positions where the column vertical velocity and terminal velocity are equal for a clast of specific size and density. The maximum range of a clast is determined in the absence of wind by the maximum width of the clast support envelope.The trajectories of clasts leaving their relevant support envelope at its maximum width have been modelled in columns from 6 to 43 km high with no wind and in a wind field. From these calculations the shapes and areas of maximum grain size contours of the air-fall deposit have been predicted. For the no wind case the theoretical isopleths show good agreement with the Fogo A plinian deposit in the Azores. A diagram has been constructed which plots, for a particular clast size, the maximum range normal to the dispersal axis against the downward range. From the diagram the column height (and hence magma discharge rate) and wind velocity can be determined. Historic plinian eruptions of Santa Maria (1902) and Mount St. Helens (1980) give maximum heights of 34 and 19 km respectively and maximum wind speeds at the tropopause of m/s and 30 m/s respectively. Both estimates are in good agreement with observations. The model has been applied to a number of other plinian deposits, including the ultraplinian phase of theA.D. 180 Taupo eruption in New Zealand which had an estimated column height of 51 km and wind velocity of 27 m/s.