分布式降雨径流物理模型的建立和应用

根据流域降雨径流的主要过程,考虑流域气象及下垫面要素的空间异质性,建立了具有物理基础的分布式降雨径流模型。模型将流域离散为栅格计算单元,并按水流特性分栅格单元为坡面单元和河网单元。在坡面单元上主要计算降雨、下渗、坡面流、壤中流等水文过程,而河网单元则主要计算河道汇流过程。模型利用空间权重插值方法将雨量站点的降雨量插值到各个计算单元,采用运动波方程来计算坡面流,将壤中流概化为垂向流和侧向流,分别用Green-Ampt公式和运动波方程来模拟,河道汇流也采用运动波方程。模型结构简单、参数的物理意义明确,大多数参数可利用DEM、土壤类型图、植被类型图直接获取,少数敏感参数通过率定确定。模型在浙江省甬江上游黄土岭流域和皎口流域进行了应用和检验,其结果令人满意。     

基于DBSCAN与网格搜索的雷电定位算法

定位精度是评价雷电定位网络的重要指标之一,定位算法直接影响雷电探测结果的精度。雷电监测系统探测数据误差不可避免,传统定位算法不具备抗误差干扰能力,迭代计算易发散,定位结果精度不高。为了满足实际应用需求,提出一种新的雷电定位算法DG-LLA（DBSCAN and grid-search lighting location algorithm）,在定位计算中引入DBSCAN（density-based spatial clustering of applications with noise）方法与网格搜索方法。通过仿真与国家雷电监测网实际定位结果对比分析定位算法性能。结果表明:到达时间差（time difference of arrival,TDOA）法和Taylor级数展开法定位误差较大,仿真区域的均方根误差分别为982 m和668 m;定位中引入DBSCAN方法后,均方根误差明显减小为406 m,引入DBSCAN方法和网格搜索方法后,均方根误差减小为349 m;在相同回击数据条件下,算法DG-LLA与国家雷电监测网相比定位数量更多,回击数据的利用率从43.4%提升到51.5%,新增定位结果周围雷达回波特征较强,定位精度更高。     

计算最小走时和射线路径的界面网全局方法

用慢度分块均匀正方形模型将介质参数化，仅在正方形单元的边界上设置计算结点，这些结点构成界面网．根据Ｈｕｖｓｅｎｓ和Ｆｅｒｍａｔ原理，由不断扩张、收缩的波前点扫描代替波前面搜索，在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似，通过比较确定最小走时和相应的次级源位置，记录在以界面网点位置为指针的３个一维数组中．借助这些数组通过向源搜索可计算任意点（包括界面网以外的点）上的全局最小走时和射线路径．这一方法不受介质慢度差异大小限制，占内存少，计算速度较快，适于走时反演和以Ｍａｓｌｏｖ射线理论为基础的波场计算．     

结合边界条件的截断高阶隐式有限差分方法(英文)

有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。     

试验模型位移场中的几种光学测量方法比较

详细地分析了几种光学测试方法在模型试验中的使用情况。自动网格法在模型表面布置网格点容易,由于开挖或裂缝等原因去掉网格点,改进后的程序仍能自动编码,必要时可通过人机对话手工编码。散斑互相关法,在模型表面布置散斑较容易,改进后的程序在开挖情况下也能自动处理,但在散斑场遭受严重破坏情况下无法比较相关性。其他光学测试技术由于受试验条件的限制,而不适合使用。     

基于ANUSPLIN的京津冀区域逐日气温格点数据集建立方法研究

为建立一个高精度、高空间分辨率的逐日气温格点数据集,满足公共气象服务对于精确信息及实时信息的需要,利用2018年6—8月京津冀区域以及临近省区共3 974个国家级及区域气象观测站质控后的逐日气温资料,采用ANUSPLIN软件对逐日气温数据进行空间内插,得到了京津冀区域逐日气温格点数据集(0.01°×0.01°),并分别利用反距离权重插值法、普通克里金插值法、样条函数法对逐日气温数据进行空间插值,采用相关系数(Corr)、平均绝对误差(MAE)、平均相对误差(MRE)等作为评估指标来检验插值精度。结果表明:1) ANUSPLIN软件满足了空间插值对精度及曲面平滑度的要求,能直观体现京津冀区域气温由北向南递增的空间分布特征; 2) 4种插值方法中,基于ANUSPLIN软件的插值结果最优,相关系数平均达0.97,其样本误差在1℃之内占比为90.59%,MAE为0.46℃,MRE为1.81%; 3)插值误差较大的区域位于冀北高原、燕山丘陵及太行山脉一带,高海拔、低站点密度等是造成插值误差的主要原因。基于ANUSPLIN插值方法建立的逐日气温格点数据集具有分辨率高、空间插值误差小的优势,ANUSPLIN对气温的空间分布具有较好的预测能力。     

基于百度地图的精细化格点预报显示

精细化格点要素预报是目前中国气象局的主推业务和未来天气预报的发展方向,如何有效地将预报产品提供给用户使用,是精细化格点预报的最终环节。本文介绍了基于百度地图API的陕西精细化格点预报显示系统,系统主要实现了:(1)降水预报,以CMORPH(NOAA Climate Prediction Center Morphing Method)卫星与全国3万个自动观测站的逐时降水量融合资料为基础,通过动态偏差订正的方法来提高格点降水的预报能力;温度预报,采用滑动线性回归的方法来改善温度预报效果。(2)任意位置的地理信息获取及对应格点240h预报时效的气象要素实时展示。(3)格点气象要素向站点转换,通过格点值提取全省98个观测站逐3h站点预报值,实时分析过去24h降水、温度预报与观测值的误差,供用户预判未来预报值的可能误差趋势;并提供未来168h逐日要素预报。(4)与以往的数据库后台支撑不同,本系统直接将3.5GB格点预报数据一次性读入内存,进行侦听,解决了数据库检索、调用效率低下的问题。     

MM5模式中不同对流参数化方案的比较试验

应用MM5中尺度模式, 在60、20和10 km模式分辨率下, 分别选用4种不同对流参数化方案 (KUO方案、GRELL方案、KAIN-FRITSCH方案和BETTS-MILLER方案, 以下简称KU、GR、KF和BM方案) 对1996年8月3~4日石家庄暴雨过程进行数值模拟试验.模拟结果的对比分析及其与观测的比较表明:主要雨带位置对参数化方案并不是十分敏感, 但随分辨率提高, 雨带分布特征的模拟更接近实况; 当分辨率提高到10 km时, 虚假的降水中心也明显增加; 模拟的暴雨中心强度随分辨率的提高而增强并随参数化方案的不同有所变化, 但均比实况偏弱.分析还发现, MM5模式的GR、KF及BM方案的次网格降水对总降水的贡献率随分辨率的提高而减小, 而KU方案的情况则呈现出不合理的缓慢增加态势.虽然4种方案下模拟的水平环流特征有较好的一致性, 但模拟的云物理特征和垂直运动特征还是存在一定差别的, 这种差别对定点、定量降水和天空状况、地面气温、湿度等要素的准确预报都会产生影响.因此, 在预报和模拟中应考虑预报和研究对象的特点来选择对流参数化方案.     

青藏高原地表分类及其东北部辐射通量的二维数值模拟

通过对青藏高原地区地表特征的格点分类及沿五道梁站的经向二维中尺度数值模拟,结果发现：青藏高原地区１３类地表特征都有,尽管其植被自西北向东南呈带状分布,但仍十分复杂;而其地表特征的这种格点分类为该地区的高分辨数值模式模拟提供了极其有用的下边界条件。地面净辐射通量及其诸分量的模拟结果与观测结果的相当接近与一致表明,该模式具有模拟青藏高原地区地面通量日变化的能力。模拟结果和观测分析还指出：夏天睛天条件下     

网格嵌套技术对一次中尺度对流系统降水过程模拟的影响

利用非静力中尺度模式MM5对2002年7月22日12:00~23日12:00(世界时,下同)长江流域的一次梅雨锋暴雨过程进行数值模拟试验,主要讨论了网格嵌套技术对降水和中尺度对流系统的影响。结果表明:三重嵌套在D1,D2域选用积云参数化方案后,模拟的雨区收缩,虚假降水中心相对减少,降水强度及分布更接近观测值。在模式非线性动力、热力及湿物理过程共同驱动下,通过嵌套网格的双向相互作用,使可分辨云尺度的细网格域D3将其信息通过嵌套边界向次细网格域D2传递,然后再通过D2域边界向粗网格域D1域传递。同样,动力、热力反馈也会反向进行。结果将有助于改进各网格域的预报效果。但对D1网格域系统位置及其发展演变过程的影响相对小些;另外,通过双向多重嵌套,可提高模式预报区域的分辨率,特别是提高模式关键预报区域的分辨率,这也就有可能改进预报的中尺度物理场,使其能够较真实地描写大气实况。