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.     

Spatial and temporal variability of stable isotopes (δ18O and δ2H) in surface waters of arid,mountainous Central Asia

Characterization of spatial and temporal variability of stable isotopes (δ¹⁸O and δ²H) of surface waters is essential to interpret hydrological processes and establish modern isotope–elevation gradients across mountainous terrains. Here, we present stable isotope data for river waters across Kyrgyzstan. River water isotopes exhibit substantial spatial heterogeneity among different watersheds in Kyrgyzstan. Higher river water isotope values were found mainly in the Issyk‐Kul Lake watershed, whereas waters in the Son‐Kul Lake watershed display lower values. Results show a close δ¹⁸O–δ²H relation between river water and the local meteoric water line, implying that river water experiences little evaporative enrichment. River water from the high‐elevation regions (e.g., Naryn and Son‐Kul Lake watershed) had the most negative isotope values, implying that river water is dominated by snowmelt. Higher deuterium excess (average d = 13.9‰) in river water probably represents the isotopic signature of combined contributions from direct precipitation and glacier melt in stream discharge across Kyrgyzstan. A significant relationship between river water δ¹⁸O and elevation was observed with a vertical lapse rate of 0.13‰/100 m. These findings provide crucial information about hydrological processes across Kyrgyzstan and contribute to a better understanding of the paleoclimate/elevation reconstruction of this region.     

Displacement analytical prediction of shallow tunnel based on unified displacement function under slope boundary

For slope condition of ground surface, the asymmetrical deformation about the vertical center line and the horizontal center line of the tunnel cross section can be formed. A unified displacement function expressed by the Fourier series is presented to express the asymmetrical deformation of the tunnel cross section. Five basic deformation modes corresponding to the expansion order 2 are a complete deformation mode to reflect deformation behaviors of the tunnel cross section under slope boundary. Such this complete displacement mode is implemented into the complex variable solution for analytically predicting tunneling-induced ground deformation under slope boundary. All of these analytical solutions are verified by good agreements of the comparison between the analytical solutions and finite element method results. A parameter study is carried out to investigate the influence of deformation modes of the tunnel cross section, geometrical conditions of the tunnel and the slope angle, and “Buoyancy effect” on the displacement field. Finally, the proposed method is consistent with measured data of the Hejie tunnel in China qualitatively. The presented solution can provide a simplified indication for evaluating the ground deformation under slope condition of ground surface.     

Temporal variation of soil erosion resistance on sloping farmland during the growth stages of maize (Zea mays L.)

Maize growth has great effects on soil properties and thus likely induces the changes in soil erosion resistance on sloping farmland. However, temporal variation of soil erosion resistance during the growth stages of maize is still unclear in the mountainous yellow soil area where maize is the dominant crop. In this study, four maize plots (MP) and four bare land plots (CK) were conducted to investigate soil erosion resistance, and multiple indicators of soil erosion resistance were measured including the total soil anti-scourability (TAS), mean weight diameter (MWD), soil erodibility K factor and soil shear strength (SH). A comprehensive soil erosion resistance index (CSERI) was employed to quantify the temporal variation of soil erosion resistance during the growth stages of maize (seedling stage, SS; jointing stage, JS; tasselling stage, TS; maturing stage, MS). The results showed that TAS, MWD, SH increased significantly with maize growth and SH decreased when at MS. But K factor decreased significantly over time. CSERI increased significantly during the growth stages of maize and the CSERI of JS, TS, MS increased on average by 74.72, 180.68 and 234.57% than that of SS. Compared to CK, CSERI of MP increased by 49.90, 66.82, 55.60 and 38.61% during the growth stages of maize. The temporal variation of soil erosion resistance was closely related to the changes in maize cover, maize roots and soil organic carbon. The findings demonstrated that it is necessary to consider the temporal variation of soil erosion resistance in the mountainous yellow soil area.     

Simulation of quasi-linear mesoscale convective systems in northern China: Lightning activities and storm structure

Two intense quasi-linear mesoscale convective systems(QLMCSs) in northern China were simulated using the WRF(Weather Research and Forecasting) model and the 3D-Var(three-dimensional variational) analysis system of the ARPS(Advanced Regional Prediction System) model.A new method in which the lightning density is calculated using both the precipitation and non-precipitation ice mass was developed to reveal the relationship between the lightning activities and QLMCS structures.Results indicate that,compared with calculating the results using two previous methods,the lightning density calculated using the new method presented in this study is in better accordance with observations.Based on the calculated lightning densities using the new method,it was found that most lightning activity was initiated on the right side and at the front of the QLMCSs,where the surface wind field converged intensely.The CAPE was much stronger ahead of the southeastward progressing QLMCS than to the back it,and their lightning events mainly occurred in regions with a large gradient of CAPE.Comparisons between lightning and non-lightning regions indicated that lightning regions featured more intense ascending motion than non-lightning regions;the vertical ranges of maximum reflectivity between lightning and non-lightning regions were very different;and the ice mixing ratio featured no significant differences between the lightning and non-lightning regions.     

高性能计算技术在气象领域的应用

高性能计算通过应用超级计算机与并行处理技术解决复杂的计算问题,是信息技术发展比较迅猛的领域之一。气象应用始终是高性能计算的重要领域,高性能计算技术有效地解决了高分辨率、高精度气象数值预报模式发展限制,在气象预报预测业务中发挥着核心支撑作用。数十年以来,由于数值天气模式研究和业务运行对计算资源的强烈需求,国内外气象领域高性能计算机及应用迅速发展起来。气象领域对高性能计算能力及系统的可靠性需求日益提升。高性能计算技术将与气象预报应用日益融合,相互影响促进,不断创新发展。为满足气象预报预测业务和科研工作需求,中国气象部门将进一步提升高性能计算能力,并致力于优化集约高性能计算系统布局,高效管理计算资源,发挥最大效益。