Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

Pixel‐oriented land use classification in energy balance modelling

Mass and energy transfer between soil, vegetation and atmosphere is the process that allows to maintain an adequate energy and water balance in the earth–atmosphere system. However, the evaluation of the energy balance components, such as the net radiation and the sensible and latent heat fluxes, is characterized by significant uncertainties related to both the dynamic nature of heat transfer processes and surfaces heterogeneity. Therefore, a detailed land use classification and an accurate evaluation of vegetation spatial distribution are required for an accurate estimation of these variables. For this purpose, in the present article, a pixel‐oriented supervised classification was applied to obtain land use maps of the Basilicata region in Southern Italy by processing three Landsat TM and ETM+ satellite images. An accuracy analysis based on the overall accuracy index and the agreement Khat of Cohen coefficient showed a good performance of the applied classification methodology and a good quality of the obtained maps. Subsequently, these maps were used in the application of a simplified two‐source energy balance model for estimating the actual evapotranspiration at a regional scale. The comparison between the simulations made by applying the simplified two‐source energy balance model and the measurements of evapotranspiration at a lysimetric station located in the study area showed the applicability and the validity of the proposed methodology. Copyright © 2012 John Wiley & Sons, Ltd.     

Deducing the spatial variability of exchange within a longitudinal channel water balance

Developing an appropriate data collection scheme to infer stream–subsurface interactions is not trivial due to the spatial and temporal variability of exchange flowpaths. Within the context of a case study, this paper presents the results from a number of common data collection techniques ranging from point to reach scales used in combination to better understand the spatial complexity of subsurface exchanges, infer the hydrologic conditions where individual influences of hyporheic and groundwater exchange components on stream water can be characterized, and determine where gaps in information arise. We start with a tracer‐based, longitudinal channel water balance to quantify hydrologic gains and losses at a sub‐reach scale nested within two consecutive reaches. Next, we look at groundwater and stream water surface levels, shallow streambed vertical head gradients, streambed and aquifer hydraulic conductivities, water chemistry, and vertical flux rates estimated from streambed temperatures to provide more spatially explicit information. As a result, a clearer spatial understanding of gains and losses was provided, but some limitations in interpreting results were identified even when combining information collected over various scales. Due to spatial variability of exchanges and areas of mixing, each technique frequently captured a combination of groundwater and hyporheic exchange components. Ultimately, this study provides information regarding technique selection, emphasizes that care must be taken when interpreting results, and identifies the need to apply or develop more advanced methods for understanding subsurface exchanges. Copyright © 2013 John Wiley & Sons, Ltd.     

Stability analysis of tunnels driven in jointed rocks by means of a homogenized limit analysis approach

Conceived as a potential alternative to the classical design methods employed for analyzing the stability of underground structures driven in jointed rocks, the homogenization approach stems from the heuristic idea that, from a macroscopic point of view, a rock mass cut by a network of joints may be perceived as a homogenized continuum. The strength properties of the latter can be theoretically obtained from the failure conditions of its individual constituents: rock matrix and joint interfaces. At the material level, the limit analysis reasoning is used in the context of homogenization to formulate the homogenized strength criterion of a jointed rock mass in the particular situation of a single set of parallel joints. As it could be expected, the obtained closed‐form expressions show the strength anisotropy induced by joint preferential orientation. The support functions (π functions) associated with the homogenized strength criterion are also determined in both plane strain and three‐dimensional cases. This criterion is then applied to the investigation of stability analysis of a tunnel excavated in a jointed rock mass. Upper bounds estimated of the stability factor are derived from the implementation of the kinematic approach directly on the homogenized underground structure. Finally, the approach is applied to analyze and discuss the collapse of the Pinheiros subway station (São Paulo, Brazil). Copyright © 2014 John Wiley & Sons, Ltd.     

Continuous measurement of whole‐tree water balance for studying urban tree transpiration

Street and garden trees in urban areas are often exposed to advection of strong vapour pressure deficit (VPD) air that can raise the whole‐tree transpiration rate (E_T), known as the oasis effect. However, urban trees tend to have small soil volume compared with natural conditions, and so they are believed to strongly regulate stomata. E_T characteristics of such urban trees have not been well understood because of a lack of reliable measurement methods. Therefore, we propose a novel weighing lysimeter method and investigate the whole‐tree water balance of an isolated container‐grown Zelkova serrata to examine (a) which biotic and abiotic factors determine E_T and (b) which spatial and temporal information is needed to predict E_T under urban conditions. Whole‐tree water balance and environmental conditions were measured from 2010 to 2012. Although leaf area substantially increased in the study period, daily E_T did not vary much. E_T increased with VPD almost linearly in 2010 but showed saturation in 2011 and 2012. Root water uptake lagged E_T by 40 min in 2012. These results suggest that the small planter box interfered with root growth and that hydraulic supply capacities did not increase sufficiently to support leaf area increase. From analysis of water balance, we believe that neglecting soil drought effects on street trees without irrigation in Japan will overestimate E_T over 4–5 sunny days at the longest. This is unlike previous studies of forest.     

Understory and small trees contribute importantly to stemflow of a lower montane cloud forest

Stemflow (Sf) measurements in tropical rain and montane forests dominated by large trees rarely include the understory and small trees. In this study, contributions of lower (1‐ to 2‐m height) and upper (>2‐m height and <5‐cm diameter at breast height [DBH]) woody understory, small trees (5 < DBH < 10 cm), and canopy trees (>10‐cm DBH) to Sf per unit ground area (Sf_a) of a Mexican lower montane cloud forest were quantified for 32 days with rainfall (P) during the 2014 wet season. Rainfall, stemflow yield (Sf_y), vegetation height, density, and basal area were measured. Subsequently, stemflow funneling ratios (SFRs) were calculated, and three common methods to scale up Sf_y from individual trees to the stand level (tree‐Sf_y correlation, P‐Sf_y correlation, and mean‐Sf_y extrapolation) were used to calculate Sf_a. Understory woody plants, small trees, and upper canopy trees represented 96%, 2%, and 2%, respectively, of the total density. Upper canopy trees had the lowest SFRs (1.6 ± 0.5 Standard Error (SE) on average), although the lower understory had the highest (36.1 ± 6.4). Small trees and upper understory presented similar SFRs (22.9 ± 5.4 and 20.2 ± 3.9, respectively). Different Sf scaling methods generally yielded similar results. Overall Sf_a during the study period was 22.7 mm (4.5% of rainfall), to which the understory contributed 70.1% (15.9 mm), small trees 10.6% (2.4 mm), and upper canopy trees 19.3% (4.4 mm). Our results strongly suggest that for humid tropical forests with dense understory of woody plants and small trees, Sf of these groups should be measured to avoid an underestimation of overall Sf at the stand level.     

Characterisation of hydroclimatological trends and variability in the Lake Naivasha basin,Kenya

Recent hydro‐climatological trends and variability characteristics were investigated for the Lake Naivasha basin with the aim of understanding the changes in water balance components and their evolution over the past 50 years. Using a Bayesian change point analysis and modified Mann–Kendall tests, time series of annual mean, maximum, minimum, and seasonal precipitation and flow, as well as annual mean lake volumes, were analysed for the period 1960–2010 to uncover possible abrupt shifts and gradual trends. Double cumulative curve analysis was used to investigate the changes in hydrological response attributable to either human influence or climatic variability. The results indicate a significant decline in lake volumes at a mean rate of 9.35 × 10⁶ m³ year^?1. Most of the river gauging stations showed no evidence of trends in the annual mean and maximum flows as well as seasonal flows. Annual minimum flows, however, showed abrupt shifts and significant (upward/downward) trends at the main outlet stations. Precipitation in the basin showed no evidence of abrupt shifts, but a few stations showed gradual decline. The observed changes in precipitation could not explain the decline in both minimum flows and lake volumes. The findings show no evidence of any impact of climate change for the Lake Naivasha basin over the past 50 years. This implies that other factors, such as changes in land cover and infrastructure development, have been responsible for the observed changes in streamflow and lake volumes. Copyright © 2015 John Wiley & Sons, Ltd.     

Unraveling complex hydrogeological processes in Andean basins in south‐central Chile: An integrated assessment to understand hydrological dissimilarity

Groundwater storage, drainage, and interbasin water exchange are common hydrological processes but often difficult to quantify due to a lack of local observations. We present a study of three volcanic mountainous watersheds located in south‐central Chile (~36.9 ° S) in the Chillán volcanic complex (Chillán, Renegado, and Diguillín river basins). These are neighboring basins that are similar with respect to the metrics normally available for characterization everywhere (e.g., precipitation, temperature, and land cover). In a hydrological sense, similar (proportional) behavior would be expected if these catchments would be characterized with this general information. However, these watersheds show dissimilar behavior when analyzed in detail. The surface water balance does not fit for any of these watersheds individually; however, the water balance of the whole system can be explained by likely interbasin water exchanges. The Renegado river basin has an average annual runoff per unit of area on the order of 60–65% less than those of the Diguillín and Chillán rivers, which is contradictory to the hydrological similarity among the basins. To understand the main processes that control streamflow generation, two analyses were performed: (a) basin metrics (land cover, geologic, topographic, and climatological maps) and hydro‐meteorological data analyses and (b) a water balance model approach. The analyses contribute to a plausible explanation for the hydrogeological processes in the system. The soils, topography, and geology of the Chillán–Renegado–Diguillín system favor the infiltration and groundwater movements from the Renegado river basin, mainly to the neighboring Diguillín basin. The interbasin water exchanges affect hydrological similarity and explain the differences observed in the hydrological processes of these three apparently similar volcanic basins. The results highlight the complexity of hydrological processes in volcanic mountainous systems and suggest that a simple watershed classification approach based on widely available data is insufficient. Simple local analyses such as specific flow analysis with a review of the geology and morphology can contribute to a better understanding of the hydrology of volcanic mountainous areas.     

The groundwater recharge regime of some slightly metamorphosed neoproterozoic sedimentary rocks: an application of natural environmental tracers

Isotope data of precipitation and groundwater in parts of the Voltaian Basin in Northern Ghana were used to explain the groundwater recharge regime in the area. Groundwater recharge is an important parameter in the development of a decision support system for the management and efficient utilization of groundwater resources in the area. It is therefore important to establish the processes and sources of groundwater recharge. δ¹⁸O and δ²H data for local precipitation suggest enrichment relative to the Global Meteoric Water Line (GMWL) and indicate that precipitation takes place at a relative humidity less than 100%. The groundwater data plot on an evaporation line with a slope of 5, suggesting a high degree of evaporative enrichment of the precipitation in the process of vertical infiltration and percolation through the unsaturated zone into the saturated zone. This finding is consistent with the observation of high evapotranspiration rates in the area and ties in with the fact that significant clay fraction in the unsaturated zone limits vertical percolation and thus exposes the percolating rainwater to the effects of high temperatures and low humidities resulting in high evapotranspiration rates. Groundwater recharge estimates from the chloride mass balance, CMB, method suggest recharge in the range of 1.8–32% of the annual average precipitation in the form of rainfall. The highest rates are associated with areas where open wells encourage significant amount of groundwater recharge from precipitation in the area. In the northern parts of the study area, groundwater recharge is lower than 12%. The recharge so computed through the application of the CMB methodology takes on a spatial distribution akin to the converse of the spatial pattern of both δ¹⁸O and δ²H in the area. As such, the locations of the highest recharge are associated with the most depleted values of the two isotopes. This observation is consistent with the assertion that low vertical hydraulic conductivities slow down vertical percolation of precipitation down to the groundwater water. The percolating precipitation water thus gets enriched in the heavier isotopes through high evapotranspiration rates. At the same time, the amount of water that finally reaches the water table is considerably reduced. Copyright © 2013 John Wiley & Sons, Ltd.     

A new standardized Palmer drought index for hydro‐meteorological use

Accepting the concept of standardization introduced by the standardized precipitation index, similar methodologies have been developed to construct some other standardized drought indices such as the standardized precipitation evapotranspiration index (SPEI). In this study, the authors provided deep insight into the SPEI and recognized potential deficiencies/limitations in relating to the climatic water balance it used. By coupling another well‐known Palmer drought severity index (PDSI), we proposed a new standardized Palmer drought index (SPDI) through a moisture departure probabilistic approach, which allows multi‐scalar calculation for accurate temporal and spatial comparison of the hydro‐meteorological conditions of different locations. Using datasets of monthly precipitation, temperature and soil available water capacity, the moisture deficit/surplus was calculated at multiple temporal scales, and a couple of techniques were adopted to adjust corresponding time series to a generalized extreme value distribution out of several candidates. Results of the historical records (1900–2012) for diverse climates by multiple indices showed that the SPDI was highly consistent and correlated with the SPEI and self‐calibrated PDSI at most analysed time scales. Furthermore, a simple experiment of hypothetical temperature and/or precipitation change scenarios also verified the effectiveness of this newly derived SPDI in response to climate change impacts. Being more robust and preferable in spatial consistency and comparability as well as combining the simplicity of calculation with sufficient accounting of the physical nature of water supply and demand relating to droughts, the SPDI is promising to serve as a competent reference and an alternative for drought assessment and monitoring. Copyright © 2013 John Wiley & Sons, Ltd.