不同混合层深度算法的全球适用性对比分析

混合层深度是研究海洋上层动力过程及海气相互作用的一个至关重要的物理量,准确估算混合层深度对上层海洋动力学和热力学的深入研究具有重要的科学意义。本文基于Argo实时观测剖面数据,分海域、分季节对比分析了目前常用的几种混合层深度算法的异同与优缺点。结果表明,理论上最大角度法的精确度最高,曲率法其次,然后是阈值法和最优线性拟合法。最大角度法和曲率法的结果比较接近,实测数据表明曲率法的时空适用性更广。阈值法、最优线性拟合法分别受梯度阈值和密度（或温度）梯度变化的制约,其计算的混合层深度相对较浅。各种算法的差异性随着季节跃层的增强而逐渐减小,且北半球的差异小于南半球。     

《长春市城市总体规划》图集的编辑制作

简要阐述了用MapGIS、CoreIDRAW两种软件结合制作《长春市城市总体规划》图集的设计、编制方法,采用图层模式有效管理地图集的数据．根据实践得出了一些经验和技巧。     

基于波形拟合的中国东海地区410 km间断面附近速度结构研究

410 km间断面是地幔转换带的顶界面,对其速度结构和起伏形态开展地震学探测有助于认识地球内部物质组成和相关的地球动力学过程.本文选取了由中国数字地震台网记录到的位于琉球俯冲区的一个中源地震P波宽频带波形资料,利用三重震相波形拟合研究了中国东海地区410 km间断面附近的精细速度结构.结果表明：中国东海地区下方410 km间断面整体表现为一尖锐的速度界面且有8~15 km的小幅抬升;该间断面之上存在52~62 km厚的低速层,其P波速度降低0.5%~1.6%;440 km深度以下存在1.0%~3.0%的P波高速异常.结合前人在该地区的层析成像结果,我们推测该高速异常体与西太平洋俯冲板片在中国东海地区地幔转换带内的滞留有关;板片内水相E分解使得转换带内水含量增加,这引发了410 km间断面的抬升;410 km间断面之上的低速层应与含水矿物脱水导致的地幔橄榄岩部分熔融有关.     

核幔耦合对地球自由核章动的激发影响

地球自由核章动（FCN）是地幔与液核相互作用的重要动力学现象,其激发机制涉及地表流体层、地幔和地核等圈层之间的耦合,此前研究多利用地表流体层角动量数据单独研究其对FCN的激发,对核幔耦合的影响考虑不足.本文基于角动量守恒理论分析了核幔耦合对FCN周期及振幅的影响,并结合多个大气及海洋角动量函数时间序列首次估算了核幔耦合在FCN激发过程中的贡献.结果表明核幔耦合对FCN周期产生的固定和时变影响对FCN激发的作用均不可忽视,尤其时变影响可达几十个微角秒,对于进一步解释FCN时变特征非常重要;核幔耦合对FCN振幅的直接影响是地表流体层的激发与实测FCN不相符的主要原因,黏滞、电磁和地形等耗散耦合的存在对地表流体的激发振幅有67%左右的减弱效果.     

基于卷积完全匹配层的非规则网格时域有限元探地雷达数值模拟

复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数m、R和κ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟.     

Spatial behaviour of Rayleigh waves in layered half‐spaces under active surface sources

In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms.     

强震记录对三个土层地震反应分析程序的检验

为了对国内外广泛使用的3个一维土层地震反应分析程序（DENSOR98、SHAKE91和LSSRLI-1）的计算精度进行检验,本文以4次地震中唐山响瞠井下三维台阵的基岩强震记录为输入地震动,利用这3个程序对此测井剖面分别进行计算,计算结果分别与同次地震中相应的台阵地表强震记录相比较,可得如下结论：在地震动不大的情况下,3个程序都可满足工程要求,其差别为SHAKE91算得的峰值偏大,DENSOR98算得的峰值偏小,LSSRLI-1得到的反应谱更接近实际,DENSOR98得到的反应谱偏低,从安全角度和经济角度综合来看,LSSRLI-1程序较好。     

黄土高填方场地地震动参数特性分析

削山造地后形成的黄土高填方场地对地震动参数特性影响较大。以实际工程项目为研究对象,构造不同填土高度的计算剖面,采用一维等效线性化方法计算土层地震动参数,分析基岩地震动输入参数和填土高度对场地地表放大效应的影响。研究表明:场地地震放大系数随填土层的高度增加呈递减趋势。在多遇地震动和基本地震动作用下,场地地震放大系数递减速度比罕遇地震动和极罕遇地震动作用下要大。当填土高度达到一定程度后放大效应趋于平稳;填土高度的变化,会改变地表加速度反应谱的形状。填土高度越大,地表反应谱长周期的频谱成分越显著,反应谱曲线向后移,反应谱峰值点均明显向长周期移动,并出现多个峰值点,反应谱特征周期值变大;下伏基岩的刚度越大,地表峰值加速度的放大效应越大。地表加速度反应谱特征值相比变小。当填土高度增大到一定程度时,下伏基岩的种类对地表地震动特性影响则不明显。该研究成果对高填方场地的地震安全评价和结构抗震设计提供参考。     

Variation in soil saturated hydraulic conductivity along the hillslope of collapsing granite gullies

Saturated hydraulic conductivity (K_s) affects the soil hydrological process and is influenced by many factors that exhibit strong spatial variations. To accurately measure K_s and its scale, spatial variability and relationship with collapsing gullies, we analysed four double-ring infiltrometer diameters in three soil layers during in situ experiments designed to measure K_s in two typical collapsing gullies (three slope sites) in Tongcheng County of China. The results showed that K_s increased with increasing inner ring diameter, but no significant difference existed between inner diameters of 30 and 40 cm. The K_s in red soil layers was higher than that in sandy soil layers, the transition layers had the lowest values. K_s also varied with slope position, gradually decreasing with distance from the gully head. The suggestion is that the spatial variation in K_s is affected not only by the intrinsic soil properties but also by the interaction with the collapsing gully.     

Analytical approach for sheet flow transport in purely acceleration-skewed oscillatory flow

An instantaneous analytical approach is developed to predict sheet flow transport in purely acceleration-skewed oscillatory flow. The approach is derived from exponential approximations of velocity and concentration profiles above a mobile seabed, and it particularly considers factors of phase lead; phase lag (i.e. phase residual and phase shift); acceleration modification; and asymmetries in shear stress, roughness height, and boundary layer development. The approach can predict net boundary layer flow above a mobile seabed, and can revert to the classical bedload model. Instantaneous and net sediment transport rates are studied using the approach. The instantaneous sediment transport rate in an onshore flow stage can be approximated by a power function of velocity in which the exponent is confirmed to range between 1 and 5 with a decrease in the phase residual. The net sediment transport rate predicted using the approach is validated using a considerable amount of measured data, and compared with existing instantaneous and half-period type models that consider the phase lag or acceleration modification. For the net sediment transport rate in purely acceleration-skewed oscillatory flow, the phase residual is less important than the acceleration-skewed boundary layer difference between onshore and offshore acceleration stages.