Auditoria à qualidade do ar interior de um edifício de serviços - RSECE QAI

O objectivo principal deste trabalho é a realização de uma auditoria, à qualidade do ar interior (QAI), a um edifício de serviços – COCIGA, SA, tendo como base o Regulamento dos Sistemas Energéticos de Climatização dos Edifícios (RSECE). A auditoria QAI implica a medição de vários parâmetros físicos, químicos, microbiológicos e também a inspecção aos componentes do sistema de climatização com a finalidade de averiguar o seu estado de limpeza e manutenção. Assim, foram seleccionados 3 espaços, para a realização de amostragens designados por Comercial - Produtos, AVAC e Mezaninne das oficinas, nos quais foi efectuada a medição de diversos parâmetros, de acordo com as imposições do RSECE, utilizando medidores portáteis ou recorrendo a métodos analíticos. Relativamente aos parâmetros físicos, registaram-se valores de temperatura, para os três espaços estudados, entre os 21 e os 24 ºC e valores médios de humidade relativa de cerca de 50 %. Outro parâmetro medido, e de grande importância para garantir o conforto dos ocupantes, foi a velocidade do ar nos postos de trabalho. De acordo com o RSECE este valor não deve ser superior a 0,2 m/s, o que se verificou em todos os pontos medidos. O último parâmetro físico medido foi a concentração de partículas (PM10) tendo-se obtido valores de cerca de 23 μg/m3ar, valor bastante inferior ao máximo permitido pelo RSECE (150 μg/m3ar). Também no que diz respeito aos parâmetros químicos, ou seja, CO2, CO, formaldeído e ozono, não se verificaram valores superiores aos regulamentares. No caso do CO2, o valor máximo encontrado, nestes três espaços, foi de 745 ppm na Mezaninne das Oficinas e para o CO, na zona AVAC com uma concentração de 0,73 ppm. A medição do formaldeído registou valores perto dos 45 μg/m3ar e o ozono apenas foi detectado, em concentração muito reduzida, na zona Comercial – Produtos. Por fim, as concentrações de bactérias e fungos, de acordo com o RSECE, não devem ultrapassar as 500 UFC/m3ar (parâmetros microbiológicos). Em qualquer dos espaços, os valores medidos foram inferiores ao máximo legal, não ultrapassando as 50 UFC/m3ar. Da avaliação do projecto AVAC, e através da medição dos caudais de insuflação/ extracção em cada zona, concluiu-se que os seus valores não estão de acordo com os valores do projecto inicial que poderá ser imputada a uma insuficiência no funcionamento do sistema detectada na altura das medições. No que diz respeito ao estado de limpeza do sistema AVAC, apenas foi possível inspeccionar as unidades de tratamento de ar, tendo-se constatado que se encontram em boas condições. Ou seja, do ponto de vista do RSECE, e referindo-nos apenas à vertente da Qualidade do Ar Interior, o edifício em causa, cumpre todos os limites impostos para as concentrações de poluentes mas, apresenta algumas deficiências no que respeita aos caudais de ar novo insuflados em cada espaço.

Susceptibilidade à Erosão Hídrica na Bacia da Ribeira Seca (Santiago, Cabo Verde)

Na Ilha de Santiago, em Cabo Verde, a erosão hídrica é o processo que afecta áreas mais extensas. A ocorrência de aguaceiros intensos e concentrados no tempo e no espaço promovem uma marcada erosividade das precipitações, com forte irregularidade regional. A grande variabilidade dos declives e das formas das vertentes, associadas à diversidade litológica, bem como à multiplicidade de ocupação do solo, permitem condições de erodibilidade muito contrastadas no espaço. O objectivo deste trabalho é obter um mapa de susceptibilidade à erosão hídrica para a bacia da Ribeira Seca (Santiago oriental) com base no modelo digital do terreno (MDT), nos mapas geológico e de ocupação do solo e na distribuição da erosividade das precipitações. Verifica-se que o sector sudeste da bacia é o mais susceptível à erosão hídrica, pois nele ocorrem a maior concentração diária das precipitações e as condições geomorfológicas e de coberto do solo de mais elevada erodibilidade.

Parametric imaging for characterizing focal liver lesions in contrast-enhanced ultrasound.

The differentiation between benign and malignant focal liver lesions plays an important role in diagnosis of liver disease and therapeutic planning of local or general disease. This differentiation, based on characterization, relies on the observation of the dynamic vascular patterns (DVP) of lesions with respect to adjacent parenchyma, and may be assessed during contrast-enhanced ultrasound imaging after a bolus injection. For instance, hemangiomas (i.e., benign lesions) exhibit hyper-enhanced signatures over time, whereas metastases (i.e., malignant lesions) frequently present hyperenhanced foci during the arterial phase and always become hypo-enhanced afterwards. The objective of this work was to develop a new parametric imaging technique, aimed at mapping the DVP signatures into a single image called a DVP parametric image, conceived as a diagnostic aid tool for characterizing lesion types. The methodology consisted in processing a time sequence of images (DICOM video data) using four consecutive steps: (1) pre-processing combining image motion correction and linearization to derive an echo-power signal, in each pixel, proportional to local contrast agent concentration over time; (2) signal modeling, by means of a curve-fitting optimization, to compute a difference signal in each pixel, as the subtraction of adjacent parenchyma kinetic from the echopower signal; (3) classification of difference signals; and (4) parametric image rendering to represent classified pixels as a support for diagnosis. DVP parametric imaging was the object of a clinical assessment on a total of 146 lesions, imaged using different medical ultrasound systems. The resulting sensitivity and specificity were 97% and 91%, respectively, which compare favorably with scores of 81 to 95% and 80 to 95% reported in medical literature for sensitivity and specificity, respectively.

Tc-99m-MAA SPECT-derived tumor-to-normal liver ratios for partition model dosimetry in selective internal radiation therapy (SIRT)

Aim: When planning SIRT using 90Y microspheres, the partition model is used to refine the activity calculated by the body surface area (BSA) method to potentially improve the safety and efficacy of treatment. For this partition model dosimetry, accurate determination of mean tumor-to-normal liver ratio (TNR) is critical since it directly impacts absorbed dose estimates. This work aimed at developing and assessing a reliable methodology for the calculation of 99mTc-MAA SPECT/CT-derived TNR ratios based on phantom studies. Materials and methods: IQ NEMA (6 hot spheres) and Kyoto liver phantoms with different hot/background activity concentration ratios were imaged on a SPECT/CT (GE Infinia Hawkeye 4). For each reconstruction with the IQ phantom, TNR quantification was assessed in terms of relative recovery coefficients (RC) and image noise was evaluated in terms of coefficient of variation (COV) in the filled background. RCs were compared using OSEM with Hann, Butterworth and Gaussian filters, as well as FBP reconstruction algorithms. Regarding OSEM, RCs were assessed by varying different parameters independently, such as the number of iterations (i) and subsets (s) and the cut-off frequency of the filter (fc). The influence of the attenuation and diffusion corrections was also investigated. Furthermore, both 2D-ROIs and 3D-VOIs contouring were compared. For this purpose, dedicated Matlab© routines were developed in-house for automatic 2D-ROI/3D-VOI determination to reduce intra-user and intra-slice variability. Best reconstruction parameters and RCs obtained with the IQ phantom were used to recover corrected TNR in case of the Kyoto phantom for arbitrary hot-lesion size. In addition, we computed TNR volume histograms to better assess uptake heterogeneityResults: The highest RCs were obtained with OSEM (i=2, s=10) coupled with the Butterworth filter (fc=0.8). Indeed, we observed a global 20% RC improvement over other OSEM settings and a 50% increase as compared to the best FBP reconstruction. In any case, both attenuation and diffusion corrections must be applied, thus improving RC while preserving good image noise (COV<10%). Both 2D-ROI and 3D-VOI analysis lead to similar results. Nevertheless, we recommend using 3D-VOI since tumor uptake regions are intrinsically 3D. RC-corrected TNR values lie within 17% around the true value, substantially improving the evaluation of small volume (<15 mL) regions. Conclusions: This study reports the multi-parameter optimization of 99mTc MAA SPECT/CT images reconstruction in planning 90Y dosimetry for SIRT. In phantoms, accurate quantification of TNR was obtained using OSEM coupled with Butterworth and RC correction.

Vesicular transport of d-serine in astrocytes

The concept of tripartite synapse suggests that astrocytes make up a functional synapse with pre- and postsynaptic neuronal elements to modulate synaptic transmission through the regulated release of neuromodulators called gliotransmitters. Release of gliotransmitters such as glutamate or D-serine has been shown to depend on Ca21-dependent exocytosis. However, the origin (cytosolic versus vesicular) of the released gliotransmitter is still a matter of debate. The existence of Ca21-regulated exocytosis in astrocytes has been questioned mostly because the nature of secretory organelles which are loaded with gliotransmitters is unknown. Here we show the existence of a population of vesicles that uptakes and stores glutamate and D-serine in astrocytes which are present in situ. Immunoisolated glial organelles expressing synaptobrevin 2 (Sb2) display morphological and biochemical features very similar to synaptic vesicles. We demonstrate that these organelles not only contain and uptake glutamate but also display a glia-specific transport activity for D-serine. Furthermore, we report that the uptake of D-serine is energized by a H1-ATPase present on the immunoisolated vesicles and that cytosolic chloride ions modulate the uptake of D-serine. Finally, we show that serine racemase (SR), the synthesizing enzyme for D-serine, is anchored to the membrane of glial organelles allowing a local and efficient concentration of the gliotransmitter to be transported. We conclude that vesicles in astrocytes do exist with the goal to store and release D-serine, glutamate and most likely other neuromodulators.

Extracellular microvesicles from astrocytes contain functional glutamate transporters: regulation by protein kinase C and cell activation.

Glutamate transport through astrocytic excitatory amino-acid transporters (EAAT)-1 and EAAT-2 is paramount for neural homeostasis. EAAT-1 has been reported in secreted extracellular microvesicles (eMV, such as exosomes) and because the protein kinase C (PKC) family controls the sub-cellular distribution of EAATs, we have explored whether PKCs drive EAATs into eMV. Using rat primary astrocytes, confocal immunofluorescence and ultracentrifugation on sucrose gradient we here report that PKC activation by phorbol myristate acetate (PMA) reorganizes EAAT-1 distribution and reduces functional [(3)H]-aspartate reuptake. Western-blots show that EAAT-1 is present in eMV from astrocyte conditioned medium, together with NaK ATPase and glutamine synthetase all being further increased after PMA treatment. However, nanoparticle tracking analysis reveals that PKC activation did not change particle concentration. Functional analysis indicates that eMV have the capacity to reuptake [(3)H]-aspartate. In vivo, we demonstrate that spinal astrocytic reaction induced by peripheral nerve lesion (spared nerve injury, SNI) is associated with a phosphorylation of PKC δ together with a shift of EAAT distribution ipsilaterally. Ex vivo, spinal explants from SNI rats release eMV with an increased content of NaK ATPase, EAAT-1 and EAAT-2. These data indicate PKC and cell activation as important regulators of EAAT-1 incorporation in eMV, and raise the possibility that microvesicular EAAT-1 may exert extracellular functions. Beyond a putative role in neuropathic pain, this phenomenon may be important for understanding neural homeostasis and a wide range of neurological diseases associated with astrocytic reaction as well as non-neurological diseases linked to eMV release.

Aerosolianalysaattorin jatkokehitys

Kiinnostus pienhiukkasia kohtaan on kasvanut voimakkaasti, koska niiden haitallisista terveysvaikutuksista on saatu uutta tietoa. Asiasta on julkaistu lukuisia tieteellisiä tutkimuksia ja viimeisimpien tietojan mukaan kohonneilla pienhiukkaspitoisuuksilla on vaikutusta jopa sydän- ja verisuonisairauksiin. Vakavien terveysvaikutusten ja kiristyvän lainsäädännön vuoksi uusille reaaliaikaisille hiukkasmittalaitteille on kova kysyntä. Tämä työ on osa suurempaa Dekati Oy:ssä toteutettua Autotest –projektia, jossa kehitetään hiukkasmittalaitteita autoteollisuudelle. Tavoitteena työssä oli kehittää hiukkasmittalaitteeseen varaaja, jossa olisi huomattavasti pienemmät pienhiukkashäviöt kuin sähköisen alipaineimpaktorin ELPI:n varaajassa. Lainsäädäntö pohjautuu nykyään pelkästään massapitoisuuden mittaamiseen eikä reaaliaikaisesti massapitoisuutta mittaavaa laitetta ole olemassa. Tässä työssä testattiin uutta hiukkasten tiheydenmääritysmenetelmää, jonka avulla on mahdollista mitata massapitoisuutta reaaliaikaisesti. Suunniteltu varaaja on helppokäyttöinen ja varaustehokkuudeltaan odotusten mukainen, mutta pienhiukkashäviöt ovat edelleen tavoiteltua suuremmat vaikkakin pienemmät kuin ELPI:ssä. Tämä johtuu osittain tilavaraushäviöistä ja osittain koronan sähkökentän vaikutuksesta näytekanavaan. Tiheysmittauksesta saatiin lupaavia tuloksia, mutta jatkokehitystä vaaditaan häiriöiden suodattamiseksi ja kuormituskestävyyden parantamiseksi.

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.

Influence of reactant absorption and dissolution in heterogeneous precipitation processes

In this research work, the aim was to investigate the volumetric mass transfer coefficient [kLa] of oxygen in stirred tank in the presence of solid particle experimentally. The kLa correlations as a function of propeller rotation speed and flow rate of gas feed were studied. The O2 and CO2 absorption in water and in solid-liquid suspensions and heterogeneous precipitation of MgCO3 were thoroughly examined. The absorption experiments of oxygen were conducted in various systems like pure water and in aqueous suspensions of quartz and calcium carbonate particles. Secondly, the precipitation kinetics of magnesium carbonate was also investigated. The experiments were performed to study the reactive crystallization with magnesium hydroxide slurry and carbon dioxide gas by varying the feed rates of carbon dioxide and rotation speeds of mixer. The results of absorption and precipitation are evaluated by titration, total carbon (TC analysis), and ionic chromatrography (IC). For calcium carbonate, the particle concentration was varied from 17.4 g to 2382 g with two size fractions: 5 µm and 45-63 µm sieves. The kLa and P/V values of 17.4 g CaCO3 with particle size of 5µm and 45-63 µm were 0.016 s-1 and 2400 W/m3. At 69.9 g concentration of CaCO3, the achieved kLa is 0.014 s-1 with particle size of 5 µm and 0.017 s-1 with particle size of 45 to 63 µm. Further increase in concentration of calcium carbonate, i.e. 870g and 2382g , does not affect volumetric mass transfer coeffienct of oxygen. It could be concluded from absorption results that maximum value of kLa is 0.016 s-1. Also particle size and concentration does affect the transfer rate to some extend. For precipitation experiments, the constant concentration of Mg(OH)2 was 100 g and the rotation speed varied from 560 to 750 rpm, whereas the used feed rates of CO2 were 1 and 9 L/min. At 560 rpm and feed rate of CO2 is 1 L/min, the maximum value of Mg ion and TC were 0.25 mol/litre and 0.12 mol/litre with the residence time of 40 min. When flow rate of CO2 increased to 9 L/min with same 560 rpm, the achieved value of Mg and TC were 0.3 mol/litre and 0.12 mol/L with shorter residence time of 30 min. It is concluded that feed rate of CO2 is dominant in precipitation experiments and it has a key role in dissociation and reaction of magnesium hydroxide in precipitation of magnesium carbonate.

STUDIES ON THE STRESS FLUCTUATIONS IN SHEARED STOKESIAN SUSPENSIONS USING CHAOS THEORY AND NONLINEAR DYNAMICS

The thesis report results obtained from a detailed analysis of the fluctuations of the rheological parameters viz. shear and normal stresses, simulated by means of the Stokesian Dynamics method, of a macroscopically homogeneous sheared suspension of neutrally buoyant non-Brownian suspension of identical spheres in the Couette gap between two parallel walls in the limit of vanishingly small Reynolds numbers using the tools of non-linear dynamics and chaos theory for a range of particle concentration and Couette gaps. The thesis used the tools of nonlinear dynamics and chaos theory viz. average mutual information, space-time separation plots, visual recurrence analysis, principal component analysis, false nearest-neighbor technique, correlation integrals, computation of Lyapunov exponents for a range of area fraction of particles and for different Couette gaps. The thesis observed that one stress component can be predicted using another stress component at the same area fraction. This implies a type of synchronization of one stress component with another stress component. This finding suggests us to further analysis of the synchronization of stress components with another stress component at the same or different area fraction of particles. The different model equations of stress components for different area fraction of particles hints at the possible existence a general formula for stress fluctuations with area fraction of particle as a parameter

Modeling the plant uptake of organic chemicals, including the soil-air-plant pathway

The soil−air−plant pathway is potentially important in the vegetative accumulation of organic pollutants from contaminated soils. While a number of qualitative frameworks exist for the prediction of plant accumulation of organic chemicals by this pathway, there are few quantitative models that incorporate this pathway. The aim of the present study was to produce a model that included this pathway and could quantify its contribution to the total plant contamination for a range of organic pollutants. A new model was developed from three submodels for the processes controlling plant contamination via this pathway: aerial deposition, soil volatilization, and systemic translocation. Using the combined model, the soil−air−plant pathway was predicted to account for a significant proportion of the total shoot contamination for those compounds with log KOA > 9 and log KAW < −3. For those pollutants with log KOA < 9 and log KAW > −3 there was a higher deposition of pollutant via the soil−air−plant pathway than for those chemicals with log KOA > 9 and log KAW < −3, but this was an insignificant proportion of the total shoot contamination because of the higher mobility of these compounds via the soil−root−shoot pathway. The incorporation of the soil−air−plant pathway into the plant uptake model did not significantly improve the prediction of the contamination of vegetation from polluted soils when compared across a range of studies. This was a result of the high variability between the experimental studies where the bioconcentration factors varied by 2 orders of magnitude at an equivalent log KOA. One potential reason for this is the background air concentration of the pollutants under study. It was found background air concentrations would dominate those from soil volatilization in many situations unless there was a soil hot spot of contamination, i.e., >100 mg kg−1.

The temperature response to stratospheric water vapour changes

This study uses an analytical model, based on the cooling-to-space approximation, and a fixed dynamical heating model to investigate the structure of the stratospheric cooling that occurs in response to a uniform increase in stratospheric water vapour (SWV). At all latitudes, the largest cooling occurs in the lower stratosphere and decreases in magnitude with height. The cooling is strongly enhanced in the Extratropics compared to the Tropics. This is markedly different to the case of an increase in CO2, which causes maximum cooling near the stratopause and a small warming in the tropical lower stratosphere. The qualitative differences in the structure of the cooling can be explained by the smaller opacity of water vapour bands in the stratosphere compared to CO2. The small opacity means that the magnitude of the initial heating rate perturbation only decreases by a factor of four between the upper and lower stratosphere for a SWV perturbation. Therefore, to balance the heating rate perturbation, the largest temperature change is required in the lower stratosphere. Increasing the background concentration of SWV causes the water vapour bands to become more opaque. For a SWV perturbation applied to a background SWV concentration ≥30 ppmv, the heating rate perturbation and temperature change structurally resemble those from an increase in CO2. In the Extratropics, the lower height of the tropopause is found to cause the enhancement in the cooling at those latitudes. By controlling the depth of atmosphere which adjusts to the SWV perturbation, the tropopause height affects the net exchange of radiation between the layers in the stratosphere as they cool. The latitudinal gradient in upwelling infrared radiation at the tropopause and variations in the background temperature are found to have only a minor effect on the structure of the stratospheric temperature response to a change in SWV.

Creep and recovery of magnetorheological fluids: experiments and simulations

A direct comparative study on the creep-recovery behavior of conventional MR fluids is carried out using magnetorheometry and particle-level simulations. Two particle concentrations are investigated (ϕ=0.05 and 0.30) at two different magnetic field strengths (53 kA•m-1 and 173 kA•m-1) in order to match the yield stresses developed in both systems for easier comparison. Simulations are mostly started with random initial structures with some additional tests of using preassembled single chains in the low concentration case. Experimental and simulation data are in good qualitative agreement. The results demonstrate three regions in the creep curves: i) In the initial viscoelastic region, the chain-like (at ϕ=0.05) or percolated three-dimensional network (at ϕ=0.30) structures fill up the gap and the average cluster size remains constant; ii) Above a critical strain of 10 %, in the retardation region, these structures begin to break and rearrange under shear. At large enough imposed stress values, they transform into thin sheet-like or thick lamellar structures, depending on the particle concentration; iii) Finally in the case of larger strain values either the viscosity diverges (at low stress values) or reaches a constant low value (at high stress values), showing a clear bifurcation behavior. For stresses below the bifurcation point the MR fluid is capable to recover the strain by a certain fraction. However, no recovery is observed for large stress values.

Many-body effect in electrorheological responses

The many-body effect in the kinetic responses of ER fluids is studied by a molecular-dynamic simulation method. The mutual polarization effects of the particles are considered by self-consistently calculating the dipole strength on each particle according to the external field and the dipole field due to all the other particles in the fluids. The many-body effect is found to increase with the enhancement of the particle concentration and the permittivity ratio between the solvent and the particles. The calculated response times are shorter than that predicted with the 'point-dipole' model and agree very well with experimental results. The many-body effect enhances the shear stresses of the fluids by several times. But they are not proportional to the many-body correction factor lambda as expected. This is due to the fact that larger interaction forces between the particles lead to coarsening of the fibers formed in the suspensions. The results show that the many-body and multipolar interaction between the particles must be treated comprehensively in the simulations in order to get more reliable results.

Holocene fires in East Amazonia (Caraja`s), new evidences, chronology and relation with paleoclimate

Past studies have evidenced the presence of charcoal in soils and lacustrine sediments of Amazonia region and suggested occurrences of widespread fires during the Middle Holocene. However, the available records do not indicate the changes in fire regime with enough time resolution. We quantified charcoal fragments in lacustrine sediments in a lake of North Carajas plateau in East Amazonia (5 degrees 50`-6 degrees 35`S and 49 degrees 30`-52 degrees 00`W). The charcoal quantification was compared to other sediment proxies, allowing a connection between paleofires and climate changes. Large variations in sediment characteristics led to distinct stages of sedimentation. From 11,800 (base of CSN 93/4) to 4750 cal yr B.P., low accumulation rates of organic matter are observed. Between 7600 cal yr B.P. (base of CSN 93/3 core) and 4750 cal yr B.P., this initial phase of sedimentation is characterized by low chlorophyll derivate accumulation rates and high accumulation rates of Botryococcus braunii, an alga resistant to episodic drought. The first phase of sedimentation would therefore correspond to, a low take level and a drier climate than today. Large biomass burning events occurred between 7450 cal yr B.P. and 4750 cal yr B.P., as indicated by the high charcoal particle concentration. From 4750 cal yr B.P. to 2800 cal yr B.P., accumulation rates of charcoal particles decreased, and the accumulation rate of chlorophyll derivate was low. From 2800 cal yr B.P. to 1300 cal yr B.P., the charcoal accumulation rates reached their lowest values in the core and a rapid increase in lacustrine production is evidenced by the increase in chlorophyll derivates and carbon accumulation rate. From 1300 cal yr B.P. to the last century, the charcoal accumulation rates increased. During the most recent period, the record is characterized by high accumulation rates of chlorophyll derivates while the charcoal particle accumulation rate decreased. This region is still unaffected by the current increase of anthropogenic fires. (c) 2007 Elsevier B.V. All rights reserved.