Some aspects of head‐variance evaluation

We compare two methods of evaluating head covariance for two‐dimensional steady‐state flow in mildly heterogeneous bounded rectangular aquifers. The quasi‐analytical approach, widely used in stochastic subsurface hydrology, is based on the Green's function representation, and involves numerical four‐fold integration. We compare this approach with a numerical solution of the two‐dimensional boundary‐value problem for head covariance. We show that the finite differences integration of this problem is computationally less expensive than numerical four‐fold integration of slowly‐convergent infinite series. This revised version was published online in July 2006 with corrections to the Cover Date.     

一种近岸区波浪破碎模型

从波浪破碎的能量关系入手,以紊流能量方程为基础,考虑破碎区内单个波在不同破碎阶段所提供的紊动能量强度的变化过程,提出了一种波浪破碎模式.通过将这一模型引入Boussinesq方程中,初步建立了一种近岸区波浪变形数学模型,并用波浪水槽实验资料对模型模拟波高和平均水位的情况进行了初步验证,得到了良好的结果.     

An approach to analysing plot scale infiltration and runoff responses to rainfall of fluctuating intensity

Simulated rainfall of fluctuating intensity was applied to runoff plots on bare dryland soils in order to explore a new method for analysing the non‐steady‐state responses of infiltration and overland flow. The rainfall events all averaged 10 mm/h but included intensity bursts of up to 70 mm/h and lasting 5–15 min, as well as periods of low intensity and intermittency of up to 25 min. Results were compared with traditional steady‐state estimates of infiltrability made under simulated rainfall sustained at a fixed intensity of 10 mm/h. Mean event infiltration rate averaged 13.6% higher under fluctuating intensities, while runoff ratios averaged only 63% of those seen under constant intensity. In order to understand the changing soil infiltrability, up to three affine Horton infiltration equations were fitted to segments of each experiment. All equations had the same final infiltrability f_c, but adjusted values for coefficients f₀ (initial infiltrability) and K_f (exponential decay constant) were fitted for periods of rainfall that followed significant hiatuses in rainfall, during which subsurface redistribution allowed near‐surface soil suction to recover. According to the fitted Horton equations, soil infiltrability recovered by up 10–24 mm/h during intra‐event rainfall hiatuses of 15 to 20‐min duration, contributing to higher overall event infiltration rates and to reduced runoff ratios. The recovery of infiltrability also reduced the size of runoff peaks following periods of low intensity rainfall, compared with the predictions based on single Horton infiltration equations, and in some cases, no runoff at all was recorded from late intensity peaks. The principal finding of this study is that, using a set of affine equations, the intra‐event time variation of soil infiltrability can be tracked through multiple intensity bursts and hiatuses, despite the lack of steady‐state conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

沉积物粒度组分空间预测的神经网络残余kriging方法

针对近海表层沉积物粒度组分空间变异的尺度差异性,提出了基于广义回归神经网络残余kriging的沉积物粒度组分空间预测方法,并以海州湾沉积物粒度数据为基础,分析了其在沉积物粒度组分空间预测和底质类型制图中的应用效果。结果表明,广义回归神经网络残余kriging方法能够获得比普通kriging方法更高的沉积物粒度组分空间预测精度,并且其底质类型的总体空间预测精度达到85%以上,相应的Kappa系数也在0.8以上,显示底质制图的预测类型与样本的实测类型具有较高的一致性。新方法对于开展定量化的沉积物粒度组分空间预测和底质类型制图具有参考价值。     

Evaluation of ASCE 7 equations for designing acceleration‐sensitive nonstructural components using data from instrumented buildings

This study uses instrumented buildings and models of code‐based designed buildings to validate the results of previous studies that highlighted the need to revise the ASCE 7 F_p equation for designing nonstructural components (NSCs) through utilizing oversimplified linear and nonlinear models. The evaluation of floor response spectra of a large number of instrumented buildings illustrates that, unlike the ASCE 7 approach, the in‐structure and the component amplification factors are a function of the ratio of NSC period to the supporting building modal periods, the ground motion intensity, and the NSC location. It is also shown that the recorded ground motions at the base of instrumented buildings in most cases are significantly lower than design earthquake (DE) ground motions. Because ASCE 7 is meant to provide demands at a DE level, for a more reliable evaluation of the F_p equation, 2 representative archetype buildings are designed based on the ASCE 7‐16 seismic provisions and exposed to various ground motion intensity levels (including those consistent with the ones experienced by instrumented buildings and the DE). Simulation results of the archetype buildings, consistent with previous numerical studies, illustrate the tendency of the ASCE 7 in‐structure amplification factor, [1 + 2(z/h)] , to significantly overestimate demands at all floor levels and the ASCE 7 limit of to in many cases underestimate the calculated NSC amplification factors. Furthermore, the product of these 2 amplification factors (that represents the normalized peak NSC acceleration) in some cases exceeds the ASCE 7 equation by a factor up to 1.50.     

透空式防波堤周围的非线性波浪传播的数值模拟

对非线形波浪在透空式防波堤周围的波浪变形进行了数值模拟,在Boussinesq波浪方程中加入与透空建筑物有关的新的耗散项,从而界定了透空建筑物引起的部分反射和透射,波浪折射衍射的传播过程通过控制方程求解。波浪控制方程通过有限差分方法求解。模型应用于模拟波浪经过具有部分反射的群桩式透空结构,结果表明透空式防波堤可以有效地衰减波浪,是重力式结构的一种替代形式。     

利用控制点三维信息标定机载双天线干涉SAR参数

采用基于敏感度方程的方法,研究基于三维重建模型下的机载双天线干涉SAR系统的干涉参数定标问题.干涉参数定标是生成高精度数字高程模型的关键.本文修正了Madsen提出的干涉SAR三维重建的视向量正交分解算法,采用电磁波波前的球面波模型,加入了载机的姿态旋转,构建一种新的干涉SAR三维重建模型.利用各干涉参数对控制点三维信息的不同的敏感性,提出分别利用地面控制点三维信息,对各干涉参数进行定标.并利用中国科学院电子学研究所自主设计、研制的机载干涉SAR系统数据,进行定标处理实验验证.     

Lagged rejuvenation of groundwater indicates internal flow structures and hydrological connectivity

Large proportions of rainwater and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time-varying processes. Chlorofluorocarbon (CFC)-based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km² forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from nine different wells with depths of 2–18 m close to the stream network. Immediately below the water table, CFC-based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age-depth-relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity.     

FluBiDi – a model to estimate flood based on runoff: validation using extreme and natural basin conditions

FluBiDi is a two-dimensional model created to simulate real events that can take days and months, as well as short events (minutes or hours) and inclusive laboratory tests. To verify the robustness of FluBiDi, it was tested using a previous study with both designed and real digital elevation models. The results highlight good agreement between the models (i.e. Mike Flood, SOBEK, ISIS 2D, and others) tested and FluBiDi (around 90% for a specific instant and 95% for the complete time simulation). In the simulated hydrographs, the discharge peak value, time to peak, and water level results were accurate, reproducing them with an error of less than 5%. The velocity differences observed in a couple of tests in FluBiDi were associated with very short periods of time (seconds). However, FluBiDi is highly accurate for simulating floods under real topographical conditions with differences of around 2 cm when water depth is around 150 cm. The average water depth and velocities are precise, and the model describes with high accuracy the pattern and extent of floods. FluBiDi has the capability to be adjusted to different types of events and only requires limited input data.     

Numerical modeling of lock-exchange gravity/turbidity currents by a high-order upwinding combined compact difference scheme

This study presents two-dimensional direct numerical simulations for sediment-laden current with higher density propagating forward through a lighter ambient water.The incompressible NavierStokes equations including the buoyancy force for the density difference between the light and heavy fluids are solved by a finite difference scheme based on a structured mesh.The concentration transport equations are used to explore such rich transport phenomena as gravity and turbidity currents.Within the framework of an Upwinding Combined Compact finite Difference(UCCD)scheme,rigorous determination of weighting coefficients underlies the modified equation analysis and the minimization of the numerical modified wavenumber.This sixth-order UCCD scheme is implemented in a four-point grid stencil to approximate advection and diffusion terms in the concentration transport equations and the first-order derivative terms in the Navier-Stokes equations,which can greatly enhance convective stability and increase dispersive accuracy at the same time.The initial discontinuous concentration field is smoothed by solving a newly proposed Heaviside function to prevent numerical instabilities and unreasonable concentration values.A two-step projection method is then applied to obtain the velocity field.The numerical algorithm shows a satisfying ability to capture the generation,development,and dissipation of the Kelvin-Helmholz instabilities and turbulent billows at the interface between the current and the ambient fluid.The simulation results also are compared with the data in published literatures and good agreements are found to prove that the present numerical model can well reproduce the propagation,particle deposition,and mixing processes of lock-exchange gravity and turbidity currents.