Modelling spatial variations in dissolved oxygen in fine‐grained streams under uncertainty

This paper presents a novel triple‐layer model, called VART DO‐3L, for simulation of spatial variations in dissolved oxygen (DO) in fine‐grained streams, characterized by a fluid mud (fluff or flocculent) layer (an advection‐dominated storage zone) as the interface between overlying stream water and relatively consolidated streambed sediment (a diffusion‐dominated storage zone). A global sensitivity analysis is conducted to investigate the sensitivity of VART DO‐3L model input parameters. Results of the sensitivity analysis indicate that the most sensitive parameter is the relative size of the advection‐dominated storage zones (A_s/A), followed by a lumped reaction term (R) for the flocculent layer, biological reaction rate (μ_o) in diffusive layer and biochemical oxygen demand concentration (L) in water column. In order to address uncertainty in model input parameters, Monte Carlo simulations are performed to sample parameter values and to produce various parameter combinations or cases. The VART DO‐3L model is applied to the Lower Amite River in Louisiana, USA, to simulate vertical and longitudinal variations in DO under the cases. In terms of longitudinal variation, the DO level decreases from 7.9 mg l at the Denham Springs station to about 2.89 mg l^?1 at the Port Vincent station. In terms of vertical variation, the DO level drops rapidly from the overlying water column to the advection‐dominated storage zone and further to the diffusive layer. The DO level (C_F) in the advective layer (flocculent layer) can reach as high as 40% of DO concentration (C) in the water column. The VART DO‐3L model may be applied to similar rivers for simulation of spatial variations in DO level. Copyright © 2014 John Wiley & Sons, Ltd.     

A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.     

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.     

Soil‐pile interaction in the pile vertical vibration considering true three‐dimensional wave effect of soil

The dynamic response of an end bearing pile embedded in a linear visco‐elastic soil layer with hysteretic type damping is theoretically investigated when the pile is subjected to a time‐harmonic vertical loading at the pile top. The soil is modeled as a three‐dimensional axisymmetric continuum in which both its radial and vertical displacements are taken into account. The pile is assumed to be vertical, elastic and of uniform circular cross section. By using two potential functions to decompose the displacements of the soil layer and utilizing the separation of variables technique, the dynamic equilibrium equation is uncoupled and solved. At the interface of soil‐pile system, the boundary conditions of displacement continuity and force equilibrium are invoked to derive a closed‐form solution of the vertical dynamic response of the pile in frequency domain. The corresponding inverted solutions in time domain for the velocity response of a pile subjected to a semi‐sine excitation force applied at the pile top are obtained by means of inverse Fourier transform and the convolution theorem. A comparison with two other simplified solutions has been performed to verify the more rigorous solutions presented in this paper. Using the developed solutions, a parametric study has also been conducted to investigate the influence of the major parameters of the soil‐pile system on the vertical vibration characteristics of the pile. Copyright © 2013 John Wiley & Sons, Ltd.     

基于Google Earth土地承包经营权工作底图的制作

对利用Google Earth影像制作农村土地承包经营权工作底图的方法和技术流程进行了全面探讨。首先在CASS 7.0软件下对集体土地所有权行政界线进行坐标转换,再使用自编程序提取界线拐点的WGS-84大地坐标,通过"地图下载器"下载Google Earth影像并拼接输出TIF影像,利用Arc Map软件的投影变换功能,将墨卡托投影的TIF影像转换为高斯投影影像,最后,在CASS 7.0软件下插入变换后的正射影像,叠加行政界线,形成完整的承包经营权外业调查工作底图。     

移动测量系统宽幅影像的几何拼接与量测

研究了移动测量系统宽幅影像的拼接与量测算法。通过建立三幅影像柱面全景投影的几何关系,能够直接拼接完成宽幅影像。利用宽幅影像和原始影像间反投影变换关系以及三幅原始影像的同、异步立体像对,可建立宽幅影像与物方空间坐标间一一对应关系,从而实现宽幅影像目标的量测算法。宽幅影像的投影拼接和量测实验表明,利用柱面全景投影拼接得到宽幅影像的方法可应用于移动测量平台宽幅影像序列的拼接;在原始影像具备同步立体像对的条件下,宽幅影像量测的绝对精度可达1 m以内,相对精度可达0.2 m以内,满足城市环境实景目标量测的精度需求。     

毫米波雷达与无线电探空对云垂直结构探测的一致性分析

在众多种针对云垂直结构的探测中,毫米波雷达可获取云底、云顶、云厚等完整的云垂直结构信息,并可以连续监测云的垂直剖面变化,是有力的探测手段之一。而无线电探空因其直接的测量优势,能直观、确切地描述大气湿度垂直结构,可将其进一步处理生成云垂直结构信息,并作为一种数据源用于与毫米波雷达云垂直结构探测结果进行比对,以评估毫米波雷达针对云宏观垂直结构的探测性能,为毫米波雷达更好地应用于云探测提供参考。通过获取位于北京南郊观象台的毫米波雷达2014年10月28日至2015年2月17日长达113天连续观测的反射率因子以及探空温、压、湿数据,设计或选取合适的方法对二者进行云边界的计算,并进行云高（包括云底高和云顶高）以及云层数的一致性比对分析。结果认为,除对于高层云的云顶高度毫米波雷达由于探测限制不能探测到10 km以上的云顶,某些时刻与探空产生较大差异外,在云底高度以及中低云的云顶高度上可以与探空观测取得很好的匹配效果,对于云的垂直分层上二者也有较强的一致性。该毫米波雷达具有较为准确并连续刻画10 km以下云垂直结构的能力。     

京津冀地区中α尺度飑线过程中大风特征分析及成因初探

应用常规天气观测资料、地面加密自动气象站资料、大风灾情报资料、京津冀地区7部多普勒天气雷达组网观测资料及VDRAS资料,从多个角度对2013年8月4日京津冀地区一次飑线过程产生的大范围大风天气过程进行了分析,结果显示:此次过程是在高空冷空气南下、低层暖湿气流北上、系统前倾及位势不稳定的有利层结条件下,由多单体风暴演变为中α尺度的强飑线所致。飑线形成于低层垂直切变加强、冷池合并之后;大风主要发生在飑线主体回波中,其次是主体回波前和中前,主体回波后很少发生。大风发生的位置取决于飑线结构中气流的性质,气流的性质与冷池前进的程度和对流的强度关系密切。大风大部分由下沉冷气流产生,少数为近地面上升暖气流导致。大风发生的范围和强度与低层风垂直切变的强度呈正比,大范围低层风垂直切变的加强增强了飑线入流和出流的强度,是大范围大风、局部强风形成的重要原因。大风发生站次与冷池的强度和范围密切相关,冷池的加强和范围的扩大加强了后侧冷入流和前侧暖入流的强度和范围,也是大范围大风形成的重要因素。     

1999年北京持续高温天气过程的诊断分析

利用北京观象台观测资料和NCEP/NCAR再分析数据,对1999年6月24日至7月2日北京一次持续性高温天气的演变和发展过程及非绝热加热作用对系统的影响进行诊断分析,结果表明：在此次高温天气发生前,欧亚大陆中高纬度环流经向度很大,欧洲北部和贝加尔湖以南为高压脊控制,中亚和我国东北地区则处于低压槽内。贝加尔湖南部的高压脊纬向延伸范围较广,在东移过程中长时间影响北京。随着贝加尔湖以南的高压脊逐渐东移,北京上空下沉增温与非绝热加热作用有所增强,北京逐渐受到高温天气影响。在高温天气发生的后半阶段,我国东北的低压槽入海后在120130°E附近维持并发展,槽前非绝热加热率很大。从垂直方向来看,加热率在500 hPa以下随高度迅速增加,根据全型涡度方程,强烈的非绝热加热率垂直分布不均作为一个明显的涡度源区,对入海低压槽的稳定维持有显著的作用。而入海低压槽的稳定维持,又阻碍了华北高压脊的东移,使其在北京地区长时间稳定少动,为北京带来多日的持续性高温天气。     

结合投影变换和相位相关的海冰运动监测

为提高对渤海海冰旋转和平移运动的监测精度,提出一种基于投影变换的相位相关跟踪方法。选取连续8景静止水色成像仪(GOCI)图像序列,根据特征图像窗口的投影变换构造辅助函数,通过寻求函数最优解得到旋转角度集合,选择修正相关系数确定最佳旋转角度,同时根据相位信息实现海冰样本间亚像素级别的平移跟踪,消除传统相位相关法中因忽略图像相频特性所造成的匹配误差。实验结果表明,以手动测量的旋转角度为基准,该方法和传统相位相关法的旋转监测均方根误差的最大值/平均值分别为0.59/0.50与1.41/0.94,跟踪速度提高了50.6%;海冰平移运动的速度矢量与辽东湾的现场实测数据及历史资料记录数据基本一致。该方法对渤海海冰旋转和平移运动具有较好的跟踪效果。