一种快速角点探测算子研究

基于SUSAN和角点几何结构分析的思想提出一种改进的快速角点探测算子.首先进行角点粗探测,然后再对粗探测的结果进行角点几何结构分析,即图像中的像元如果为角点,它必须满足角点的圈层结构.通常认为角点处各个方向上灰度差变化都较大,所以角点响应函数采用反映角点突出性的响应函数,使提取的角点含有更高的信息量.实验证明,本算法提取的角点定位准确,快速有效,可以适用于实时图像处理.     

On the scavenging of SO2 by cloud and raindrops: II. An experimental study of SO2 absorption and desorption for water drops in air

For the purpose of testing our previously described theory of SO₂ scavenging a laboratory investigation was carried out in the UCLA 33 m long rainshaft. Drops with radii between 250 and 2500 m were allowed to come to terminal velocity, after which they passed through a chamber of variable length filled with various SO₂ concentrations in air. After falling through a gas separating chamber consisting of a fluorocarbon gas the drops were collected and analyzed for their total S content in order to determine the rate of SO ₂ absorption.The SO₂ concentration in air studied ranged between 1 and 60% (v). Such relatively large concentrations were necessary due to the short times the drops were exposed to SO₂ in the present setup. The present experimental results were therefore not used to simulate atmospheric conditions but rather to test our previously derived theory which is applicable to any laboratory or atmospheric condition. Comparison of our studies with the results from our theory applied to our laboratory conditions led to predicted values for the S concentration in the drops which agreed well with those observed if the drops had radii smaller than 500 m. In order to obtain agreement between predicted and observed S concentrations in larger drops, an empirically derived eddy diffusivity for SO₂ in water had to be included in the theory to take into account the effect of turbulent mixing inside such large drops.In a subsequent set of experiments, drops initially saturated with S (IV) were allowed to fall through S-free air to determine the rate of SO ₂ desorption. The results of these studies also agreed well with the results of our theoretical model, thus justifying the reversibility assumption made in our theoretical models.In a final set of experiments, the effects of oxidation on SO₂ absorption was studied by means of drops containing various amounts of H₂O₂. For comparable exposure times to SO₂, the S concentration in drops with H₂O₂ was found to be up to 10 times higher than the concentration in drops in which no oxidation occurred.     

The characteristics of shear wave splitting in the source region of the April 20, 2013 Lushan earthquake

Using seismic data of the aftershocks sequence of the April 20, 2013 Lushan earthquake recorded by seismic temporary and permanent stations in the source region, with the visual inspection of particle motion diagrams, this paper preliminarily contains the polarization directions of fast shear wave and the time-delays of split shear waves at every station, and analyzes the crustal anisotropic characteristics in the source region. In the study area, the polarization direc- tions at stations BAX, TQU, L 132, L 133, L 134, and L 135 are northeast, which is consistent with the strike of Dachuan- Shuangshi fault. There are two polarization directions at MDS and L131, which are northeast and southeast. The scatter of polarization directions suggests the complex stress field around these two stations where two faults intersect. For the normalized time-delays at every station, the range is 1.02-8.64 ms/km. The largest time-delay is from L134 which is closest to the mainshock, and the smallest one is from L133. The variations in time-delays show the decreasing at stations BAX, L134, and L135 because of the stress-relaxation after earthquake.     

缓变地形下β效应的Rossby代数孤立波

正压流体中,从准地转位涡方程出发采用摄动方法和时空伸长变换推导了在缓变地形下β效应的Rossby代数孤立波方程,得到Rossby波振幅满足带有缓变地形非齐次Benjamin-Davis-Ono(BDO)方程的结论.通过分析孤立Rossby波振幅的演变,指出了β效应、地形效应是诱导Rossby孤立波产生的重要因素,说明了在缓变地形强迫效应和非线性作用相平衡的假定下,Rossby孤立波振幅的演变满足非齐次BDO方程,给出在切变基本气流下缓变地形和正压流体中Rossby波的相互作用.     

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second-order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite-element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.     

地震波动方程的局部间断有限元方法数值模拟

近年来,随着地震波数值模拟对计算精度和效率的要求越来越高,间断有限元方法开始受到越来越多的关注.本文中,针对具有吸收边界条件的二维地震声波波动方程,作者提出了一种基于局部间断有限元方法的数值模拟算法.该算法在空间上使用局部间断有限元方法进行离散,在时间上采用了显式蛙跳格式.在这种时空离散的组合方式下,每个时间步上,此算法在空间剖分的每个单元上的求解计算是相互独立的,因而具有极高的并行性.通过数值算例,我们将该算法与连续有限元方法进行了比较.结果表明,本算法不仅具有对起伏构造的良好适应性,而且在计算效率和计算精度等方面,都具有优越性.     

Spatio-temporal variations of reference crop evapotranspiration and pan evaporation in the West Songnen Plain of China

Abstract

The spatio-temporal variations of reference crop evapotranspiration (ET_ref) reflect the combined effects of meteorological variables, primarily wind speed, relative humidity, net radiation and air temperature. This study investigated the spatial distribution and temporal trends of ET_ref (calculated by the FAO-56 Penman-Monteith equation), pan evaporation (E_pan) and pan coefficient (K_p) in a 140?×?10³ km² semi-humid to semi-arid area in China. The results show that: (i) although the spatial distributions of ET_ref and E_pan are roughly similar and their spatial correlation is high over the growing season, K_p varied considerably in space due to high humidity in the east of the region and low humidity in the southwest; (ii) the monthly variations of ET_ref and E_pan are similar to that of net radiation and opposite to that of relative humidity, while the monthly variation of K_p is similar to that of relative humidity and opposite to that of wind speed, and the long-term trend is slightly increasing for ET_ref and E_pan, while significantly (10% significance level) increasing for K_p; and (iii) generally, the time series of ET_ref and E_pan from 1951 to 2001 could be divided into three phases due to variations of meteorological variables.

Citation Liang, L.-Q., Li, L.-J. & Liu, Q. (2011) Spatio-temporal variations of reference crop evapotranspiration and pan evaporation in the West Songnen Plain of China. Hydrol. Sci. J. 56(7), 1300–1313.     

Fission of a weakly nonlinear interfacial solitary wave at a step

The transformation of a weakly nonlinear interfacial solitary wave in an ideal two-layer flow over a step is studied. In the vicinity of the step the wave transformation is described in the framework of the linear theory of long interfacial waves, and the coefficients of wave reflection and transmission are calculated. A strong transformation arises for propagation into shallower water, but a weak transformation for propagation into deeper water. Far from the step, the wave dynamics is described by the Korteweg-de Vries equation which is fully integrable. In the vicinity of the step, the reflected and transmitted waves have soliton-like shapes, but their parameters do not satisfy the steady-state soliton solutions. Using the inverse scattering technique it is shown that the reflected wave evolves into a single soliton and dispersing radiation if the wave propagates from deep to shallow water, and only dispersing radiation if the wave propagates from shallow to deep water. The dynamics of the transmitted wave is more complicated. In particular, if the coefficient of the nonlinear quadratic term in the Korteweg-de Vries equation is not changed in sign in the region after the step, the transmitted wave evolves into a group of solitons and radiation, a process similar to soliton fission for surface gravity waves at a step. But if the coefficient of the nonlinear term changes sign, the soliton is destroyed completely and transforms into radiation. The effects of cubic nonlinearity are studied in the framework of the extended Korteweg-de Vries (Gardner) equation which is also integrable. The higher-order nonlinear effects influence the amplitudes of the generated solitons if the amplitude of the transformed wave is comparable with the thickness of lower layer, but otherwise the process of soliton fission is qualitatively the same as in the framework of the Korteweg-de Vries equation.     

On the Seddon speed formula

Abstract

The Seddon speed formula expressed mathematically as c = dQ/dA (or alternatively, as c = (1 + m)u_o; where Q is the discharge, A is the area of cross-section, c is the wave speed, u_o is the normal flow velocity, and m is a dimensionless parameter) is revisited in the context of elasticity and thermodynamics. Its link with the linearized solution of St Venant's equations for wave celerity, which does not appear to have been reported in the hydrological literature, is established. The rating exponent m is shown to be equivalent to the dimensionless relative celerity and is found to be the ratio of two specific heats, viz. c_p and c_v which are the specific heats at constant pressure and volume, respectively. The use of the parameter m as a complex variable helps describe shallow wave characteristics, the damping capacity of a wave, and the mechanism of occurrence of the hysteretic phenomenon. The damping capacity is found to describe the magnitude of wave subsidence, whereas the hysteresis also describes the speed of subsidence.