A probabilistic approach for computing water retention of particulate systems from statistics of grain size and tessellated pore network

The paper offers an analytical determination of the hydraulic properties of an unsaturated soil with reference to its retention curve, which describes the relationship between the volumetric water content and capillarity through matric suction. The analysis combines a particulate approach focused on the physics at the pore scale, including microstructural aspects, with a probabilistic approach where the void space and grain size are considered as random variables. In the end, the soil water characteristic curve of an unsaturated granular medium along a drying path can be derived analytically based on the sole information of particle size distribution. The analysis hinges on the tessellation of a wet granular system into an assemblage of tetrahedral unit cells revealing a pore network upon which capillary physics are computed with respect to pore throat invasion by a non-wetting fluid with evolving pendular capillary bridges. The crux of the paper is to pass from particle size probability distribution to a matching void space distribution to eventually reveal key information such as void cell and solid volume statistics. Making reasonable statistically based assumptions to render calculations tractable, the water retention curve can be readily constructed. Model predictions compare quite favourably with experimental data available for actual soils, especially in the high saturation range. Having a sound scientific basis, the model can be made amenable to address a variety of soils with a wider range of particle sizes.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

Suspended sediment and total dissolved solid yield patterns at the headwaters of Urumqi River,northwestern China: a comparison between glacial and non‐glacial catchments

In order to understand the differences in the suspended sediment and total dissolved solid (TDS) yield patterns between the glacial and non‐glacial catchments at the headwaters of Urumqi River, northwestern China, water samples were collected from a glacier catchment and an empty cirque catchment within the region, during three melting seasons from 2006 to 2008. These samples were analyzed to estimate suspended sediment and TDS concentrations, fluxes and erosion rates in the two adjoining catchments. There were remarked differences in suspended sediment and TDS yield patterns between the two catchments. Suspended sediment concentrations were controlled mainly by the sediment source, whereas TDS concentrations were primarily related to the hydrologic interaction with soil minerals. Generally, the glacial catchment had much higher suspended sediment and TDS yields, together with higher denudation rates, than the non‐glacial catchment. Overall, glacial catchment was mainly dominated by physical denudation process, whereas the non‐glacial catchment was jointly influenced by physical and chemical denudation processes. The observed differences in material delivery patterns were mainly controlled by the runoff source and the glacial processes. The melting periods of glacier and snow were typically the most important time for the suspended sediment and TDS yields. Meanwhile, episodic precipitation events could generate disproportionately large yields. Subglacial hydrology dynamics, glaciers pluck and grind processes could affect erodibility, and the large quantities of dust stored on the glacier surface provided additional sources for suspended sediment transport in the glacial catchment. These mechanisms imply that, in response to climate change, the catchment behaviour will be modified significantly in this region, in terms of material flux. Copyright © 2013 John Wiley & Sons, Ltd.     

The response in floodplain respiration of an alpine river to experimental inundation under different temperature regimes

The respiratory potential [i.e. electron transport system activity (ETSA)] of soils and sediments from five floodplain habitats (channel, gravel, islands, riparian forest and grassland) of the Urbach River, Switzerland, and actual respiration rate (R) of the same samples exposed to experimental inundation were measured. Measurements were carried out at three incubation temperatures (4°C, 12°C and 20°C), and ETSA/R ratios (i.e. exploitation of the overall metabolic capacity) were investigated to better understand the effects of temperature and inundation on floodplain functional heterogeneity. Furthermore, ETSA/R ratios obtained during experimental inundation were compared with ETSA/R ratios from field measurements to investigate the exploitation in total metabolic potential at different conditions. Lowest ETSA and R were measured in samples from channel and gravel habitats, followed by those from islands. Substantially higher values were measured in soils from riparian forest and grassland. Both ETSA and R increased with increasing temperature in samples from all habitats, while the ETSA/R ratio decreased because of a rapid response in microbial community respiration to higher temperatures. The metabolic capacity exploitation (i.e. ETSA/R) during experimental inundation was lowest in predominantly terrestrial samples (riparian forest and grassland), indicating the weakest response to wetted conditions. Comparison of experimentally inundated and field conditions revealed that in rarely flooded soils, the metabolic capacity was less exploited during inundation than during non‐flooded conditions. The results suggest high sensitivity in floodplain respiration to changes in temperature and hydrological regime. ETSA/R ratios are considered good indicators of changes in metabolic activity of floodplain soils and sediments, and thus useful to estimate the impact of changes in hydrological regime or to evaluate success of floodplain restoration actions. Copyright © 2015 John Wiley & Sons, Ltd.     

Submarine groundwater discharge revealed by aerial thermal infrared imagery: a case study on Jeju Island,Korea

Submarine groundwater discharge (SGD) is a global phenomenon that carries large volumes of groundwater and dissolved chemical species such as nutrient, metals, and organic compounds to coastal zones. We report the influence of SGD on the coastal waters of Jeju Island, Korea, using high‐resolution aerial thermal infrared (TIR) mapping techniques and field investigations. An aircraft‐based system was implemented using a cost‐effective TIR camera for aerial TIR mapping. Ground‐based calibrations and system integration with GPS/IMU (global positioning system/inertial measurement unit) were performed for the aerial systems. The aerial surveys showed distinct low‐temperature signatures of SGD along the coasts of Jeju Island, revealing large groundwater inputs from the coastal aquifers to the ocean. Multiple aerial surveys over a range of seasons and tidal stages revealed that SGD rates dynamically affect the sea surface temperature (SST) of the coastal zone. The in‐situ measurements supported that SGD has a substantial influence on the coastal water chemistry as well as SST. Our observations highlight the extent to which aerial‐based TIR mapping can serve as a powerful tool for studying SGD and other coastal processes. Copyright © 2016 John Wiley & Sons, Ltd.     

Trends and multi‐annual variability of water temperatures in the river Danube,Serbia

The relationship between air (Ta) and water temperature (Tw) is very important because it shows how the temperature of a water body might respond to future changes in surface Ta. Mean monthly Tw records of three gauging stations (Bezdan, Bogojevo i Veliko Gradi?te) were analysed alongside mean monthly discharge (Q) for the same stations. Additionally, Ta series from two meteorological stations (Sombor and Veliko Gradi?te) were correlated with Tw variations over the period 1950–2012. Locally weighted scatter point smoothing (LOWESS) was used to investigate long‐term trends in the raw data, alongside the Mann–Kendall (MK) trend test. Trend significance was established using Yue–Pilon's pre‐whitening approaches to determine trends in climate data. Also, the rescaled adjusted partial sums (RAPS) method was used to detect dates of possible changes in the time series. Statistically significant warming trends were observed for annual and seasonal minimum and maximum Tw at all investigated sites. The strongest warming was observed at Bogojevo gauging station for seasonal maximum Tw, with +0.05 °C per year on average. RAPS established that the trend began in the 1980s. This behaviour is linked to climate patterns in the North and East Atlantic which determine the amount of heat advected onto mainland Europe. Statistically significant correlations were found for all Tw on an annual basis. Overall, the strongest correlations (p < 0.01) between Tw residuals and the North Atlantic Oscillation (NAO) were recorded for the winter period. These findings suggest possible predictability of Tw over seasonal time‐scales. Copyright © 2016 John Wiley & Sons, Ltd.     

Dissolved inorganic carbon evolution in neutral discharge from mine tailings piles

We measured the concentrations of dissolved inorganic carbon (DIC) and major ions and the stable carbon isotope ratios of DIC (δ¹³C_DIC) in two creeks discharging from carbonate‐rich sulphide‐containing mine tailings piles. Our aim was to assess downstream carbon evolution of the tailings discharge as it interacted with the atmosphere. The discharge had pH of 6.5–8.1 and was saturated with respect to carbonates. Over the reach of one creek, the DIC concentrations decreased by 1.1 mmol C/l and δ¹³C_DIC increased by ~4.0‰ 200 m from the seep source. The decrease in the DIC concentrations was concomitant with decreases in the partial pressure of CO_2(aq) because of the loss of excess CO_2(aq) from the discharge. The corresponding enrichment in the δ¹³C_DIC is because of kinetic isotope fractionation accompanying the loss of CO_2(g). Over the reach of the other creek, there was no significant decrease in the DIC concentrations or notable changes in the δ¹³C_DIC. The insignificant change in the DIC concentrations and the δ¹³C_DIC is because the first water sample was collected 160 m away from the discharge seep, not accessible during this research. In this case, most of the excess CO_2(aq) was lost before our first sampling station. Our results indicate that neutral discharges from tailings piles quickly lose excess CO_2(aq) to the atmosphere and the DIC becomes enrich in ¹³C. We suggest that a significant amount of carbon cycling in neutral discharges from tailings piles occur close to the locations where the discharge seeps to the surface. Copyright © 2015 John Wiley & Sons, Ltd.     

Root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate

The root‐zone moisture replenishment mechanisms are key unknowns required to understand soil hydrological processes and water sources used by plants. Temporal patterns of root‐zone moisture replenishment reflect wetting events that contribute to plant growth and survival and to catchment water yield. In this study, stable oxygen and hydrogen isotopes of twigs and throughfall were continuously monitored to characterize the seasonal variations of the root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate in South Australia. The two studied hillslopes (the north‐facing slope [NFS] and the south‐facing slope [SFS]) had different environmental conditions with opposite aspects. The twig and throughfall samples were collected every ~20 days over 1 year on both hillslopes. The root‐zone moisture replenishment, defined as percentage of newly replenished root‐zone moisture as a complement to antecedent moisture for plant use, calculated by an isotope balance model, was about zero (±25% for the NFS and ± 15% for the SFS) at the end of the wet season (October), increased to almost 100% (±26% for the NFS and ± 29% for the SFS) after the dry season (April and May), then decreased close to zero (±24% for the NFS and ± 28% for the SFS) in the middle of the following wet season (August). This seasonal pattern of root‐zone moisture replenishment suggests that the very first rainfall events of the wet season were significant for soil moisture replenishment and supported the plants over wet and subsequent dry seasons, and that NFS completed replenishment over a longer time than SFS in the wet season and depleted the root zone moisture quicker in the dry season. The stable oxygen isotope composition of the intraevent samples and twigs further confirms that rain water in the late wet season contributed little to root‐zone moisture. This study highlights the significant role of the very first rain events in the early wet season for ecosystem and provides insights to understanding ecohydrological separation, catchment water yield, and vegetation response to climate changes.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.