Testing a statistical curve‐fitting procedure for quantifying sediment populations within multi‐modal particle‐size distributions

This paper describes a series of tests designed to evaluate the capacity of a personal computer (PC) based statistical curve‐fitting program called MIX to quantify composite populations within multi‐modal particle‐size distributions. Three natural soil samples were analysed by a Coulter Multisizer, and their particle‐size distributions analysed using MIX software to identify the modes, standard deviations and proportions of their composite populations. The particle‐size distributions of the three natural soil samples were then numerically combined in equal proportions using a spreadsheet program to create synthetic particle‐size distributions of known populations. MIX was then tested on the synthetic particle‐size distributions to see if the modes and proportions it identified were similar to those modes and proportions known to characterize the synthetic particle‐size distributions. The main outcome is that MIX can very accurately describe the modal particle size and proportions of the major composite populations within a particle‐size distribution. However MIX has difficulty in identifying small populations (those contributing <10 per cent of a total particle‐size distribution), particularly when they are located in the central sections of particle‐size distributions, overlain by larger populations, or when positioned in the fine tails of distributions. Despite these minor shortcomings, MIX is a valuable tool for the examination and interpretation of particle‐size data. Copyright © 2005 John Wiley & Sons, Ltd.     

Geostatistical simulation of random fields with bivariate isofactorial distributions by adding mosaic models

This work deals with the geostatistical simulation of a family of stationary random field models with bivariate isofactorial distributions. Such models are defined as the sum of independent random fields with mosaic-type bivariate distributions and infinitely divisible univariate distributions. For practical applications, dead leaf tessellations are used since they provide a wide range of models and allow conditioning the realizations to a set of data via an iterative procedure (simulated annealing). The model parameters can be determined by comparing the data variogram and madogram, and enable to control the spatial connectivity of the extreme values in the realizations. An illustration to a forest dataset is presented, for which a negative binomial model is used to characterize the distribution of coniferous trees over a wooded area.     

Scattering of plane SV waves by cylindrical canyons in saturated porous medium

On the basis of Biot dynamic theory, an analytic solution of two-dimensional scattering and diffraction of plane SV waves by circular cylindrical canyons in a half space of saturated porous media is presented in this paper for the first time. The solution is obtained by employing the Fourier–Bessel series expansion technique. Parametric studies had been carried out, which includes: the angle of incidence, the frequency of the incident SV wave, the porosity of saturated porous medium and the stiffness and Poisson's ratio of the solid-skeleton. All the outcomes are useful for the seismic analysis of the surface topography conditions.     

On size and boundary effects in scaled model blasts spatial problems

This paper addresses size and boundary effects on wave propagation, fracture pattern development and fragmentation in small scale laboratory-size specimens for model blasting. Small block type specimens are centre-line loaded by linear explosive charges and supersonically detonated. Using elastic wave propagation theory and fracture mechanics it is shown that the type of boundary conditions which prevail at the outer boundary of the cylinder control the extension of bore-hole cracking and fragmentation within the body of the cylinder. In the case of a composite block where a cylindrical core of different material is embedded, the level of fracturing and fragmentation is controlled by the separation of the interface which in turn depends on the relative dimensions of the core and the block. The most important parameter is the ratio between the length of the pulse (space-wise or time-wise) and the characteristic dimensions of the models, i.e. in this case the dimensions of the core and the mantel. Stress wave superposition effects occur in the corner sections of the mantel. Theoretical results are in good agreement with recent experimental findings.     

兰州市2001年沙尘气溶胶质量浓度的特征分析

 分析了2001年沙尘暴期间兰州与靠近沙尘源区的武威的沙尘浓度和粒径分布特征,并运用对数正态分布规律拟合了沙尘粒径的分布。通过对比武威、皋兰和兰州沙尘暴期间沙尘浓度的变化以及武威与兰州的沙尘粒径分布特征,揭示了河西走廊沙漠对兰州市沙尘暴的影响。     

Stress‐dependent elastic wave velocity of microfractured sandstone

A model for the stress‐dependent elastic wave velocity response of fractured rock mass is proposed based on experimental evidence of stress‐dependent fracture normal and shear stiffness. Previously proposed models and previous experimental studies on stress‐dependent fracture stiffness have been reviewed to provide a basis for the new model. Most of the existing stress‐dependent elastic wave velocity models are empirical, with model parameters that do not have clear physical meanings. To propose the new model, the rock mass is assumed to have randomly oriented microscopic fractures. In addition, the characteristic length of microfractures is assumed to be sufficiently short compared to the rock mass dimensions. The macroscopic stress‐dependent elastic wave velocity response is assumed to be attributed to the stress dependency of fracture stiffness. The stress‐dependent fracture normal stiffness is defined as a generalized power law function of effective normal stress, which is a modification of the Goodman's model. On the other hand, the stress dependency of fracture shear stiffness is modeled as a linear function of normal stress based on experimental data. Ultrasonic wave velocity responses of a dry core sample of Berea sandstone were tested at effective stresses ranging from 2 to 55 MPa. Visual observation of thin sections obtained from the Berea sandstone confirms that the assumptions made for microstructure of rock mass model are appropriate. It is shown that the model can describe the stress‐dependent ultrasonic wave velocity responses of dry Berea sandstone with a set of reasonable material parameter values. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

Modeling nonlinear wave–wave interactions with the elliptic mild slope equation

An approach is developed to simulate wave–wave interactions using nonlinear elliptic mild-slope equation in domains where wave reflection, refraction, diffraction and breaking effects must also be considered. This involves the construction of an efficient solution procedure including effective boundary treatment, modification of the nonlinear equation to resolve convergence issues, and validation of the overall approach. For solving the second-order boundary-value problem, the Alternating Direction Implicit (ADI) scheme is employed, and the use of approximate boundary conditions is supplemented, for improved accuracy, with internal wave generation method and dissipative sponge layers. The performance of the nonlinear model is investigated for a range of practical wave conditions involving reflection, diffraction and shoaling in the presence of nonlinear wave–wave interactions. In addition, the transformation of a wave spectrum due to nonlinear shoaling and breaking, and nonlinear resonance inside a rectangular harbor are simulated. Numerical calculations are compared with the results from other relevant nonlinear models and experimental data available in literature. Results show that the approach developed here performs reasonably well, and has thus improved the applicability of this class of wave transformation models.     

Frequency analysis of nonstationary annual maximum flood series using the time‐varying two‐component mixture distributions

The most popular practice for analysing nonstationarity of flood series is to use a fixed single‐type probability distribution incorporated with the time‐varying moments. However, the type of probability distribution could be both complex because of distinct flood populations and time‐varying under changing environments. To allow the investigation of this complex nature, the time‐varying two‐component mixture distributions (TTMD) method is proposed in this study by considering the time variations of not only the moments of its component distributions but also the weighting coefficients. Having identified the existence of mixed flood populations based on circular statistics, the proposed TTMD was applied to model the annual maximum flood series of two stations in the Weihe River basin, with the model parameters calibrated by the meta‐heuristic maximum likelihood method. The performance of TTMD was evaluated by different diagnostic plots and indexes and compared with stationary single‐type distributions, stationary mixture distributions and time‐varying single‐type distributions. The results highlighted the advantages of TTMD with physically‐based covariates for both stations. Besides, the optimal TTMD models were considered to be capable of settling the issue of nonstationarity and capturing the mixed flood populations satisfactorily. Copyright © 2016 John Wiley & Sons, Ltd.     

极化干涉SAR森林树高提取算法的对比分析

在对极化干涉SAR森林树高反演的DEM差值算法、相干相位-幅度综合反演算法进行分析的基础上,对基于极化干涉相干优化方法的改进算法进行了探讨。利用黑龙江大兴安岭地区的一对ALOS全极化干涉SAR数据进行实验,并对比分析各算法的反演结果。结果表明,在使用改进的算法进行森林树高反演时可以获取精度较高的反演结果,并且在一定程度上提高了森林树高反演的稳定性,为森林树高反演工作的业务化运行提供一定的依据。