Analysis of snowpack dynamics during the spring melt season for a sub‐alpine site using point measurements and numerical modeling

Snowpack dynamics through October 2014–June 2017 were described for a forested, sub‐alpine field site in southeastern Wyoming. Point measurements of wetness and density were combined with numerical modeling and continuous time series of snow depth, snow temperature, and snowpack outflow to identify 5 major classes of distinct snowpack conditions. Class (i) is characterized by no snowpack outflow and variable average snowpack temperature and density. Class (ii) is characterized by short durations of liquid water in the upper snowpack, snowpack outflow values of 0.0008–0.005 cm hr^?1, an increase in snowpack temperature, and average snow density between 0.25–0.35 g cm^?3. Class (iii) is characterized by a partially saturated wetness profile, snowpack outflow values of 0.005–0.25 cm hr^?1, snowpack temperature near 0 °C, and average snow density between 0.25–0.40 g cm^?3. Class (iv) is characterized by strong diurnal snowpack outflow pattern with values as high as 0.75 cm hr^?1, stable snowpack temperature near 0 °C, and stable average snow density between 0.35–0.45 g cm^?3. Class (v) occurs intermittently between Classes (ii)–(iv) and displays low snowpack outflow values between 0.0008–0.04 cm hr^?1, a slight decrease in temperature relative to the preceding class, and similar densities to the preceding class. Numerical modeling of snowpack properties with SNOWPACK using both the Storage Threshold scheme and Richards' equation was used to quantify the effect of snowpack capillarity on predictions of snowpack outflow and other snowpack properties. Results indicate that both simulations are able to predict snow depth, snow temperature, and snow density reasonably well with little difference between the 2 water transport schemes. Richards' equation more accurately simulates the timing of snowpack outflow over the Storage Threshold scheme, especially early in the melt season and at diurnal timescales.     

An original semi‐discrete approach to assess gas conductivity of concrete structures

For civil engineering structures with a tightness role, structural permeability is a key issue. In this context, this paper presents a new proposition of a numerical modelling of leakage rate through a cracked concrete structure undergoing mode I cracking. The mechanical state of the material, considered in the framework of continuum mechanics based on finite element modelling, is described by means of the stress‐based nonlocal damage model which takes into account the stress state and provides realistic local mechanical fields. A semi‐discrete method based on the strong discontinuity approach to estimate crack opening is then considered in the post‐treatment phase. Using a Poiseuille's like relation, the coupling between the mechanical state of the material and its dry gas conductivity is performed. For validation purposes, an original experimental campaign is conducted on a dry concrete disc loaded in a splitting setup. During the loading, gas conductivity and digital image correlation analysis are performed. The comparison with the 3D experimental mechanical global response highlights the performance of the mechanical model. The comparison between crack openings measured by digital image correlation and estimated by the strong discontinuity method shows a good agreement. Finally, the results of the semi‐discrete approach coupled with the gas conductivity compared with experimental data show a good estimation of the structural conductivity. Consequently, if the mechanical problem is well modelled at the global scale, then the proposed approach provides good estimation of gas conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Efficient DEM‐based overland flow routing using integrated recursive algorithms

Two primary concerns in performing watershed overland flow routing are the numerical instability and computational efficiency. The stability of executing an explicit scheme has to be maintained by observing the Courant–Friedrich–Lewy criterion, which is adopted to confirm that the numerical marching speed is larger than the wave celerity. Moreover, there is another criterion of time step devised in previous studies to avoid back‐and‐forth refluxing between adjacent grids. The situation of refluxing usually occurs on flat regions. In light of this, the selection of a small time increment to honor both restrictions simultaneously is believed to decrease the computational efficiency in performing overland flow routing. This study aims at creating a robust algorithm to relax both restrictions. The proposed algorithm was first implemented on a one‐dimensional overland plane to evaluate the accuracy of the numerical result by comparing it with an analytical solution. Then, the algorithm was further applied to a watershed for 2D runoff simulations. The results show that the proposed integrated algorithm can provide an accurate runoff simulation and achieve satisfactory performance in terms of computational speed.     

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.     

An analytical study on three‐dimensional versus two‐dimensional water table‐induced flow patterns in a Tóthian basin

Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.     

Changes in potential evapotranspiration and surface runoff in 1981–2010 and the driving factors in Upper Heihe River Basin in Northwest China

Changes in potential evapotranspiration and surface runoff can have profound implications for hydrological processes in arid and semiarid regions. In this study, we investigated the response of hydrological processes to climate change in Upper Heihe River Basin in Northwest China for the period from 1981 to 2010. We used agronomic, climatic and hydrological data to drive the Soil and Water Assessment Tool model for changes in potential evapotranspiration (ET₀) and surface runoff and the driving factors in the study area. The results showed that increasing autumn temperature increased snow melt, resulting in increased surface runoff, especially in September and October. The spatial distribution of annual runoff was different from that of seasonal runoff, with the highest runoff in Yeniugou River, followed by Babaohe River and then the tributaries in the northern of the basin. There was no evaporation paradox at annual and seasonal time scales, and annual ET₀ was driven mainly by wind speed. ET₀ was driven by relative humidity in spring, sunshine hour duration in autumn and both sunshine hour duration and relative humility in summer. Surface runoff was controlled by temperature in spring and winter and by precipitation in summer (flood season). Although surface runoff increased in autumn with increasing temperature, it depended on rainfall in September and on temperature in October and November. Copyright © 2016 John Wiley & Sons, Ltd.     

Review and enhancement of 3D concrete models for large‐scale numerical simulations of concrete structures

The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd.     

Introducing a new conceptual instantaneous unit hydrograph model based on a hydraulic approach

ABSTRACT

A novel approach is introduced for simulation of instantaneous unit hydrographs (IUHs). The model consists of a series of linear reservoirs that are connected to each other, and is referred to as the inter-connected linear reservoir model (ICLRM). By assuming the flow between two reservoirs is a linear function of the difference between the water levels in the reservoirs, a system of first-order linear differential equations is obtained as the ICLRM governing equation. By solving the equations, the discharge from the last reservoir is considered as an IUH. A small-scale laboratory device was constructed for the simulation of IUHs using the model. By studying four hydrographs extracted from the literature, and simulating them using both the ICLRM and the Nash model, it is concluded that the ICLRM can predict these hydrographs more accurately than the Nash model. Due to the simplicity of the construction and operation of the ICLRM and, more importantly, its visual aspect, the ICLRM may be considered as an effective educational tool for studying IUHs.     

Soil‐pile interaction in the pile vertical vibration considering true three‐dimensional wave effect of soil

The dynamic response of an end bearing pile embedded in a linear visco‐elastic soil layer with hysteretic type damping is theoretically investigated when the pile is subjected to a time‐harmonic vertical loading at the pile top. The soil is modeled as a three‐dimensional axisymmetric continuum in which both its radial and vertical displacements are taken into account. The pile is assumed to be vertical, elastic and of uniform circular cross section. By using two potential functions to decompose the displacements of the soil layer and utilizing the separation of variables technique, the dynamic equilibrium equation is uncoupled and solved. At the interface of soil‐pile system, the boundary conditions of displacement continuity and force equilibrium are invoked to derive a closed‐form solution of the vertical dynamic response of the pile in frequency domain. The corresponding inverted solutions in time domain for the velocity response of a pile subjected to a semi‐sine excitation force applied at the pile top are obtained by means of inverse Fourier transform and the convolution theorem. A comparison with two other simplified solutions has been performed to verify the more rigorous solutions presented in this paper. Using the developed solutions, a parametric study has also been conducted to investigate the influence of the major parameters of the soil‐pile system on the vertical vibration characteristics of the pile. Copyright © 2013 John Wiley & Sons, Ltd.