Spatial distribution of runoff depth from a watershed located in a cuesta relief area of São Paulo State, Brazil

This study was undertaken in a 1566 ha drainage basin situated in an area with cuesta relief in the state of São Paulo, Brazil. The objectives were: 1) to map the maximum potential soil water retention capacity, and 2) to simulate the depth of surface runoff in each geographical position of the area based on a typical rainfall event. The database required for the development of this research was generated in the environment of the geographical information system ArcInfo v.10.1. Undeformed soil samples were collected at 69 points. The ordinary kriging method was used in the interpolation of the values of soil density and maximum potential soil water retention capacity. The spherical model allowed for better adjustment of the semivariograms corresponding to the two soil attributes for the depth of 0 to 20 cm, while the Gaussian model enabled a better fit of the spatial behavior of the two variables for the depth of 20 to 40 cm. The simulation of the spatial distribution revealed a gradual increase in the depth of surface runoff for the rainfall event taken as example (25 mm) from the reverse to the peripheral depression of the cuesta (from west to east). There is a positive aspect observed in the gradient, since the sites of highest declivity, especially those at the front of the cuesta, are closer to the western boundary of the watershed where the lowest depths of runoff occur. This behavior, in conjunction with certain values of erodibility and depending on the land use and cover, can help mitigate the soil erosion processes in these areas.

Hydrological processes obtained on the plot scale under four simulated rainfall tests during the cycle of different crop systems

The cropping system influences the interception of water by plants, water storage in depressions on the soil surface, water infiltration into the soil and runoff. The aim of this study was to quantify some hydrological processes under no tillage cropping systems at the edge of a slope, in 2009 and 2010, in a Humic Dystrudept soil, with the following treatments: corn, soybeans, and common beans alone; and intercropped corn and common bean. Treatments consisted of four simulated rainfall tests at different times, with a planned intensity of 64 mm h-1 and 90 min duration. The first test was applied 18 days after sowing, and the others at 39, 75 and 120 days after the first test. Different times of the simulated rainfall and stages of the crop cycle affected soil water content prior to the rain, and the time runoff began and its peak flow and, thus, the surface hydrological processes. The depth of the runoff and the depth of the water intercepted by the crop + soil infiltration + soil surface storage were affected by the crop systems and the rainfall applied at different times. The corn crop was the most effective treatment for controlling runoff, with a water loss ratio of 0.38, equivalent to 75 % of the water loss ratio exhibited by common bean (0.51), the least effective treatment in relation to the others. Total water loss by runoff decreased linearly with an increase in the time that runoff began, regardless of the treatment; however, soil water content on the gravimetric basis increased linearly from the beginning to the end of the rainfall.

Glacier mass balances of Stubacher Sonnblickkees, Hohe Tauern Range, Eastern Alps, Austria, 1958/1959 to 2012/2013

Stubacher Sonnblickkees (SSK) is located in the Hohe Tauern Range (Eastern Alps) in the south of Salzburg Province (Austria) in the region of Oberpinzgau in the upper Stubach Valley. The glacier is situated at the main Alpine crest and faces east, starting at elevations close to 3050 m and in the 1980s terminated at 2500 m a.s.l. It had an area of 1.7 km² at that time, compared with 1 km² in 2013. The glacier type can be classified as a slope glacier, i.e. the relief is covered by a relatively thin ice sheet and there is no regular glacier tongue. The rough subglacial topography makes for a complex shape in the surface topography, with various concave and convex patterns. The main reason for selecting this glacier for mass balance observations (as early as 1963) was to verify on a complex glacier how the mass balance methods and the conclusions - derived during the more or less pioneer phase of glaciological investigations in the 1950s and 1960s - could be applied to the SSK glacier. The decision was influenced by the fact that close to the SSK there was the Rudolfshütte, a hostel of the Austrian Alpine Club (OeAV), newly constructed in the 1950s to replace the old hut dating from 1874. The new Alpenhotel Rudolfshütte, which was run by the Slupetzky family from 1958 to 1970, was the base station for the long-term observation; the cable car to Rudolfshütte, operated by the Austrian Federal Railways (ÖBB), was a logistic advantage. Another factor for choosing SSK as a glaciological research site was the availability of discharge records of the catchment area from the Austrian Federal Railways who had turned the nearby lake Weißsee ('White Lake') - a former natural lake - into a reservoir for their hydroelectric power plants. In terms of regional climatic differences between the Central Alps in Tyrol and those of the Hohe Tauern, the latter experienced significantly higher precipitation , so one could expect new insights in the different response of the two glaciers SSK and Hintereisferner (Ötztal Alps) - where a mass balance series went back to 1952. In 1966 another mass balance series with an additional focus on runoff recordings was initiated at Vernagtfener, near Hintereisferner, by the Commission of the Bavarian Academy of Sciences in Munich. The usual and necessary link to climate and climate change was given by a newly founded weather station (by Heinz and Werner Slupetzky) at the Rudolfshütte in 1961, which ran until 1967. Along with an extension and enlargement to the so-called Alpine Center Rudolfshütte of the OeAV, a climate observatory (suggested by Heinz Slupetzky) has been operating without interruption since 1980 under the responsibility of ZAMG and the Hydrological Service of Salzburg, providing long-term met observations. The weather station is supported by the Berghotel Rudolfshütte (in 2004 the OeAV sold the hotel to a private owner) with accommodation and facilities. Direct yearly mass balance measurements were started in 1963, first for 3 years as part of a thesis project. In 1965 the project was incorporated into the Austrian glacier measurement sites within the International Hydrological Decade (IHD) 1965 - 1974 and was afterwards extended via the International Hydrological Program (IHP) 1975 - 1981. During both periods the main financial support came from the Hydrological Survey of Austria. After 1981 funds were provided by the Hydrological Service of the Federal Government of Salzburg. The research was conducted from 1965 onwards by Heinz Slupetzky from the (former) Department of Geography of the University of Salzburg. These activities received better recognition when the High Alpine Research Station of the University of Salzburg was founded in 1982 and brought in additional funding from the University. With recent changes concerning Rudolfshütte, however, it became unfeasible to keep the research station going. Fortunately, at least the weather station at Rudolfshütte is still operating. In the pioneer years of the mass balance recordings at SSK, the main goal was to understand the influence of the complicated topography on the ablation and accumulation processes. With frequent strong southerly winds (foehn) on the one hand, and precipitation coming in with storms from the north to northwest, the snow drift is an important factor on the undulating glacier surface. This results in less snow cover in convex zones and in more or a maximum accumulation in concave or flat areas. As a consequence of the accentuated topography, certain characteristic ablation and accumulation patterns can be observed during the summer season every year, which have been regularly observed for many decades . The process of snow depletion (Ausaperung) runs through a series of stages (described by the AAR) every year. The sequence of stages until the end of the ablation season depends on the weather conditions in a balance year. One needs a strong negative mass balance year at the beginning of glacier measurements to find out the regularities; 1965, the second year of observation resulted in a very positive mass balance with very little ablation but heavy accumulation. To date it is the year with the absolute maximum positive balance in the entire mass balance series since 1959, probably since 1950. The highly complex ablation patterns required a high number of ablation stakes at the beginning of the research and it took several years to develop a clearer idea of the necessary density of measurement points to ensure high accuracy. A great number of snow pits and probing profiles (and additional measurements at crevasses) were necessary to map the accumulation area/patterns. Mapping the snow depletion, especially at the end of the ablation season, which coincides with the equilibrium line, is one of the main basic data for drawing contour lines of mass balance and to calculate the total mass balance (on a regular-shaped valley glacier there might be an equilibrium line following a contour line of elevation separating the accumulation area and the ablation area, but not at SSK). - An example: in 1969/70, 54 ablation stakes and 22 snow pits were used on the 1.77 km² glacier surface. In the course of the study the consistency of the accumulation and ablation patterns could be used to reduce the number of measurement points. - At the SSK the stratigraphic system, i.e. the natural balance year, is used instead the usual hydrological year. From 1964 to 1981, the yearly mass balance was calculated by direct measurements. Based on these records of 17 years, a regression analysis between the specific net mass balance and the ratio of ablation area to total area (AAR) has been used since then. The basic requirement was mapping the maximum snow depletion at the end of each balance year. There was the advantage of Heinz Slupetzky's detailed local and long-term experience, which ensured homogeneity of the series on individual influences of the mass balance calculations. Verifications took place as often as possible by means of independent geodetic methods, i.e. monoplotting , aerial and terrestrial photogrammetry, more recently also the application of PHOTOMODELLER and laser scans. The semi-direct mass balance determinations used at SSK were tentatively compared with data from periods of mass/volume change, resulting in promising first results on the reliability of the method. In recent years re-analyses of the mass balance series have been conducted by the World Glacier Monitoring Service and will be done at SSK too. - The methods developed at SSK also add to another objective, much discussed in the 1960s within the community, namely to achieve time- and labour-saving methods to ensure continuation of long-term mass balance series. The regression relations were used to extrapolate the mass balance series back to 1959, the maximum depletion could be reconstructed by means of photographs for those years. R. Günther (1982) calculated the mass balance series of SSK back to 1950 by analysing the correlation between meteorological data and the mass balance; he found a high statistical relation between measured and determined mass balance figures for SSK. In spite of the complex glacier topography, interesting empirical experiences were gained from the mass balance data sets, giving a better understanding of the characteristics of the glacier type, mass balance and mass exchange. It turned out that there are distinct relations between the specific net balance, net accumulation (defined as Bc/S) and net ablation (Ba/S) to the AAR, resulting in characteristic so-called 'turnover curves'. The diagram of SSK represents the type of a glacier without a glacier tongue. Between 1964 and 1966, a basic method was developed, starting from the idea that instead of measuring years to cover the range between extreme positive and extreme negative yearly balances one could record the AAR/snow depletion/Ausaperung during one or two summers. The new method was applied on Cathedral Massif Glacier, a cirque glacier with the same area as the Stubacher Sonnblickkees, in British Columbia, Canada. during the summers of 1977 and 1978. It returned exactly the expected relations, e.g. mass turnover curves, as found on SSK. The SSK was mapped several times on a scale of 1:5000 to 1:10000. Length variations have been measured since 1960 within the OeAV glacier length measurement programme. Between 1965 and 1981, there was a mass gain of 10 million cubic metres. With a time lag of 10 years, this resulted in an advance until the mid-1980s. Since 1982 there has been a distinct mass loss of 35 million cubic metres by 2013. In recent years, the glacier has disintegrated faster, forced by the formation of a periglacial lake at the glacier terminus and also by the outcrops of rocks (typical for the slope glacier type), which have accelerated the meltdown. The formation of this lake is well documented. The glacier has retreated by some 600 m since 1981. - Since August 2002, a runoff gauge installed by the Hydrographical Service of Salzburg has recorded the discharge of the main part of SSK at the outlet of the new Unterer Eisboden See. The annual reports - submitted from 1982 on as a contractual obligation to the Hydrological Service of Salzburg - document the ongoing processes on the one hand, and emphasize the mass balance of SSK and outline the climatological reasons, mainly based on the met-data of the observatory Rudolfshütte, on the other. There is an additional focus on estimating the annual water balance in the catchment area of the lake. There are certain preconditions for the water balance equation in the area. Runoff is recorded by the ÖBB power stations, the mass balance of the now approx. 20% glaciated area (mainly the Sonnblickkees) is measured andthe change of the snow and firn patches/the water content is estimated as well as possible. (Nowadays laserscanning and ground radar are available to measure the snow pack). There is a net of three precipitation gauges plus the recordings at Rudolfshütte. The evaporation is of minor importance. The long-term annual mean runoff depth in the catchment area is around 3.000 mm/year. The precipitation gauges have measured deficits between 10% and 35%, on average probably 25% to 30%. That means that the real precipitation in the catchment area Weißsee (at elevations between 2,250 and 3,000 m) is in an order of 3,200 to 3,400 mm a year. The mass balance record of SSK was the first one established in the Hohe Tauern region (and now since the Hohe Tauern National Park was founded in 1983 in Salzburg) and is one of the longest measurement series worldwide. Great efforts are under way to continue the series, to safeguard against interruption and to guarantee a long-term monitoring of the mass balance and volume change of SSK (until the glacier is completely gone, which seems to be realistic in the near future as a result of the ongoing global warming). Heinz Slupetzky, March 2014

Runoff and sediment transport in a degraded area

Gully erosion occurs by the combined action of splash, sheetwash and rill-wash (interrill and rill erosion). These erosion processes have a great capacity for both sediment production and sediment transport. The objectives of this experiment were to evaluate hydrological and sediment transport in a degraded area, severely dissected by gullies; to assess the hydraulic flow characteristics and their aggregate transport capacity; and to measure the initial splash erosion rate. In the study area in Guarapuava, State of Paraná, Brazil (lat 25º 24' S; long 51º24' W; 1034 m asl), the soil was classified as Cambissolo Húmico alumínico, with the following particle-size composition: sand 0.116 kg kg-1; silt 0.180 kg kg-1; and clay 0.704 kg kg-1. The approach of this research was based on microcatchments formed in the ground, to study the hydrological response and sediment transport. A total of eight rill systems were simulated with dry and wet soil. An average rainfall of 33.7 ± 4.0 mm was produced for 35 to 54 min by a rainfall simulator. The equipment was installed, and a trough was placed at the end of the rill to collect sediments and water. During the simulation, the following variables were measured: time to runoff, time to ponding, time of recession, flow velocity, depth, ratio of the initial splash and grain size. The rainsplash of dry topsoil was more than twice as high as under moist conditions (5 g m-2 min-1 and 2 g m-2 min-1, respectively). The characteristics of the flow hydraulics indicate transition from laminar to turbulent flow [Re (Reynolds number) 1000-2000]. In addition, it was observed that a flow velocity of 0.12 m s-1 was the threshold for turbulent flow (Re > 2000), especially at the end of the rainfall simulation. The rill flow tended to be subcritical [Fr (Froude Number) < 1.0]. The variation in hydrological attributes (infiltration and runoff) was lower, while the sediment yield was variable. The erosion in the rill systems was characterized as limited transport, although the degraded area generated an average of 394 g m-2 of sediment in each simulation.

Available content, surface runoff and leaching of phosphorus forms in a typic hapludalf treated with organic and mineral nutrient sources

The application of animal manure to soil can increase phosphorus availability to plants and enhance transfer of the nutrient solution drained from the soil surface or leached into the soil profile. The aim of this study was to evaluate the effect of successive applications of organic and mineral nutrient sources on the available content, surface runoff and leaching of P forms in a Typic Hapludalf in no-tillage systems. Experiment 1 was set up in 2004 in the experimental area of UFSM, in Santa Maria (RS, Brazil). The treatments consisted of: control (without nutrient application) and application of pig slurry (PS), pig deep-litter (PL), cattle slurry (CS), and mineral fertilizers (NPK). The rates were determined to meet the N crop requirements of no-tillage black oat and maize, grown in the 2010/2011 growing season. The soil solution was collected after each event (rain + runoff or leaching) and the soluble, particulate and total P contents were measured. In November 2008, soil was collected in 2 cm intervals to a depth of 20 cm, in 5 cm intervals to a depth of 40 cm, and in 10 cm intervals to a depth of 70 cm. The soil was dried and ground, and P determined after extraction by anion exchange resin (AER). In experiment 2, samples collected from the Typic Hapludalf near experiment 1 were incubated for 20, 35, 58, 73 and 123 days after applying the following treatments: soil, soil + PS, soil + PL, soil + CS and soil + NPK. Thereafter, the soil was sampled and P was analyzed by AER. The applications of nutrient sources over the years led to an increase in available P and its migration in the soil profile. This led to P transfer via surface runoff and leaching, with the largest transfer being observed in PS and PL treatments, in which most P was applied. The soil available P and P transfer via surface runoff were correlated with the amounts applied, regardless of the P source. However, P transfer by leaching was not correlated with the applied nutrient amount, but rather with the solution amount leached in the soil profile.

Extension And Testing Of A 2D Hydrodynamic Model For Direct Rainfall Runoff Simulation

While the simulation of flood risks originating from the overtopping of river banks is well covered within continuously evaluated programs to improve flood protection measures, flash flooding is not. Flash floods are triggered by short, local thunderstorm cells with high precipitation intensities. Small catchments have short response times and flow paths and convective thunder cells may result in potential flooding of endangered settlements. Assessing local flooding and pathways of flood requires a detailed hydraulic simulation of the surface runoff. Hydrological models usually do not incorporate surface runoff at this detailedness but rather empirical equations are applied for runoff detention. In return 2D hydrodynamic models usually do not allow distributed rainfall as input nor are any types of soil/surface interaction implemented as in hydrological models. Considering several cases of local flash flooding during the last years the issue emerged for practical reasons but as well as research topics to closing the model gap between distributed rainfall and distributed runoff formation. Therefore, a 2D hydrodynamic model, depth-averaged flow equations using the finite volume discretization, was extended to accept direct rainfall enabling to simulate the associated runoff formation. The model itself is used as numerical engine, rainfall is introduced via the modification of waterlevels at fixed time intervals. The paper not only deals with the general application of the software, but intends to test the numerical stability and reliability of simulation results. The performed tests are made using different artificial as well as measured rainfall series as input. Key parameters of the simulation such as losses, roughness or time intervals for water level manipulations are tested regarding their impact on the stability.

Tracking water pathways in steep hillslopes by d18O depth profiles of soil water

Assessing temporal variations in soil water flow is important, especially at the hillslope scale, to identify mechanisms of runoff and flood generation and pathways for nutrients and pollutants in soils. While surface processes are well considered and parameterized, the assessment of subsurface processes at the hillslope scale is still challenging since measurement of hydrological pathways is connected to high efforts in time, money and personnel work. The latter might not even be possible in alpine environments with harsh winter processes. Soil water stable isotope profiles may offer a time-integrating fingerprint of subsurface water pathways. In this study, we investigated the suitability of soil water stable isotope (d18O) depth profiles to identify water flow paths along two transects of steep subalpine hillslopes in the Swiss Alps. We applied a one-dimensional advection–dispersion model using d18O values of precipitation (ranging from _24.7 to _2.9‰) as input data to simulate the d18O profiles of soil water. The variability of d18O values with depth within each soil profile and a comparison of the simulated and measured d18O profiles were used to infer information about subsurface hydrological pathways. The temporal pattern of d18O in precipitation was found in several profiles, ranging from _14.5 to _4.0‰. This suggests that vertical percolation plays an important role even at slope angles of up to 46_. Lateral subsurface flow and/or mixing of soil water at lower slope angles might occur in deeper soil layers and at sites near a small stream. The difference between several observed and simulated d18O profiles revealed spatially highly variable infiltration patterns during the snowmelt periods: The d18O value of snow (_17.7 ± 1.9‰) was absent in several measured d18O profiles but present in the respective simulated d18O profiles. This indicated overland flow and/or preferential flow through the soil profile during the melt period. The applied methods proved to be a fast and promising tool to obtain time-integrated information on soil water flow paths at the hillslope scale in steep subalpine slopes.

Multifilter rotating shadowband radiometer calibration for spectral aerosol optical depth retrievals over Sao Paulo City, Brazil

Multifilter rotating shadowband radiometer (MFRSR) calibration values for aerosol optical depth (AOD) retrievals were determined by means of the general method formulated by Forgan [Appl. Opt. 33, 4841 (1994)] at a polluted urban site. The obtained precision is comparable with the classical method, the Langley plot, applied on clean mountaintops distant of pollution sources. The AOD retrieved over Sao Paulo City with both calibration procedures is compared with the Aerosol Robotic Network data. The observed results are similar, and, except for the shortest wavelength (415 nm), the MFRSR`s AOD is systematically overestimated by similar to 0.03. (c) 2008 Optical Society of America.

Depth-dependent local structures in thin films unraveled by grazing-incidence X-ray absorption spectroscopy

A method of using X-ray absorption spectroscopy together with resolved grazing-incidence geometry for depth profiling of atomic, electronic or chemical local structures in thin films is presented. The quantitative deconvolution of thickness-dependent spectral features is performed by fully considering both scattering and absorption formalisms. Surface oxidation and local structural depth profiles in nanometric FePt films are determined, exemplifying the application of the method.

Measurement of the depth of maximum of extensive air showers above 10(18) eV

We describe the measurement of the depth of maximum, X(max), of the longitudinal development of air showers induced by cosmic rays. Almost 4000 events above 10(18) eV observed by the fluorescence detector of the Pierre Auger Observatory in coincidence with at least one surface detector station are selected for the analysis. The average shower maximum was found to evolve with energy at a rate of (106 +/- 35-21) g/cm(2)/decade below 10(18.24) +/- (0.05) eV, and d24 +/- 3 g/cm(2)/ecade above this energy. The measured shower-to-shower fluctuations decrease from about 55 to 26 g/cm(2). The interpretation of these results in terms of the cosmic ray mass composition is briefly discussed.

Runoff sources and land cover change in the Amazon: an end-member mixing analysis from small watersheds

The flowpaths by which water moves from watersheds to streams has important consequences for the runoff dynamics and biogeochemistry of surface waters in the Amazon Basin. The clearing of Amazon forest to cattle pasture has the potential to change runoff sources to streams by shifting runoff to more surficial flow pathways. We applied end-member mixing analysis (EMMA) to 10 small watersheds throughout the Amazon in which solute composition of streamwater and groundwater, overland flow, soil solution, throughfall and rainwater were measured, largely as part of the Large-Scale Biosphere-Atmosphere Experiment in Amazonia. We found a range in the extent to which streamwater samples fell within the mixing space determined by potential flowpath end-members, suggesting that some water sources to streams were not sampled. The contribution of overland flow as a source of stream flow was greater in pasture watersheds than in forest watersheds of comparable size. Increases in overland flow contribution to pasture streams ranged in some cases from 0% in forest to 27-28% in pasture and were broadly consistent with results from hydrometric sampling of Amazon forest and pasture watersheds that indicate 17- to 18-fold increase in the overland flow contribution to stream flow in pastures. In forest, overland flow was an important contribution to stream flow (45-57%) in ephemeral streams where flows were dominated by stormflow. Overland flow contribution to stream flow decreased in importance with increasing watershed area, from 21 to 57% in forest and 60-89% in pasture watersheds of less than 10 ha to 0% in forest and 27-28% in pastures in watersheds greater than 100 ha. Soil solution contributions to stream flow were similar across watershed area and groundwater inputs generally increased in proportion to decreases in overland flow. Application of EMMA across multiple watersheds indicated patterns across gradients of stream size and land cover that were consistent with patterns determined by detailed hydrometric sampling.

Node-Depth Encoding and Multiobjective Evolutionary Algorithm Applied to Large-Scale Distribution System Reconfiguration

The power loss reduction in distribution systems (DSs) is a nonlinear and multiobjective problem. Service restoration in DSs is even computationally hard since it additionally requires a solution in real-time. Both DS problems are computationally complex. For large-scale networks, the usual problem formulation has thousands of constraint equations. The node-depth encoding (NDE) enables a modeling of DSs problems that eliminates several constraint equations from the usual formulation, making the problem solution simpler. On the other hand, a multiobjective evolutionary algorithm (EA) based on subpopulation tables adequately models several objectives and constraints, enabling a better exploration of the search space. The combination of the multiobjective EA with NDE (MEAN) results in the proposed approach for solving DSs problems for large-scale networks. Simulation results have shown the MEAN is able to find adequate restoration plans for a real DS with 3860 buses and 632 switches in a running time of 0.68 s. Moreover, the MEAN has shown a sublinear running time in function of the system size. Tests with networks ranging from 632 to 5166 switches indicate that the MEAN can find network configurations corresponding to a power loss reduction of 27.64% for very large networks requiring relatively low running time.

Energy dissipation in depth-sensing indentation as a characteristic of the nanoscratch behavior of coatings

Wear behavior of coatings has usually been described in terms of mechanical properties such as hardness (H) and effective elastic modulus (E*). Alternatively, an energy approach appears as a promising analysis taking into account the influence of those properties. In a nanoindentation test, the dissipated energy depends not only on the hardness and elastic modulus, but also on the elastic recovery (W(e)). This work aims to establish a relation between plastic deformation energy (E(p)) during depth-sensing indentation method and the grooving resistance of coatings in nanoscratch tests. An energy dissipation coefficient (K(d)) was defined, calculated as the ratio of the plastic to the total deformation energy (E(p)/E(t)), which represents the energy dissipation of materials. Reactive depositions using titanium as the target and nitrogen and methane as reactive gases were obtained by triode magnetron sputtering, in order to assess wear and nanoindentation data. A topographical, chemical and microstructural characterization has been conducted using X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), wave dispersion spectroscopy (WDS), scanning electron (SEM) and atomic force microscopy (AFM) techniques. Nanoscratch results showed that the groove depth was well correlated to the energy dissipation coefficient of the coatings. On the other hand, a reduction in the coefficient was found when the elastic recovery was increased. (C) 2009 Elsevier B.V. All rights reserved.