Effect of Sorption Heterogeneity on Moments of Solute Residence Time in Convergent Flows

In this paper, we analyze the impact of physical and chemical heterogeneity on solute travel time to a pumping well. We consider a solute undergoing reversible linear instantaneous equilibrium sorption. Both the distribution coefficient, K _d , and the transmissivity field, T, are considered spatially variable, and are modeled as partially correlated spatial random functions. Groundwater flow and solute transport are then solved within the context of a numerical Monte Carlo framework. The results are analyzed on the basis of dimensional analysis techniques. Simple and compact expressions characterizing the dependence of the target travel time moments on relevant dimensionless groups are proposed. The functional form of these expressions is inspired by, and is consistent with, the previous works of Sanchez-Vila and Rubin (Water Resour. Res. 39(4):1086, 2003) and Riva et al. (J. Contam. Hydrol. 82:23–43, 2006) A key result is that the effects of the chemical and physical heterogeneities on the mean travel time can be decoupled consistently with existing analytical results. The relative role of physical and geochemical heterogeneities in travel time variance is more complex, and such a decoupling is not observed. Potential uses of this work include the assessment of aquifer reclamation time by means of a single pumping well.     

Comment on: ‘High-temperature dehydration melting and decompressive P–T path in a granulite complex from the Eastern Ghats India’ by S. Bhattacharya and R. Kar

Thermodynamic and chemographic modelling of complex reaction textures observed in Mg-Al-rich pelitic granulites is an important tool to unravel the P–T evolutionary history of high-grade rocks. In the Eastern Ghats Belt, India, several studies have been carried out on these fascinating aluminous granulites, and the results of these studies have revealed complex P–T–t histories (Dasgupta and Sengupta 1995; Sengupta et al. 1999; Rickers et al. 2001a, 2001b; Gupta et al. 1999; Dobmeier and Simmat 2002; Dobmeier and Raith 2003). In recent communication, Bhattacharya and Kar (2002) reported reaction textures from a suite of Mg-Al granulites from the Paderu area of the Eastern Ghats Belt. Combining the textural relations and thermodynamic calibration of some construed reactions, the authors have put forward a single phase metamorphic evolution of the area along a clockwise pressure–temperature trajectory. Combining the petrological features from the Paderu area with those reported from the Chilka Lake complex, the authors proposed a general tectonic model for the entire Eastern Ghats Belt. Incidentally, the rocks in and around Paderu have been studied in some detail by several other workers (Lal et al. 1987; Mohan et al. 1997; Sengupta et al. 1997). The purpose of this comment is to demonstrate that the conclusions made in the paper are inconsistent with the petrological features described in the text. Further, the thermodynamic treatment used in the paper has serious errors in many places, and hence, is often in complete disagreement with the existing experimental data and theoretical analyses on the Mg-Al-rich assemblages. There are also significant problems arising from the poor quality of the analytical database. Unfortunately, the authors cite only a few published works (mostly their own) ignoring many other relevant studies from this belt (cited above). Our observations are organised according to the sections of the paper.Editorial responsibility: T.L. Grove     

Residual and Fully Softened Strength Evaluation of Soils using Artificial Neural Networks

A backpropagation artificial neural network (ANN) model is developed to predict the secant friction angle of residual and fully softened soils, using data reported by Stark et al. (J Geotech Geoenviron Eng ASCE 131:575–588, 2005). In the ANN model, index properties such as liquid limit, plastic limit, activity, clay fraction and effective normal stress are used as input variables while secant residual friction angle is used as output variable. The model is verified using data that were not used for model training and testing. The results also indicate that the secant residual friction angle of cohesive soils can be predicted quite accurately using liquid limit, clay fraction and effective normal stress as input variables with R ² = 0.93. The sensitivity analysis results indicate that plastic limit and activity have no appreciable effect on ANN predicted secant friction angles. The secant friction angle predictions of the ANN model were also compared with those of Stark’s et al. (2005) curves and the empirical formulas suggested for the same data sets by Wright (Evaluation of soil shear strengths for slope and retaining wall stability with emphasis on high plasticity clays, 2005). The comparison shows that the ANN model predictions are very close to those suggested by the Stark et al. (2005) curves but much better than the prediction of Wright’s (2005) empirical equations. The results also show that ANN is an alternative powerful tool to predict the secant friction angle of soils.     

An integrated model simulating the initiation and motion of earthquake and rain induced rapid landslides and its application to the 2006 Leyte landslide

A gigantic rapid landslide claiming over 1,000 fatalities was triggered by rainfalls and a small nearby earthquake in the Leyte Island, Philippines in 2006. The disaster presented the necessity of a new modeling technology for disaster risk preparedness which simulates initiation and motion. This paper presents a new computer simulation integrating the initiation process triggered by rainfalls and/or earthquakes and the development process to a rapid motion due to strength reduction and the entrainment of deposits in the runout path. This simulation model LS-RAPID was developed from the geotechnical model for the motion of landslides (Sassa 1988) and its improved simulation model (Sassa et al. 2004b) and new knowledge obtained from a new dynamic loading ring shear apparatus (Sassa et al. 2004a). The examination of performance of each process in a simple imaginary slope addressed that the simulation model well simulated the process of progressive failure, and development to a rapid landslide. The initiation process was compared to conventional limit equilibrium stability analyses by changing pore pressure ratio. The simulation model started to move in a smaller pore pressure ratio than the limit equilibrium stability analyses because of progressive failure. However, when a larger shear deformation is set as the threshold for the start of strength reduction, the onset of landslide motion by the simulation agrees with the cases where the factor of safety estimated by the limit equilibrium stability analyses equals to a unity. The field investigation and the undrained dynamic loading ring shear tests on the 2006 Leyte landslide suggested that this landslide was triggered by the combined effect of pore water pressure due to rains and a very small earthquake. The application of this simulation model could well reproduce the initiation and the rapid long runout motion of the Leyte landslide.     

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.     

GIS-based determination of the mean long-term groundwater recharge in Lower Saxony

This paper presents the findings of a collaborative research project of the Geological Survey of Lower Saxony (NLfB) and the Programme Group Systems Analysis and Technology Evaluation (STE) of Research Centre Jülich on the GIS-based determination of the mean long-term groundwater recharge in Lower Saxony using high-resolution digital data (Dörhöfer and others 2001). The model calculations were performed on the basis of the water-balance model GROWA (Kunkel and Wendland 2002) with a spatial resolution of 100x100 m². The accuracy of the calculated groundwater recharge values for the period 1961–1990 was verified on the basis of data from gauging stations and displayed a good agreement between observed runoff values and model results.     

Latent Variable Modeling for Integrating Output from Multiple Climate Models

Numerical models of atmosphere–ocean circulation are widely used to understand past climate and to project future climate change. Although the same laws of physics, chemistry, and fluid dynamics govern any general circulation model, each model’s formulations and parameterizations are different, yielding different projections. Notwithstanding, models within an ensemble will have varying degrees of similarity for different outputs of interest. Multi-model ensembles have been used to increase forecast skill by using simple or weighted averages where weights have been obtained by considering factors such as estimated model bias and consensus with other models (Giorgi and Mearns, J. Clim. 15:1141–1158, 2002, Geophys. Res. Lett. 30:1629–1632, 2003; Tebaldi et al., Geophys. Res. Lett. 31:L24213, 2004, J. Clim. 18:1524–1540, 2005). This paper considers an alternative view of multi-model ensembles. For use with the North American Regional Climate Change Assessment Program (NARCCAP), multivariate statistical models are employed to characterize modes of similarity within the members of an ensemble. Specifically, we propose a spatially-correlated latent variable model which facilitates the exploration of when, where, and how regional climate models are similar, and what factors best predict observed locations of model convergence.     

The inclusive and simplified forms of Bayesian interpolation for general and monotonic models using Gaussian and Generalized Beta distributions with application to Monte Carlo simulations

A recently developed Bayesian interpolation method (BI) and its application to safety assessment of a flood defense structure are described in this paper. We use a one-dimensional Bayesian Monte Carlo method (BMC) that has been proposed in (Rajabalinejad 2009) to develop a weighted logical dependence between neighboring points. The concept of global uncertainty is adequately explained and different uncertainty association models (UAMs) are presented for linking the local and global uncertainty. Based on the global uncertainty, a simplified approach is introduced. By applying the global uncertainty, we apply the Guassian error estimation to general models and the Generalized Beta (GB) distribution to monotonic models. Our main objective in this research is to simplify the newly developed BMC method and demonstrate that it can dramatically improve the simulation efficiency by using prior information from outcomes of the preceding simulations. We provide theory and numerical algorithms for the BI method geared to multi-dimensional problems, integrate it with a probabilistic finite element model, and apply the coupled models to the reliability assessment of a flood defense for the 17th Street Flood Wall system in New Orleans.     

Model or measurements? A discussion of the key issue in Chapman and Pollack’s critique of Hamza et al.’s re-evaluation of oceanic heat flux and the global power

Chapman and Pollack (C and P)[2007, Int J Earth Sci] criticize Hamza et al. [2007, Int J Earth Sci] for using actual heat flux measurements in young oceanic crust instead of values from 1-D cooling models. The rationalization of C and P and previous authors is that hydrothermal circulation causes the discrepancy between model and measurement. However, the discrepancy between model values and measured heat flux exists over the entire ocean floor and is opposite to the perturbations that hydrothermal circulation would superimpose on a conductive system [Hofmeister and Criss (2005) Tectonophysics 409:199–203]. The error lies in force-fitting a 1-D cooling model to the 3-D oceanic crust [Hofmeister and Criss (2005) Tectonophysics 395:159–177]. Shortcomings of the 1-D model include mathematical errors, such as use of volumetric rather than linear thermal expansivity to describe contraction which, by assumption, is limited only to the Z -direction [Hofmeister and Criss (2006) Tectonophysics]. This 3× error, traceable to McKenzie and Sclater [1969, Bull Vocanol 33–1:101–118], accidentally provides good agreement of model values with globally averaged seafloor depths for young, but not old ages, and is the sole rationale for using the simplistic cooling model. There is no justification for selective substitution of erroneous 1-D model values for measurements only for the younger half of the 3-D oceanic crust, as stridently and arbitrarily promoted by C and P. Hamza et al. [2007, Int J Earth Sci], in contrast, use the scientific method, which calls for discarding models that do not well describe physical phenomena. The remainder of this report summarizes the shortcomings of cooling models, particularly the half-space cooling (HSC) model touted by C and P, and explains how hydrothermal circulation affects heat flux. We focus on the basics, as these have been misunderstood. With the key issues of C and P being erroneous, it is not necessary to address their remaining comments, many of which enumerate the vote for an imagined, gargantuan circulation of hot fluid through oceanic basins that is somehow warmed without removing heat from the rocks. The use of “consensus” to belittle valid challenge is the enemy of the scientific method.     

Seismic liquefaction analysis of a reservoir dam foundation in the South–North Water Diversion project in China. Part I: Liquefaction potential assessment

In this paper, the dynamic characteristics of a liquefiable silt substratum within the foundation soil of a reservoir dam in the Tianjin area are investigated by means of standard penetration resistance and dynamic triaxial tests. Properties including N-values, factors influencing liquefaction as a cyclic stress, consolidation pressure, structure, and particle composition are considered in this research. Parameters used to evaluate liquefaction potential are obtained through testing. A comprehensive program based on the Chinese code and standard for geological investigation (Ministry of Water Resources of China 1999a; Ministry of Construction of China 2001a) and Seed’s simplified method (Seed and Idriss 1971; in J Geotech Eng Div ASCE 109(3): 458–482, 1983) was carried out to evaluate the potential of liquefaction within the reservoir dam foundation. Liquefaction potentials were also assessed in response to the Chinese codes for seismic design (Hydropower Research Institute of China 2000; Ministry of Construction of China 2001b). The evaluation shows that saturated surface silt in the reservoir dam foundations is vulnerable to liquefaction at seismic intensities of VII and above. The two assessment methods are in good agreement with each other, and the research results can provide useful information for the safe construction and normal operation of the reservoir.     

Environmental hazards posed by the Los Angeles Basin urban oilfields: an historical perspective of lessons learned

Urban encroachment into areas historically reserved for oil and gas field operations is an ever-present problem within the Los Angeles Basin. The recent frenzy in real estate development has only intensified what can be characterized as a conflict in land usage. Subsurface mineral rights are severed from surface ownership, often resulting in developments being approved without adequate consideration of the underlying oil and gas field consequences. Also, surface operations are frequently co-located within residential areas without consideration of the health and safety consequences of emissions of toxics to air. This paper presents a review of the environmental, health and safety hazards posed by urban oilfield operations, with an emphasis upon the lessons learned from the L.A. Basin: Original Urban Oilfield Legend (see Castle and Yerkes 1976; Denton and others 2001; Endres and others 2002; Kouznetsov and others 1994; Katz and others 1994; Schumacher and Abrams 1994; and Schoell 1983). The Los Angeles Basin has provided the authors with one of the largest natural laboratories in the world for studying the consequences of these issues. The results presented are part of a long-term research program based upon the application of geoscience and petroleum engineering principles in obtaining a fundamental understanding of the root causes of the environmental hazards posed. Topics addressed include: (1) vertical migration of gas to the surface along faults and improperly completed or abandoned wellbores (e.g., due to poor cementing practices), (2) subsidence caused by the fluid production and declining reservoir pressures, (3) soil and groundwater contamination resulting from historic oil and gas field operations, and (4) air toxics resulting from surface operations. A number of case histories are discussed that illustrate the seriousness of the problem. A clear case is made for the urgent need for closer coordination and education by the petroleum industry of the local government planning departments. These departments have the principal role in determining land use policies, acting as the lead agency in performing environmental site assessments (e.g., under the California Environmental Quality Act), and in establishing mitigation measures for dealing with the long-term environmental hazards. This paper establishes prudent practices on the part of oilfield operators for the monitoring and mitigation of these hazards.     

Stochastic methods for safety assessment of the flood defense system in the Scheldt Estuary of the Netherlands

In this study, we address the effective method to apply a novel reliability method integrated with finite element models to the safety assessment of pilot site Scheldt in the Netherlands. This site was considered as one of the three main pilot sites in Europe to assess the application of newly suggested techniques in order to reduce and manage the flood risk in the Floodsite project. , 2004–2009). The novel method of dynamic bounds (DB) is applied to this site after a successful experience in (Rajabalinejad in Reliability methods for finite element models, 1 edn. IOS Press, Amsterdam, 2009). In this study, the bi-functional response of the finite element model is considered, and the dimensional uncertainty is defined presenting the expected uncertainty for a certain dimension in the DB method. The uncertainty is used as a judgment tool to choose the dimension for the DB method for the desired accuracy. The results obtained by applying this technique are presented in this paper.     

Extended paraelasticity and its application to a boundary value problem

Paraelasticity (PE) is a hysteretic, fully reversible model with rate-independent damping and with variable secant stiffness. PE has been proposed in the companion paper (Niemunis et al. in Acta Geotech, 2011) for the 3D case. Here, paraelasticity is extended for demanding application in geomechanics by implementation of a reversible dilatancy–contractancy. Moreover, an implicit FE implementation for Abaqus and calculation of a boundary value problem are presented.     

Thermochemistry of sulfide liquids III: Ni-bearing liquids at 1 bar

Experiments were performed in the system O–S–Fe–Ni designed to extend our understanding of the chemistry of sulfide liquids. Results indicate that adding nickel to Fe-rich sulfide liquids in equilibrium with silicate liquids extends their stability field to much higher oxygen fugacities and lower sulfur fugacities. Increasing Ni/Fe at a given temperature and sulfur and oxygen fugacity is accompanied by a significant decrease in the oxygen content of the sulfide liquid. Results of these experiments are combined with data from the literature to calibrate an associated regular solution model for O–S–Fe–Ni liquids. This model represents a complete refit of the associated regular solution model of Kress (Contrib Mineral Petrol 139:316–325, 2000). The resulting model is combined with the olivine solution model of Hirschmann (Am Mineral 76:1232–1248, 1991) to explore the effect of variations in oxygen and sulfur fugacities on the distribution of Fe and Ni between olivine and sulfide liquid. Predicted olivine–sulfide distribution trends parallel those observed by Gaetani and Grove (Geochim Cosmochim Acta 61:1829–1846, 1997), Gaetani and Grove (Earth Planet Sci Lett 169:147–163, 1999), Brenan and Caciagli (Geochim Cosmochim Acta 64:307–320, 2000) and Brenan (Geochim Cosmochim Acta 67:2663–2681, 2003), but are systematically offset toward lower predicted Ni in the sulfide. Nevertheless our results are consistent with the assertion that low K _D ^os values in magmatic ore deposits such as the J-M Reef reflect high iron contents in the sulfides combined with relatively high oxygen fugacities.

---

Experimental Validation of an Method for In Situ Estimation of Hydraulic Conductivity of Water Treatment Granular Materials

Domestic wastewater-treatment system is today a widely held technique. In order to evaluate the efficiency of these systems, it is necessary to determinate the in situ hydraulic conductivity of the water-treatment granular materials constituting them. The in situ measurement of the hydraulic conductivity of soils has proven to be imprecise, take a long time. Empirical equations based on the physical properties of soil have been proposed to overcome these difficulties. In this context, the originality of this paper is to propose an in situ method to obtain reliable input parameters for the predictive equations of Chapuis (Can Geotech J 41(5): 787–795, 2004) and Carrier (J Geotech Geoenviron Eng 129(11): 1054–1056, 2003) by using both a penetrometer and a geoendoscope. This method is described in the first part of this paper. Then, the validation of the method based on laboratory tests performed under controlled conditions for three kinds of soils is presented. Hydraulic conductivity obtained when applying the method is compared to that measured by a Darcy permeameter. The difference between these two hydraulic conductivities is less than 25%. Finally, the precision of the results is discussed.     

How to cut a sediment core for <Superscript>210</Superscript>Pb geochronology: a supplement

To obtain the profile of excess ²¹⁰Pb in sediment column, it is necessary to section the collected core to date. Sectioning thickness is generally not explicit. An approach when sedimentation rate is constant has been suggested by Lu and Matsumoto (Environ Geol 47:804–810, 2005). This paper addresses the case when sedimentation rate is variable. Actually, it is possible to assume that sedimentation rate is constant within the sectioning thickness and based on this concept, a tentative approach for determining reasonable sectioning thickness is proposed.     

Oxygen isotope heterogeneity of the mantle beneath the Canary Islands: a discussion of the paper of Gurenko et al.

Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) report laser-assisted fluorination (LF) and secondary ionization mass spectrometry (SIMS) ¹⁸O/¹⁶O datasets for olivine grains from the Canary Islands of Gran Canaria, Tenerife, La Gomera, La Palma and El Hierro. As with prior studies of oxygen isotopes in Canary Island lavas (e.g. Thirlwall et al. Chem Geol 135:233–262, 1997; Day et al. Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010), these authors find variations in δ¹⁸O_ol (~4.6–6.0 ‰) beyond that measured for mantle peridotite olivine (Mattey et al. Earth Planet Sci Lett 128:231–241, 1994) and interpret this variation to reflect contributions from pyroxenite-peridotite mantle sources. Furthermore, Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) speculate that δ¹⁸O_ol values for La Palma olivine grains measured by LF (Day et al. Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010) may be biased to low values due to the presence of altered silicate, possibly serpentine. The range in δ¹⁸O_ol values for Canary Island lavas are of importance for constraining their origin. Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) took a subset (39 SIMS analyses from 13 grains from a single El Hierro lava; EH4) of a more extensive dataset (321 SIMS analyses from 110 grains from 16 Canary Island lavas) to suggest that δ¹⁸O_ol is weakly correlated (R ² = 0.291) with the parameter used by Gurenko et al. (Earth Planet Sci Lett 277:514–524, 2009) to describe the estimated weight fraction of pyroxenite-derived melt (Xpx). With this relationship, end-member δ¹⁸O values for HIMU-peridotite (δ¹⁸O = 5.3 ± 0.3 ‰) and depleted pyroxenite (δ¹⁸O = 5.9 ± 0.3 ‰) were defined. Although the model proposed by Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) implicates similar pyroxenite-peridotite mantle sources to those proposed by Day et al. (Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010) and Day and Hilton (Earth Planet Sci Lett 305:226–234, 2011), there are significant differences in the predicted δ¹⁸O values of end member components in the two models. In particular, Day et al. (Geochim Cosmochim Acta 74:6565–6589, 2010) proposed a mantle source for La Palma lavas with low-δ¹⁸O (<5 ‰), rather than higher-δ¹⁸O (c.f. the HIMU-peridotite composition of Gurenko et al. in Contrib Mineral Petrol 162:349–363, 2011). Here we question the approach of using weakly correlated variations in δ¹⁸O_ol and the Xpx parameter to define mantle source oxygen isotope compositions, and provide examples of why this approach appears flawed. We also provide reasons why the LF datasets previously published for Canary Island lavas remain robust and discuss why LF and SIMS data may provide complementary information on oxygen isotope variations in ocean island basalts (OIB), despite unresolved small-scale uncertainties associated with both techniques.     

Climate change,environmental degradation and migration

Climate change will have a progressively increasing impact on environmental degradation and environmentally dependent socio-economic systems with potential to cause substantial population displacement. The key concerns in Less Developed Countries (LDCs) will include serious threats to food security and health, considerable economic decline, inundation of coastal areas, and degradation of land and fresh water resources (Reuveny in Polit Geogr, 2007). The relationship between environmental change and potential humanitarian crises has been captured by: McGregor (Geography and refugees: patterns and processes of change, Belhaven Press, London, pp 159–70, 1993), Kibreab (Environment and Population Change, International Union for the Scientific Study of Population, Liège, 1994), Kibreab (Disasters 21(1):20–38, 1997), Myers (Bioscience 43:752–761, 1993), Myers and Kent (Environmental exodus: an emergent crisis in the global arena, Climate Institute, Washington, DC, 1995), Black (New Issues in Refugee Research, Working Paper no. 34, 2001), Lee (Environmental matters: conflict, refugees and international relations, World Human Development Institute Press, Seoul and Tokyo, 2001), Castles (Environmental Change and Induced Migration: Making Sense of the Debate Working Paper No. 70, 2002), Christian Aid (Human tide: the real migration crisis, Christian Aid, London, 2007), and Massey et al. (, 2007). However, we know little about the interplay between environmental change and stresses on ecological systems, resulting socio-economic vulnerability and potential outcomes in terms of population displacement or induced migration. So far these relationships are poorly conceptualized, lack systematic investigation, and are reduced to simplistic causal explanations. This leads to misleading conclusions that deny the complex multivariate processes—environmental, political, social, and economic— which are the root causes of environmentally induced migration and/or conflict. When people are faced with severe environmental degradation they have one of three options: (1) stay and adapt to mitigate the effects; (2) stay, do nothing and accept a lower quality of life; or (3) leave the affected area. The process of movement and migration is usually subject to a complex set of push and pull forces, where push forces relate to the source area while pull factors relate to the destination. These forces are in constant flux, as much as environmental change, and interact with socio-economic and political conditions including state or government decision making powers, which can tip the balance at any point by either denying movement or the right to settle elsewhere. The paper focuses on how environmental change and environmental hazards contribute to the migration by exploring the mechanisms through which vulnerability and migration are linked—via livelihoods, relocation policies, and other factors. The paper begins by outlining important definitions of what is environmentally induced migration. The paper also considers the question of whether migration is a process that reduces or increases vulnerability. The paper draws on multidisciplinary literature including ecology, environment, and climate change; sociology of migration; anthropology of displacement; and economics; but also on preliminary from various case studies in Egypt, Vietnam, and Mozambique.     

Soil/Geotextile Interface Behaviour in Direct Shear and Pullout Movements

This paper deals with soil/reinforced geotextile interface behaviour in direct shear and pullout movements. A soil/geosynthetic direct shear apparatus, developed in accordance to EN ISO 12957-1 (2004) and ASTM D5321-92 (1992) standards, is presented. Some details of the test (like: specimen fixation, influence of the vertical stress on the registered horizontal force, type of test and measurement of the vertical displacement) are discussed and modifications in the test procedures are adopted. Then, the reinforced geotextile and the residual soil of granite used in the research are described. The behaviour of soil/geosynthetic interface in direct shear is characterized based on modified direct shear tests and maximum interface shear stress is determined at peak and residual for a confining stress of 50 kPa. The modification in the pullout test apparatus described by Lopes and Ladeira (1996a; 1996b) and by Lopes and Lopes (1999) are noticed. The behaviour of soil/geosynthetic interface in pullout is characterized based on pullout tests performed, in accordance EN 13738 (2004), and interface shear stress at maximum pullout force is defined for a confining stress of 50 kPa. Finally, the values of interface coefficient at soil/geosynthetic interface are obtained in direct shear and in pullout and then compared. The main conclusions that can be outlined from the present study are the following: modifications should be made to EN ISO 12957-1 (2004) standard, namely in what concerns the dependence of the measured horizontal force from the vertical stress, the difficulties to perform constant area direct shear test with the lower half box filled with soil and on the measurement of the vertical displacement of the load plate; on the contrary to which is normally accepted the characteristics of the behaviour of soil/geosynthetic interface in pullout, when the geosynthetic has a full plane contact area with the soil, are not able to be obtained based on results of direct shear tests, as in this type of tests the contribution of the geosynthetic deformation on the characteristics of the interface in pullout is not considered.