Simulation of the direct runoff hydrograph at basin outlet

A simple model describing the transformation of effective rainfall to direct runoff through the overland flow mechanism is presented. The model is based on the classical representation of a watershed by a combination of planes and channels. The dynamics of overland flow in each plane is simulated by the non-linear kinematic wave, but the outflow from a given plane is concentrated in the middle of the corresponding drainage channel. The water routing in the channels is carried out by a piece-wise linearized formulation in space of the kinematic wave approximation. Using synthetic events on 10 watersheds, the model was tested by comparing it with results obtained by applying the non-linear kinematic wave to all the elements of the watershed. The model was found to be adequate, even in a form that simplifies the geometric features of the planes through an averaging procedure based on the Horton–Strahler ordering scheme of the watershed. © 1998 John Wiley & Sons, Ltd.     

浅埋地下结构和土层在动荷载作用下的反应分析

在日本阪神大地震中,地下结构遭到了严重破坏。正确分析地下结构的地震反应,提出合理的抗震设计和安全性评价方法十分必要。用FLAC v4．0对某地的一个地铁车站进行了动力分析,发现随着埋深变浅,土体的加速度时程曲线表现出以某一确定的频率振动,该振动频率与模型的基频相近;位移在正方向和负方向的振幅都随埋深的增加而变小：在地表处达到最大值。一些认识和结论对于正确的进行抗震设计有指导意义。     

Model of a strong volcanic blast and a method of estimating the mass ejected

A hydrodynamic model of a strong volcanic blast associated with large (caldera-forming) ash–gas-flow eruptions is presented. We show that the mass erupted from a vent influences considerably the damping of the shock wave. If the ejected mass flux is constant, the shock amplitude as a function of distance will decrease slower (in accordance with the inverse 4/3 power law) than that of an atomic explosion (the inverse cube power law). From the modelling, we extract a relation connecting the effective energy of a blast with the mass ejected during the blast phase. The mass released during the largest explosion of the 1883 eruption of Krakatau, and the lateral blast energy at Mount St. Helens on 1980 May 18 are estimated on the basis of the equation obtained.     

结合被动微波和光学遥感的“超级寒潮”广东 地表温度特征及其对马铃薯生长的影响

利用被动微波可穿透云层的优势,基于GCOM-W1/AMSR2（Global Change Observation Mission-Water/ Advanced Microwave Scanning Radiometer 2）星载被动微波传感器的多通道亮度温度数据,反演了2016年1月下旬“超级寒潮”过程的广东省每日地表温度,分析了寒潮过程中广东省地表温度的时空变化特征。结果表明：在寒潮前后,广东省地表温度呈现出明显的先下降而后上升的趋势,寒潮过程带来的地表最低温度在广东北部大部分地区达到了277 K（4℃）,西南大部分地区为280~282 K（7~9℃）,其余大部分地区278~279 K（5~6℃）;降温幅度在广东南部地区达到了8~12 K（8~12℃）,北部大部分地区为5~7 K（5~7℃）。此外,利用归一化植被指数（NDVI）对植被生长的指示作用,基于HJ-1 A/B卫星的CCD传感器的多通道反射率数据,计算了近3年冬季广东省典型马铃薯种植区域马铃薯叶片的NDVI值,结果表明：马铃薯叶片的NDVI值在此次寒潮过程后出现了明显的下降,大部分地区（55.4%）马铃薯NDVI降低了0.1~0.2,部分地区（17.2%）马铃薯NDVI下降了0.2~0.3;进一步对比此次“超级寒潮”后同时段的前2年数据,发现寒潮后马铃薯叶片NDVI相对于正常年份也低了0.2,表明马铃薯叶片NDVI下降的主要原因是寒潮带来的低温,而不只是马铃薯的正常老化现象。NDVI的降低表明大部分地区马铃薯生长受到了抑制,马铃薯生长状况的这一变化也通过现场调查得到了证实。     

Generalised observations of wave characteristics on near‐horizontal shore platforms: Synthesis of six case studies from the North Island,New Zealand

Most field studies of wave processes on shore platforms in front of eroding cliffs focus on a single site, revealing aspects of wave dynamics at that location. Here, we analyse data from six platforms around northeastern New Zealand and describe the fundamental control of shore platform width, gradient and elevation on wave processes, including greater attenuation of short‐period waves at lower tidal stages and increases in longer period wave energy towards the cliff toe. These data suggest that empirical formulae developed from coral‐reef environments provide better predictions of wave height on platforms than formulae currently used in shore platform models.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.