Evaluating the effect of land development on sediment transport using a probability density function

An important problem in sedimentation analysis is the development of a channel section that preserves, as best as possible, the current sedimentation regime even though the flood frequency tendencies have been altered due to land development within the catchment. In order to accomplish this task, a methodology is needed that estimates sediment transport capacity for various channel configurations. Such a procedure is described which allows the computation of the total sediment transport capacity for each of several T-year return frequency runoff hydrographs. This information is used to obtain an approximate probability distribution for the total sediment transport capacity, and the mean and standard deviation of this distribution are computed.Comparing the results for the catchment in its present state with a future developed state, using a selection of new channel parameters, indicates how to improve the channel to control changes in sedimentation due to development. The analysis procedure provides a basis for estimating a new channel configuration such that the new flow conditions retain, as best as possible, the existing condition sedimentation effects, and hence retain the natural sediment supply and transport trends even though runoff flow rates have changed due to land development within the catchment.The results of Wilson Creek are typical of the several sites examined, see Table 3 below. The T=2, T=5, T=25, and T=100 year values for total sediment transport capacity, in kilotons, are 6.9, 39.4, 61.3, and 96.7 with a mean of 17.1 and standard deviation of 19.3. After development with no change in the channel the respective values increase to: 17.9, 84.6, 128.1, and 258.0 with a mean of 39.1 and standard deviation of 44.3. A new channel can be constructed which will reduce these sediment transport capacity values, after development, to 5.2, 41.0, 62.0, and 124.8 with a mean of 17.4 and standard deviation of 22.0.     

新地震区划图的几个基本问题

我国新的地震区划图(1990年版)是采用地震危险性慨率分析方法编制的。该图给出的是场点地震烈度值,该值在50年内被突破的概率为0.1。人们普遍关注该图与我国曾经编制的地震区划图(1957年版,1977年版)的区别,该图超越概率概念的内含和外延以及超越概率水平为什么采用50年超越概率0.1。本文围绕这些问题进行了讨论。分析结果表明,前两张地震区划图编图的基本着眼点都是地震预测,而新的地震区划着眼于场点的地震动预测。新的地震区划图是按场点地震危险性分析方法给出的,它所表示的地震危险性只能针对具体的场点,不能完全反映区域的地震危险性特征。而弄清场点地震危险性和区域地震危险性的差异是正确进行区域防灾对策的基础。作者希望这些讨论能对正确使用新的地震区划图有所裨益。     

第23周太阳活动峰年北京天文台太阳射电观测

An introduction to the instrumonts to be used for solar radio observations in Beijing in Solar Cycle 23 is made in this paper. They are 10cm solar radio telescope which has been used for a long time,and 1.0～2.0GHz,2.6～3.8GHz and 5.2～7.6GHz spectrometers. The former two spectrometers has passed a test,showing high guality,and the 5.2～7.6GHz one will be in operation from 1999. It is believed that highl gualified data should be obtained in the next solar cycle.     

吉林大阳岔上寒武统凤山组—下奥陶统冶里组层序地层和化学地层研究

根据岩相序列、沉积、成岩作用特征,将吉林省大阳岔地区晚寒武世凤山期至早奥陶世冶里期层序地层分为６个三级层序,代表了６个三级海侵—海退旋回;与之相应,该时期碳酸盐岩δ１３Ｃ值变化曲线出现５个低谷,其中凤山期至冶里早期泥晶灰岩δ１３Ｃ值出现４次幅度较大的上升、下降飘移,而冶里中、晚期泥晶灰岩的δ１３Ｃ值是在负值范围内波动。δ１３Ｃ值低谷与层序界面和海退事件基本拟合。牙形类Ｃ．ｉｎｔｅｒｍｅｄｉｕｓ带下部海绿石质凝缩段是寒武纪末—奥陶纪初期最大海泛事件沉积,也是寒武纪与奥陶纪之交等时性的年代地层对比标志。     

Short Term and Short Range Seismicity Patterns in Different Seismic Areas of the World

The aim of this work is to quantitatively set up a simple hypothesis for occurrence of earthquakes conditioned by prior events, on the basis of a previously existing model and the use of recent instrumental observations. A simple procedure is presented in order to determine the conditional probability of pairs of events (foreshock-mainshock, mainshock-aftershock) with short time and space separation. The first event of a pair should not be an aftershock, i.e., it must not be related to a stronger previous event. The Italian earthquake catalog of the Istituto Nazionale di Geofisica (ING) (1975–1995, M 3.4), the earthquake catalog of the Japan Meteorological Agency (JMA) (1983–1994, M 3.0) and that of the National Observatory of Athens (NOA) (1982–1994, M 3.8) were analyzed. The number of observed pairs depends on several parameters: the size of the space-time quiescence volume defining nonaftershocks, the inter event time, the minimum magnitude of the two events, and the spatial dimension of the alarm volume after the first event. The Akaike information criterion has been adopted to assess the optimum set of space-time parameters used in the definition of the pairs, assuming that the occurrence rate of subsequent events may be modeled by two Poisson processes with different rates: the higher rate refers to the space-time volume defined by the alarms and the lower one simulates earthquakes that occur in the nonalarm space-time volume. On the basis of the tests carried out on the seismic catalog of Italy, the occurrence rate of M 3.8 earthquakes followed by a M 3.8 mainshock within 10 km and 10 days (validity) is 0.459. We have observed, for all three catalogs, that the occurrence rate density for the second event of a couple (mainshock or aftershock) of magnitude M₂ subsequent to a nonaftershock of magnitude M₁ in the time range T can be modeled by the following relationship: (T, M₂) = 10^{a + b(M1 - M2)} with b varying from 0.74 (Japan) to 1.09 (Greece). The decrease of the occurrence rate in time for a mainshock after a foreshock or for large aftershocks after a mainshock, for all three databases, obeys the Omori's law with p changing from 0.94 (Italy) to 2.0 (Greece).     

副高异常偏西背景下西北一次漏报的雷暴大风中尺度分析*

摘要：利用西安多普勒天气雷达、L波段风廓线雷达和加密自动站探测资料，结合天气实况，对2018年7月26日发生在陕西一次副热带高压影响下的的强风暴过程进行了中尺度分析。结果表明：（1）本次强风暴伴随的阵风锋共维持了4h，其中有3h出现7级以上的大风，且最大风力10级。（2）在副热带高压影响下，陕西处于高温、高湿气团中，大气层结极不稳定。（3）此次强风暴在高的对流有效位能环境下（CAPE）下，抬升触发的关键因子是关中地区中尺度辐合线，当初生的对流云团下山后，中尺度辐合线触发的对流风暴形成小范围冷池出流与环境风场形成新的辐合线，加强对流风暴发展。（4）当阵风锋移动过程中遇到前方的对流云团时，将低层暖湿空气抬升，并随着上升气流输送到主体对流风暴中，迅速补充了主体风暴的能量，使得主体风暴再次强烈发展，延长了阵风锋的生命史。     

Topological properties of disjoint channel networks within enclosed regions

Enclosure of some portion of one or more natural stream-drainage basins by superposition of a rectangle on a map of drainage network results in fragmentation of the natural basins into a set of disjoint channel networks. Each of these may have some channel links and forks of the natural network plus truncated links intersected by the enclosure boundary. The topological properties of the network elements in the enclosure are used to set up a model of random network patterns, in which the number of disjoint channel networks is expressed as a function of the number of links and forks in the enclosures. This function is shown to be a multiplicative constant times the square root of the number of links or forks. Empirical data on square and rectangular enclosures of several sizes from the Inez (Kentucky)Quadrangle map showed that the predicted multiplicative constants do not agree with observation, but that the square-root relation seems to hold at least to a first approximation. The models thus can be used as a norm against which deviations of real-world enclosures from network pattern randomness can be studied.     

On solution of the inverse problem for confined aquifer flow via maximum likelihood

Joint estimation of transmissivity (T) and storativity (S) in a confined aquifer is done via maximum likelihood (ML). The differential equation of groundwater flow is discretized by the finite-element method, leading to equation _t+x _t=u _t. Elements of matrices and , as well as estimated covariance matrix of noise termu _t, are functions of T and S. By minimizing the negative loglikelihood function corresponding to discretized groundwater flow equation with respect to T and S, ML estimators are obtained. The ML approach is found to yield accurate estimates of T and S (within 9 and 10% of their actual values, respectively) and showed quadratic convergence in Newton's search technique. Prediction of aquifer response, using ML estimators, results in estimated piezometric heads accurate to ±0.5 m from their actual, exact values. Statistical properties of ML estimators are derived and some basic results for statistical inference are given.     

A computer method to calculate the convolution of statistical distributions

It is becoming increasingly important to determine probability distributions of combinations of random variables. Convolution is a technique by which the distribution of a sum of random variables can be determined. This paper presents some simplifications in order to reduce the numerical integrations and computer time. In addition, the method may be used with empirical nonanalytic distributions. While Monte Carlo methods are also appropriate for calculating the distribution, convolution can give at least as much accuracy as Monte Carlo methods with a reduction in computation. Two applications are presented: one approximates the distribution of percent sand in an area, and the other indicates a method of determining sample size when using the distribution of means to approximate normality.