高能离子束流引起低频电磁波动的特性

理论和数值模拟研究了日冕区高速低密离子束流与等离子体相互作用产生的低频电磁波动特征. 结果表明,在线性阶段,系统波动主要表现为短波长的非共振模的激发;在非线性阶段,长波模占优,短波长波逐渐衰减,此时系统波动表现为具有Alfvén波动特征,波的能谱为双幂律谱. 这些结果可用来解释太阳风中的Alfvén湍动现象.     

垂直切变风场中声重波的谱传播

本文研究声重波在垂直切变风场中传播时谱结构的改变。根据可压缩情况下的流体力学方程组,导出垂直切变风场中大气声重波的谱方程。并指出,这个谱方程是Taylor-Goldstein方程在可压缩和任意垂直波长情况下的推广。由谱方程出发,考虑低边界上的声重波强迫为高斯谱型,采用差分方法,我们计算了ω-K_H二维空间中的高层谱响应。结果表明,在具有垂直切变的背景风场作用下,声重波的谱结构发生了改变:1.在原强迫谱所覆盖的谱区段内,一些谱成份被挖除,形成了尖锐的谱峰群;2.最大谱峰发生了偏移;3.波谱变得狭窄。这些变化是垂直切变风场对声重波谱滤波的结果。文中还将上述结论与电离层观测结果进行比较,解释了某些观测现象。     

空间尘埃等离子体中的重力波特性

建立尘埃等离子体中重力波的基本方程,推导尘埃等离子体中重力波的色散关系,分析地球极区中间层顶处尘埃等离子体层中的重力波特性,研究了重力波在电子密度垂直分层的尘埃等离子体中的反射. 结果表明尘埃等离子体改变了通常大气中的重力内波的色散关系,限制了小水平波数重力内波的传播,改变了波的能量特性,减小了重力波在不均匀大气中垂直向上传播时振幅的增长;在尘埃等离子体中传播时重力波可被电子密度垂直分层的结构反射而导致波能量的集中, 它产生的湍动所导致的空间电子密度的不均匀性分布是极区上空PMSEs的可能机制.     

P波和S波接收函数的贝叶斯联合反演

S波接收函数对于研究岩石圈速度结构具有重要价值. 本文利用合成地震图技术研究了S波接收函数的动力学特征. 在接收函数非线性复谱比反演方法的基础上,发展了基于贝叶斯理论的P波和S波接收函数的非线性联合反演方法. 结果表明:(1)适用于S波接收函数反演的震中距范围约为55°~80°,S波接收函数反演要求所用远震事件的震级大于5级; (2)与陡变的岩石圈底部界面(LAB)相比,梯度带类型LAB上生成的S_LP转换波相对较弱,台站下方的沉积盖层有助于相对增强S_LP震相; (3)由于S波接收函数径向分量不符合δ脉冲,不依赖于等效震源假定的三分量接收函数多道最大或然性反褶积方法更适合S波接收函数的估计;(4)数值检验的结果表明,在初始模型速度参数偏离真实模型20%的情况下,本文的方法能够预测300 km深度范围内的P波和S波速度结构;(5)观测数据的反演结果表明,由于P波接收函数低频分量相对不足,本文的联合反演方法对于大于100 km深度上地幔的S波速度结构约束相对较弱.     

用数字测高仪漂移测量研究电离层声重波扰动

数字测高仪Digisonde中的漂移测量,常用来研究小尺度电离层扰动,如电离层小不均匀体的漂移。文中提出一种新的分析方法,利用漂移测量数据中多普勒频移和到达角参量的最大熵动态功率谱,估算声重波一类大尺度电离层扰动的水平传播速度和传播方向。作为实例,研究了Millstone Hill测高仪站的漂移观测资料,并对处理结果进行了初步分析。分析结果表明,从数字测高仪漂移测量数据中,可有效地提取声重波一类大尺度电离层扰动的传播参量,在电离层动力过程的研究中很有意义。     

断续节理岩体中地震波响应的数值与理论分析

为研究地震产生的应力波在断续节理岩体中的传播规律和应力分布趋势,首先,采用数值模拟方法分析应力波通过贯通节理的传播规律,并与已有理论研究结果进行对比,验证数值分析的准确性和适用性;然后,对应力波在断续节理岩体中的传播进行数值模拟,分析透射系数在水平方向的分布趋势以及不同节理连续性对波传播的影响,并结合波的衍射原理,给出定性的理论解释。结果表明:应力波通过断续节理时,节理的透射作用会使应力波振幅减小,引起波的衰减,岩桥的衍射作用则会使波阵面由平面变为曲面,波的传播方向发生改变,从而导致应力波振幅在水平方向的分布发生变化;应力波通过断续节理的透射系数与岩桥尺寸L_r和衍射角μ相关,当衍射角比较小时,透射系数主要受岩桥尺寸L_r的影响,当衍射角较大时,岩桥尺寸L_r和衍射角μ共同影响应力波在岩体中的传播。      

强震动记录H/V谱比法计算处理的若干关键环节

强震动记录H/V谱比法自提出以来,已广泛应用于地震工程各研究领域中。但对数据处理的2个关键环节——Taper预处理及傅氏谱平滑未有详细研究。因此,本文针对强震动记录H/V谱比法计算涉及的重要数据处理过程,对四川地区19个强震动台站在汶川地震余震中获取的642组强震动记录开展研究。研究实例表明:S波H/V谱比振幅在周期 < 1s时高于全时程,论证计算H/V谱比时截取S波窗口的必要性;S波窗口的截断会引起傅里叶振幅谱的边瓣效应,若不加以处理,将显著影响低频部分H/V谱比结果,Taper预处理对于消除这种截断误差具有良好效果;在兼顾平滑效果及卓越周期识别准确度的基础上,通过不同带宽的Parzen窗试算,认为0.5Hz带宽宜用于傅里叶H/V谱比曲线平滑。     

平面P波入射下深埋圆形复合式衬砌隧道抗减震机理研究

基于Fourier-Bessel级数展开法,研究深埋圆形三层复合式衬砌洞室在平面P波入射下的动应力集中问题,并给出三层衬砌洞室动应力集中系数级数解析解;依托某IX度地震区管道隧道实际工程,分析不同衬砌刚度组合和厚度组合对洞室动应力集中系数的影响。研究表明:注浆加固洞室围岩和设置减震层都可以降低二次衬砌动应力集中系数;增大围岩注浆区弹性模量和厚度,有利于减小衬砌动应力集中系数,最优围岩注浆区厚度为1倍洞室净空半径;减震层弹性模量降低,减震层厚度增大,二次衬砌动应力集中系数变小,减震层弹性模量宜低于围岩弹模1/20,最优减震层厚度宜取1/50的洞室内净空半径。最后针对实际管道隧道抗减震技术,考虑围岩稳定性,提出"围岩-加固圈-减震层-衬砌"新型减震结构,分析结果表明:对比其他三种抗减震措施,新型减震结构的减震效果最好。     

风积土高铁路基震动波的人工合成

地震荷载作用形成的震动波具有非平稳性、随机性等特点。在室内动三轴试验中不能准确将实际的震动波施加给土体试样,也就无法准确研究风积土土体在地震荷载作用时的动力特性。文中首先基于对风积土路基周边地震资料的分析,以与抗震规范相容的功率谱为目标谱,确定需要控制的反应谱的坐标点和反应谱的容许误差,利用三角级数法和快速傅里叶变换法合成了人工地震波;然后利用M ATLAB和震动持时、频率和幅值等实测数据合成了具有一定频谱特性的人工波,从而极大地方便和简化了对震动荷载作用下风积土的动力特性的研究。     

初始波法线角对中低纬非导管哨声波传播的影响

本文对初始波法线角δ_s不同的中低纬非导管哨声进行了射线追踪研究,结果表明,当δ_s≤0°时,哨声波可到达较大的L值,在另一半球较高纬度处形成“聚焦区”,这组哨声色散值较大;当δ>0°时,随以的增加,哨声射线内移,δ>10°的哨声射线汇聚于另一半球较低纬度处,这组哨声色散值较小。据此提出解释中低纬哨声色散值连续变化的一种可能机制:电离层电子浓度水平梯度的连续变化引起初始波法线角以的连续变化,从而导致射线路径、群时延和色散值的变化。     

A note on Stokes production of turbulence kinetic energy in the oceanic mixed layer: observations in the Baltic Sea

Shear- and convection-driven turbulence coexists with wind-generated surface gravity waves in the upper ocean. The turbulent Reynolds stresses in the oceanic mixed layer can therefore interact with the shear of the wave-generated Stokes drift velocity to extract energy from the surface waves and inject it into turbulence, thus augmenting the mean shear-driven turbulence. Stokes production of turbulence kinetic energy (TKE) is difficult to measure in the field, since it requires simultaneous measurement of the turbulent stress and the Stokes drift profiles in the water column. However, it is readily inferred using second moment closure models of the oceanic mixed layer provided: (1) wave properties are available, along with the usual water mass properties, and radiative and air–sea fluxes needed to drive the mixed layer model and (2) the model skill can be assessed by comparing the model results against the observed dissipation rates of TKE. Comprehensive measurements made during the Reynolds 2002 campaign in the Baltic Sea have made the estimation of Stokes production possible, and in this paper, we report on the effort and the conclusions reached. Measurements of air–sea exchange parameters and water mass properties during the campaign allowed a mixed layer model to be run and the turbulent stress in the water column to be inferred. Simultaneous wave spectrum measurements enabled Stokes drift profile to be deduced and wave breaking to be included in the model run, and the Stokes production of TKE in the water column estimated. Direct measurements of the TKE dissipation rate from an upward traversing microstructure profiler were used to assure that the model could reproduce the turbulent dissipation rate in the water column. The model results indicate that the Stokes production of TKE in the mixed layer is of the same order of magnitude as the shear production and must therefore be included in mixed layer models.     

利用卫星和GPS观测研究低纬局部等离子体浓度增强现象

本文同时利用DMSP、ROCSAT-1卫星数据和地基的GPS观测数据,研究一种与低纬等离子体泡相伴随的局部等离子体浓度增强现象.地基GPS的观测表明电离层总电子含量(TEC)也能反映这种等离子体浓度增强.通过4个观测事例的详细分析表明:这种等离子体浓度增强主要出现在磁纬±10°～±20°的局部区域,有时在近磁赤道区和中纬地区的电离层顶部也能观测到;与等离子体泡的出现规律相似,这种等离子体浓度增强主要出现在地方时21∶00以后,并在午夜后也能观测到.当等离子体浓度增强和等离子体泡发生时,在午夜前一般对应着背景垂直速度明显向上的扰动,在午夜后一般处于等离子体垂直速度下降至反向前的时间段,表明东向电场对于低纬不规则体的产生有非常重要的作用.     

大气重力波触发的Rayleigh－Taylor不稳定性的时空演变

对大气重力波触发的赤道电离层Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性的时空演变进行了数值模拟．结果表明，重力波能在Ｆ区底部触发Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性，这种等离子体不稳定扰动穿过Ｆ峰，到达Ｆ区顶部，形成等离子体泡结构．等离子体泡出现向西倾斜和分岔等特征．当夜Ｆ区较高时，产生等离子体泡的时间仅１９００ｓ．数值模拟结果证实了重力波触发．Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性的理论，解释了大量电离层观测现象．     

大气重力波触发的Rayleigh－Taylor不稳定性的时空演变

对大气重力波触发的赤道电离层Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性的时空演变进行了数值模拟．结果表明，重力波能在Ｆ区底部触发Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性，这种等离子体不稳定扰动穿过Ｆ峰，到达Ｆ区顶部，形成等离子体泡结构．等离子体泡出现向西倾斜和分岔等特征．当夜Ｆ区较高时，产生等离子体泡的时间仅１９００ｓ．数值模拟结果证实了重力波触发．Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性的理论，解释了大量电离层观测现象．     

Measurement of long waves at the harbor of Marina di Carrara,Italy

This paper presents the analysis of measurements of long waves at the harbor of Marina di Carrara, Italy. Nine pressure gauges (eight in the inner harbor and one outside the harbor) were used to record continuously the water surface elevation, for about 3 years; incoming short waves were measured using a directional buoy at about 500 m offshore the breakwaters. The analysis is carried out splitting the energy of the long waves into two bands, i.e., very long waves (VLW; f < 0.003 Hz) and long waves (LW; 0.003 < f < 0.030 Hz); it is found that LW energy is strongly correlated with the energy of the incoming short waves while no correlation (or very poor) is found for VLW. During energetic sea states, surface elevation spectra of VLW and LW, normalized against the wave energy of each frequency band, appear to be autosimilar; this suggests that the spectra of the incoming long waves have a mostly constant shape. Simple numerical computations of harbor resonance, carried out using a mild-slope equation linearized model, indicate that the amplification of the height of the incoming LW generated by short waves (calculated using empirical formulae from previous researches by Melito et al. 2007) inside the harbor can be reasonably predicted using such simple approach, assuming that the LW spectra are rectangular. The shape of the amplification diagram calculated by the numerical model shows reasonable (but not perfect) agreement with the shape of the normalized measured spectra inside the harbor.     

L'EGALISATION SPECTRALE: UN MOYEN D'AMELIORER LA QUALITE DES DONNEES SISMIQUES*

In land seismic surveys spectrum equalization can increase the quality of seismic data in a selected frequency band. The power of lower frequencies in the spectrum of input traces is generally greater than that of higher frequencies, particularly in land seismic surveys because of ground roll. In order to improve the quality of seismic data it is necessary to raise the energy of higher frequencies to the same level as that of lower frequencies, without alteration of the phases. The first step of the method is to compute the amplitude spectrum of each input trace to determine a weighting function which is then applied to the amplitude spectrum in order to balance it. The function is the inverse of the short wavelength variation of the amplitude spectrum. The short wavelength variation can be obtained by interpolation between average values of the modulus of the amplitude spectrum computed in narrow bands within a selected band of frequencies. Another way of obtaining the short wavelength variation is to apply a low-pass filter to the amplitude spectrum. The calculations are readily performed in the frequency domain by the Fourier transform. Spectrum equalization is automatically adjusted to each trace and does not modify the average amplitude in the time domain. However, as the frequency band and energy of the ground roll both vary according to the distance from the shot, spectrum equalization tends to make the spectrum of output traces independent of the offset distance. The use of spectrum equalization before any two-dimensional filtering improves ground roll elimination. Continuity and resolution of horizons are also increased by spectrum equalization before CDP stack. Several examples of applications of spectrum equalization to seismic land and marine surveys are shown.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.

     

青藏高原上空UTLS区域一次地形重力波过程中的物质上传

利用美国航空航天局MERRA(Modern-Era Retrospective Analysis for Research and Applications)再分析资料和MODIS(Moderate-Resolution Imaging Spectroradiometer)卫星资料以及欧洲气象中心ECMWF-Interim(European Centre for Medium-Range Weather Forecasts)再分析资料,分析了发生于青藏高原北侧上空的一次地形重力波事件,并使用中尺度预报模式WRF-ARW.V3.0(Weather Research and Forecasting model,V3.0)对其进行了数值模拟.在此基础上,诊断分析了此次地形重力波在UTLS(Upper Troposphere and Lower Stratosphere)区域造成的物质和能量垂直传输特征.分析结果表明这一中尺度地形重力波信号的水平波长约为600km,与地形扰动水平尺度接近,重力波在对流层中传播的垂直波长约为3km,在垂直方向上随着高度的增加呈现出由东向西倾斜的结构特征.此次地形重力波上传进入平流层并在150hPa附近破碎,波破碎后动量通量在短时间内发生了强烈的衰减,重力波携带的能量在破碎高度附近释放.重力波破碎的同时垂直方向湍流混合变得异常强烈,湍流交换系数可在短时间内增加到背景值的8倍以上,剧烈湍流混合过程导致了对流层上层的空气进入平流层,使下平流层空气出现了位势涡度和臭氧的低值区,在浮力频率的垂直剖面中也可以看到由于地形重力波过程造成的平流层下层浮力频率异常低值区.     

A consistent derivation of the wave-energy equation from basic hydrodynamic principles

Based on a decomposition of the velocity into mean flow, turbulent and wave components, momentum and hereafter a wave-energy equation is derived. It contains a turbulent energy dissipation term which is closed by applying a wave-related mixing length model and linear wave theory solutions. This closure produces a non-linear turbulent wave-energy dissipation including the wave energy in a 5/2 power law. The theory is able to predict correctly the shape of deep-water wave spectra according to Phillips' similarity law.Responsible Editor: Hans Burchard