中国地区集合预报产品自适应递减平均偏差订正法的改进研究

基于中国地区T213集合预报产品2 m温度预报数据,采用卡尔曼滤波类型的自适应递减平均法进行偏差订正处理,原方案在剧烈降温天气订正效果表现不理想。通过对递减平均参数w的重新构建得到改进的订正方案w(i,p)(i为站点信息,p为天气过程信息),在此基础上进一步优化对历史信息的有效提取,得到改进的方案w(i,p)相似法和w(i,p)统计法,并进行效果检验。结果表明:改进为包含空间和天气过程信息的函数w(i,p)后方案的订正效果得到不同程度的提高,其中24 h剧烈降温预报各成员预报均方根误差平均减小了0. 15℃;而进一步改进的w(i,p)统计法在当前几种剧烈降温预报中订正效果最优,其集合平均偏差与w(i,p)方案相比减小2. 54℃。     

Evaluation of grid-point atmospheric model of IAP LASG version 2 (GAMIL2)

The Grid-point Atmospheric Model of IAP LASG version 2 (GAMIL2) has been developed through upgrading the deep convection parameterization, cumulus cloud fraction and two-moment cloud microphysical scheme, as well as changing some of the large uncertain parameters. In this paper, its performance is evaluated, and the results suggest that there are some significant improvements in GAMIL2 compared to the previous version GAMIL1, for example, the components of the energy budget at the top of atmosphere (TOA) and surface; the geographic distribution of shortwave cloud radiative forcing (SWCF); the ratio of stratiform versus total rainfall; the response of atmospheric circulation to the tropical ocean; and the eastward propagation and spatiotemporal structures of the Madden Julian Oscillation (MJO). Furthermore, the indirect aerosols effect (IAE) is -0.94 W m-2, within the range of 0 to -2 W m-2 given by the IPCC 4th Assessment Report (2007). The influence of uncertain parameters on the MJO and radiation fluxes is also discussed.     

基于CMIP6强迫模拟分析人为气溶胶的气候效应(一)——介绍NUIST模式评估结果

为减少不同气候模式评估气溶胶气候效应的差异,第六次耦合模式比较计划(Coupled Model Intercomparison Project Phase 6,CMIP6)直接给定了人为气溶胶强迫数据。因此,有必要基于此强迫数据重新评估气溶胶气候效应。本研究首先将CMIP6给出的描述人为气溶胶强迫的模块引入南京信息工程大学(Nanjing University of Information Science and Technology,NUIST)的地球系统模式(The NUIST Earth System Model,NESM)。之后,利用NESM模式评估地球辐射收支平衡对此人为气溶胶强迫的响应,并分析模式模拟结果的不确定性。评估给出的人为气溶胶有效辐射强迫为-0. 45(±0. 28) W·m~(-2)。其中,气溶胶直接辐射效应为-0. 34(±0. 01) W·m~(-2),与第二次气溶胶比较计划(The second phase of Aerosol Comparisons between Observations and M odels,Aero ComⅡ)的评估结果基本一致;气溶胶对云辐射强迫的影响(包括半直接效应和间接效应)为-0. 10(±0. 30) W·m~(-2),明显受到模式内部变率的干扰,具有较大的不确定性。     

Constraints on rupture speed of the 2001 <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript> 8.1 West Kunlun Mountain Pass earthquake by co-seismic surface rupture slip displacements based on the slip-weakening mechanism with frictional undershoot involved

With co-seismic surface rupture slip displacements provided by the field observation for the 2001 M_S8.1 West Kunlun Mountain Pass earthquake, this paper estimates the rupture speed on the main faulting segment with a long straight fault trace on the surface based on a simple slip-weakening rupture model, in which the frictional overshoot or undershoot are involved in consideration of energy partition during the earthquake faulting. In contrast to the study of Bouchon and Vallée, in which the rupture propagation along the main fault could exceed the local shear-wave speed, perhaps reach the P-wave speed on a certain section of fault, our results show that, under a slip-weakening assumption combined with a frictional undershoot (partial stress drop model), average rupture speed should be equal to or less than the Rayleigh wave speed with a high seismic radiation efficiency, which is consistent with the result derived by waveform inversion and the result estimated from source stress field. Associated with the surface rupture mechanism, such as partial stress drop (frictional undershoot) associated with the apparent stress, an alternative rupture mechanism based on the slip-weakening model has also been discussed.     

1976年MS7.8唐山地震断层动态破裂及近断层强地面运动特征

采用美国南加州地震委员会（SCEC）Steven Day博士提供的三维有限差分断层瞬态破裂动力学模型（3D-FDM）,以1976年唐山M_S7.8地震为例,从简化的断层双侧破裂模式出发,对该地震发震断层的动态破裂过程及近断层地表运动特征进行了仿真模拟和计算.研究区域为围绕发震断层200 km×140 km×40 km（深度）的长方形块体组成,模拟计算的空间分辨率和时间分辨率分别为200 m和0.012 s,形成的空间网格节点数为1051×701×201.在DELL小型工作站上,我们实现了对源程序的移植和并行计算.同时,通过引进计算机可视化技术,对模拟数据进行了3D/4D解释分析.另外,在对源程序修改过程中,实现了对京津唐地区三维地壳速度结构的嵌入,在一定程度上增强了对地震波传播以及地面运动模拟的真实性,并讨论了地震破裂的方向性对近断层地表运动的影响.最后根据初步研究结果结合京津唐地区活动断层构造特征,对唐山M_S7.8级主震后随之而来的1976滦县M_S7.1级余震及宁河M_S6.9级余震的动态触发机制提出了新的解释.由于受主震破裂方向性作用的影响,使得主震对后续两个较大余震产生的动态应力变化的峰值在断层的走滑方向上较大,为2~3 MPa,在逆冲方向上较小,为0.1~0.2 MPa.即唐山主震的发生使得其周边的应力场有一个瞬态的应力调整,唐山主震对后续余震的发生有促发作用.     

沈阳市供水系统抗震功能可靠性分析

城市供水管网的抗震研究能够服务于供水管网的规划、设计和改造。以沈阳市供水管网系统为背景,介绍了一类求解大型供水网络节点及系统功能可靠度的方法。对沈阳市主干供水网络进行了抗震功能可靠性分析。研究结果证明,建议方法是合理评价供水管网地震后工作性能的工具,可以为大型供水管网的抗震可靠度分析及抗震设计服务提供基础。     

基于华北区域地震活动性分布的地震危险性评价模型

采用了分布式地震活动性模型. 该模型无需潜在震源区划分,同时简化了地震危险性概率分析方法. 根据破坏性地震目录建立了3个地震活动性模型,利用高斯光滑函数获得了华北区域内的a值分布特征,使用3种典型的衰减模型,分别计算了50年内超越概率10%, 5%和2%的地震动峰值加速度分布. 其分析结果显示了峰值加速度分布特征与我国第四代区划图大体一致,特定地震活动区（太原、 石家庄等地区）的峰值加速度略高于第四代区划图的结果,而这种峰值加速度分布特征与该地区较高的地震活动性特征是一致的. 概率危险性曲线结果表明,唐山、太原和北京等地区的潜在地震危险性比华北区域内其它城市高.      

自维持断层蠕滑过程和余震发生率的衰减特征—以汶川MW7.9地震序列为例

基于Dieterich地震活动性理论,本文推导出计算余震发生率和余震累积次数的一般表达式,其中主震后发震断层内部的剪切应力随时间的演化过程遵从Jeffreys-Lomnitz蠕变模型,且与修正Omori定律直接相关。修正Omori定律中的p值与震后断层的短时应力加卸载过程正相关。采用Rubin和Ampuero 给出的震后断层自维持蠕滑模型本文得出计算余震发生率的近似表达式,并对2008年汶川地震序列进行拟合。结果表明,p值的大小直接对应了速率-状态摩擦定律中摩擦参量b/a,而修正Omori定律中的c值则与速率-状态摩擦定律中的临界滑移D_c相关。对于汶川余震序列而言,拟合结果显示b/a约为1.13,D_c约为2—3 cm。Rubin-Ampuero震后自维持蠕滑描述了震后孕震层内部短暂的速率变化特征,是孕震断层演化过程不可缺少的环节。      

低频电场变化探测阵列建设及其初步运行结果

在低频(LF)和甚低频(VLF)频段探测全闪活动是近年来闪电探测技术发展的一个重要方向.中国气象科学研究院于2014~2015年在广州地区建设了由9个快电场变化测量仪构成的低频闪电电场探测阵列(Low-frequency E-field Detection Array,LFEDA).本文介绍了LFEDA构成以及针对闪电低频电场信号建立的闪电定位算法,基于蒙特卡罗法和人工触发闪电试验对LFEDA的定位精度和探测效率进行了理论模拟和客观评估.前者表明站网中,具有较长基线的定位站有助于改善垂直于该基线方向的水平定位精度和沿着该基线方向的高度定位精度.后者表明LFEDA对触发闪电事件和回击的探测效率分别为100%和95%,回击平面定位误差平均值为102m.定位结果与雷暴结构的对比分析表明,闪电脉冲放电事件与反射率柱所表征的对流特征较强区域有较好对应,脉冲放电事件在垂直方向上呈现合理的分布.此外,通过单次云闪、地闪的定位结果发现,LFEDA能够描绘闪电发展形态,基于三维定位结果获得的云闪初始阶段垂直发展速度以及地闪先导垂直发展速度变化趋势与之前观测结果类似,研究结果表明,LFEDA具备了全闪电三维定位能力,为闪电发展特征研究以及雷暴电学研究提供了新的技术手段.     

近断层应力变化快速计算方法及对强余震空间分布的指示意义:以M_w7.9汶川地震为例

构造地震一般由断层摩擦失稳所致.断层内部及周边所累积的剪切形变则通过同震滑动位移得到局部释放.因此,震后断层内部及近断层周边的静态剪切应力变化量的空间分布可通过断层面上的滑动位移分布计算得到.本文采用傅氏变换方法(FTM)计算单一有限断层同震滑移场所形成的静态剪切应力变化分布,近断层两侧的应力变化计算可由波数域内应力近似算法获得.结果表明,FTM快速有效、易于实现,有效地避免了常规应力计算中奇异值的出现.以2008年Mw7.9中国汶川大地震为例,采用前人所得有限断层滑动位移模型,得到了断层面和近断层周边准3D剪切应力分布解,并同主震后中强余震的空间分布特征作了比较.结果表明,大部分的中强余震震源位置处于剪切应力变化值为正的区域,由同震位移所产生的局部应力降峰值和均值大小同应力变化的正值大致相当,从而表明了快速且有效地计算断层内部及近断层附近的应力变化分布可以为主震后强余震发生的潜在区域提供指示意义.需要强调的是,应力变化空间分布特征的获取强烈地依赖于有限断层滑移模型解.有关滑动位移反演解的多解性对应力变化计算结果的影响,本文作了必要的讨论.