Frequency analysis of nonstationary annual maximum flood series using the time‐varying two‐component mixture distributions

The most popular practice for analysing nonstationarity of flood series is to use a fixed single‐type probability distribution incorporated with the time‐varying moments. However, the type of probability distribution could be both complex because of distinct flood populations and time‐varying under changing environments. To allow the investigation of this complex nature, the time‐varying two‐component mixture distributions (TTMD) method is proposed in this study by considering the time variations of not only the moments of its component distributions but also the weighting coefficients. Having identified the existence of mixed flood populations based on circular statistics, the proposed TTMD was applied to model the annual maximum flood series of two stations in the Weihe River basin, with the model parameters calibrated by the meta‐heuristic maximum likelihood method. The performance of TTMD was evaluated by different diagnostic plots and indexes and compared with stationary single‐type distributions, stationary mixture distributions and time‐varying single‐type distributions. The results highlighted the advantages of TTMD with physically‐based covariates for both stations. Besides, the optimal TTMD models were considered to be capable of settling the issue of nonstationarity and capturing the mixed flood populations satisfactorily. Copyright © 2016 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.     

Groundwater recharge in an arid grassland as indicated by soil chloride profile and multiple tracers

Previous studies have shown that shallow groundwater in arid regions is often not in equilibrium with near‐surface boundary conditions due to human activities and climate change. This is especially the case where the unsaturated zone is thick and recharge rate is limited. Under this nonequilibrium condition, the unsaturated zone solute profile plays an important role in estimating recent diffuse recharge in arid environments. This paper combines evaluation of the thick unsaturated zone with the saturated zone to investigate the groundwater recharge of a grassland in the arid western Ordos Basin, NW China, using the soil chloride profiles and multiple tracers (²H, ¹⁸O, ¹³C, ¹⁴C, and water chemistry) of groundwater. Whereas conventional water balance and Darcy flux measurements usually involve large errors in recharge estimations for arid areas, chloride mass balance has been widely and generally successfully used. The results show that the present diffuse recharge beneath the grassland is 0.11–0.32 mm/year, based on the chloride mass balance of seven soil profiles. The chloride accumulation age is approximately 2,500 years at a depth of 13 m in the unsaturated zone. The average Cl content in soil moisture in the upper 13 m of the unsaturated zone ranges from 2,842 to 7,856 mg/L, whereas the shallow groundwater Cl content ranges from 95 to 351 mg/L. The corrected ¹⁴C age of shallow groundwater ranges from 4,327 to 29,708 years. Stable isotopes show that the shallow groundwater is unrelated to modern precipitation. The shallow groundwater was recharged during the cold and wet phases of the Late Pleistocene and Holocene humid phase based on palaeoclimate, and consequently, the groundwater resources are nonrenewable. Due to the limited recharge rate and thick unsaturated zone, the present shallow groundwater has not been in hydraulic equilibrium with near‐surface boundary conditions in the past 2,500 years.     

Extended finite element framework for fault rupture dynamics including bulk plasticity

We present an explicit extended finite element framework for fault rupture dynamics accommodating bulk plasticity near the fault. The technique is more robust than the standard split‐node method because it can accommodate a fault propagating freely through the interior of finite elements. To fully exploit the explicit algorithmic framework, we perform mass lumping on the enriched finite elements that preserve the kinetic energy of the rigid body and enrichment modes. We show that with this technique, the extended FE solution reproduces the standard split‐node solution, but with the added advantage that it can also accommodate randomly propagating faults. We use different elastoplastic constitutive models appropriate for geomaterials, including the Mohr–Coulomb, Drucker–Prager, modified Cam‐Clay, and a conical plasticity model with a compression cap, to capture off‐fault bulk plasticity. More specifically, the cap model adds robustness to the framework because it can accommodate various modes of deformation, including compaction, dilatation, and shearing. Copyright © 2013 John Wiley & Sons, Ltd.