Granular element method for three‐dimensional discrete element calculations

This paper endows the recently‐proposed granular element method (GEM) with the ability to perform 3D discrete element calculations. By using non‐uniform rational B‐Splines to accurately represent complex grain geometries, we proposed an alternative approach to clustering‐based and polyhedra‐based discrete element methods whereby the need for complicated and ad hoc approaches to construct 3D grain geometries is entirely bypassed. We demonstrate the ability of GEM in capturing arbitrary‐shaped 3D grains with great ease, flexibility, and without excessive geometric information. Furthermore, the applicability of GEM is enhanced by its tight integration with existing non‐uniform rational B‐Splines modeling tools and ability to provide a seamless transition from binary images of real grain shapes (e.g., from 3D X‐ray CT) to modeling and discrete mechanics computations.Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Long‐term streamflow trends in Hawai'i and implications for native stream fauna

Climate change has fundamentally altered the water cycle in tropical islands, which is a critical driver of freshwater ecosystems. To examine how changes in streamflow regime have impacted habitat quality for native migratory aquatic species, we present a 50‐year (1967–2016) analysis of hydrologic records in 23 unregulated streams across the five largest Hawaiian Islands. For each stream, flow was separated into direct run‐off and baseflow and high‐ and low‐flow statistics (i.e., Q10 and Q90) with ecologically important hydrologic indices (e.g., frequency of flooding and low flow duration) derived. Using Mann–Kendall tests with a running trend analysis, we determined the persistence of streamflow trends through time. We analysed native stream fauna from ~400 sites, sampled from 1992 to 2007, to assess species richness among islands and streams. Declines in streamflow metrics indicated a general drying across the islands. In particular, significant declines in low flow conditions (baseflows), were experienced in 57% of streams, compared with a significant decline in storm flow conditions for 22% of streams. The running trend analysis indicated that many of the significant downward trends were not persistent through time but were only significant if recent decades (1987–2016) were included, with an average decline in baseflow and run‐off of 10.90% and 8.28% per decade, respectively. Streams that supported higher native species diversity were associated with moderate discharge and baseflow index, short duration of low flows, and negligible downward trends in flow. A significant decline in dry season flows (May–October) has led to an increase in the number of no‐flow days in drier areas, indicating that more streams may become intermittent, which has important implications for mauka to makai (mountain to ocean) hydrological connectivity and management of Hawai'i's native migratory freshwater fauna.     

A multidirectional semi‐analytical method for analysis of laterally loaded pile groups in multi‐layered elastic strata

A semi‐analytical method for calculating the response of single piles and pile groups subjected to lateral loading is developed in this paper. Displacements anywhere in the soil domain are tied to the displacements of the piles through decay functions. The principle of virtual work and the calculus of variations are used to derive the governing differential equations that describe the response of the piles and soil. The eigenvalue method and the finite difference technique are used to solve the system of coupled differential equations for the piles and soil, respectively. The proposed method takes into account the soil surface displacement along and perpendicular to the loading direction and produces displacement fields that are very close to those produced by the finite element method but at lower computational effort. Compared with the previous method that considered only the soil displacement along the loading direction, accounting for the multi‐directional soil displacement field produces responses for the piles and soil that are closer to those approximated by the finite element method. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling landslides in unsaturated slopes subjected to rainfall infiltration using material point method

This paper presents a dynamic fully coupled formulation for saturated and unsaturated soils that undergo large deformations based on material point method. Governing equations are applied to porous material while considering it as a continuum in which the pores of the solid skeleton are filled with water and air. The accuracy of the developed method is tested with available experimental and numerical results. The developed method has been applied to investigate the failure and post‐failure behaviour of rapid landslides in unsaturated slopes subjected to rainfall infiltration using two different bedrock geometries that lie below the top soil. The models show different failure and post‐failure mechanisms depending on the bedrock geometry and highlight the negative effects of continuous rain infiltrations. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessing uncertainty in best management practice effectiveness under future climate scenarios

Uncertainty of best management practice (BMP) performance in future climates is an important consideration for water resources managers. The objective of this study was to quantify the level of uncertainty in performance of seven agricultural BMPs due to climate change in reducing sediment, total nitrogen, and total phosphorus loads. The Soil and Water Assessment Tool coupled with mid‐21^st century climate data from the Community Climate System Model were used to develop climate change scenarios for the Tuttle Creek Lake Watershed of Kansas and Nebraska. Uncertainty level of each BMP was determined using Latin Hypercube Sampling, a constrained Monte Carlo sampling technique. Samples were taken from distributions of several variables (monthly precipitation, temperature, CO₂, and BMP implementation parameters). Cumulative distribution functions were constructed for each BMP, pollutant, and climate scenario combination. Results demonstrated that BMP performance uncertainty is amplified in the extreme climate scenario. Among BMPs, native grass replacement generally had higher uncertainty level but also had the greatest reductions. This study highlights the importance of incorporating uncertainty analysis into mitigation strategies aiming to reduce negative impacts of climate change on water resources. Copyright © 2013 John Wiley & Sons, Ltd.     

Phase equilibria modelling and LASS monazite petrochronology: P–T–t constraints on the evolution of the Priest River core complex,northern Idaho

Phase equilibria modelling, laser‐ablation split‐stream (LASS)‐ICP‐MS petrochronology and garnet trace‐element geochemistry are integrated to constrain the P–T–t history of the footwall of the Priest River metamorphic core complex, northern Idaho. Metapelitic, migmatitic gneisses of the Hauser Lake Gneiss contain the peak assemblage garnet + sillimanite + biotite ± muscovite + plagioclase + K‐feldspar ± rutile ± ilmenite + quartz. Interpreted P–T paths predict maximum pressures and peak metamorphic temperatures of ~9.6–10.3 kbar and ~785–790 °C. Monazite and xenotime ²⁰⁸Pb/²³²Th dates from porphyroblast inclusions indicate that metamorphism occurred at c. 74–54 Ma. Dates from HREE‐depleted monazite formed during prograde growth constrain peak metamorphism at c. 64 Ma near the centre of the complex, while dates from HREE‐enriched monazite constrain the timing of garnet breakdown during near‐isothermal decompression at c. 60–57 Ma. Near‐isothermal decompression to ~5.0–4.4 kbar was followed by cooling and further decompression. The youngest, HREE‐enriched monazite records leucosome crystallization at mid‐crustal levels c. 54–44 Ma. The northernmost sample records regional metamorphism during the emplacement of the Selkirk igneous complex (c. 94–81 Ma), Cretaceous–Tertiary metamorphism and limited Eocene exhumation. Similarities between the Priest River complex and other complexes of the northern North American Cordillera suggest shared regional metamorphic and exhumation histories; however, in contrast to complexes to the north, the Priest River contains less partial melt and no evidence for diapiric exhumation. Improved constraints on metamorphism, deformation, anatexis and exhumation provide greater insight into the initiation and evolution of metamorphic core complexes in the northern Cordillera, and in similar tectonic settings elsewhere.     

Multi‐temporal scale changes of streamflow and sediment load in a loess hilly watershed of China

Global climate change and diverse human activities have resulted in distinct temporal–spatial variability of watershed hydrological regimes, especially in water‐limited areas. This study presented a comprehensive investigation of streamflow and sediment load changes on multi‐temporal scales (annual, flood season, monthly and daily scales) during 1952–2011 in the Yanhe watershed, Loess Plateau. The results indicated that the decreasing trend of precipitation and increasing trend of potential evapotranspiration and aridity index were not significant. Significant decreasing trends (p < 0.01) were detected for both the annual and flood season streamflow, sediment load, sediment concentration and sediment coefficient. The runoff coefficient exhibited a significantly negative trend (p < 0.01) on the flood season scale, whereas the decreasing trend on the annual scale was not significant. The streamflow and sediment load during July–August contributed 46.7% and 86.2% to the annual total, respectively. The maximum daily streamflow and sediment load had the median occurrence date of July 31, and they accounted for 9.7% and 29.2% of the annual total, respectively. All of these monthly and daily hydrological characteristics exhibited remarkable decreasing trends (p < 0.01). However, the contribution of the maximum daily streamflow to the annual total progressively decreased (?0.07% year^?1), while that of maximum daily sediment load increased over the last 60 years (0.08% year^?1). The transfer of sloping cropland for afforestation and construction of check‐dams represented the dominant causes of streamflow and sediment load reductions, which also made the sediment grain finer. Copyright © 2015 John Wiley & Sons, Ltd.     

Sensitivity of a hydraulic model to channel erosion uncertainty during extreme flooding

Recent research into flood modelling has primarily concentrated on the simulation of inundation flow without considering the influences of channel morphology. River channels are often represented by a simplified geometry that is implicitly assumed to remain unchanged during flood simulations. However, field evidence demonstrates that significant morphological changes can occur during floods to mobilize the boundary sediments. Despite this, the effect of channel morphology on model results has been largely unexplored. To address this issue, the impact of channel cross‐section geometry and channel long‐profile variability on flood dynamics is examined using an ensemble of a 1D–2D hydraulic model (LISFLOOD‐FP) of the ~1 : 2000 year recurrence interval floods in Cockermouth, UK, within an uncertainty framework. A series of simulated scenarios of channel erosional changes were constructed on the basis of a simple velocity‐based model of critical entrainment. A Monte‐Carlo simulation framework was used to quantify the effects of this channel morphology together with variations in the channel and floodplain roughness coefficients, grain size characteristics and critical shear stress on measures of flood inundation. The results showed that the bed elevation modifications generated by the simplistic equations reflected an approximation of the observed patterns of spatial erosion that enveloped observed erosion depths. The effect of uncertainty on channel long‐profile variability only affected the local flood dynamics and did not significantly affect the friction sensitivity and flood inundation mapping. The results imply that hydraulic models generally do not need to account for within event morphodynamic changes of the type and magnitude of event modelled, as these have a negligible impact that is smaller than other uncertainties, e.g. boundary conditions. Instead, morphodynamic change needs to happen over a series of events to become large enough to change the hydrodynamics of floods in supply limited gravel‐bed rivers such as the one used in this research. Copyright © 2014 John Wiley & Sons, Ltd.     

Integrating isolated and riparian wetland modules in the PHYSITEL/HYDROTEL modelling platform: model performance and diagnosis

Mathematical modelling is a well‐accepted framework to evaluate the effects of wetlands on stream flow and watershed hydrology in general. Although the integration of wetland modules into a distributed hydrological model represents a cost‐effective way to make this assessment, the added value brought by landscape‐specific modules to a model's ability to replicate basic hydrograph characteristics remains unclear. The objectives of this paper were the following: (i) to present the adaptation of PHYSITEL (a geographic information system) to parameterize isolated and riparian wetlands; (ii) to describe the integration of specific isolated wetland and riparian wetland modules into HYDROTEL, a distributed hydrological model; and (iii) to evaluate the performance of the updated modelling platform with respect to the capacity of replicating various hydrograph characteristics. To achieve this, two sets of simulations were performed (with and without wetland modules), and the added value was assessed at three river segments of the Becancour River watershed, Quebec, Canada, using six general goodness‐of‐fit indicators and 14 water flow criteria. A sensitivity analysis of the wetland module parameters was performed to characterize their impact on stream flows of the modelled watershed. Results of this study indicate the following: (i) integration of specific wetland modules can slightly increase the capacity of HYDROTEL to replicate basic hydrograph characteristics; and (ii) the updated modelling platform allows for the explicit assessment of the impact of wetlands (e.g. typology and location) on watershed hydrology. Copyright © 2015 John Wiley & Sons, Ltd.