声波方程频率域有限元参数反演

推导出频率域有限元声波正演方程,为了消除边界反射,将Clayton-Engquist旁轴波动方程吸收边界条件引入频率域,并对有限元刚度矩阵和质量矩阵进行压缩存储,利用广义共轭梯度法求解有限元方程获得正演解.在此基础上,推导出在某一频率下波场数据残差δU与单元物性参数修改量δλ之间关系的Jacobi矩阵,反演方法允许利用地面二维炮集全波场资料与给出初始模型参数的正演值的差值δU,迭代求得δλ.由于计算机内存的限制,方法计算不允许有过多数目的未知数个数,因此还提出了对同一介质物性单元的Jacobi矩阵元素进行压缩组装的措施,从而使反演的未知量个数减少,结合采用共轭梯度迭代法,使得只需利用有效波频段的少数一些频率即可进行迭代反演.正演和反演理论模型的数值模拟结果表明方法是有效的.     

各向异性半空间上方垂直磁偶极子电磁场的球面波展开式

文中推导出单轴各向异性媒质半空间上方垂直磁偶极源电磁场的精确解析表达式. 应用圆柱波函数的球面波展开式和超几何函数理论,场分量中的索末菲尔德型积分被表示成快速、绝对收敛的球面波函数系展开式;展开系数是以物性参数为复宗量的勒让德多项式. 该展开式数学物理意义明显,并且不受场点和源点的位置、媒质的物性参数和频率等条件的限制. 利用本文的结果可十分方便地计算和分析任意场点处的电磁场分布.     

Fully normalized spherical cap harmonics: application to the analysis of sea-level data from TOPEX/POSEIDON and ERS-1

We introduce two sets of fully normalized harmonics for the spectral analysis of functions defined on a spherical cap. The harmonics are the products of Fourier functions and the fully normalized associated Legendre functions of non-integer degree. Using Sturm-Liouville theory for boundary-value problems, we present two convenient and stable formulae for computing the zeros of the associated Legendre functions that form two sets of orthogonal functions. Formulae for the stable numerical evaluation of the fully normalized associated Legendre functions of non-integer degree that avoid the gamma function are also derived. The result from the expansions of sea-level anomaly from altimetry into Set 2 fully normalized cap harmonics shows fast convergence of the series, and the degree variances decay rapidly without aliasing effects. The zero-degree coefficients (Set 2) of sea-level anomaly from TOPEX/ POSEIDON (T/P) and ERS-1 indicate an El Niño event during 1993 January-1993 July, and a La Niño event during 1993 November-1994 July, although the ERS-1 result is less obvious. Ocean circulations over the South China Sea and the Kuroshio area are clearly identified with the low-degree expansions of sea-surface topography (SST) from T/P and ERS-1. A cold-core eddy of 4° in diameter centred at 17.5°N, 118 E was detected with the expansion of SST from T/P cycle 47, and a property of the cap harmonics is used to compute this eddy's kinetic energy. The kinetic energy is at a low in winter and high in summer, and its variation seems to be periodic with an amplitude of 0.4 m² S^-2.     

A discrete global grid system for earth system modeling

To support Earth system modeling, we propose a discrete global grid system that expresses multi-resolution spatial data. Specifically, a unified coding model that expresses a grid of nodes, edges, and cells is constructed for a triangular discrete global grid system. To fulfill the requirements of practical applications, we design a code-based topological query method for this grid system and an algorithm to transform between grid codes and geographic coordinates. We evaluate the Global Finite Volume Community Ocean Model (Global-FVCOM) on the triangular discrete global grid system in the proposed uniform coding model. The ocean tidal waves simulated by the Global-FVCOM running on the coded grid are then compared with results obtained using a traditional irregular spherical grid system, and the results display comparable accuracy. The uniform coding model proposed in this paper provides a triangular discrete global grid system that can represent multi-resolution spatial data and can be used in Earth system models. This unified coding model can also be applied to the geographic coordinate system made up of latitudes and longitudes, as well as diamond and hexagonal grids.     

500hPa气候异常高度场强度及谱结构的季节变化与半球际差异

用半球气候异常场强度及球函数分析改进方案，分析了北、南半球500hPa气候异常高度场的环流特征，得到如下主要结论：1)半球异常冬强于夏,北半球强于南半球。2)500hPa气候异常位势高度场具有低阶、低维的特征，它们主要由超长波及长波波段(0≤m、k≤6)的球函数构成。3)北半球500hPa气候异常位势高度场集的球函数谱结构较南半球复杂．夏季较冬季复杂。仅用超长波及长波波段的球函数拟合半球500hPa异常位势高度场，就可保证其有足够的精确度。     

小波检验方法在国内业务模式中的最新应用研究进展

本文主要介绍小波检验方法在国内区域业务模式中的最新应用研究进展,该方法被应用于模式降水预报的个例检验和统计检验,并针对两种解决二元域限制问题的技术方案进行阐述和讨论。结果表明：D8小波可以提供一种从不同降水率阈值和空间尺度对模式定量降水预报进行评估的方案,该方案首先对预报场和分析场进行阈值获得二进制误差场,再对二进制误差场进行尺度分解,从而可以评估模式在各个降水率阈值和空间尺度上的预报技巧和偏差。统计检验结果显示模式对较大尺度降水事件表现出较高的预报技巧评分,较小尺度降水事件和强降水的预报技巧评分相对较低,并且总体上呈现出预报降水事件多于观测降水事件的特征。Padding方案和Tiling方案对应的模式统计检验结果并无明显差异,由于后者不对原始降水场做任何改变,因此是一种更为稳健和可靠的方案。     

全球扰动的结构、演变及其球面Rossby波包表示

利用按谱函数展开的方法,计算和分析了正压球面Rossby波的连续谱和离散谱的波包结构、演变和传播特征,发现无论是结构、演变还是传播,连续谱波包和离散谱波包都很不相同.连续谱波包的传播速度不仅大小而且方向都随纬度变化;离散谱波包的传播速度只有大小随纬度变化,其方向在任何纬度上都是向西的.连续谱的波包和离散谱的波包在初始时刻虽然都表现为子午圈上不规则的波动,但随着时间的增加,二者的变化规律不同,连续谱的波包具有局部特征,主要表现在两个半球的中高纬地区,且随时间振幅迅速衰减;而离散谱的波包具有全球性质.将二者与整个扰动波包进行比较,发现:在中高纬度地区,连续谱波包的结构与扰动波包的结构相似,而在赤道地区,离散谱的波包与扰动波包相似.