Numerical simulation of seawater intrusion from estuary into river using a coupled modeling system

One of the most important causes of the freshwater shortage in estuarine area is the increasing seawater intrusion into the river. To simulate seawater intrusion properly, two important factors should be considered. One is the bidirectional and time-dependent coupling effects between river discharges and tidal forces. The other is the three-dimensional and stratified structure of dynamic processes involved. However, these two factors have rarely been investigated simultaneously, or they were often simplified in previous researches, especially for the estuary connected with an upstream river network through multiple outlets such as the Pearl River Estuary (PRE). In order to consider these two factors, a numerical modeling system, which couples a one-dimensional river network model with a three-dimensional unstructured-grid Finite-Volume Coastal Ocean Model (FVCOM), has been developed and successfully applied to the simulation of seawater intrusion into rivers emptying into the PRE. By treating the river network with a one-dimensional model, computational efficiency has been improved. With coupling 1D and 3D models, the specification of up-stream boundary conditions becomes more convenient. Simulated results are compared with field measured data. Good agreement indicates that the modeling system may correctly capture the physical processes of seawater intrusion into rivers.     

ADI barotropic ocean model for simulation of Kuroshio intrusion into China southeastern waters

In this numerical model for simulating the Kuroshio intrusion into the East and South China Seas, vertically averaged marine hydrodynamic equations governing ocean currents and long-period waves are approximated by a set of two-time-level semi-implicit fimite difference equations. The major terms including the local acoeleration, sea-surface slope, Coriolis force and the bottom friction are approximated with the Crank-Nicholson scheme, which is of second order accuracy. The advection terms are approximated with the Leith scheme. The difference equations are split into two sets of alternating direction implicit equations, each of which has a tridiagonal matrix and can be easily solved. The model reproduces a major Kuroshio intrusion north of Luzon Island, one north of Taiwan Island, and one west of the Tokara Strait. The model shows a current system running from the Luzon Strait to the coast of Vietnam and Hainan Island, through the Taiwan Strait and then into the Tsushima Strait. The summer and winter monsoons generate several eddies in the South China Sea. Project supported by the National Natural Science Foundation of China.     

Evaluation of a regional climate model for atmospheric simulation over Arctic river basins

Evaluation on a regional climate model was made with five-month atmospheric simulations over the Arctic river basins. The simulations were performed with a modified mesoscale model, Polar MM5 coupled to the NCAR Land Surface Model （LSM） to illustrate the skill of the coupled model （Polar MM5＋LSM） in simulating atmospheric circulation over the Arctic river basins. Near-surface and upper-air observations were used to verify the simulations. Sensitivity studies between the Polar MM5 and Polar MM5＋LSM simulations revealed that the coupled model could improve the forecast skill for surface variables at some sites. In addition, the extended evaluations of the coupled model simulations on the North American Arctic domain during December 15, 2002 to May 15, 2003 were carried out. The time series plots and statistics of the observations and Polar MM5＋LSM simulations at six stations for near-surface and vertical profiles at 850 hPa and 500 hPa were analyzed. The model was found capable of reproducing the observed atmospheric behavior in both magnitude and variability, especially for temperature and near-surface wind direction.     

基于滑面分区段力学模型的高速滑坡运动过程能量转化研究

高速滑坡具有运动速度快、波及范围广的灾害特征,因此对滑坡启动、加速和静止整个运动过程进行研究很有必要。基于滑面力学特性将滑面分为弹性介质区和应变弱化区,构建了高速滑坡二维力学模型,提出了滑坡启动动能计算公式和滑坡运动过程的能量计算公式;以千将坪滑坡为例,采用启动动能计算公式得出滑坡的启动速度为2.35 m/s;依据滑面形态将其运动轨迹划分为快速加速、平稳加速、平稳减速、急剧减速4个阶段,进行运动过程分析,得出滑坡最大速度为16.8 m/s,以滑坡前缘高程所在平面为势能基准面,分析不同能量与总能量占比的变化情况,在滑坡的4个运动阶段中,动能占比分别为9.1%,25.6%,15.1%,0%;摩擦损耗能量占比分别为：5.5%,58.8%,81.7%,95.5%;势能占比分别为：85.2%,14.2%,0%,0%;其他阻力能耗占比分别为：0.2%,1.4%,3.2%,4.5%。研究结论对高速滑坡致灾机制和风险分析具有重要意义。     

Effects of river flow velocity on the formation of landslide dams

Natural dams are formed when landslides are triggered by heavy rainfall during extreme weather events in the mountainous areas of Taiwan.During landslide debris movement, two processes occur simultaneously: the movement of landslide debris from a slope onto the riverbed and the erosion of the debris under the action of high-velocity river flow. When the rate of landslide deposition in a river channel is higher than the rate of landslide debris erosion by the river flow, the landslide forms a natural dam by blocking the river channel. In this study, the effects of the rates of river flow erosion and landslide deposition(termed the erosive capacity and depositional capacity, respectively) on the formation of natural dams are quantified using a physics-based approach and are tested using a scaled physical model.We define a dimensionless velocity index vde as the ratio between the depositional capacity of landslide debris(vd) and the erosive capacity of water flow(ve).The experimental test results show that a landslidedam forms when landslide debris moves at high velocity into a river channel where the river-flow velocity is low, that is, the dimensionless velocity index vde 54. Landslide debris will not have sufficient depositional capacity to block stream flow when the dimensionless velocity index vde 47. The depositional capacity of a landslide can be determined from the slope angle and the friction of the sliding surface, while the erosive capacity of a dam can be determined using river flow velocity and rainfall conditions. The methodology described in this paper was applied to seven landslide dams that formed in Taiwan on 8 August 2009 during Typhoon Morakot,the Tangjiashan landslide dam case, and the YingxiuWolong highway K24 landslide case. The dimensionless velocity index presented in this paper can be used before a rainstorm event occurs to determine if the formation of a landslide dam is possible.     

Damming of large river by debris flow: Dynamic process and particle composition

The frequency and extent of debris flows have increased tremendously due to the extreme weather and the Wenchuan earthquake on May 12, 2008. Previous studies focused on the debris flow from gullies damming the mountain streams. In this paper, an equation for the run-out distance of debris flow in the main river is proposed based on the dynamic equation of debris flow at different slopes given by Takahashi. By undertaking field investigations and flume experiments, a new calculation method of the volume of debris flow damming large river is obtained. Using the percolation theory and the renormalization group theory it was deduced that the large particles should comprise more than 50% for forming a stable debris flow dam. Hence, the criteria of damming large river by debris flow is presented in terms of run-out distance and grain composition which was then validated through the event of damming river by debris flow at Gaojia gully, the upper reaches of the Minjiang River, Sichuan, China, on July 3, 2011.     

Dynamic process analysis of the Niumiangou landslide triggered by the Wenchuan earthquake using the finite difference method and a modified discontinuous deformation analysis

The Niumiangou landslide was the largest landslide triggered by the 2008 Wenchuan earthquake, which was significantly affected by the amplification effect of seismic acceleration. The ringshear experiments indicated that the materials in the source area of the Niumiangou landslide were subjected to friction degradation under a big shear displacement, which may result in rapid movement of the landslide. In order to better understand the landslide movement and study the effect of the friction degradation on movement mechanisms, the dynamic process of Niumiangou landslide was simulated with a new numerical method, which combines the finite difference method(FDM) and the discontinuous deformation analysis(DDA). First, the FDM was used to study the initiation time, amplification effect and velocity of the landslide. Afterwards, these initiation velocities were applied to the blocks in the DDA model by corresponding coordination in the FDM model. A displacementdependent friction model of the sliding surface was incorporated into DDA code to further understand the kinetic behavior of the landslide. The results show that the displacement-dependent friction strongly decreases the friction coefficient of sliding surface under a big displacement, which can obviously promote the run-out and velocity of landslide. The model output well matches the topographic map formed by the landslide. This implies that the proposed model can be applied to the simulation of earthquake-induced landslides with amplification effect, and the friction degradation model is important to clarify the movement mechanism of high-speed and long-distance landslides.     

考虑真实气体的低渗气藏动态预测模型

针对气藏的生产特征,基于稳定渗流理论,以气体高速非达西渗流为基础,考虑真实气体PVT参数动态变化,结合气藏物质平衡方程,推导建立同时考虑滑脱效应、启动压力梯度、应力敏感等因素共同影响的低渗气藏动态预测模型,分析各因素对气井产量变化的影响.结果表明:在考虑高速非达西流的基础上,考虑滑脱效应将导致渗流阻力降低而引起气井稳产年限增加和地层压力下降变快,考虑启动压力梯度和应力敏感效应时将导致渗流阻力增加而引起气井稳产年限降低和地层压力下降变慢;启动压力梯度和渗透率变形系数越大,气井稳产年限下降幅度越大,地层压力下降越慢,而滑脱因子越大,气井稳产年限增加幅度越大,地层压力下降越快.     

多功能高速远程滑坡运动堆积过程物理模型试验装置设计与应用

高速远程滑坡运动堆积过程影响因素众多,在物理模型试验装置研发过程中应满足多因素变化需求,从而实现多功能目的.该试验装置初始状态下(垂直角度=20°,水平角度=0°)长3.40m,宽0.56m,高1.35m.设计有滑体体积调节、上滑槽坡度调节和下滑槽水平角度调节三大模块.材料上采用了3mm厚不锈钢钢板与8mm厚钢化玻璃2种材料,使用耐久性较好,并且采用分部件组装形式,安装有万向轮,便于实验仪器的搬运.基于砂子与卵石颗粒材料,应用所研发装置,初步开展了高速远程滑坡运动堆积过程物理模型试验,简要分析了滑体坡度、水平角度、滑体高度及基底材质等参数对滑动距离的影响规律.     

Observations of Kuroshio intrusion into the South China Sea

To discuss the intrusion of the Kuroshio into the SCS, we examined the mixing between the North Pacific and South China Sea (SCS) waters based on in-situ CTD data collected in August and September 2008 and the moored ADCP data taken from mid September 2008 to early July 2009. The CTD survey included four meridional sections from 119°E to 122°E around the Luzon Strait, during which pressure, temperature, and salinity were measured. The CTD data show that the isopycnal surface tilted from the SCS to the North Pacific; and it was steeper in the lower layers than in the upper ones. Meanwhile, we found strong vertical mixing taken place in the areas near 121°E. The Kuroshio in high temperature and salinity intruded westward through Luzon Strait. The frequency of buoyancy was one order of magnitude greater than that of the common ones in the ocean, suggesting stronger stratification in the northeastern SCS. On the other hand, the long-term ADCP data show that before late October 2008, the direction of water flow in the SCS was eastward, and from November 2008 to late February 2009, it turned northwestward in the layers shallower than 150 m, while remained unchanged in deep layers from 200 to 450 m. From March to June 2009, the direction shifted with increasing depth from northward to southward, akin to the Ekman spiral. EOF analysis of the current time series revealed dominant empirical modes: the first mode corresponded to the mean current and showed that the Kuroshio intrusion occurred in the upper layers only from late December to early March. The temporal coefficient of the first and the second mode indicated clearly a dominant signal in a quasi-seasonal cycle.     

Numerical simulation of a high-speed landslide in Chenjiaba,Beichuan, China

High-speed landslide is a catastrophic geological disaster in the mountainous area of southwest China. To predict the movement process of landslide reactivation in Chenjiaba town, Beichuan county, Sichuan province, China, we simulated the movement process of two landslide failures in Chenjiaba via rapid mass movement simulation and unmanned aerial vehicle images(UAV), and obtained the movement characteristic parameters of the landslides. According to a back analysis, the most remarkable fitting rheological parameters were friction coefficient(μ=0.18) and turbulence(). The parameter of landslide pressure was applied as the zoning index of landslide hazard to obtain the influence zone and hazard zoning map of the Chenjiaba landslide. Results show that the Duba River was blocked quickly with a landslide accumulation at the maximum height of 44.14 mwhen the Chenjiaba deposits lost stability. The hazard zoning map indicated that the landslide hazard degree is positively correlated with the slope.This landslide assessment is a quantitative hazard assessment method based on a landslide movement process and is suitable for high-speed landslide. Such method can provide a scientific basis for urban construction and planning in the landslide hazard area to avoid hazards effectively.     

Dynamic model for oil slick dispersion into a water column—A wind-driven wave tank experiment

A dynamic experiment for oil dispersion into a water column was performed with a 21 m long, 0.5 m wide, and 1 m high wind-driven wave tank. At wind velocity between 6–12 m/s and with the oil slick kept constant (about 1 μm), the rate of the oil content increase in the water column could be approximated from the difference between the dispersion rate (R) of the oil slick and the coagulation rate (R′) of the dispersed oil slick. Assuming the coagulation rate is directly proportional to the concentration of the water dispersed oil slick (i. e.R′-K ^*C), the integral form of the dynamic model can be expressed asC=R ^* [1−exp(−K ^*t)]/K and parametersR andK can be regressed with a computer. The relative deviation of model results from the experimental data was mainly less than 10%. The oil slick dispersion rate (R) had exponential relationship with the wind velocity (U), and can be fitted with a formulaR=A ^* (U+1) ^B . The fitted constant of the coagulation rate,K(0.8–3.0^* 10⁻³ min⁻¹) did not have significant relationship with the wind velocity.     

Numerical simulation on the process of saltwater intrusion and its impact on the suspended sediment concentration in the Changjiang (Yangtze) estuary

To study the relationship between sediment transportation and saltwater intrusion in the Changjiang (Yangtze) estuary, a three-dimensional numerical model for temperature, salinity, velocity field, and suspended sediment concentration was established based on the ECOMSED model. Using this model, sediment transportation in the flood season of 2005 was simulated for the Changjiang estuary. A comparison between simulated results and observation data for the tidal level, flow velocity and direction, salinity and suspended sediment concentration indicated that they were consistent in overall. Based on model verification, the simulation of saltwater intrusion and its effect on sediment in the Changjiang estuary was analyzed in detail. The saltwater intrusion in the estuary including the formation, evolution, and disappearance of saltwater wedge and the induced vertical circulation were reproduced, and the crucial impact of the wedge on cohesive and non-cohesive suspended sediment distribution and transportation were successfully simulated. The result shows that near the salinity front, the simulated concentrations of both cohesive and non-cohesive suspended sediment at the surface layer had a strong relationship with the simulated velocity, especially when considering a 1-hour lag. However, in the bottom layer, there was no obvious correlation between them, because the saltwater wedge and its inducing vertical circulation may have resuspended loose sediment on the bed, thus forming a high-concentration area near the bottom even if the velocity near the bottom was very low during the transition phase from flood to ebb.     

Numerical simulation on the process of saltwater intrusion and its impact on the suspended sediment concentration in the Changjiang(Yangtze)estuary

To study the relationship between sediment transportation and saltwater intrusion in the Changjiang (Yangtze) estuary, a three-dimensional numerical model for temperature, salinity, velocity field, and suspended sediment concentration was established based on the ECOMSED model. Using this model, sediment transportation in the flood season of 2005 was simulated for the Changjiang estuary. A comparison between simulated results and observation data for the tidal level, flow velocity and direction, salinity and suspended sediment concentration indicated that they were consistent in overall. Based on model verification, the simulation of saltwater intrusion and its effect on sediment in the Changjiang estuary was analyzed in detail. The saltwater intrusion in the estuary including the formation, evolution, and disappearance of saltwater wedge and the induced vertical circulation were reproduced, and the crucial impact of the wedge on cohesive and non-cohesive suspended sediment distribution and transportation were successfully simulated. The result shows that near the salinity front, the simulated concentrations of both cohesive and non-cohesive suspended sediment at the surface layer had a strong relationship with the simulated velocity, especially when considering a 1-hour lag. However, in the bottom layer, there was no obvious correlation between them, because the saltwater wedge and its inducing vertical circulation may have resuspended loose sediment on the bed, thus forming a high-concentration area near the bottom even if the velocity near the bottom was very low during the transition phase from flood to ebb.     

Rapid prediction models for 3D geometry of landslide dam considering the damming process

The geometry of a landslide dam plays a critical role in its stability and failure mode, and is influenced by the damming process. However, there is a lack of understanding of the factors that affect the 3D geometry of a landslide dam. To address this gap,we conducted a study using the smoothed particle hydrodynamics numerical method to investigate the evolution of landslide dams. Our study included 17 numerical simulations to examine the effects of several factors on the geometry of landslide d...     

远源细粒辫状河心滩坝演化与河流分叉的交互沉积过程:现代沉积启示与数值模拟分析

远源细粒辫状河广泛发育于河流的中下游，其形成的致密储层是深入挖潜油气面临的挑战。结合现代远源细粒辫状河演变的遥感记录影像，运用基于泥沙水动力学的沉积过程数值模拟软件Delft3D正演远源细粒辫状河的生长演化过程，并探讨物源供给对远源细粒辫状河形态的影响，深入剖析和解释远源细粒辫状河心滩坝演化与河流分叉的交互作用和沉积过程，并进行了现代沉积验证。研究表明:①河道沿下游方向由窄变宽是触发远源细粒辫状河形成的主要地形条件。河流流速和河岸牢固程度是决定河流宽度的重要因素。②远源细粒辫状河心滩坝演化与河流分叉的交互沉积过程是心滩坝泥沙与河流流水交互作用的沉积演化过程。心滩坝泥沙与河流流水交互作用具体包括泥沙在河流中的沉积作用和河流对泥沙的侵蚀作用。③远源细粒辫状河心滩坝演化与河流分叉的交互沉积过程包括3种:河道内心滩坝加积导致河流分叉、心滩坝与河岸分离导致河流分叉、心滩坝的局部冲裂分解导致河流分叉。该研究丰富了远源细粒辫状河沉积理论，并且有助于分析远源细粒辫状河沉积相的储层空间展布规律并应用于油气的深入勘探开发。      

Method for determining the position of landslide slip-surface with a typical inclinometric curves

In the actual monitoring of deep hole displacement, the identification of slip surfaces is primarily based on abrupt changes observed in the inclinometric curve. In conventional identification methods, inclinometric curves exhibiting indications of sliding can be categorized into three types: B-type, D-type, and r-type. The position of the slip surface is typically determined by identifying the depth corresponding to the point of maximum displacement mutation. However, this method is sensitive t...     

Snow cover variation and streamflow simulation in a snow-fed river basin of the Northwest Himalaya

Snowmelt is an important component of any snow-fed river system.The Jhelum River is one such transnational mountain river flowing through India and Pakistan.The basin is minimally glacierized and its discharge is largely governed by seasonal snow cover and snowmelt.Therefore,accurate estimation of seasonal snow cover dynamics and snowmeltinduced runoff is important for sustainable water resource management in the region.The present study looks into spatio-temporal variations of snow cover for past decade and stream flow simulation in the Jhelum River basin.Snow cover extent(SCE) was estimated using MODIS(Moderate Resolution Imaging Spectrometer) sensor imageries.Normalized Difference Snow Index(NDSI) algorithm was used to generate multi-temporal time series snow cover maps.The results indicate large variation in snow cover distribution pattern and decreasing trend in different sub-basins of the Jhelum River.The relationship between SCE-temperature,SCE-discharge and discharge-precipitation was analyzed for different seasons and shows strong correlation.For streamflow simulation of the entire Jhelum basin Snow melt Runoff Model(SRM) used.A good correlation was observed between simulated stream flow and in-situ discharge.The monthly discharge contribution from different sub-basins to the total discharge of the Jhelum River was estimated using a modified version of runoff model based on temperature-index approach developed for small watersheds.Stream power - an indicator of the erosive capability of streams was also calculated for different sub-basins.     

三峡库区麻柳林滑坡变形特征及演化模拟

在库水位波动和降雨作用的共同影响下，库岸滑坡的变形规律往往更为复杂。以三峡库区麻柳林滑坡为例，基于野外调查、钻探编录、深部位移监测以及数值模拟等手段，分析了库水位波动和降雨作用下滑坡变形特征及演化规律。结果表明:麻柳林滑坡在粉质黏土层和块石层交界处发育一个次级滑带，目前该滑坡主要沿次级滑带运动，导致次级滑动的原因与坡体物质的差异性有关；S_i(S_f)指标分析法揭示滑坡的滑带还未完全破坏，滑坡仍处于蠕变状态；根据三峡水库水位调度规律，将一个完整水文年划分为6个阶段，数值模拟结果表明滑坡在库水位缓慢下降阶段变形速率较小、在快速下降阶段和低水位阶段变形速率持续增大、在快速上升阶段和缓慢上升阶段以及高水位阶段变形速率则保持平稳。其中，降雨的直接影响和降雨导致库水位波动进而对滑坡变形造成的间接影响，使得麻柳林滑坡在低水位阶段的变形显著增加、稳定性最差，应加强该时段内滑坡的监测和预警。