DISCUSSION: ‘HYDROLOGICAL MODEL OF PEAT-MOUND FORM WITH VERTICALLY VARYING HYDRAULIC CONDUCTIVITY’ BY ADRIAN C. ARMSTRONG

In this comment we argue that the premise on which the peat mound model developed by Armstrong (Earth Surface Process and Landforms, 1995, 20 , 473–477) is based, that hydraulic conductivity shows an exponential decline with depth in bog peats, is unsound. Empirical evidence in the literature for such an exponential decline is less sound than Armstrong suggests. In addition, Armstrong's suggestion that the hypothesis of Baird and Gaffney (Earth Surface Processes and Landforms, 1995, 20 , 561–566) supports an exponential decline is shown to be erroneous.     

VELOCITY–DISCHARGE RELATIONSHIPS DERIVED FROM DYE TRACER EXPERIMENTS IN GLACIAL MELTWATERS: IMPLICATIONS FOR SUBGLACIAL FLOW CONDITIONS

Repeated dye tracer tests were undertaken from individual moulins at Haut Glacier d'Arolla, Switzerland, over a number of diurnal discharge cycles during the summers of 1989–1991. It was hoped to use the concepts of at-a-station hydraulic geometry to infer flow conditions in subglacial channels from the form of the velocity–discharge relationships derived from these tests. The results, however, displayed both clockwise and anticlockwise velocity–discharge hysteresis, in addition to the simple power function relationship assumed in the hydraulic geometry approach. Clockwise hysteresis seems to indicate that a moulin drains into a small tributary channel rather than directly into an arterial channel, and that discharges in the two channels vary out of phase with each other. Anticlockwise hysteresis is accompanied by strong diurnal variations in the value of dispersivity derived from the dye breakthrough curve, and is best explained by hydraulic damming of moulins or sub/englacial passageways. Despite the complex velocity–discharge relationships observed, some indication of subglacial flow conditions may be obtained if tributary channels comprise only a small fraction of the drainage path and power function velocity–discharge relationships are derived from dye injections conducted during periods when the supraglacial discharge entering the moulin and the bulk discharge vary in phase. Analyses based on this premise suggest that both open and closed channel flow occur beneath Haut Glacier d'Arolla, and that flow conditions are highly variable at and between sites.     

地心和月心引力常数及月球形心与质心的确定

指出了地心引力常数GMe、月心引力常数GMl及月球形状中心和月球质量中心的意义和重要性;介绍了用空间探测器观测数据测定GMe与GMl及月球形状中心和月球质量中心的原理和方法;综合给出了利用空间探测器测定的数值。     

平面点集凸壳的一种快速算法

提出一种计算平面点集凸壳的快速算法———八方向极值快速凸壳算法。该算法首先对平面点集进行一次扫描,从而快速查找到东、南、西、北、东南、西南、东北、西北8个方向上的极值点,构造出一个更接近凸壳的初始凸壳,从而在后续的点集扫描中可以排除更多的内点,使该算法计算效率更高。该算法的空间复杂度为O(N);其时间复杂度虽然无法突破最坏情况下O(NlogN)的理论下限,但其期望时间复杂度已达到线性水平,并且可以容易地扩展到三维和高维空间。     

Implementation of the semi-Lagrangian advection scheme on a quasi-uniform overset grid on a sphere

The semi-Lagrangian advection scheme is implemented on a new quasi-uniform overset （Yin-Yang） grid on the sphere. The Yin-Yang grid is a newly developed grid system in spherical geometry with two perpendicularly-oriented latitude-longitude grid components （called Yin and Yang respectively） that overlapp each other, and this effectively avoids the coordinate singularity and the grid convergence near the poles. In this overset grid, the way of transferring data between the Yin and Yang components is the key to maintaining the accuracy and robustness in numerical solutions. A numerical interpolation for boundary data exchange, which maintains the accuracy of the original advection scheme and is computationally efficient, is given in this paper. A standard test of the solid-body advection proposed by Williamson is carried out on the Yin-Yang grid. Numerical results show that the quasi-uniform Yin-Yang grid can get around the problems near the poles, and the numerical accuracy in the original semi-Lagrangian scheme is effectively maintained in the Yin-Yang grid.     

利用毛管压力曲线分形分维方法研究流动单元

利用取心井铸体薄片获得的图像资料和毛管压力曲线,通过图像分形几何学方法以分维数的形式定量地表征出了复杂的微观孔隙喉道结构特征,发现能够很好地划分和评价孔隙岩石中油、气、水的渗流差异,可以用于储层微观流动单元表征。文中阐述了岩石微观孔隙喉道结构分形的理论基础、计算方法和应用于表征流动单元的依据。建立了中国西部砾岩低渗透油藏微观孔隙喉道分维数与孔隙度、渗透率之间计算图版,据此在油藏中利用常规测井资料获得的孔隙度、渗透率参数计算微观孔隙喉道分维数,开展全油藏流动单元划分与评价,取得了良好的效果。研究结果表明,利用毛管压力曲线分形分维方法研究储层微观流动单元是一种很有效的途径。     

长江三峡巴东复杂斜坡系统成因研究

三峡库区巴东县新城区所在的扇形大斜坡是一个复杂斜坡系统。斜坡后缘边界的挤压变形带是斜坡整体沿其底界面（T1j^3-T2b^1）之间岩层界面向长江发生重力蠕滑作用的结果,不存在“巴东断裂”问题。斜坡系统内部未发现成规模的褶皱,也不存在数量众多的内动力成因的断裂带。巴东斜坡区浅表生变形破坏类型可概括为侧向拉裂、“雪崩式”垮塌、表层风化、侧向滑移张剪破裂、岩层折断变位、软岩膝折剪破裂和压扭破裂又重胶结7种重力卸荷破坏力学模式。斜坡系统内部发育分期分区（块）滑动（滑坡）现象,但不存在构造地质意义上的“断层”形迹。在地貌形态上,巴东大斜坡被4条冲沟分割成5个斜坡单元。在斜坡地质结构、成因类型与空间分布上,巴东大斜坡可划分为3个层次,即表层崩塌滑坡成因为主的堆积层（第一层次）,冲沟分割且浅表生地质形迹发育的层状顺倾的中间基岩层（第二层次）和整体连续顺倾的深层基岩层（第三层次）。分形几何计算证明,巴东大斜坡的地貌形态尚处于侵蚀发育的青（幼）年期。FLAC^3D数值模拟发现,在长江侵蚀下切的不同阶段,巴东斜坡体前缘和后缘接近底界面位置塑性变形区分布集中,但不具备沿深层界面发生整体滑动剪出的可能性。基本认识是,“巴东复杂斜坡系统”是在官渡口-东瀼口向斜南翼（单斜山）的地质背景下,持续经受长江快速侵蚀下切外动力作用,河谷岸坡快速临空导致其自身重力产生强烈的侧向卸荷与滑移等浅表生地质改造作用过程而形成的,可简单地概括为“重力成因论”。     

Simulating cross-sectional geometry of the main channel in response to changes in water and sediment in Lower Yellow River

Journal of Geographical Sciences - To understand the non-equilibrium morphological adjustment of a river in response to environmental changes, it is essential to (i) accurately identify how past...     

URANS investigation of the interaction between the free surface and a shallowly submerged underwater vehicle at steady drift

An axisymmetric underwater vehicle (UV) at a steady drift angle experiences the complex three-dimensional crossflow separation. This separation arises from the unfavorable circumferential pressure gradient developed from the windward side toward the leeward side. As is well known, the separated flow in the leeward side gives rise to the formation of a pair of vortices, which affects considerably the forces and moments acting on the UV. In this regard, the main purpose of the present study is to evaluate the role of the leeward vortical flow structure in the hydrodynamic behavior of a shallowly submerged UV at a moderate drift angle traveling beneath the free surface. Accordingly, the static drift tests are performed on the SUBOFF UV model using URANS equations coupled with a Reynolds stress turbulence model. The simulations are carried out in the commercial code STARCCM+ at a constant advance velocity based on Froude number equal to F_n = 0.512 over submergence depths and drift angles ranging from h = 1.1D to h = ∞ and from β = 0 to β = 18.11°, respectively. The validation of the numerical model is partially conducted by using the existing experimental data of the forces and moment acting on the totally submerged bare hull model. Significant interaction between the low-pressure region created by the leeward vortical flow structure and the free surface is observed. As a result of this interaction, the leeward vortical flow structure appears to be largely responsible for the behavior of the forces and moments exerted on a shallowly submerged UV at steady drift.     

Temporal variation of modulated-wave-train geometries and their influence on vertical bending moments of a container ship

A towing experiment was conducted using a modulated wave train to investigate the vertical bending responses of a hydro-structural container ship model. In the experiment, a spatially periodic modulated wave train, as a model of a freak wave in successive high waves mimicking the so-called three sisters, was generated by the recently established higher-order spectral method wave generation (HOSM-WG) method. HOSM-WG enables us to control the location and timing of the maximum crest height in a wave tank. With precise control of the towing carriage, an experiment was conducted in which the timing of the encounters between the ship model and the modulated wave train was accurately determined. The maximum sagging moment (SM) was found to increase in proportion with the encounter wave height. However, because of differences in the relative depth of the fore and aft troughs, the maximum SM is highly variable for a given wave height. The temporal wave-geometry evolution caused the relative trough-depth to vary significantly within a wave period in the vicinity of the maximum crest height. As a result, depending on the encounter timing, the SM varied considerably for a given wave height. The temporal variation of the wave geometry is a robust feature of a modulated wave train and is common between the spatially periodic and temporally periodic modulated wave trains.