中国土壤温度的空间插值方法比较

利用中国698个气象站点1971~2000年的地面气候资料,采用三种不同方法预测中国0cm、20cm和40cm深度年均土壤温度的空间分布,其中普通克里格和泛克里格法直接以年均土壤温度数据为源数据、回归克里格法以中国年均气温数据和中国DEM数据为源数据进行预测。预测结果的准确性通过平均绝对误差(MAE)和均方根误差(RMSE)值来评价。结果表明回归克里格法预测的MAE值和RMSE值均为最小,说明其预测结果的准确性最好、预测的极端误差也最小;其次为泛克里格法;普通克里格法预测的效果最差。回归克里格法预测结果由于采用了中国DEM数据进行修正,在空间特征表达方面能够更好地表达复杂地形地区的局部变异,其平滑效应明显小于泛克里格法和普通克里格法的预测结果。     

华北地区烈度衰减模型建立及其用于震中区域和震级的定量估算.

利用地震烈度资料定量估算历史地震震中区域和震级的方法是由Bakun和Wentworth于1997年首先提出的．该方法定量程度高,对烈度数据较少或发生在近海的历史地震的定位和震级估算尤为有效．按照Baku和Wentworth给出的思路,笔者尝试了对我国华北地区的历史和近代地震的震中和震级进行初步分析．首先,选取20世纪以来发生在该地区的10次有仪器记录的地震（5.3le;MSle;7.8）,对该地区的烈度-震级-震中距衰减关系进行标定并给出了烈度衰减模型,表明华北地区烈度随震中距增大而衰减的速率明显小于美国加州地区（约50%）．在此衰减模型的基础上,提出了确定震级和震中区域的网格搜索试算方法（GSTSL）,并给出了适用于华北地区的圈定震中区域和烈度震级的等值线置信值．最后,讨论了计算震中区域等值线时所引进的权因子Wi及其中参数b对震中区域等值线圈闭性的影响．利用该方法,对发生在1668年莒县——郯城地震,1679年三河——平谷地震,以及1966年隆尧地震和1969年渤海地震进行了分析．需要指出的是,该方法也可推广应用于我国其它历史地震资料丰富的地区．      

A uniform representation in the calculation of seismic radiation efficiency and rupture speed

According to energy partition principle,we have derived a generalized formula used to calculate the radiation effi-ciency.It overcomes the numerical errors involved in the estimation of radiation efficiency caused by the uncer-tainties owe to the different frictional mechanisms which is used to describe the earthquake faulting.Using real seismic data published in the previous studies,we have recalculated the radiation efficiency with associated to the different frictional mechanisms,and compared with previous results lack of consideration of frictional mechanism.The results indicate that,based on the new formula,the radiation efficiency we have estimated is always less than 1,and the unphysical conclusion of radiation efficiency greater than 1 is avoided completely.     

Constraints on rupture speed of the 2001 M S 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 MS8.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.     

土坡有限元稳定分析若干问题探讨

详细探讨了用有限元方法分析土坡稳定所遇到的问题及处理方法 ,主要包括: 有限元前后处理 ;有限元计算范围及网格划分对计算结果的影响 ;安全系数的定义 ;最危险滑动圆心位置及临界滑动面的形状等问题。     

亚像素细分在光学自准直系统观测中的应用

介绍一种基于CCD的光学自准直系统观测方法及相关的数字图像处理方法。关于图像处理,首先采用自动阈值分割提取十字丝的像素级边缘,然后使用积分法找到十字丝的像素级边缘位置,接着运用三次多项式拟合出图像边缘实现十字丝边缘的亚像素定位,最后对十字丝亚像素边缘进行直线拟合得到边缘直线方程,进而得到十字丝中心。实验结果表明,亚像素细分技术能够较好地提取十字丝亚像素边缘。     

天津地区一次大到暴雪天气的诊断分析

应用常规探测资料和NCEP再分析资料, 对2011年2月下旬典型华北回流形势下天津地区一次大到暴雪天气进行了诊断分析。结果表明：回流降雪过程中,华北上空西风环流以纬向型为主,冷空气主体偏北,主要影响系统为华北回流冷高压和低压倒槽。同时,回流降雪中有浅薄的冷空气垫,其上有暖湿气流在爬升,爬升高度大约为650 hPa。回流降雪期间有来自西南和东北两个方向的水汽在天津地区交绥,西南方向的水汽较为暖湿,东北方向的水汽相对干冷,低空和超低空为一致的东北气流,900 hPa附近有超低空急流,700 hPa以上为西南暖湿气流。降雪过程中对流层低层到高层均为一致的强上升运动,上升高度可达200 hPa,对应于低空和超低空有强的辐合。降雪开始前天津及其周边地区有较强的对流不稳定能量和对称不稳定性,有利于对流的发展。     

基于经验模型和物理模型研究2013 M_S7.0芦山地震余震序列

基于中国地震台网中心2013 MS7.0芦山地震余震数据我们首先确定了余震空间分布范围并根据G-R关系计算了主震后半小时内的完备震级Mc=3.5,并且得到了ML≥3.5和ML≥3.0的地震在2001年至芦山地震前的背景场地震发生率.通过Omori-Ustu经验定律和两种Dieterich模型对芦山地震余震发生率的拟合,我们发现阶梯型Dieterich模型只能模拟p=1的情况,从而造成了模拟曲线与观测数据的差别;前人研究表明震后滑移同样是产生余震的原因,如果假设余震序列由主震静态剪应力Δτ和震后滑移共同作用所产生,我们数值模拟得到的对数型Dieterich模型能够较好地推断余震发生率R随时间t增加而衰减的趋势,能够从物理机制上解释MS7.0芦山地震余震序列衰减指数大于1这一现象.通过对数型Dieterich模型的拟合并结合Andrews的方法,我们还得到MS7.0芦山地震Aσ约为0.155 MPa,ta约为8.4年,这一值与前人研究结果十分接近.     

与特征断层结合的地震活动空间光滑模型在衡水市活断层地震危险性评价中的应用

在衡水市活断层地震危险性评价中,应用研究区域的地震活动性数据建立恰当的分布式地震活动性空间光滑模型,并结合断层对周边地震危险性的影响,计算前磨头断裂和衡水断裂在不同超越概率下上限震级分别为6.5、6.0级的区域基岩峰值加速度。研究发现,研究区域内50年超越概率为63%的基岩峰值加速度略有增加;50年超越概率为10%、5%的基岩峰值加速度在断层附近有所增加。研究表明在较低超越概率的情况下,断层周围如深州市、冀州市和衡水市辖区的地震危险性大于其他地区;利用仪器记录资料得到的地震活动性模型并结合断层资料计算出的地震危险性结果能够反映衡水地区现今的地震活动水平和危险性水平。     

FY2G卫星新一代高能带电粒子探测器观测数据分析

风云二号系列卫星是我国开展动态空间天气事件和空间环境监测及预警业务的重要观测平台,各系列星上均安装有高能带电粒子探测仪器开展卫星轨道空间带电粒子辐射环境连续实时的动态监测.FY2G卫星于2015年1月发射,星上采用了全新的高能粒子探测器,包括:一台高能电子探测器可监测200keV-4 MeV的高能电子,一台高能质子重离子探测器可监测4~300 MeV的高能质子,从而实现对带电粒子更宽、更精细能谱的监测.本文给出了FY2G高能带电粒子探测器在2015年1月至2015年10月期间几起典型的带电粒子动态观测结果,结合太阳和地磁活动相关参数,对高能带电粒子通量在亚暴、磁暴和太阳爆发等扰动影响下细节变化过程和特征作出了较为详细的分析描述,展现了FY2G卫星高能带电粒子探测器对轨道空间粒子环境动态变化的准确响应能力,表明观测数据可开展更加精细的轨道粒子环境评估.针对FY2G高能带电粒子探测结果进一步开展了与GOES系列卫星同期观测的比对分析,结果反映出在较小的扰动条件下多星观测到的带电粒子响应和通量变化可基本趋于一致或保持相对稳定的偏差,而扰动条件的显著变化会加大多星观测带电粒子响应和通量变化的差异,这些结果可为今后开展多星数据同化应用提供参考,也为发展磁层对扰动响应的更加复杂的图像提供了新的可能.