Preliminary regional rockfall hazard mapping using LIDAR-based slope frequency distribution and CONEFALL modeling

A factor limiting preliminary rockfall hazard mapping at regional scale is often the lack of knowledge of potential source areas. Nowadays, high resolution topographic data (LiDAR) can account for realistic landscape details even at large scale. With such fine-scale morphological variability, quantitative geomorphometric analyses become a relevant approach for delineating potential rockfall instabilities. Using digital elevation model (DEM)-based ?slope families? concept over areas of similar lithology and cliffs and screes zones available from the 1:25,000 topographic map, a susceptibility rockfall hazard map was drawn up in the canton of Vaud, Switzerland, in order to provide a relevant hazard overview. Slope surfaces over morphometrically-defined thresholds angles were considered as rockfall source zones. 3D modelling (CONEFALL) was then applied on each of the estimated source zones in order to assess the maximum runout length. Comparison with known events and other rockfall hazard assessments are in good agreement, showing that it is possible to assess rockfall activities over large areas from DEM-based parameters and topographical elements.

Debris flow hazard modelling on medium scale: Valtellina di Tirano, Italy

Debris flow hazard modelling at medium (regional) scale has been subject of various studies in recent years. In this study, hazard zonation was carried out, incorporating information about debris flow initiation probability (spatial and temporal), and the delimitation of the potential runout areas. Debris flow hazard zonation was carried out in the area of the Consortium of Mountain Municipalities of Valtellina di Tirano (Central Alps, Italy). The complexity of the phenomenon, the scale of the study, the variability of local conditioning factors, and the lacking data limited the use of process-based models for the runout zone delimitation. Firstly, a map of hazard initiation probabilities was prepared for the study area, based on the available susceptibility zoning information, and the analysis of two sets of aerial photographs for the temporal probability estimation. Afterwards, the hazard initiation map was used as one of the inputs for an empirical GIS-based model (Flow-R), developed at the University of Lausanne (Switzerland). An estimation of the debris flow magnitude was neglected as the main aim of the analysis was to prepare a debris flow hazard map at medium scale. A digital elevation model, with a 10 m resolution, was used together with landuse, geology and debris flow hazard initiation maps as inputs of the Flow-R model to restrict potential areas within each hazard initiation probability class to locations where debris flows are most likely to initiate. Afterwards, runout areas were calculated using multiple flow direction and energy based algorithms. Maximum probable runout zones were calibrated using documented past events and aerial photographs. Finally, two debris flow hazard maps were prepared. The first simply delimits five hazard zones, while the second incorporates the information about debris flow spreading direction probabilities, showing areas more likely to be affected by future debris flows. Limitations of the modelling arise mainly from the models applied and analysis scale, which are neglecting local controlling factors of debris flow hazard. The presented approach of debris flow hazard analysis, associating automatic detection of the source areas and a simple assessment of the debris flow spreading, provided results for consequent hazard and risk studies. However, for the validation and transferability of the parameters and results to other study areas, more testing is needed.

Flow-R, a model for susceptibility mapping of debris flows and other gravitational hazards at a regional scale: some case studies

Every year, debris flows cause huge damage in mountainous areas. Due to population pressure in hazardous zones, the socio-economic impact is much higher than in the past. Therefore, the development of indicative susceptibility hazard maps is of primary importance, particularly in developing countries. However, the complexity of the phenomenon and the variability of local controlling factors limit the use of processbased models for a first assessment. A debris flow model has been developed for regional susceptibility assessments using digital elevation model (DEM) with a GIS-based approach.. The automatic identification of source areas and the estimation of debris flow spreading, based on GIS tools, provide a substantial basis for a preliminary susceptibility assessment at a regional scale. One of the main advantages of this model is its workability. In fact, everything is open to the user, from the data choice to the selection of the algorithms and their parameters. The Flow-R model was tested in three different contexts: two in Switzerland and one in Pakistan, for indicative susceptibility hazard mapping. It was shown that the quality of the DEM is the most important parameter to obtain reliable results for propagation, but also to identify the potential debris flows sources.

Debris flow susceptibility mapping at a regional scale along the National Road N7, Argentina

A first assessment of debris flow susceptibility at a large scale was performed along the National Road N7, Argentina. Numerous catchments are prone to debris flows and likely to endanger the road-users. A 1:50,000 susceptibility map was created. The use of a DEM (grid 30 m) associated to three complementary criteria (slope, contributing area, curvature) allowed the identification of potential source areas. The debris flow spreading was estimated using a process- and GISbased model (Flow-R) based on basic probabilistic and energy calculations. The best-fit values for the coefficient of friction and the mass-to-drag ratio of the PCM model were found to be ? = 0.02 and M/D = 180 and the resulting propagation on one of the calibration site was validated using the Coulomb friction model. The results are realistic and will be useful to determine which areas need to be prioritized for detailed studies.

Combining palaeodistribution modelling and phylogeographical approaches for identifying glacial refugia in Alpine Primula

Aim We investigated the late Quaternary history of two closely related and partly sympatric species of Primula from the south-western European Alps, P. latifolia Lapeyr. and P. marginata Curtis, by combining phylogeographical and palaeodistribution modelling approaches. In particular, we were interested in whether the two approaches were congruent and identified the same glacial refugia. Location South-western European Alps. Methods For the phylogeographical analysis we included 353 individuals from 28 populations of P. marginata and 172 individuals from 15 populations of P. latifolia and used amplified fragment length polymorphisms (AFLPs). For palaeodistribution modelling, species distribution models (SDMs) were based on extant species occurrences and then projected to climate models (CCSM, MIROC) of the Last Glacial Maximum (LGM), approximately 21 ka. Results The locations of the modelled LGM refugia were confirmed by various indices of genetic variation. The refugia of the two species were largely geographically isolated, overlapping only 6% to 11% of the species' total LGM distribution. This overlap decreased when the position of the glacial ice sheet and the differential elevational and edaphic distributions of the two species were considered. Main conclusions The combination of phylogeography and palaeodistribution modelling proved useful in locating putative glacial refugia of two alpine species of Primula. The phylogeographical data allowed us to identify those parts of the modelled LGM refugial area that were likely source areas for recolonization. The use of SDMs predicted LGM refugial areas substantially larger and geographically more divergent than could have been predicted by phylogeographical data alone

Technical note: Preliminary estimation of rockfall runout zones

Rockfall propagation areas can be determined using a simple geometric rule known as shadow angle or energy line method based on a simple Coulomb frictional model implemented in the CONEFALL computer program. Runout zones are estimated from a digital terrain model (DTM) and a grid file containing the cells representing rockfall potential source areas. The cells of the DTM that are lowest in altitude and located within a cone centered on a rockfall source cell belong to the potential propagation area associated with that grid cell. In addition, the CONEFALL method allows estimation of mean and maximum velocities and energies of blocks in the rockfall propagation areas. Previous studies indicate that the slope angle cone ranges from 27° to 37° depending on the assumptions made, i.e. slope morphology, probability of reaching a point, maximum run-out, field observations. Different solutions based on previous work and an example of an actual rockfall event are presented here.

Structural characterization of rockfall sources in Yosemite Valley from remote sensing data and field surveys

Yosemite Valley poses significant rockfall hazard and related risk due to its glacially steepened walls and approximately 4 million visitors annually. To assess rockfall hazard, it is necessary to evaluate the geologic structure that contributes to the destabilization of rockfall sources and locate the most probable future source areas. Coupling new remote sensing techniques (Terrestrial Laser Scanning, Aerial Laser Scanning) and traditional field surveys, we investigated the regional geologic and structural setting, the orientation of the primary discontinuity sets for large areas of Yosemite Valley, and the specific discontinuity sets present at active rockfall sources. This information, combined with better understanding of the geologic processes that contribute to the progressive destabilization and triggering of granitic rock slabs, contributes to a more accurate rockfall susceptibility assessment for Yosemite Valley and elsewhere.

Brief communication: Report on the impact of the 27 February 2010 earthquake (Chile, Mw 8.8) on rockfalls in the Las Cuevas valley, Argentina

Numerous rockfalls were detected in the Las Cuevas valley, Argentina, after the 27 February 2010 earthquake in Chile. Live rockfalls were observed during aftershocks of 11 March 2010. Many rockfall source areas coincide with known thrust fault and some areas presented a rockfall activity even after the tremors. Some rockfalls crossed the National Road 7 but no damages to houses or vehicles were reported. This study illustrates how the 27 February 2010 earthquake impacted on unstable slopes in a valley far from the earthquakes epicentre. It is an interesting addition to previous studies on landslides caused by earthquakes because of the high magnitude of the event and of its aftershocks.

Budgeting rockfall and modeling sedimentary delivery in torrent systems

This thesis is a compilation of projects to study sediment processes recharging debris flow channels. These works, conducted during my stay at the University of Lausanne, focus in the geological and morphological implications of torrent catchments to characterize debris supply, a fundamental element to predict debris flows. Other aspects of sediment dynamics are considered, e.g. the coupling headwaters - torrent, as well as the development of a modeling software that simulates sediment transfer in torrent systems. The sediment activity at Manival, an active torrent system of the northern French Alps, was investigated using terrestrial laser scanning and supplemented with geostructural investigations and a survey of sediment transferred in the main torrent. A full year of sediment flux could be observed, which coincided with two debris flows and several bedload transport events. This study revealed that both debris flows generated in the torrent and were preceded in time by recharge of material from the headwaters. Debris production occurred mostly during winter - early spring time and was caused by large slope failures. Sediment transfers were more puzzling, occurring almost exclusively in early spring subordinated to runoffconditions and in autumn during long rainfall. Intense rainstorms in summer did not affect debris storage that seems to rely on the stability of debris deposits. The morpho-geological implication in debris supply was evaluated using DEM and field surveys. A slope angle-based classification of topography could characterize the mode of debris production and transfer. A slope stability analysis derived from the structures in rock mass could assess susceptibility to failure. The modeled rockfall source areas included more than 97% of the recorded events and the sediment budgets appeared to be correlated to the density of potential slope failure. This work showed that the analysis of process-related terrain morphology and of susceptibility to slope failure document the sediment dynamics to quantitatively assess erosion zones leading to debris flow activity. The development of erosional landforms was evaluated by analyzing their geometry with the orientations of potential rock slope failure and with the direction of the maximum joint frequency. Structure in rock mass, but in particular wedge failure and the dominant discontinuities, appear as a first-order control of erosional mechanisms affecting bedrock- dominated catchment. They represent some weaknesses that are exploited primarily by mass wasting processes and erosion, promoting not only the initiation of rock couloirs and gullies, but also their propagation. Incorporating the geological control in geomorphic processes contributes to better understand the landscape evolution of active catchments. A sediment flux algorithm was implemented in a sediment cascade model that discretizes the torrent catchment in channel reaches and individual process-response systems. Each conceptual element includes in simple manner geomorphological and sediment flux information derived from GIS complemented with field mapping. This tool enables to simulate sediment transfers in channels considering evolving debris supply and conveyance, and helps reducing the uncertainty inherent to sediment budget prediction in torrent systems. Cette thèse est un recueil de projets d'études des processus de recharges sédimentaires des chenaux torrentiels. Ces travaux, réalisés lorsque j'étais employé à l'Université de Lausanne, se concentrent sur les implications géologiques et morphologiques des bassins dans l'apport de sédiments, élément fondamental dans la prédiction de laves torrentielles. D'autres aspects de dynamique sédimentaire ont été abordés, p. ex. le couplage torrent - bassin, ainsi qu'un modèle de simulation du transfert sédimentaire en milieu torrentiel. L'activité sédimentaire du Manival, un système torrentiel actif des Alpes françaises, a été étudiée par relevés au laser scanner terrestre et complétée par une étude géostructurale ainsi qu'un suivi du transfert en sédiments du torrent. Une année de flux sédimentaire a pu être observée, coïncidant avec deux laves torrentielles et plusieurs phénomènes de charriages. Cette étude a révélé que les laves s'étaient générées dans le torrent et étaient précédées par une recharge de débris depuis les versants. La production de débris s'est passée principalement en l'hiver - début du printemps, causée par de grandes ruptures de pentes. Le transfert était plus étrange, se produisant presque exclusivement au début du printemps subordonné aux conditions d'écoulement et en automne lors de longues pluies. Les orages d'été n'affectèrent guère les dépôts, qui semblent dépendre de leur stabilité. Les implications morpho-géologiques dans l'apport sédimentaire ont été évaluées à l'aide de MNT et études de terrain. Une classification de la topographie basée sur la pente a permis de charactériser le mode de production et transfert. Une analyse de stabilité de pente à partir des structures de roches a permis d'estimer la susceptibilité à la rupture. Les zones sources modélisées comprennent plus de 97% des chutes de blocs observées et les bilans sédimentaires sont corrélés à la densité de ruptures potentielles. Ce travail d'analyses des morphologies du terrain et de susceptibilité à la rupture documente la dynamique sédimentaire pour l'estimation quantitative des zones érosives induisant l'activité torrentielle. Le développement des formes d'érosion a été évalué par l'analyse de leur géométrie avec celle des ruptures potentielles et avec la direction de la fréquence maximale des joints. Les structures de roches, mais en particulier les dièdres et les discontinuités dominantes, semblent être très influents dans les mécanismes d'érosion affectant les bassins rocheux. Ils représentent des zones de faiblesse exploitées en priorité par les processus de démantèlement et d'érosion, encourageant l'initiation de ravines et couloirs, mais aussi leur propagation. L'incorporation du control géologique dans les processus de surface contribue à une meilleure compréhension de l'évolution topographique de bassins actifs. Un algorithme de flux sédimentaire a été implémenté dans un modèle en cascade, lequel divise le bassin en biefs et en systèmes individuels répondant aux processus. Chaque unité inclut de façon simple les informations géomorpologiques et celles du flux sédimentaire dérivées à partir de SIG et de cartographie de terrain. Cet outil permet la simulation des transferts de masse dans les chenaux, considérants la variabilité de l'apport et son transport, et aide à réduire l'incertitude liée à la prédiction de bilans sédimentaires torrentiels. Ce travail vise très humblement d'éclairer quelques aspects de la dynamique sédimentaire en milieu torrentiel.

Debris flow modeling for susceptibility mapping at regional to national scale in Norway

Debris flows and related landslide processes occur in many regions all over Norway and pose a significant hazard to inhabited areas. Within the framework of the development of a national debris flows susceptibility map, we are working on a modeling approach suitable for Norway with a nationwide coverage. The discrimination of source areas is based on an index approach, which includes topographic parameters and hydrological settings. For the runout modeling, we use the Flow-R model (IGAR, University of Lausanne), which is based on combined probabilistic and energetic algorithms for the assessment of the spreading of the flow and maximum runout distances. First results for different test areas have shown that runout distances can be modeled reliably. For the selection of source areas, however, additional factors have to be considered, such as the lithological and quaternary geological setting, in order to accommodate the strong variation in debris flow activity in the different geological, geomorphological and climate regions of Norway.

08004-001 Summit Lake Final Report

The Summit Lake Watershed Improvement Project is a watershed-based sediment control project designed to greatly reduce to nearly eliminate sedimentation of an existing lake that is being renovated for use as a water source in southern Iowa. Summit Lake is owned by the City of Creston and was once a water source lake until around 1984. The watershed improvements will include lakeshore stabilization and erosion control practices as a precursor for related improvements to the lake and overall 4,900-acre watershed. Best practices included in this phase are the implementation of riprap, a rain garden, grade stabilization structures, grassed waterways, terraces, basins, water use and access ordinances, education and outreach, water monitoring, and other stream bank improvements. These improvements, along with leveraged work to be done by strategic partners, will enable the lake to be used for local and regional water supplies by sustaining the lake for many years to come. Without the lake rehabilitation, the lake will likely be filled with sedimentation to the point that it will have no recreational value. Key partners are the City of Creston, IDNR, Southern Iowa Rural Water Association, Union County, the Union County NRCS office, Southwestern Community College, and the Summit Lake Association, which is a non-profit group of landowners working to protect the lake. The project will address WIRB targets: a) streambank stabilization, b) livestock runoff, c) agricultural runoff and drainage, d) stormwater runoff, and e) a section of inadequately sewered community.

Geological structure, recharge processes and underground drainage of a glacierised karst aquifer system, Tsanfleuron-Sanetsch, Swiss Alps

The relationships between stratigraphic and tectonic setting, recharge processes and underground drainage of the glacierised karst aquifer system `Tsanfleuron-Sanetsch' in the Swiss Alps have been studied by means of various methods, particularly tracer tests (19 injections). The area belongs to the Helvetic nappes and consists of Jurassic to Palaeogene sedimentary rocks. Strata are folded and form a regional anticlinorium. Cretaceous Urgonian limestone constitutes the main karst aquifer, overlain by a retreating glacier in its upper part. Polished limestone surfaces are exposed between the glacier front and the end moraine of 1855/1860 (Little Ice Age); typical alpine karrenfields can be observed further below. Results show that (1) large parts of the area are drained by the Glarey spring, which is used as a drinking water source, while marginal parts belong to the catchments of other springs; (2) groundwater flow towards the Glarey spring occurs in the main aquifer, parallel to stratification, while flow towards another spring crosses the entire stratigraphic sequence, consisting of about 800 m of marl and limestone, along deep faults that were probably enlarged by mass movements; (3) the variability of glacial meltwater production influences the shape of the tracer breakthrough curves and, consequently, flow and transport in the aquifer.

Climate and sea-level variations along the northwestern Tethyan margin during the Valanginian C-isotope excursion: Mineralogical evidence from the Vocontian Basin (SE France)

A high resolution mineralogical study (bulk-rock and clay-fraction) was carried out upon the hemipelagic strata of the Angles section (Vocontian Basin, SE France) in which the Valanginian positive C-isotope excursion occurs. To investigate sea-level fluctuations and climate change respectively, a Detrital Index (DI: (phyllosilicates and quartz)/calcite) and a Weathering Index (WI: kaolinite/(illite + chlorite)) were established and compared to second-order sea-level fluctuations. In addition, the mineralogical data were compared with the High Nutrient Index (HNI, based on calcareous nannofossil taxa) data obtained by Duchamp-Alphonse et al. (2007), in order to assess the link between the hydrolysis conditions recorded on the surrounding continents and the trophic conditions inferred for the Vocontian Basin. It appears that the mineralogical distribution along the northwestern Tethyan margin is mainly influenced by sea-level changes during the Early Valanginian (Pertransiens to Stephanophorus ammonite Zones) and by climate variations from the late Early Valanginian to the base of the Hauterivian (top of the Stephanophorus to the Radiatus ammonite Zones). The sea-level fall observed in the Pertransiens ammonite Zone (Early Valanginian) is well expressed by an increase in detrital inputs (an increase in the DI) associated with a more proximal source and a shallower marine environment, whereas the sea-level rise recorded in the Stephanophorus ammonite Zone corresponds to a decrease in detrital influx (a decrease in the DI) as the source becomes more distal and the environment deeper. Interpretation of both DI and WI, indicates that the positive C-isotope excursion (top of the Stephanophorus to the Verrucosum ammonite Zones) is associated with an increase of detrital inputs under a stable, warm and humid climate, probably related to greenhouse conditions, the strongest hydrolysis conditions being reached at the maximum of the positive C-isotope excursion. From the Verrucosum ammonite Zone to the base of the Hauterivian (Radiatus ammonite Zone) climatic conditions evolved from weak hydrolysis conditions and, most likely, a cooler climate (resulting in a decrease in detrital inputs) to a seasonal climate in which more humid seasons alternated with more arid ones. The comparison of the WI to the HNI shows that the nutrification recorded al: the Angles section from the top of the Stephanophorus to the Radiatus ammonite Zones (including the positive C-isotope shift), is associated with climatic changes in the source areas. At that time, increased nutrient inputs were generally triggered by increased weathering processes in the source areas due to acceleration in the hydrological cycle under greenhouse conditions This scenario accords with the widely questioned palaeoenvironmental model proposed by Lini et al., (1992) and suggests that increasing greenhouse conditions are the main factor that drove the palaeoenvironmental changes observed in the hemipelagic realm of the Vocontian Basin, during the Valanginian positive C-isotope shift. This high-resolution mineralogical study highlights short-term climatic changes during the Valanginian, probably associated to rapid changes in the C-cycle. Coeval Massive Parana-Etendeka flood basalt eruptions may explain such rapid perturbations. (C) 2011 Elsevier B.V. All rights reserved.

Growth and demise of Permian biogenic chert along northwest Pangea: evidence for end-Permian collapse of thermohaline circulation

The Permian Chert Event (PCE) was a 30 Ma long episode of unusual chert accumulation along the northwest margin of Pangea, and possibly worldwide. The onset of the PCE occurred at about the Sakmarian-Artinskian boundary in the Sverdrup Basin, Canadian Arctic, where it coincides with a maximum flooding event, the ending of high-frequency/high-amplitude shelf cyclicity, the onset of massive biogenic chert deposition in deep-water distal areas, and a long-term shift from warm- to cool-water carbonate sedimentation in shallow-water proximal areas. A similar and coeval shift is observed from the Barents Sea to the northwestern USA. A landward and southward expansion of silica factories occurred during the Middle and Late Permian at which time warm-water carbonate producers disappeared completely from the northwest margin of Pangea. Biotically impoverished and increasingly narrow cold-water carbonate factories (characterised by non-cemented bioclasts of sponges, bryozoans, echinoderms and brachiopods) were then progressively replaced by silica factories. By Late Permian time, little carbonate sediments accumulated in the Barents Sea and in the Sverdrup Basin. where the deep- to shallow-water sedimentary spectrum was occupied by siliceous sponge spicules. By that time, biogenic silica sedimentation was common throughout the world. Silica factories collapsed in the Late Permian, abruptly bringing the PCE to an end. In northwest Pangea, the end- Permian collapse of the PCE was associated with a major transgression and with a return to much warmer oceanic and continental climatic conditions. Chert deposition resumed in the distal oceanic areas during the early Middle Triassic (Anisian) after a 8-10 Ma interruption (Early Triassic Chert Gap). The conditions necessary for the onset, expansion and zenith of the PCE were provided by the thermohaline circulation of nutrient-rich cold waters along the northwestern and western margin of Pangea, and possibly throughout the world oceans. These conditions provided an efficient transportation mechanism that constantly replenished the supply of silica in the area, created a nutrient- and oxygen-rich environment favouring siliceous biogenic productivity. established cold sea-floor conditions, hindering silica dissolution, while increasing calcium carbonate solubility, and provided conditions adverse to organic and inorganic carbonate production, The northwest margin of Pangea was, for nearly 30 Ma. bathed by cold waters presumably derived from the seasonal melting of northern sea ice, the assumed engine for thermohaline circulation. This process started near the Sakmarian-Artinskian boundary. intensified throughout Middle and Late Permian time and ceased suddenly in latest Permian time, It led to oceanic conditions much colder than normally expected from the palaeolatitudes. and the influence of cold northerly-derived water was felt as far south southern Nevada. The demise of silica factories was caused by the rapid breakdown of these conditions and the establishment of a much warmer marine environment accompanied by sluggish circulation and perhaps a reduced input of dissolved silica to the ocean. Complete thawing of northern sea ice would have ended thermohaline circulation and led to warm and sluggish oceanic conditions inimical to the production. accumulation and preservation of biogenic silica.

Community replacement of neritic carbonate organisms during the late Valanginian platform demise: A new record from the Provence Platform

The Valanginian is marked by a major platform demise inducing a hiatus in the northern Tethyan neritic carbonate record from the top of the lower Valanginian to the lower Hauterivian. New biostratigraphic and chemostratigraphic data from the Ollioules section (Provence Platform, southern France) are presented here, demonstrating that a large part of the upper Valanginian is preserved in an inner platform environment. The thick, upper Valanginian, aggrading carbonate succession is observed in an aborted rift domain, implying relatively low subsidence. In this context, a relatively long-term sea-level rise was required to sustain a keep-up style of carbonate production. Like the Apulian Platform, the remarkable preservation of the Provence Platform may have been favored by its remoteness from terrigenous source areas, as suggested by the low clastic inputs and low P-accumulation rates. Two main biotic community replacements are observed in Ollioules. The first saw the development of abundant microbialites and algae at the onset of the late Valanginian. A Tubiphytes concentration occurred during the coolest climatic conditions and the transition towards arid conditions, whereas the subsequent Lithocodium-Bacinella and orbitolinids assemblages developed under low nutrient conditions during a warmer interval. Both assemblages may have been triggered by increased alkalinity. The second community replacement saw the installation of coral- and rudist-dominated communities during the latest Valanginian to early Hauterivian. They indicate a change to oligotrophic, open marine conditions. Six medium-scale sequences have been defined in Ollioules, indicating short-term transgressive-regressive trends superimposed on a long-term transgression. Low nutrient inputs and relatively low subsidence in an aggradational context may explain the survival of the isolated Provence Carbonate Platform during a time of widespread drowning episodes and platform demise in the northern Tethyan domain. (c) 2012 Elsevier B.V. All rights reserved.