Rainfall intensity affects runoff responses in a semi-arid catchment

The source and hydrochemical makeup of a stream reflects the connectivity between rainfall, groundwater, the stream, and is reflected to water quantity and quality of the catchment. However, in a semi-arid, thick, loess covered catchment, temporal variation of stream source and event associated behaviours are lesser known. Thus, the isotopic and chemical hydrographs in a widely distributed, deep loess, semi-arid catchment of the northern Chinese Loess Plateau were characterized to determine the source and hydrochemical behaviours of the stream during intra-rainfall events. Rainfall and streamflow were sampled during six hydrologic events coupled with measurements of stream baseflow and groundwater. The deuterium isotope (²H), major ions (Cl⁻, SO₄²⁻, NO₃⁻, Ca²⁺, K⁺, Mg²⁺, and Na⁺) were evaluated in water samples obtained during rainfall events. Temporal variation of ²H and Cl⁻ measured in the groundwater and stream baseflow prior to rainfall was similar; however, the isotope compositions of the streamflow fluctuated significantly and responded quickly to rainfall events, likely due to an infiltration excess, overland dominated surface runoff during torrential rainfall events. Time source separation using ²H demonstrated greater than 72% on average, the stream composition was event water during torrential rainfall events, with the proportion increasing with rainfall intensity. Solutes concentrations in the stream had loglinear relationships with stream discharge, with an outling anomaly with an example of an intra-rainfall event on Oct. 24, 2015. Stream Cl⁻ behaved nonconservative during rainfall events, temporal variation of Cl⁻ indicated a flush and washout at the onset of small rainfall events, a dilution but still high concentration pattern in high discharge and old water dominated in regression flow period. This study indicates rainfall intensity affects runoff responses in a semi-arid catchment, and the stored water in the thick, loess covered areas was less connected with stream runoff. Solute transport may threaten water quality in the area, requiring further analysis of the performance of the eco-restoration project.     

A New Wavelet Method for Identification of Eddies and Assessment of Incidents on Islands of the Eastern Caribbean

Using two dimensional continuous wavelet transforms, a novel method for identification of mesoscale eddies is presented to facilitate extraction of characteristics for area, amplitude, type, and location from maps of sea level anomalies. In comparison with the previously established growing method for eddy identification, it is found that the wavelet method identifies more than twice the number of eddies and is particularly better at resolving small eddies down to the 0.25 degree resolution of the data. Such research into eddy identification and tracking is significant to the assessment of eddies with potential to impact on coastlines of small islands. The method is applied to the identification of eddies on tracks towards islands of the Eastern Caribbean over 23?years. Spatial and temporal variation in rate of occurrence and magnitude is established. For Barbados there is an average of 9 anticyclonic incidents a year with maximum amplitude of typically 0.22?m in the dry seasons and 0.16?m in the wet seasons. Seasonal variation is reversed for the other islands with twice the number of anticyclonic incidents having maximum amplitudes of about 0.20?m annually.     

复杂地形景区三维导览图的设计与实现

为满足复杂地形景区对三维地图导览的需求,构建更为清晰的三维地图模型,提供良好的三维导览地图设计方案显得尤为重要。为了更好地与卫星影像进行贴合,DEM数据的采样密度要与卫星影像的分辨率一致。本文通过对比当下常用的四种空间插值方法的适用范围与运用特点,设置相关阈值及权重,直观比较了四种插值方法产生的插值结果,选择更适合复杂景区DEM插值的插值方法并对其进行精度提升,以此得到符合精度要求的DEM数据。最后利用相应地区的卫星影像进行地图投影及影像贴图,两种数据结合构建复杂景区的三维地图模型,给游客带来更为精确直观的定位信息和空间要素信息。     

Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

沉积物中硬质粘土层的研究进展

大量的钻孔资料表明,平原地区及其海域沉积物中普遍有硬质粘土层的存在。硬质粘土层在工程地质、地理学研究等方面具有重要的价值,它不仅是高层建筑的持力层,而且是研究全球环境变化事件、古气候演变的良好载体。该文阐述了硬质粘土层的一般特征、成因以及在粒度、地球化学元素、年代学、沉积环境等方面的研究进展,并展望了硬质粘土层在年代、物质来源和沉积环境等方面的研究方向,认为黄河三角洲硬质粘土层研究较少,需要开展大范围的、系统深入的、多学科的宏观和微观综合分析,挖掘硬质粘土层蕴含的古环境意义和应用价值。     

Smectite formation in metalliferous sediments near the East Pacific Rise at 13°N

A 43 cm long E271 sediment core collected near the East Pacific Rise(EPR) at 13°N were studied to investigate the origin of smectite for understanding better the geochemical behavior of hydrothermal material after deposition.E271 sediments are typical metalliferous sediments. After removal of organic matter, carbonate, biogenic opal,and Fe-Mn oxide by a series of chemical procedures, clay minerals(2 μm) were investigated by X-ray diffraction,chemical analysis and Si isotope analysis. Due to the influence of seafloor hydrothermal activity and close to continent, the sources of clay minerals are complex. Illite, chlorite and kaolinite are suggested to be transported from either North or Central America by rivers or winds, but smectite is authigenic. It is enriched in iron, and its contents are highest in clay minerals. Data show that smectite is most likely formed by the reaction of hydrothermal Fe-oxyhydroxide with silica and seawater in metalliferous sediments. The Si that participates in this reaction may be derived from siliceous microfossils(diatoms or radiolarians), hydrothermal fluids, or detrital mineral phases. And their δ30 Si values are higher than those of authigenic smectites, which implies that a Si isotope fractionation occurs during the formation because of the selective absorption of light Si isotopes onto Feoxyhydroxides. Sm/Fe mass ratios(a proxy for overall REE/Fe ratio) in E271 clay minerals are lower than those in metalliferous sediments, as well as distal hydrothermal plume particles and terrigenous clay minerals. This result suggests that some REE are lost during the smectite formation, perhaps because their large ionic radii of REE scavenged by Fe-oxyhydroxides preclude substitution in either tetrahedral or octahedral lattice sites of this mineral structure, which decreases the value of metalliferous sediments as a potential resource for REE.     

杭州湾中部实测波浪特性分析

利用杭州湾中部一年实测波浪资料,分析该地区的波参数统计变化特性,采用最小二乘法拟合得出相关参数之间的关系式,并统计分析了频谱特性。结果表明:杭州湾中部以小浪和轻浪为主,在各月分布较为均匀,年平均周期为2.97 s;最大一次波浪过程为冷空气影响所致,影响强度大于台风莫拉克;常浪向分布在东北至东南向,其中又以小浪出现较多,对应的波周期以2~4 s为主;强浪向主要分布在西北和东北方向,对应的波周期主要分布在3~5 s;多数特征波参数之间相关关系较好;由谱分析得知风浪占多数,且以单峰为主。研究结果可为相关工程和理论研究提供基础资料。