Uncertainty,entropy, scaling and hydrological stochastics. 1. Marginal distributional properties of hydrological processes and state scaling / Incertitude,entropie, effet d'échelle et propriétés stochastiques hydrologiques. 1. Propriétés distributionnelles marginales des processus hydrologiques et échelle d'état

Abstract

The well-established physical and mathematical principle of maximum entropy (ME), is used to explain the distributional and autocorrelation properties of hydrological processes, including the scaling behaviour both in state and in time. In this context, maximum entropy is interpreted as maximum uncertainty. The conditions used for the maximization of entropy are as simple as possible, i.e. that hydrological processes are non-negative with specified coefficients of variation (CV) and lag one autocorrelation. In this first part of the study, the marginal distributional properties of hydrological variables and the state scaling behaviour are investigated. Application of the ME principle under these very simple conditions results in the truncated normal distribution for small values of CV and in a nonexponential type (Pareto) distribution for high values of CV. In addition, the normal and the exponential distributions appear as limiting cases of these two distributions. Testing of these theoretical results with numerous hydrological data sets on several scales validates the applicability of the ME principle, thus emphasizing the dominance of uncertainty in hydrological processes. Both theoretical and empirical results show that the state scaling is only an approximation for the high return periods, which is merely valid when processes have high variation on small time scales. In other cases the normal distributional behaviour, which does not have state scaling properties, is a more appropriate approximation. Interestingly however, as discussed in the second part of the study, the normal distribution combined with positive autocorrelation of a process, results in time scaling behaviour due to the ME principle.     

Preliminary results of seismic hazard assessment in the Sannio-Matese area of Southern Italy

The seismic hazard in the Sannio-Matese area has been worked out by a modification of the McGuire (1976) computing programme, taking into account the influence of nine potential seismic source zones.The method uses truncated-quadratic intensity-frequency distribution and azimuth-dependent intensity attenuation derived from isoseismal maps for each of the seismogenetic sources. A new modification has been introduced to take into account different decay of the intensity in the near (to VIII degree) and far (from VIII degree) field.Different assumptions about maximum possible intensities and truncation of intensity-frequency laws are used to evaluate the effects of the uncertainties on the computed hazard at high intensities. Intensities associated with different level of annual probability are computed for five test sites in the considered area. Maps displaying the expected intensity for a mean return period of 500 years (pa 0.002) are presented and compared with observed intensities.Presented at the XXIst General Assembly of the European Seismological Commission, Symposium on Methods of Seismic Hazard Assessment in Europe, Sofia, 23–27 August 1988.     

Statistical estimation of a shock center: Slate Islands astrobleme

Statistical models are proposed for the distribution of directions in three dimensions that are thought to point toward a single source. These models are based on the Fisher distribution. The method of maximum likelihood is used to obtain an estimate of position of the source and to provide corre-sponding confidence regions. When applied to shatter cone data from the Slate Islands, northern Lake Superior, the method yields estimates comparable with those obtained by Stesky and Halls (1982), as well as statistically valid confidence regions.     

重磁解释中健全的极大似然法及三维反演

本文介绍了JOAO.B.C.SILVA提出的重磁解释中健全的(Robust)极大似然反演方法的基本原理,讨论了重力二度体非线性几何参数的反演情况,并将该方法扩展到三维重、磁非线性几何参数反演中。重点讨论了场值中有较大随机干扰、有不同类型的地质噪音源影响下的反演结果。理论反演结果表明,该方法有比最小二乘反演方法更接近实际的解,从而提高了解的可靠性,对于重磁局部异常或孤立异常有较好的反演效果。     

山西地震带历史地震序列的统计分析

在工程地震工作中需要知道某地震带或某地区未来100年内的地震活动水平,以便适当选择历史地震的统计时段,从而得到和地震活动水平相适应的不同潜在震源区的各级地震的年平均发生率。 本文利用极值理论、最大熵原理、马尔可夫模型的方法、莫尔纳方法和伯努里模型方法计算了山西地震带未来100年内中强震的年平均发生率、平均重现期大于等于某级地震的个数和发震概率,计算结果表明,山西地震带未来百年内将发生≥6.0级地震3次,发震概率为0.94左右;将发生≥6.5级地震1—2次,发震概率为0.84左右;有可能发生≥7.0级地震,发震概率为0.53左右;发生≥7.5级地震的可能性不大,发震概率为0.32左右;发生≥8.0级地震的可能性很小,发震概率为0.15左右。     

澜沧—耿马7.6级地震前中小地震活动信息熵和分维的研究

从信息论原理和分形几何学的观点出发,采用多种方法,分析研究1988年澜沧—耿马7.6级地震前地震活动的信息、信息熵和分维的变化特征表明,地震前1—2年内,中、小地震活动的信息熵、缺信量和容量维、信息维出现明显的低值变化过程,反映了地震前兆的减熵、降维特点,从一个侧面反映了大震前,一定区域内中、小地震活动的时空变化从无序向有序的演化过程。     

Seismic recognition of local inhomogeneities

A simple method to contour local inhomogeneities using seismic data is proposed. It formalizes an approximate inversion method which is based on the interpretation of local inhomogeneities as making the differences between an actual seismic data set and a previous reference model. It uses the optimal statistical criteria of parameter estimation and recognition and the ray representation of the waves spreading. Any combination of direct, reflected and/or other types of waves may be used as the database. Inhomogeneities, having a size two times above the wavelength of the seismic waves, can be resolved. Laboratory experiments, using ultrasonic waves and analysis of data from field experiments, confirmed the theoretical results. The method can be used to search for ore bodies, kimberlite cubes, oiltraps, etc.     

A SEARCH TECHNIQUE FOR CHOOSING OPTIMUM PARAMETERS IN NUMERICAL WEATHER PREDICTION

The steepest descent(or ascent)search is employed for finding optimum diffusion coefficients in T42L9Gmodel,with a view to improving the model's computational stability or prediction accuracy.The method of thesteepest descent search is first described,in which the golden section search is chosen as the fundamental one-dimensional search used in the multi-dimentional steepest descent search,and then the optimization of the dif-fusion coefficients is described.     

Entropy estimation of hydrological extremes

The method of Relative Entropy with Fractile constraints (REF method) is explained and applied to model extreme compound hydrological phenomena, such as extreme sea levels under storm conditions. Also presented is a simple method of Tail Entropy Approximation (TEA), which amounts to a correction of traditional statistical estimates for extreme observations.Distribution assumptions are necessary but downplayed in the REF method, relegating the prior distribution to the role of an extrapolation function. The estimates are objective in an information-theoretical sense. They also satisfy a strict requirement of self-consistency that is generally not satisfied by standard statistical methods: invariance under monotonic transformations of the random variable.Historical records of storm surge levels in the Netherlands and annual maximum tidal heights for Sheerness, UK, are used as examples. Comparison is made with distributions obtained using other methods.It is concluded that the tail entropy approximation provides simple, objective estimates of extremes in the tail beyond the range of observations.     

Research of pattern dynamics parameters of crustal deformation field in seismogenic process

ＲｅｓｅａｒｃｈｏｆｐａｔｔｅｒｎｄｙｎａｍｉｃｓｐａｒａｍｅｔｅｒｓｏｆｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄｉｎｓｅｉｓｍｏｇｅｎｉｃｐｒｏｃｅｓｓＳｈｕｏ－ＹｕＺＨＯＵ（周硕愚），ＹｕｎＷＵ（吴云），Ｒｕｏ－ＢｏＷＡＮＧ（王若柏）ａｎｄ...