变分同化结合广义变分最佳分析对微波散射计资料进行海面风场反演

针对微波散射计资料的海面风场反演,本文首先利用QuikSCAT卫星散射计的近实时资料,借助荷兰皇家气象学会的SeaWinds散射计数据处理模式,在不同风速和海域的情况下,采用多解方案结合二维变分同化方法进行风向的模糊去除,达到消除解的不唯一性的目的,且对4种反演结果（低速风、中速风、高速风和近海岸风）与美国国家海洋和大气管理局卫星资料中心的近实时风场数据的对比分析,说明了该技术的有效性.在此基础上,对反演得到的中速风场利用广义变分最佳分析方法进行了调整,调整后的风场结构有了进一步的改善.此方法可为我国散射计资料反演海面风场的研究提供新的思路和方法.     

合成孔径雷达资料反演海面风场的正则化方法研究

针对基于散射计地球物理模型函数的合成孔径雷达资料反演海面风场中背景风向存在很大误差的问题,提出了利用正则化方法反演海面风场,并通过L曲线准则对正则化参数进行最优选择,然后对该方法进行仿真试验.试验结果显示,正则化方法反演海面风场的风速和风向的反演精度均优于变分同化方法,分析不同真实风向情况下的反演精度表明,真实风向在45°和135°附近时,风速反演的精度高,真实风向在0°和90°时,风向反演的精度高,这种现象与地球物理模型函数中后向散射截面相对于风向的导数相关. 关键词： 合成孔径雷达 正则化 L曲线准则 变分同化     

散射计资料对台风海平面气压场的反演和定位的新方法研究

文章首先利用QuikSCAT散射计获取的海面风场资料,利用华盛顿大学行星边界层(UWPBL)模式对台风海平面气压场进行反演(称为直接反演),并对其强度进行了分析;然后对获取的海面风场,采用变分方法对风场进行分解,分解出无旋和有旋风场,对分解出的有旋风场,再利用UWPBL模式对魔羯台风的海平面气压场进行了反演,对所得结果与通常直接反演、NCEP海平面气压场及台风报文资料进行一一对比,结果表明:新的反演定位的方法对台风海平面气压场中心气压的定位精度有了进一步的提高,此种定位的方法为散射计资料在台风中的实际应用开辟了一条新的思路. 关键词： 散射计风场 变分方法 台风 海平面气压场     

多参数背景场误差模型在散射计资料台风风场反演中的应用

利用散射计资料反演海面风场时, 台风区域普遍存在降雨使得风场反演误差很大, 引入降雨地球物理模型函数 (GMF+Rain) 及多解方案 (MSS), 结合二维变分(2DVAR) 模糊去除思想风速反演误差很大程度减小, 但风向反演误差仍有待进一步改善, 如何进一步减小风向反演误差有待进一步研究. 文章介绍了2DVAR模糊去 除方法的基本思想, 针对背景场误差较大时, 2DVAR模糊去除风向误差较大, 引入包含若干参数的背景场误差模型. 基于台风个例数值试验结果, 着重从理论分析角度讨论各参数关于2DVAR模糊去除效果的敏感性, 进而提出最优参数设置方案以改善风向模糊去除效果. 2006年“摩羯”台风QuikSCAT数据风场反演数值试验结果结合理论分析表明: 引入多参数误差模型, 通过设置粗糙误差概率等于0, 2DVAR风向模糊去除效果明显改善; 同时, 背景场的影响可通过增大背景场误差方差, 减小背景场误差相关尺度和减小粗糙误差概率而减小, 进而减小在背景场误差较大情况下的风向反演误差. 关键词： 台风风场反演 二维变分 多参数误差模型 散射计资料     

合成孔径雷达反演海面风场变分模型分析

为考察合成孔径雷达反演海面风场变分模型精度,开展了误差分析试验.在背景场误差为极大值条件下分析场误差低于背景场误差,且随背景场风向增大呈周期性变化;在背景场误差逐渐变大条件下分析场误差逐渐增大,误差偏离方向与背景场误差偏离方向一致;在特定背景场条件下与直接反演模型相比,低风速时反演误差略高,中高风速时反演误差明显较低.总体来讲,变分模型误差小于背景场误差,风速反演误差小于1.60 m/s,风向误差小于17.15°.     

凑合反推法──流体力学变分原理建立的一条新途径

凑合反推法是刘商联系统方法[1]的进一步发展，应用这种方法可以方便地构造各种亚广义变分原理及广义变分原理，并可以消除临界变分现象．对于任何二维守恒型流体力学方程，作者推导得到了其广义变分通用公式．几个实例证明这种方法是有效的、简单的，并具有普遍的意义．     

GPS掩星一维弯曲角算子在四维变分资料同化系统中的实现方法研究

GPS无线电掩星资料的良好特性使其作为十分重要的观测资料在欧洲 及美国数值天气预报中效果明显. 本文设计了一维弯曲角算子在四维变分资料同化系统中的实现方法和质量控制策略. 详细阐述了一维弯曲角正演算子及其切线性/伴随算子在全球四维变分资料同化系统中的实现步骤和过程, 其中折射率计算中采用Rüeger的折射率最佳平均因子的三项表达式, 计算出更加精确的折射率值. 对2012年8月份的GPS无线电掩星资料进行了观测误差统计, 引入了更加符合四维变分资料同化系统实际情况的一维弯曲角观测误差模型和质量控制策略. 同时, 对比了观测误差统计试验中同化GPS无线电掩星资料的模式背景场和未同化GPS无线电掩星资料的背景场的观测误差统计结果, 发现同化GPS无线电掩星资料之后的观测误差统计结果在全球区域正作用十分明显, 尤其是在中低纬度地区, 高度在10 km 以上的区域观测误差减小十分明显. 关键词： GPS无线电掩星 一维弯曲角算子 四维变分同化系统 质量控制     

优化模式物理参数的扩展四维变分同化方法

数值模拟的一个重要误差来源是模式物理参数,为提高模拟准确率,如何改进模式物理参数是亟需解决的问题.本文对经典四维变分同化技术进行了改进,提出了一种新的利用观测资料来同时优化模式初始场和物理参数的扩展四维变分同化方法,并以Ekman边界层模式和Lorenz模式为例进行了数值试验.结果表明,利用本文提出的新方法,通过对观测资料的变分同化,可以在实现对模式初始场进行优化的同时,纠正了模式物理参数中的误差,从而有效提高了模式的模拟准确率.该方法对于改进数值模式物理参数有着重要的促进意义.     

C波段机载合成孔径雷达海面风场反演新方法

针对基于散射计地球物理模型函数的机载合成孔径雷达(SAR)海面风场反演中存在的风向获取依赖于图像风条纹或数值预报、 散射计数据和浮标等背景场资料, 风向与SAR图像时空分辨率不匹配, 进而影响机载SAR海面风场反演精度等问题, 本文根据机载SAR对海探测特点, 研究一种适用于C波段机载SAR的海面风场反演新方法. 利用SAR图像距离向不同入射角的后向散射系数, 依据地球物理模型构造最小代价函数, 通过代价函数的求解直接从机载SAR数据同时反演出海面风速和风向. 利用论文提出的海面风场反演方法分别对仿真SAR数据和实测C波段机载SAR数据进行风向、 风速的反演误差分析及试验验证.研究结果表明, 该方法适用于机载SAR海面风场反演, 可不依赖背景风向直接反演出精度较高的风速和风向; 雷达后向散射系数误差是决定海面风速、风向反演精度的关键因素, 辐射定标精度越高则反演误差越小; 海面风速反演误差随着风速的提高而增大, 当海面风速大于18 m/s时, 风速反演误差显著增加, 而海面风向的反演误差与风速无明显关系. 关键词： 机载合成孔径雷达 海面风场 多入射角     

基于区域逐步分析的集合变分资料同化方法

集合变分数据同化方法的同化效果对集合样本容量具有很强的依赖性,研究发现此问题的出现是因为其计算过程中分析增量被表示为集合扰动向量或其展开正交基向量的线性组合.这样的处理方法虽然避免了计算梯度而引入伴随模式,但是因为物理控制变量个数远大于集合样本容量,就会导致物理量的同化分析值对集合样本容量很敏感.根据此原因,提出了区域逐步分析方法,减小了同化分析区域内物理变量个数与集合样本容量数之间的比值,使问题得到解决.利用浅水方程模式进行资料同化数值试验表明,基于区域逐步分析的集合变分资料同化方法可以得到较好的结果,能明显提高同化的精度.     

基于星载激光雷达数据的海面风速探测

使用新版(Version 3.01)星载激光雷达CALIPSO Level 1剖面数据和Level 2的大气气溶胶光学厚度(AOD)数据反演全球海面风速。将AOD数据用于大气双程透射率的校正。对CALIPSO所采集的2007年1月、4月、7月和10月532nm波长无云激光测量数据进行应用计算,采用532nm退偏比经验修正海面白帽及水下次表层水体散射的影响,同时选用准同步的AMSR-E海面风速测量值为真值进行对比。对应2007年的4个月,由CALIPSO星载激光雷达532nm单脉冲测量数据反演海面风速的标准偏差分别为1.24、1.24、1.24、1.20m/s,由5km滚动平均数据反演海面风速的标准偏差分别为0.98、1.02、0.98、0.94m/s。结果表明,采用的数据反演方法可行且数据利用率得到提高。     

Polarimetric scattering from a two-dimensional improved sea fractal surface

This paper studies the influence of wind parameters and fractal dimension from an improved two-dimensional sea fractal surface on the polarimetric scattering by using facet integration.A two-dimensional improved sea surface simulated is discretized into three matrices of sea surface facets including a height matrix and two slope matrices on orthogonal directions.Based on the Kirchhoff approximation,the polarimetric scattered field is derived in the Cartesian coordinate system by integration of three matrices mentioned above.Finally,the fully polarised radar cross section is numerically simulated and the dependence of the polarimetric scattering on the sea fractal surface,such as the wind speed,the wind direction,as well as the fractal dimension,is discussed in detail.     

高度计风速与辐射计风速的变分融合研究

为实现高度计风速与辐射计风速的融合, 采用Kriging插值方法将辐射计风速插值到高度计路径得到高度计风速观测算子, 建立高度计风速与辐射计风速融合的代价函数; 利用变分方法求解分析风速, 得出融合结果. 分别在仅高度计风速有误差、 仅辐射计风速有误差和两者都有误差的情况下开展了模拟试验, 结果表明变分融合结果更接近于理想风速, 尤其在高度计路径上更为明显. 选取Jason-1高度计第241周期风速资料和与其时空匹配的国防气象卫星计划 F17上所搭载的专用传感器微波成像仪辐射计风速资料开展了实例试验, 结果表明融合后的结果更加接近于浮标观测结果, 证实了变分融合方法的有效性. 经统计, 上述高度计和辐射计资料有60%时空匹配, 开展风速融合有重要的理论及应用价值. 该变分融合方法可应用于我国海洋二号卫星所载高度计与辐射计风速反演结果的融合. 关键词： 高度计 辐射计 风速 融合     

Independent Plaquette Trial Action for 4-Dimensional Lattice Gauge Theory

Based on the explicit expressions of the plaquette formulations, the independent plaquette trial action for 4-dimensional lattice gauge theory is introduced. As an example, the mean plaquette energy EP for the SU(2) lattice gauge theory is calculated by using action variational approach with the independent trial action. The results are in good agreement with the Monte Carlo results in the strong coupling and the crossover region, and the curve is smooth in the whole region, which show that 4-dimensional SU(2) theory has only a single, confining phase. The unwanted discontinuity of EP given by the single link trial action, which is used in the earlier variational calculations has been avoided.     

Two-dimensional electron systems in inversion layers of p-type Hg0.8Zn0.2Te metal-insulator-semiconductor structures

Strong oscillations on capacitance and conductance have been observed in p-type Hg0.8Zn0.2Te metal-insulator-semiconductor structures, made by using a recent process for the interface passivation. This behaviour is attributed to a two-dimensional electron gas in the n-inversion layer and the variation of the conductance maximums with temperature indicates that the dominant perpendicular transport mechanism for electrons is an incoherent two-step tunnelling through deep levels in the gap. Three models have been used to describe the quantum confinement: the simple variational method, the triangular potential approximation and the propagation matrix method. The later approach takes into account the non parabolicity of the conduction band structure and uses a finite height barrier at the insulator-semiconductor interface. A very good agreement between experimental and calculated values for the two lowest subband energy is obtained. Received 9 February 1999     

聚合物脊形光波导的变分有效折射率法分析

聚合物脊形光波导是聚合物集成光电子器件的重要构成单元.利用有效折射率法计算聚合物脊形光波导的横向折射率分布及有效折射率,将各区域中的光场分布近似用分段函数表达.基于导模满足的标量波动方程,利用变分法确定变分参量,以求得准确的横向光场分布.对聚合物脊形多模光波导基模和高阶模的色散特性与横向场分布进行分析,研究了波导结构参数对色散特性的影响,计算出TM基模和高阶模的光场分布,得出了聚合物脊形光波导的单模传输条件.研究表明,该方法计算量小、精度高,对聚合物光电子器件中脊形光波导的理论分析与设计优化提供了简单高效的方法.     

降雨对C波段散射计测风的影响及其校正

传统观点认为C波段散射计工作波长大于雨滴直径,受降雨散射衰减的影响很小, 因此往往忽略降雨对C波段散射计测风的影响.本文基于降雨引起C波段散射计信号的衰减、 后向体散射及雨滴落入海面后的扰动作用,推导了降雨条件下的雷达方程, 构建了2010年全年的ASCAT散射计、降雨雷达和欧洲中期天气预报中心数值预报的匹配数据集, 定量分析了降雨对C波段散射计测风的影响,发现其信号衰减随降雨强度和入射角的增大而增强; 后向体散射和雨表面扰动作用随降雨强度的增大而增强、随入射角的增大而减小, 其中雨表面扰动作用对散射计测风的影响大于后向体散射.另外, 利用降雨条件下的雷达方程和匹配数据集,本文建立了降雨条件下的C波段主动微波辐射传输模型, 实验表明,该模型能够改善降雨条件下C波段散射计测风的精度.