Characterizing rainfall–groundwater dynamics in a hard‐rock aquifer system using time series,geographic information system and geostatistical modelling

The aim of this study was to investigate rainfall–groundwater dynamics over space and annual time scales in a hard‐rock aquifer system of India by employing time series, geographic information system and geostatistical modelling techniques. Trends in 43‐year (1965–2007) annual rainfall time series of ten rainfall stations and 16‐year (1991–2006) pre‐monsoon and post‐monsoon groundwater levels of 140 sites were identified by using Mann–Kendall, Spearman rank order correlation and Kendall rank correlation tests. Trends were quantified by Kendall slope method. Furthermore, the study involves novelty of examining homogeneity of pre‐monsoon and post‐monsoon groundwater levels, for the first time, by applying seven tests. Regression analysis between rainfall and post‐monsoon groundwater levels was performed. The pre‐monsoon and post‐monsoon groundwater levels for four periods – 1991–1994, 1995–1998, 1999–2002 and 2003–2006 – were subjected to geographic information system‐based geostatistical modelling. The rainfall showed considerable spatiotemporal variations, with a declining trend at the Mavli rainfall station (p‐value < 0.05). The Levene's tests revealed spatial homogeneity of rainfall at α = 0.05. Regression analyses indicated significant relationships (r² > 0.5) between groundwater level and rainfall for eight rainfall stations. Non‐homogeneity and declining trends in the groundwater level, attributed to anthropogenic and hydrologic factors, were found at 5–61 more sites in pre‐monsoon compared with post‐monsoon season. The groundwater declining rates in phyllite–schist, gneiss, schist and granite formations were found to be 0.18, 0.26, 0.21 and 0.14 m year^?1 and 0.13, 0.19, 0.16 and 0.02 m year^?1 during the pre‐monsoon and post‐monsoon seasons, respectively. The geostatistical analyses for four time periods revealed linkages between the rainfall and groundwater levels. Copyright © 2013 John Wiley & Sons, Ltd.     

Characterizing and classifying urban watersheds with compositional and structural attributes

Current land-use classifications used to assess urbanization effects on stream water quality date back to the 1980s when limited information was available to characterize watershed attributes that mediate non-point source pollution. With high resolution remote sensing and widely used GIS tools, there has been a vast increase in the availability and precision of geospatial data of built environments. In this study, we leverage geospatial data to expand the characterization of developed landscapes and create a typology that allows us to better understand the impact of complex developed landscapes across the rural to urban gradient. We assess the ability of the developed landscape typology to reveal patterns in stream water chemistry previously undetected by traditional land-cover based classification. We examine the distribution of land-cover, infrastructure, topography and geology across 3876 National Hydrography Dataset Plus catchments in the Piedmont region of North Carolina, USA. From this dataset, we generate metrics to evaluate the abundance, density and position of landscape features relative to streams, catchment outlets and topographic wetness metrics. While impervious surfaces are a key distinguishing feature of the urban landscape, sanitary infrastructure, population density and geology are better predictors of baseflow stream water chemistry. Unsupervised clustering was used to generate a distinct developed landscape typology based on the expanded, high-resolution landscape feature information. Using stream chemistry data from 37 developed headwater catchments, we compared the baseflow water chemistry grouped by traditional land-cover based classes of urbanization (rural, low, medium and high density) to our composition and structure-based classification (a nine-class typology). The typology based on 22 metrics of developed landscape composition and structure explained over 50% of the variation in NO₃⁻-N, TDN, DOC, Cl⁻, and Br⁻ concentration, while the ISC-based classification only significantly explained 23% of the variation in TDN. These results demonstrate the importance of infrastructure, population and geology in defining developed landscapes and improving discrete classes for water management.     

Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data

Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources.     

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.     

Quantifying flows along hydrological pathways by applying a new filtering algorithm in conjunction with master recession curve analysis

Quantifying the proportion of the river hydrograph derived from the different hydrological pathways is essential for understanding the behaviour of a catchment. This paper describes a new approach using the output from master recession curve analysis to inform a new algorithm based on the Lyne and Hollick ‘one‐parameter’ signal analysis filtering algorithm. This approach was applied to six catchments (including two subcatchments of these) in Ireland. The conceptual model for each catchment consists of four main flow pathways: overland flow, interflow, shallow groundwater and deep groundwater. The results were compared with those of the master recession curve analysis, a recharge coefficient approach developed in Ireland and the semi‐distributed, lumped and deterministic hydrological model Nedbør‐Afstrømings‐Model. The new algorithm removes the ‘free variable’ aspect that is typically associated with filtering algorithms and provides a means of estimating the contribution of each pathway that is consistent with the results of hydrograph separation in catchments that are dominated by quick response pathways. These types of catchments are underlain by poorly productive aquifers that are not capable of providing large baseflows in the river. Such aquifers underlie over 73% of Ireland, ensuring that this new algorithm is applicable in the majority of catchments in Ireland and potentially in those catchments internationally that are strongly influenced by the quick‐responding hydrological pathways. Copyright © 2013 John Wiley & Sons, Ltd.     

智能卡和出租车轨迹数据中蕴含城市人群活动模式的差异性分析

当前采用交通流数据量化城市人群活动模式研究已经取得了丰硕的研究成果,但是对于同一区域、同一时段不同类型交通流数据反映城市人群活动模式的共性与差异性仍然知之甚少,直接影响了城市人群活动模式挖掘结果的可解释性与实际应用效果。为此,本文旨在对目前广泛采用的智能卡数据（公交和地铁刷卡）和出租车轨迹数据2种重要的交通流数据,从时空分布模式的差异性、行程距离及距离衰减效应的差异性、空间社团结构的差异性3个方面,探索二者反映城市人群活动模式的差异性：① 采用北京市六环以内区域2016年5月9日至15日的智能卡和出租车轨迹数据进行实验分析,研究发现：① 2种交通流反映出行需求的空间分布呈现出高度相关性,但是在同一空间单元上,2种交通流反映出行需求的时间相关性较低;② 2种交通流的使用率在不同空间位置存在明显差异,仅在城市中心区域使用率较为均衡;③ 2种交通流反映人群行程距离的空间分布、距离衰减效应存在明显差异,公共交通对于促进长距离出行更为重要;④ 从2种交通流发现的空间社团结构都显示了城市的多中心结构特征,但是二者发现社团结构存在的差异性表明两种交通方式对城市空间交互起着不同的作用。本研究有助于深入理解多源交通流反映城市人群活动的内在机理,提升城市人群活动模式在城市规划、交通管理等领域的应用效果。     

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.     

Three‐dimensional shakedown solutions for anisotropic cohesive‐frictional materials under moving surface loads

Previous work on three‐dimensional shakedown analysis of cohesive‐frictional materials under moving surface loads has been entirely for isotropic materials. As a result, the effects of anisotropy, both elastic and plastic, of soil and pavement materials are ignored. This paper will, for the first time, develop three‐dimensional shakedown solutions to allow for the variation of elastic and plastic material properties with direction. Melan's lower‐bound shakedown theorem is used to derive shakedown solutions. In particular, a generalised, anisotropic Mohr–Coulomb yield criterion and cross‐anisotropic elastic stress fields are utilised to develop anisotropic shakedown solutions. It is found that shakedown solutions for anisotropic materials are dominated by Young's modulus ratio for the cases of subsurface failure and by shear modulus ratio for the cases of surface failure. Plastic anisotropy is mainly controlled by material cohesion ratio, the rise of which increases the shakedown limit until a maximum value is reached. The anisotropic shakedown limit varies with frictional coefficient, and the peak value may not occur for the case of normal loading only. Copyright © 2013 John Wiley & Sons, Ltd.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Mechanical transitional behavior of binary mixtures via DEM: Effect of differences in contact-type friction coefficients

This report presents a numerical investigation of the shear behavior of binary mixtures via a two-dimensional discrete element method. The effect of the coarse particle content on the peak shear strength of mixtures is investigated for different contact-type friction coefficients. A detailed analysis on the anisotropies enables us to understand the microscopic mechanisms that result in the dependency of the peak shear strengths on the coarse particle content. The contributions of different contact types to the peak shear strength are quantified. Lastly, the inter-particles structures are examined when the binary mixtures become coarse-particle-supported structures.