A thermo‐mechanical model for variably saturated soils based on hypoplasticity

The paper presents a mechanical model for non‐isothermal behaviour of unsaturated soils. The model is based on an incrementally non‐linear hypoplastic model for saturated clays and can therefore tackle the non‐linear behaviour of overconsolidated soils. A hypoplastic model for non‐isothermal behaviour of saturated soils was developed and combined with the existing hypoplastic model for unsaturated soils based on the effective stress principle. Features of the soil behaviour that are included into the model, and those that are not, are clearly distinguished. The number of model parameters is kept to a minimum, and they all have a clear physical interpretation, to facilitate the model usefulness for practical applications. The step‐by‐step procedure used for the parameter calibration is described. The model is finally evaluated using a comprehensive set of experimental data for the thermo‐mechanical behaviour of an unsaturated compacted silt. Copyright © 2011 John Wiley & Sons, Ltd.     

Flood frequency analysis of Turkey using L‐moments method

The index flood procedure coupled with the L‐moments method is applied to the annual flood peaks data taken at all stream‐gauging stations in Turkey having at least 15‐year‐long records. First, screening of the data is done based on the discordancy measure (D_i) in terms of the L‐moments. Homogeneity of the total geographical area of Turkey is tested using the L‐moments based heterogeneity measure, H, computed on 500 simulations generated using the four parameter Kappa distribution. The L‐moments analysis of the recorded annual flood peaks data at 543 gauged sites indicates that Turkey as a whole is hydrologically heterogeneous, and 45 of 543 gauged sites are discordant which are discarded from further analyses. The catchment areas of these 543 sites vary from 9·9 to 75121 km² and their mean annual peak floods vary from 1·72 to 3739·5 m³ s^?1. The probability distributions used in the analyses, whose parameters are computed by the L‐moments method are the general extreme values (GEV), generalized logistic (GLO), generalized normal (GNO), Pearson type III (PE3), generalized Pareto (GPA), and five‐parameter Wakeby (WAK). Based on the L‐moment ratio diagrams and the |Z_dist|‐statistic criteria, the GEV distribution is identified as the robust distribution for the study area (498 gauged sites). Hence, for estimation of flood magnitudes of various return periods in Turkey, a regional flood frequency relationship is developed using the GEV distribution. Next, the quantiles computed at all of 543 gauged sites by the GEV and the Wakeby distributions are compared with the observed values of the same probability based on two criteria, mean absolute relative error and determination coefficient. Results of these comparisons indicate that both distributions of GEV and Wakeby, whose parameters are computed by the L‐moments method, are adequate in predicting quantile estimates. Copyright © 2011 John Wiley & Sons, Ltd.     

A discrete particle representation of hillslope hydrology: hypothesis testing in reproducing a tracer experiment at Gårdsjön,Sweden

Despite the long history of the continuum equation approach in hydrology, it is not a necessary approach to the formulation of a physically based representation of hillslope hydrology. The Multiple Interacting Pathways (MIPs) model is a discrete realization that allows hillslope response and transport to be simultaneously explored in a way that reflects the potential occurrence of preferential flows and lengths of pathways. The MIPs model uses random particle tracking methods to represent the flow of water within the subsurface alongside velocity distributions that acknowledge preferential flows and transition probability matrices, which control flow pathways. An initial realization of this model is presented here in application to a tracer experiment carried out in Gårdsjön, Sweden. The model is used as an exploratory tool, testing several hypotheses in relation to this experiment. Copyright © 2011 John Wiley & Sons, Ltd.     

Process identification at a slow‐moving landslide in the Vorarlberg Alps

A fine‐grained slope that exhibits slow movement rates was investigated to understand how geohydrological processes contribute to a consecutive development of mass movements in the Vorarlberg Alps, Austria. For that purpose intensive hydrometeorological, hydrogeological and geotechnical observations as well as surveying of surface movement rates were conducted during 1998–2001. Subsurface water dynamics at the creeping slope turned out to be dominated by a three‐dimensional pressure system. The pressure reaction is triggered by fast infiltration of surface water and subsequent lateral water flow in the south‐western part of the hillslope. The related pressure signal was shown to propagate further downhill, causing fast reactions of the piezometric head at 5·5 m depth on a daily time scale. The observed pressure reactions might belong to a temporary hillslope water body that extends further downhill. The related buoyancy forces could be one of the driving forces for the mass movement. A physically based hydrological model was adopted to model simultaneously surface and subsurface water dynamics including evapotranspiration and runoff production. It was possible to reproduce surface runoff and observed pressure reactions in principle. However, as soil hydraulic functions were only estimated on pedotransfer functions, a quantitative comparison between observed and simulated subsurface dynamics is not feasible. Nevertheless, the results suggest that it is possible to reconstruct important spatial structures based on sparse observations in the field which allow reasonable simulations with a physically based hydrological model. Copyright © 2004 John Wiley & Sons, Ltd.     

Least‐square support vector machine applied to settlement of shallow foundations on cohesionless soils

This paper examines the potential of least‐square support vector machine (LSVVM) in the prediction of settlement of shallow foundation on cohesionless soil. In LSSVM, Vapnik's ε‐insensitive loss function has been replaced by a cost function that corresponds to a form of ridge regression. The LSSVM involves equality instead of inequality constraints and works with a least‐squares cost function. The five input variables used for the LSSVM for the prediction of settlement are footing width (B), footing length (L), footing net applied pressure (P), average standard penetration test value (N) and footing embedment depth (d). Comparison between LSSVM and some of the traditional interpretation methods are also presented. LSSVM has been used to compute error bar. The results presented in this paper clearly highlight that the LSSVM is a robust tool for prediction of settlement of shallow foundation on cohesionless soil. Copyright © 2008 John Wiley & Sons, Ltd.     

Assessment of δ13C and C/N ratios in bulk organic matter as palaeosalinity indicators in Holocene and Lateglacial isolation basin sediments,northwest Scotland

Carbon isotopes (δ¹³C) and C/N ratios from bulk organic matter have recently been used as alternative proxies for relative sea‐level (RSL) reconstruction where there are problems associated with conventional biological indictors. A previous study on a single isolation basin (Upper Loch nan Eala) in northwest Scotland has shown a clear relationship between δ¹³C, C/N ratios and palaeosalinity from Younger Dryas and Holocene aged sediments. In this paper we present results of δ¹³C and C/N ratio analyses from other isolation basins in northwest Scotland over the Holocene and the Lateglacial period in order to validate this technique. The results from the Holocene sequences support the earlier findings that this technique can be used to identify RSL change from isolation basins over the Holocene in this region. The relationship between δ¹³C, C/N ratios and RSL change is not apparent in sediments of Lateglacial age. Other environmental variables such as atmospheric CO₂ concentration, poor vegetation development and temperature influence δ¹³C values during this period. Copyright © 2007 John Wiley & Sons, Ltd.     

Holocene coastal environmental changes on the periphery of an area of glacio‐isostatic uplift: an example from Scapa Bay,Orkney, UK

The first detailed investigation of a deep, coastal, sedimentary basin in Orkney reveals a complex Holocene history of back‐barrier morphodynamics. At Scapa Bay, the sea flooded a freshwater marsh after ca. 9400 yr BP at ca. ?5.4 m OD. Before ca. 7800 BP, abundant sediment from nearby cliffs was mobilised inland into a series of gravel barriers across the valley mouth. By ca. 7500 BP, direct marine influence was restricted in the back‐barrier area, although saltmarsh persisted until ca. 5900 BP. By then, at least four gravel ridges had enclosed the backing lagoon, where freshwater inputs became dominant. As terrestrial sediments filled the basin, another freshwater marsh developed. The multiple barrier complex demonstrates progradation resulting from continuous sediment supply in a sheltered embayment. The progressively rising height of the barrier crests seawards probably resulted from a combination of factors such as barrier morphodynamics, increased storminess and long‐term rising relative sea levels. The dominant vegetation surrounding Scapa Bay changed from open grassland to scrub ca. 9400 BP, then to deciduous woodland ca. 7800 BP, and to dwarf‐shrub heath ca. 2600 BP, the latter probably a response to a combination of climate change and human activity. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatial scale effect on seasonal streamflows in permafrost catchments on the Qinghai–Tibet Plateau

The scale issue is of central concern in hydrological processes to understand the potential upscaling or downscaling methodologies, and to develop models for scaling the dominant processes at different scales and in different environments. In this study, a typical permafrost watershed in the Qinghai‐Tibet Plateau was selected. Its hydrological processes were monitored for 4 years from 2004 to 2008, measuring the effects of freezing and thawing depth of active soil layers on runoff processes. To identify the nature and cause of variation in the runoff response in different size catchments, catchments ranging from 1·07 to 112 km² were identified in the watershed. The results indicated that the variation of runoff coefficients showed a ‘V’ shape with increasing catchment size during the spring and autumn seasons, when the active soil was subjected to thawing or freezing processes. A two‐stage method was proposed to create runoff scaling models to indicate the effects of scale on runoff processes. In summer, the scaling transition model followed an exponential function for mean daily discharge, whereas the scaling model for flood flow exhibited a linear function. In autumn, the runoff process transition across multiple scales followed an exponential function with air temperature as the driving factor. These scaling models demonstrate relatively high simulation efficiency and precision, and provide a practical way for upscaling or downscaling runoff processes in a medium‐size permafrost watershed. For permafrost catchments of this scale, the results show that the synergistic effect of scale and vegetation cover is an important driving factor in the runoff response. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of fully saturated porous materials

Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd.     

Glacier fluctuations in the southern Peruvian Andes during the late‐glacial period,constrained with cosmogenic 3He

The occurrence of pronounced climate reversals during the last glacial termination has long been recognised in palaeoclimate records from both hemispheres and from high to low latitudes. Accurate constraint of both the timing and magnitude of events, such as the Younger Dryas and Antarctic Cold Reversal, is vital in order to test different hypotheses for the causes and propagation of abrupt climate change. However, in contrast to higher‐latitude regions, well‐dated records from the Tropics are rare and the structure of late‐glacial tropical climate remains uncertain. As a step toward addressing this problem, we present an in situ cosmogenic ³He surface exposure chronology from Nevado Coropuna, southern Peru, documenting a significant fluctuation of the ice margin during the late‐glacial period. Ten tightly clustered ages from a pair of moraines located halfway between the modern glacier and the Last Glacial Maximum terminus range from 11.9 to 13.9 ka and give an arithmetic mean age of 12.8 ± 0.7 ka (1σ). These data constitute direct evidence for a readvance, or prolonged stillstand, of glaciers in the arid Andes of southwestern Peru. Copyright © 2011 John Wiley & Sons, Ltd.