Application of bivariate extreme value distribution to flood frequency analysis: a case study of Northwestern Mexico

In Mexico, poverty has forced people to live almost on the water of rivers. This situation along with the occurrence of floods is a serious problem for the local governments. In order to protect their lives and goods, it is very important to account with a mathematical tool that may reduce the uncertainties in computing the design events for different return periods. In this paper, the Logistic model for bivariate extreme value distribution with Weibull-2 and Mixed Weibull marginals is proposed for the case of flood frequency analysis. A procedure to estimate their parameters based on the maximum likelihood method is developed. A region in Northwestern Mexico with 16 gauging stations has been selected to apply the model and regional at-site quantiles were estimated. A significant improvement occurs, measured through the use of a goodness-of-fit test, when parameters are estimated using the bivariate distribution instead of its univariate counterpart. Results suggest that it is very important to consider the Mixed Weibull distribution and its bivariate option when analyzing floods generated by a␣mixture of two populations.     

Two-dimensional modelling of some CFC replacement compounds

The Cambridge 2-D Eulerian model has been used to study the potential atmospheric distributions and lifetimes of a number of CFC replacement compounds and their degradation products. The study has focused on HFC 134a and HCFCs 123, 141b and 142b and the major products formed by their atmospheric degradation. The loss of these compounds and their products by hydroxyl radical attack, photolysis and in-cloud hydrolysis have been investigated. The study has shown that HCFCs 141b and 142b have sufficiently long lifetimes to enter the stratosphere in significant quantities, where degradation leads to an increase in the total stratospheric chlorine concentration. The study has also highlighted areas where further experimental work would be valuable, in particular characterisation of the product channels for the degradation reactions and determination of the removal rates of the products in the aqueous phase.     

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.     

Estimation of Scalar Source/Sink Distributions in Plant Canopies Using Lagrangian Dispersion Analysis: Corrections for Atmospheric Stability and Comparison with a Multilayer Canopy Model

Source/sink distributions of heat, water vapour andCO₂ within a rice canopy were inferred using aninverse Lagrangian dispersion analysis and measuredmean profiles of temperature, specific humidity andCO₂ mixing ratio. Monin–Obukhov similarity theorywas used to account for the effects of atmosphericstability on _w(z), the standard deviation ofvertical velocity and _L(z), the Lagrangian timescale of the turbulence. Classical surface layer scaling was applied in the inertial sublayer (z > z_ruf)using the similarity parameter = (z - d)/L, where z is height above ground, d is the zero plane displacementheight for momentum, L is the Obukhov length,and z_ruf 2.3h_c, where h_c iscanopy height. A single length scale h_c, was usedfor the stability parameter 3 = h_c/L in the height range 0.25 < z/h_c < 2.5. This choice is justified by mixing layer theory, which shows that within the roughness sublayer there is one dominant turbulence length scaledetermined by the degree of inflection in the windprofile at the canopy top. In the absence of theoretical or experimental evidence for guidance,standard Monin–Obukhov similarity functions, with = h_c/L, were used to calculate the stabilitydependence of _w(z) and _L(z) in the roughness sublayer. For z/h_c < 0.25 the turbulence length and time scales are influenced by the presence of the lowersurface, and stability effects are minimal. With theseassumptions there was excellent agreement between eddycovariance flux measurements and deductions from theinverse Lagrangian analysis. Stability correctionswere particularly necessary for night time fluxes whenthe atmosphere was stably stratified.The inverse Lagrangian analysis provides a useful toolfor testing and refining multilayer canopy models usedto predict radiation absorption, energy partitioningand CO₂ exchanges within the canopy and at thesoil surface. Comparison of model predictions withsource strengths deduced from the inverse analysisgave good results. Observed discrepancies may be dueto incorrect specification of the turbulent timescales and vertical velocity fluctuations close to theground. Further investigation of turbulencecharacteristics within plant canopies is required toresolve these issues.     

我国北方地区雷电活动的时空特征

利用地闪定位仪（ＤＦ）探测了我国大兴安岭（东北地区）和北京地区的雷电活动。文中详细分析了地闪的空间分布,并讨论了不同地区各雷电参量的时间变化规律。两地区雷电活动的季节变化与日变化均有较大差异。这些差异的产生可能与本地区的天气背景及地理特征有密切的联系。要较深入地了解上述差异产生的可能原因,将雷电活动与天气因子的给合讨论是必需的。     

海南省澄迈县雷暴气候特征及其灾害防御

通过统计分析海南省澄迈县1959-2004年逐日雷暴资料,找出雷暴发生时空分布特征,统计1995-2004年10a雷暴发生的影响天气系统,计算出各天气系统影响下雷暴发生的概率。分析结果表明：澄迈县雷暴日46a平均雷暴频次437次,年际变化呈波动减少趋势,5-8月是雷暴发生的高发期,占全年雷暴的53%,16-17时是一天中发生雷暴的最高期,西南方向发生的雷暴略多于其他方向,澄迈县受西南低压槽、南海低压槽和副热带高压等天气系统影响时发生雷暴的概率较大。     

Fractals and the accuracy of geographical measures

The problems of estimating line length, area, and point characteristics in the earth sciences have generated substantial but independent literatures. All three problems are of increasing concern given the current interest in digital capture, processing, and the storage of geographically referenced data. In the case of qualitative maps, all three are shown to be related to Mandelbrot's fractional dimension D (Mandelbrot, 1977) which allows the dependence of each on sampling density to be predicted. The general results are confirmed by simulation on surfaces of constant D. They also imply that certain improvements can be made in a number of previously proposed methods.     

BIB-SEM characterization of pore space morphology and distribution in postmature to overmature samples from the Haynesville and Bossier Shales

Four Haynesville Shale and four Bossier Shale samples were investigated using a combination of Scanning Electron Microscopy (SEM) and Broad Ion Beam (BIB) polishing. This approach enables the microstructure and porosity to be studied down to the mesopore size (<50 nm) in representative areas at the scale of the BIB cross-sections. The samples vary in mineralogy, grain size and TOC and the maturity ranges from 2.42 to 2.58 VR_r in the Haynesville Shale to 1.79–2.26 VR_r in the Bossier Shale. This variety within the samples enabled us to study controls on the porosity distribution in these shales. Visible pores exist as intraparticle pores mainly in carbonate grains and pyrite framboids and as interparticle pores, mainly in the clay-rich matrix. Pores in organic matter show a characteristic porosity with respect to the type of organic matter, which mainly consists of mixtures of amorphous organic matter and minerals, organic laminae and discrete macerals. A clear positive trend of organic-matter porosity with maturity was found. Pore sizes are power law distributed in the range of 4.4 μm to at least 36 nm in equivalent diameter. The differences in power law exponents suggest that a more grain supported, coarse-grained matrix may prevent pores from mechanical compaction. Porosities measured in the BIB cross-sections were significantly lower in comparison to porosities obtained by Mercury Intrusion Porosimetry (MIP). This difference is mainly attributed to the different resolution achieved with BIB-SEM and MIP and type of pore network. Extrapolation of pore size distributions (PSDs) enables the BIB-SEM porosity to be estimated down to the resolution of the MIP and thus to upscale microstructural observation at the confined space of the BIB-SEM method to bulk porosity measurement. These inferred porosities are in good agreement with the MIP determined porosities, which underpins the assumption that pores segmented in BIB-SEM mosaics are representative of the MIP methodology.