An alternative Boussinesq equation considering the effect of hysteresis on coastal groundwater waves

This paper presents an alternative Boussinesq equation considering hysteresis effect via a third‐order derivative term. By introducing an improved moisture–pressure retention function, this equation describes, with reasonable precision, groundwater propagation in coastal aquifers subject to Dirichlet boundary condition of different oscillation frequencies. Test results confirmed that it is necessary to consider horizontal and vertical flows in unsaturated zone, because of their variable influences on hysteresis. Hysteresis in unsaturated zone can affect the water table wave number of groundwater wave motion, such as wave damping rate and phase lag. Oscillations with different periods exert different hysteresis effect on wave propagation. Truncation/shrinkage of unsaturated zones also affects the strength of hysteresis. These impacts can be reflected in the alternative Boussinesq equation by adjusting the parameter representing the variation rate of moisture associated with pressure change, as opposed to traditional computationally expensive hysteresis algorithms. The present Boussinesq equation is simple to use and can provide feasible basis for future coupling of groundwater and surface water models. Copyright © 2016 John Wiley & Sons, Ltd.     

Hydrologic response of engineered media in living roofs and bioretention to large rainfalls: experiments and modeling

The hydrologic response of engineered media plays an important role in determining a stormwater control measure's ability to reduce runoff volume, flow rate, timing, and pollutant loads. Five engineered media, typical of living roof and bioretention stormwater control measures, were investigated in laboratory column experiments for their hydrologic responses to steady, large inflow rates. The inflow, medium water content response, and outflow were all measured. The water flow mechanism (uniform flow vs. preferential flow) was investigated by analyzing medium water content response in terms of timing, magnitude, and sequence with depth. Modeling the hydrologic process was conducted in the HYDRUS‐1D software, applying the Richards equation for uniform flow modeling, and a mobile–immobile model for preferential flow modeling. Uniform flow existed in most cases, including all initially dry living roof media with bimodal pore size distributions and one bioretention medium with unimodal pore size distribution. The Richards equation can predict the outflow hydrograph reasonably well for uniform flow conditions when medium hydraulic properties are adequately represented by appropriate functions. Preferential flow was found in two media with bimodal pore size distributions. The occurrence of preferential flow is more likely due to the interaction between the bimodal pore structure and the initial water content rather than the large inflow rate.     

Modeling nonlinear wave–wave interactions with the elliptic mild slope equation

An approach is developed to simulate wave–wave interactions using nonlinear elliptic mild-slope equation in domains where wave reflection, refraction, diffraction and breaking effects must also be considered. This involves the construction of an efficient solution procedure including effective boundary treatment, modification of the nonlinear equation to resolve convergence issues, and validation of the overall approach. For solving the second-order boundary-value problem, the Alternating Direction Implicit (ADI) scheme is employed, and the use of approximate boundary conditions is supplemented, for improved accuracy, with internal wave generation method and dissipative sponge layers. The performance of the nonlinear model is investigated for a range of practical wave conditions involving reflection, diffraction and shoaling in the presence of nonlinear wave–wave interactions. In addition, the transformation of a wave spectrum due to nonlinear shoaling and breaking, and nonlinear resonance inside a rectangular harbor are simulated. Numerical calculations are compared with the results from other relevant nonlinear models and experimental data available in literature. Results show that the approach developed here performs reasonably well, and has thus improved the applicability of this class of wave transformation models.     

U–Pb dating of silicic lavas,sills and syneruptive resedimented volcaniclastic deposits of the Lower Devonian Crudine Group,Hill End Trough,New South Wales

The Hill End Trough of central‐western New South Wales was an elongate deep marine basin that existed in the Lachlan Fold Belt from the early Late Silurian to late Early Devonian. It is represented by a regionally extensive, unfossiliferous sequence of interbedded turbidites and hemipelagites of substantially silicic volcanic derivation, which passes laterally into contemporaneous shallow‐water sedimentary rocks. The Turondale and Merrions Formations of the Lower Devonian Crudine Group are two prominent volcanogenic formations in the predominantly sedimentary trough sequence. They contain a range of primary and resedimented volcanic facies suitable for U–Pb dating. These include widespread subaqueous silicic lavas and/or lava cryptodomes, and thick sequences of crystal‐rich volcaniclastic sandstone emplaced by a succession of mass‐flows that were generated by interaction between contemporaneous subaerial pyroclastic flows and the sea. Ion microprobe dating of the two volcanogenic formations by means of the commonly used SL 13 zircon standard yields ages ranging between 411.3 ± 5.1 and 404.8 ± 4.8 Ma. Normalising the data against a different zircon standard (QGNG) yields preferred slightly older mean ages that range between 413.4 ± 6.6 and 407.1 ± 6.9 Ma. These ages broadly approximate the Early Devonian age that has been historically associated with the Crudine Group. However, the biostratigraphically inferred late Lochkovian ‐ early Emsian (mid‐Early Devonian) age for the Merrions Formation is inconsistent with the current Australian Phanerozoic Timescale, which assigns an age of 410 Ma to the Silurian‐Devonian boundary, and ages of 404.5 Ma and 395.5 Ma to the base and top of the Pragian, respectively. There is, however, good agreement if the new ages are compared with the most recently published revision of the Devonian time‐scale. This suggests that the Early Devonian stage boundaries of the Australian Phanerozoic Timescale need to be revised downward. The new ages for the Merrions Formation could also provide a time point on this time‐scale for the Pragian to early Emsian, for which no data are presently available.     

Universal hydrofracturing algorithm for shear-thinning fluids: Particle velocity based simulation

A universal particle velocity based algorithm for simulating hydraulic fractures, previously proposed for Newtonian fluids, is extended to the class of shear-thinning fluids. The scheme is not limited to any particular elasticity operator or crack propagation regime. The computations are based on two dependent variables: the crack opening and the reduced particle velocity. The application of the latter facilitates utilization of the local condition of Stefan type (speed equation) to trace the fracture front. The condition is given in a general explicit form which relates the crack propagation speed (and the crack length) to the solution tip asymptotics. The utilization of a modular structure, and the adaptive character of its basic blocks, result in a flexible numerical scheme. The computational accuracy of the proposed algorithm is validated against a number of analytical benchmark solutions.     

天山南麓雷暴日短期预报研究

利用阿克苏及邻近地区12个气象站1980—2013年雷暴资料,以及同期高空资料,统计了各站年均雷暴日数,对发生区域雷暴天气的环流形势进行分类,归纳出各型的入型指标。通过逐步回归法,建立阿克苏及邻近地区区域雷暴概率回归预报模型,并对2013年进行试预报。结果表明:(1)阿克苏及邻近地区区域雷暴的影响系统主要分为4类:巴湖低槽型、急流型、西北气流型和温度槽型。(2)对2002—2012年5—9月(共1683 d)历史资料进行判别,满足入型条件的样本数为876 d,消空率为48%;对2013年5—9月(共153 d)历史资料进行判别,入型样本数为80 d,消空率为48%。(3)对2002—2012年5—9月所有入型样本进行回代检验,平均准确率为72.0%(平均TS评分为30.1%);对2013年5—9月所有入型样本进行试预报,平均准确率为63.2%(平均TS评分为28.2%)。     

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

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

中国东部城市群发展对南海夏季风爆发影响的模拟研究

本文利用NCAR开发的CAM5.1(Community Atmosphere Model Version 5.1)模式,针对我国东部大规模城市下垫面发展对南海夏季风爆发的影响进行了数值模拟研究。结果表明我国东部大规模城市群发展可能使得南海夏季风提前1候爆发;机理分析表明:在南海夏季风爆发之前,中国东部城市群发展引起的陆面增温,使得南海及其附近地区南北温差提前逆转、中国东部区域海平面气压降低,导致中南半岛到南海地区西南气流加强,中南半岛到南海地区降水增加,而凝结潜热垂直变化强迫出的异常环流,促进了南亚高压的加强及提前北跳,相伴随的高层抽吸作用有助于季风对流的建立和西太平洋副高的减弱东撤,从而形成了有利于南海夏季风爆发的高低层环流条件,导致南海夏季风提前爆发。另外,观测结果表明1993年之后南海夏季风爆发的日期相对上一个年代明显提前约2候,城市化快速发展阶段与南海夏季风爆发的年代际变化存在时间段的吻合,表明城市下垫面发展可能是南海夏季风提前爆发的原因之一。     

利用压密方程计算分析西安市地面沉降特征

利用西安市地下水位监测资料,基于水文地质三维结构模型,在不考虑粘土层滞后压缩变形的理想情况下,采用压密方程计算获得西安市抽取地下水可引起的理论地面沉降量及沉降分布特征,结果表明:计算所得沉降显著区位于西安市西南部的高新区及东南部的曲江新区,沉降量呈由北向南递增的特征,与InSAR监测结果整体趋势具有较好的一致性;地面沉降分布特征受到活动地裂缝影响,沉降曲线呈近NE向偏转展布。研究结果可为沉降灾害预防研究提供宏观的参考信息。     

一种裸露土壤湿度反演方法

针对目前土壤湿度反演方法研究较少且缺少实时性的现状,该文提出一种土壤湿度反演方法——最小二乘支持向量机技术。以积分方程模型为正向算法,数值模拟不同雷达参数(频率、入射角及极化)下后向散射系数随土壤含水量和地表粗糙度的变化情况。经过数据敏感性分析,选取C-波段和X-波段、小入射角下的同极化后向散射系数作为支持向量回归的训练样本信息;经过适当的训练,利用支持向量回归技术对土壤含水量进行了反演研究;并考虑通过多频率、多极化、多入射角数据的组合,消除地表粗糙度的影响,提高反演精度。模拟结果表明,该方法反演土壤湿度具有较高的精度和较好的实时性;同时,与人工神经网络方法的结果比较,证明了该方法的有效性,为土壤湿度的反演研究提供了一种方法。