Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

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.     

Towards a method for postcolonial development geography? Possibilities and challenges

In this paper we explore the contours of a 'method' for postcolonial development geography, which makes it possible to imagine another 'world-picturing'. We suggest three steps towards such a method. First, we propose that a postcolonial method involves thinking about why we are doing research in the south in the first place; how we come to and produce our questions; and how we analyze and represent our findings based on our subject positionings. Second, that we need to recognize theorization as an inherent part of method, rethink how we currently theorize and reconfigure our methods of theorization to address wider political aims. Problematizing theorization helps challenge the universalism of Eurocentric theories, thus enabling development geography to move towards more decolonized versions and visions. Finally, that this must be accompanied by firmer recognition of our multiple investments – personal, institutional and geopolitical – and how they frame the possibilities for change. These are some possible steps that we think can reconfigure the 'scholarly track' that postcolonial development geographers traverse.     

Dissolved trace elements in river water: spatial distribution and the influencing factor, a study for the Pearl River Delta Economic Zone, China

Twenty-nine water samples were collected from different river channels of the Pearl River Delta Economic Zone, China. An inductively coupled plasma-mass spectromonitor (ICP-MS) was used to measure concentrations of the trace elements in these samples. The results suggest that the average concentrations of rare earth elements in river water show an increasing trend from the West River, the North River, the rivers of the Pearl River Delta, and the Shenzhen River to the East River. Relatively high concentrations of heavy metals appear in the East River, the rivers of the Pearl River Delta and the Shenzhen River, while the West River and the North River have relatively low heavy metal concentrations. Trace element concentrations in samples collected near urban or industrial areas are much higher than those of samples collected from distant areas, away from urban and industrial areas. After natural conditions, human activities have significant influence on the trace element concentrations in river water. This trace element concentration’s spatial distribution in the river water from the Pearl River Delta Economic Zone is actually an integrated effect of natural conditions and human activity.     

云南元-磨公路某高堑坡灾害发生机制及防治对策

论文根据现场地质调查,应用地质力学分析方法,分析了元-磨公路K235 160~600段高堑坡地质条件、病害发生机制和诱发因素。堑坡所在区域位于红河深大断裂与哀牢山深大断裂所夹持的哀牢山褶皱带内,属强烈侵蚀切割的低中山区,河流下切严重,河流阶地不发育。南溪河及其支流小曼萨河、叉河依附于不同时期生成的构造发育而成,堑坡位于小曼萨河注入南溪河的交汇部位小曼萨河的左岸。小曼萨河及南溪河所夹持的山体,尤其是堑坡所在的部位,已严重破碎、松弛,岩石似大块石堆垒状,坡面上的树木均明显地向小曼萨河方向倾倒。经地质力学配套分析,堑坡所依附的山梁形成于红河深大断裂生成期(γ4),堑坡岩体受多期构造活动的作用,其中NW65/°SW64°及NW63/°SW35°两组倾向临空的构造面控制坡体的稳定。堑坡具备发生变形破坏的地质条件,分析认为由于组成坡体的岩土强度不足而引起既有构造面的大型坍塌。连续降雨及温差是诱发因素。基于目前国内边坡防治技术、施工能力,提出了以“治”或“避”为主的防治对策。     

基于有限差分法的抗滑桩计算机辅助设计

基于地基系数“m-m”法、“m-k”法、“k-k”法的原理,考虑桩顶和桩底边界条件以及桩在滑动面处位移、转角、弯矩和剪力的连续条件,可解得桩身各节点的位移和内力,提出了进行抗滑桩全桩内力计算的有限差分法。根据差分方程并用VB6.0编制了实用的计算程序,既可避免繁琐的查表计算,提高计算速度,又可提高计算精度,直观生动,真正实现了人机交互,在界面的引导下,设计人员可完成全部计算,并绘出内力图形和抗滑桩截面配筋图,使设计更方便快捷,该软件可以极大地提高生产效率,降低工程造价,从而实现抗滑桩的优化设计。最后采用上述方法对某滑坡的悬臂抗滑桩进行了设计与计算。     

福建平潭虎潮山滑坡成因机理分析与评价

虎潮山滑坡位于平宏公路12k 704~12k 884段北侧山坡体上,是平潭海岛通往陆地的唯一公路通道。论文通过定性的工程地质分析和定量计算评价了虎潮山坡体稳定性的变化趋势,提出改建公路超挖深切坡脚和持续长时间强降雨导致的孔隙水压力作用是边坡发生滑动破坏的主要因素。