Trimmed L-moments (1,0) for the generalized Pareto distribution

Abstract

Statistical analysis of extremes is often used for predicting the higher return-period events. In this paper, the trimmed L-moments with one smallest value trimmed—TL-moments (1,0)—are introduced as an alternative way to estimate floods for high return periods. The TL-moments (1,0) have an ability to reduce the undesirable influence that a small value in the statistical sample might have on a large return period. The main objective of this study is to derive the TL-moments (1,0) for the generalized Pareto (GPA) distribution. The performance of the TL-moments (1,0) was compared with L-moments through Monte Carlo simulation based on the streamflow data of northern Peninsular Malaysia. The result shows that, for some cases, the use of TL-moments (1,0) is a better option as compared to L-moments in modelling those series.

Citation Ahmad, U.N., Shabri, A. & Zakaria, Z.A. (2011) Trimmed L-moments (1,0) for the generalized Pareto distribution. Hydrol.Sci. J. 56(6), 1053–1060.     

雷达测雨误差及其对淮河流域径流模拟的影响

将雷达测雨数据与分布式水文模型相耦合进行径流过程模拟,分析雷达测雨误差及其径流过程模拟效果,研究雷达测雨误差对径流过程模拟的影响效应.在对淮河流域气象中心业务化的5种淮河流域雷达测雨数据进行误差分析的基础上,采用雷达测雨数据驱动HEC-HMS水文模型,模拟分析淮河息县水文站以上流域2007年7月1-10日强降雨集中期的径流过程.结果表明:利用雷达测雨数据的径流模拟结果与实测资料的模拟结果基本吻合,各种雷达测雨数据误差经过HEC-HMS水文模型传递后,误差明显减小.联合校准法对应的模拟效果最好,过程流量相对误差NBs'和洪峰流量相对误差Z'分别为-20.2%和-13.3%.     

华北地区强震前的信号震及其预测意义

根据对华北地区1970年以来M_S≥6地震以前中小地震活动空间图像变化特征的分析, 研究了“信号震”发生的时空特征及其地震活动背景, 由“场-源”关系特征对一般地震进行严格的筛选识别, 从而得出信号震的有关预测指标。 信号震一般发生在强震前的2年之内, 多数发生在15个月内; 信号震与强震的距离不超过200 km, 多数在100 km之内; 震级强度一般为M_L4.0～5.3。 信号震一般发生在局部的M_L≥4.0地震平静区内, 一般发生在中小地震条带上或条带附近, 在其周围或附近存在中小地震空区。 检验结果表明, 信号震发生后的9个月之内, 其预测概率P_t即超过0.5, R_t值达到0.27; 预测区域半径在距信号震震中100 km之内时, 其发生概率P_d可以达到0.73; 预测强震震级一般为M_S≥6.0。 研究表明, 信号震的环境应力值τ₀明显高于其他地震, 显示了高应力背景的异常显著性, 它所辐射的地震波中很可能含有未来强震孕震区的大量的本质性信息。     

大庆地震台网近震震级偏差分析与校正

选取大庆地震台网2002-2009年近震报告资料,以《黑龙江省地震目录》给出的震级为标准,进行震级偏差统计.按地区分段校正后,震级精度明显提高.     

基于损伤部位向量法的钢框架损伤识别研究

对损伤部位向量(DLV)法作了简单介绍,并用该方法对钢框架进行了损伤识别和损伤定位。该方法假定结构损伤前后为线性,对结构损伤前后柔度矩阵差进行奇异值分解,将奇异值为零所对应的向量,作为静荷载施加在无损结构的测点位置,则应力为零的单元为可能损伤的单元。对3种不同工况的钢框架进行了振动模态试验,用前3阶模态参数构造框架的柔度矩阵,按照DLV法对其进行了损伤识别,识别结果与已知损伤情况相一致。从测试自由度不完备、噪声和振型质量归一化系数这3个方面对识别效果进行了分析,结果表明:当损伤使结构动力特性有微小改变时,使用该方法不易定位损伤,应结合局部损伤识别方法进行判定;当损伤使结构动力特性有较大改变时,该方法能有效识别损伤的单元。DLV方法概念简单,理论明确,不受结构类型的限制,不需要结构的数学模型和模型缩聚或扩展技术,只需获得结构损伤前后的前几个低阶模态参数,即可识别结构一处或多处损伤,实际应用时可操作性强。     

新疆天山地区的波速比异常分析

Based on the seismic observation report data provided by the Xinjiang Digital Seismic Network from 2009 to 2014,we calculate the wave velocity ratio and its background value for medium and small earthquakes by using the multi-station method in Tianshan,Xinjiang.This paper analyzes the variation of the wave velocity ratio disturbance value to highlight the abnormal,and also back-traces 7 moderate earthquakes at the research area.The results show that:(1)the background value of the wave velocity ratio is almost 1.70,the wave velocity ratio obviously decreases in the middle-eastern part of Tianshan and the region near the Puchang fault;(2)the wave velocity ratio disturbance value is mostly low in the epicenter before four earthquakes of M≥5.0 from 2011 to 2013 in the study area;(3)before 7 moderate strong earthquakes,the earthquake events with low value of the wave velocity ratio account for over 60% of corresponding total events near the epicenters,and the low value of the wave velocity ratio is relatively obvious before moderate earthquakes.     

CT及MRI检查对胰腺囊性病变的诊断价值

目的:通过病理对照分析,探讨CT、MRI对胰腺囊性病变的诊断价值。方法:选取我院2013年1月至2016年12月收治的胰腺囊性病变患者60例,所有患者均分别采用CT、MRI检查,依据检查方法分为CT组和MRI组,对两组患者检查结果和病理结果比较、检查时间进行统计分析。结果:MRI诊断假性囊肿、真性囊肿、囊腺瘤、囊腺癌、胰腺导管内乳头状粘液瘤、癌变胰岛细胞瘤、实性假乳头状瘤与病理符合率分别为100.0%(13/13)、100.0%(16/16)、91.7%(11/12)、100.0%(7/7)、100.0%(6/6)、100.0%(3/3)、100.0%(3/3),CT诊断与病理符合率分别为84.6%(11/13)、93.8%(15/16)、100.0%(12/12)、85.7%(6/7)、83.3%(5/6)、66.7%(2/3)、100.0%(3/3);MRI诊断与病理符合率略高于CT诊断,差异不显著(P>0.05)。MRI组患者的检查时间显著长于CT组(P<0.01)。结论:CT与MRI对胰腺囊性病变的诊断率均较高,虽然MRI具有较长的检查时间,但其能够多方位成像,同时能够将更多的信息提供给临床对囊性成分的分析,特别是较为敏感的出血及囊性病变的大小、数量、范围、境界、成份、性质等,因此二者均具有一定特征性,建议临床结合应用。     

Evaporation modelling by heuristic regression approaches using only temperature data

Accurate estimation of pan evaporation (E_pan) is very important in water resources management, irrigation scheduling and water budget of lakes. This study investigates the accuracy of two heuristic regression approaches, multivariate adaptive regression splines (MARS) and M5 model tree (M5Tree) in estimating pan evaporation using only temperature data as input. Monthly minimum temperature, maximum temperature and E_pan data from three Turkish stations were used, with month number (periodicity information) added as input to see its effect on estimation accuracy. The models were compared with the calibrated Hargreaves-Samani (CHS), Stephens-Stewart (SS) and multiple linear regression methods. Three different train-test splitting strategies (50%–50%, 60%–40% and 75%–25%) were employed for better evaluation of the applied methods. The results show that the MARS method generally estimated monthly E_pan with higher accuracy compared to the M5Tree, CHS and SS methods. When extraterrestrial radiation, calculated from Julian date and latitude information, was used as input to the SS instead of solar radiation, satisfactory estimates were obtained. A positive effect on model accuracy was observed when involving periodicity information in inputs and increasing training data length.     

GOCE Gravity Field Processing Strategy

A processing strategy and the corresponding software architecture for the processing of GOCE (Gravity field and steady-state Ocean Circulation Explorer) observables is presented and described, with the major objective to compute a high-accuracy, high-resolution spherical harmonic model of the Earth's gravity field. The combination of two numerical solution strategies, i.e. the rigorous solution of the corresponding large normal equation systems applying parallel processing (on a PC cluster) as the core solver, and the fast semianalytic approach as a quick-look gravity field analysis (QL-GFA) tool, is proposed. Such a method fusion benefits from the advantages of the individual components: the rigorous inversion of the system providing also the full variance-covariance information, and the quickness enabling the consecutive production of intermediate gravity field solutions, for the purpose to analyse partial and incomplete data sets and to derive a diagnosis of the performance of the GOCE measurement system. The functionality and operability of the individual components are demonstrated in the framework of a closed loop simulation, which is based on a realistic mission scenario both in terms of the orbit configuration and the coloured measuring noise. Special concern is given to the accuracy of the recovered coefficients, the numerical behaviour, the required computing time, and the particular role of the individual modules within the processing chain. In the case of the core solver, it is demonstrated that the assembling and rigorous solution of large normal equation systems can be handled by using Beowulf clusters within a reasonable computing time. The application of the quick-look tool to partial data sets with short-term data gaps is demonstrated on the basis of several case studies. Additionally, the spectral analysis of the residuals of the adjustment is presented as a valuable tool for the verification of the noise characteristics of the GOCE gradiometer.     

On the Incorporation of the Effect of Crustal Structure into Empirical Strong Ground Motion Estimation

This article has two purposes. Firstly, a validation exercise of the modal summation technique for the computation of synthetic strong-motion records is performed for two regions of Europe (Umbria-Marche and south Iceland), using a variety of region specific crustal structure models, by comparing the predicted ground motion amplitudes with observed motions. It is found that the rate of decay of ground motions is well predicted by the theoretical decay curves but that the absolute size of the ground motions is underpredicted by the synthetic time-histories. This is thought to be due to the presence of low-velocity surface layers that amplify the ground motions but are not included in the crustal structure models used to compute the synthetic time-histories. Secondly, a new distance metric based on the computed theoretical decay curves is introduced which should have the ability to model the complex decay of strong ground motions. The ability of this new distance metric to reduce the associated scatter in empirically derived equations for the estimation of strong ground motions is tested. It is found that it does not lead to a reduction in the scatter but this is thought to be due to the use of crustal structure models that are not accurate or detailed enough for the regions studied. This revised version was published online in July 2006 with corrections to the Cover Date.