Vehicle emissions and PM2.5 mass concentrations in six Brazilian cities

In Brazil, the principal source of air pollution is the combustion of fuels (ethanol, gasohol, and diesel). In this study, we quantify the contributions that vehicle emissions make to the urban fine particulate matter (PM2.5) mass in six state capitals in Brazil, collecting data for use in a larger project evaluating the impact of air pollution on human health. From winter 2007 to winter 2008, we collected 24-h PM2.5 samples, employing gravimetry to determine PM2.5 mass concentrations; reflectance to quantify black carbon concentrations; X-ray fluorescence to characterize elemental composition; and ion chromatography to determine the composition and concentrations of anions and cations. Mean PM2.5 concentrations in the cities of Sao Paulo, Rio de Janeiro, Belo Horizonte, Curitiba, Porto Alegre, and Recife were 28, 17.2, 14.7, 14.4, 13.4, and 7.3 mu g/m(3), respectively. In Sao Paulo and Rio de Janeiro, black carbon explained approximately 30% of the PM2.5 mass. We used receptor models to identify distinct source-related PM2.5 fractions and correlate those fractions with daily mortality rates. Using specific rotation factor analysis, we identified the following principal contributing factors: soil and crustal material; vehicle emissions and biomass burning (black carbon factor); and fuel oil combustion in industries (sulfur factor). In all six cities, vehicle emissions explained at least 40% of the PM2.5 mass. Elemental composition determination with receptor modeling proved an adequate strategy to identify air pollution sources and to evaluate their short- and long-term effects on human health. Our data could inform decisions regarding environmental policies vis-a-vis health care costs.

Mapping the Energy Distribution of SERRS Hot Spots from Anti-Stokes to Stokes Intensity Ratios

The anomalies in the anti-Stokes to Stokes intensity ratios in single-molecule surface-enhanced resonance Raman scattering were investigated. Brilliant green and crystal violet dyes were the molecular probes, and the experiments were carried out on an electrochemically activated Ag surface. The results allowed new insights into the origin of these anomalies and led to a new method to confirm the single-molecule regime in surface-enhanced Raman scattering. Moreover, a methodology to estimate the distribution of resonance energies that contributed to the imbalance in the anti-Stokes to Stokes intensity ratios at the electromagnetic hot spots was proposed. This method allowed the local plasmonic resonance energies on the metallic surface to be spatially mapped.

LC-DAD/ESI-MS/MS study of phenolic compounds in ash (Fraxinus excelsior L. and F. americana L.) heartwood. Effect of toasting intensity at cooperage

The phenolic composition of heartwood extracts from Fraxinus excelsior L. and F. americana L., both before and after toasting in cooperage, was studied using LC-DAD/ESI-MS/MS. Low-molecular weight (LMW) phenolic compounds, secoiridoids, phenylethanoid glycosides, dilignols and oligolignols compounds were detected, and 48 were identified, or tentatively characterized, on the basis of their retention time, UV/Vis and MS spectra, and MS fragmentation patterns. Some LMW phenolic compounds like protocatechuic acid and aldehyde, hydroxytyrosol and tyrosol, were unlike to those for oak wood, while ellagic and gallic acid were not found. The toasting of wood resulted in a progressive increase in lignin degradation products with regard to toasting intensity. The levels of some of these compounds in medium-toasted ash woods were much higher than those normally detected in toasted oak, highlighting vanillin levels, thus a more pronounced vanilla character can be expected when using toasted ash wood in the aging wines. Moreover, in seasoned wood, we found a great variety of phenolic compounds which had not been found in oak wood, especially oleuropein, ligstroside and olivil, along with verbascoside and isoverbascoside in F. excelsior, and oleoside in F. americana. Toasting mainly provoked their degradation, thus in medium-toasted wood, only four of them were detected. This resulted in a minor differentiation between toasted ash and oak woods. The absence of tannins in ash wood, which are very important in oak wood, is another peculiar characteristic that should be taken into account when considering its use in cooperage. Copyright (C) 2012 John Wiley & Sons, Ltd.

Metabolic, hemodynamic and structural adjustments to low intensity exercise training in a metabolic syndrome model

Background The increase in fructose consumption is paralleled by a higher incidence of metabolic syndrome, and consequently, cardiovascular disease mortality. We examined the effects of 8 weeks of low intensity exercise training (LET) on metabolic, hemodynamic, ventricular and vascular morphological changes induced by fructose drinking in male rats. Methods Male Wistar rats were divided into (n = 8 each) control (C), sedentary fructose (F) and ET fructose (FT) groups. Fructose-drinking rats received D-fructose (100 g/l). FT rats were assigned to a treadmill training protocol at low intensity (30% of maximal running speed) during 1 h/day, 5 days/week for 8 weeks. Measurements of triglyceride concentrations, white adipose tissue (WAT) and glycemia were carried out together with insulin tolerance test to evaluate metabolic profile. Arterial pressure (AP) signals were directly recorded. Baroreflex sensitivity (BS) was evaluated by the tachycardic and bradycardic responses. Right atria, left ventricle (LV) and ascending aorta were prepared to morphoquantitative analysis. Results LET reduced WAT (−37.7%), triglyceride levels (−33%), systolic AP (−6%), heart weight/body weight (−20.5%), LV (−36%) and aortic (−76%) collagen fibers, aortic intima-media thickness and circumferential wall tension in FT when compared to F rats. Additionally, FT group presented improve of BS, numerical density of atrial natriuretic peptide granules (+42%) and LV capillaries (+25%), as well as the number of elastic lamellae in aorta compared with F group. Conclusions Our data suggest that LET, a widely recommended practice, seems to be particularly effective for preventing metabolic, hemodynamic and morphological disorders triggered by MS.

Relations between soil surface roughness, tortuosity, tillage treatments, rainfall intensity and soil and water losses from a red yellow latosol

The soil surface roughness increases water retention and infiltration, reduces the runoff volume and speed and influences soil losses by water erosion. Similarly to other parameters, soil roughness is affected by the tillage system and rainfall volume. Based on these assumptions, the main purpose of this study was to evaluate the effect of tillage treatments on soil surface roughness (RR) and tortuosity (T) and to investigate the relationship with soil and water losses in a series of simulated rainfall events. The field study was carried out at the experimental station of EMBRAPA Southeastern Cattle Research Center in São Carlos (Fazenda Canchim), in São Paulo State, Brazil. Experimental plots of 33 m² were treated with two tillage practices in three replications, consisting of: untilled (no-tillage) soil (NTS) and conventionally tilled (plowing plus double disking) soil (CTS). Three successive simulated rain tests were applied in 24 h intervals. The three tests consisted of a first rain of 30 mm/h, a second of 30 mm/h and a third rain of 70 mm/h. Immediately after tilling and each rain simulation test, the surface roughness was measured, using a laser profile meter. The tillage treatments induced significant changes in soil surface roughness and tortuosity, demonstrating the importance of the tillage system for the physical surface conditions, favoring water retention and infiltration in the soil. The increase in surface roughness by the tillage treatments was considerably greater than its reduction by rain action. The surface roughness and tortuosity had more influence on the soil volume lost by surface runoff than in the conventional treatment. Possibly, other variables influenced soil and water losses from the no-tillage treatments, e.g., soil type, declivity, slope length, among others not analyzed in this study.

Effects of high-intensity intermittent training on carnitine palmitoyl transferase activity in the gastrocnemius muscle of rats

We examined the capacity of high-intensity intermittent training (HI-IT) to facilitate the delivery of lipids to enzymes responsible for oxidation, a task performed by the carnitine palmitoyl transferase (CPT) system in the rat gastrocnemius muscle. Male adult Wistar rats (160-250 g) were randomly distributed into 3 groups: sedentary (Sed, N = 5), HI-IT (N = 10), and moderate-intensity continuous training (MI-CT, N = 10). The trained groups were exercised for 8 weeks with a 10% (HI-IT) and a 5% (MI-CT) overload. The HI-IT group presented 11.8% decreased weight gain compared to the Sed group. The maximal activities of CPT-I, CPT-II, and citrate synthase were all increased in the HI-IT group compared to the Sed group (P < 0.01), as also was gene expression, measured by RT-PCR, of fatty acid binding protein (FABP; P < 0.01) and lipoprotein lipase (LPL; P < 0.05). Lactate dehydrogenase also presented a higher maximal activity (nmol·min-1·mg protein-1) in HI-IT (around 83%). We suggest that 8 weeks of HI-IT enhance mitochondrial lipid transport capacity thus facilitating the oxidation process in the gastrocnemius muscle. This adaptation may also be associated with the decrease in weight gain observed in the animals and was concomitant to a higher gene expression of both FABP and LPL in HI-IT, suggesting that intermittent exercise is a "time-efficient" strategy inducing metabolic adaptation.

Emissions generated by sugarcane burning promote genotoxicity in rural workers: a case study in Barretos, Brazil

Abstract Background To determine the possible genotoxic effect of exposure to the smoke generated by biomass burning on workers involved in manual sugar cane harvesting. Methods The frequency of micronuclei in exfoliated buccal cells and peripheral blood lymphocytes was determined in sugarcane workers in the Barretos region of Brazil, during the harvest season and compared to a control population, comprised of administrative employees of Barretos Cancer Hospital. Results The frequency of micronuclei was higher in the sugar cane workers. The mean frequency in blood lymphocytes (micronuclei/1000 cells) in the test group was 8.22 versus 1.27 in the control group. The same effect was observed when exfoliated buccal cells were considered (22.75 and 9.70 micronuclei/1000 cells for sugar cane workers and controls, respectively). Conclusion Exposure to emissions produced by the burning of sugar cane during harvesting induces genomic instability in workers, indicating the necessity of adopting more advanced techniques of harvesting sugar cane to preserve human health.

Muscular and pulmonary O2 uptake kinetics during moderate- and high-intensity sub-maximal knee-extensor exercise in humans

[EN] The purpose of this investigation was to determine the contribution of muscle O(2) consumption (mVO2) to pulmonary O(2) uptake (pVO2) during both low-intensity (LI) and high-intensity (HI) knee-extension exercise, and during subsequent recovery, in humans. Seven healthy male subjects (age 20-25 years) completed a series of LI and HI square-wave exercise tests in which mVO2 (direct Fick technique) and pVO2 (indirect calorimetry) were measured simultaneously. The mean blood transit time from the muscle capillaries to the lung (MTTc-l) was also estimated (based on measured blood transit times from femoral artery to vein and vein to artery). The kinetics of mVO2 and pVO2 were modelled using non-linear regression. The time constant (tau) describing the phase II pVO2 kinetics following the onset of exercise was not significantly different from the mean response time (initial time delay + tau) for mVO2 kinetics for LI (30 +/- 3 vs 30 +/- 3 s) but was slightly higher (P < 0.05) for HI (32 +/- 3 vs 29 +/- 4 s); the responses were closely correlated (r = 0.95 and r = 0.95; P < 0.01) for both intensities. In recovery, agreement between the responses was more limited both for LI (36 +/- 4 vs 18 +/- 4 s, P < 0.05; r = -0.01) and HI (33 +/- 3 vs 27 +/- 3 s, P > 0.05; r = -0.40). MTTc-l was approximately 17 s just before exercise and decreased to 12 and 10 s after 5 s of exercise for LI and HI, respectively. These data indicate that the phase II pVO2 kinetics reflect mVO2 kinetics during exercise but not during recovery where caution in data interpretation is advised. Increased mVO2 probably makes a small contribution to during the first 15-20 s of exercise.

Effects of recovery mode on performance, O2 uptake, and O2 deficit during high-intensity intermittent exercise

[EN] The aim of this study was to determine the influence of activity performed during the recovery period on the aerobic and anaerobic energy yield, as well as on performance, during high-intensity intermittent exercise (HIT). Ten physical education students participated in the study. First they underwent an incremental exercise test to assess their maximal power output (Wmax) and VO2max. On subsequent days they performed three different HITs. Each HIT consisted of four cycling bouts until exhaustion at 110% Wmax. Recovery periods of 5 min were allowed between bouts. HITs differed in the kind of activity performed during the recovery periods: pedaling at 20% VO2max (HITA), stretching exercises, or lying supine. Performance was 3-4% and aerobic energy yield was 6-8% (both p < 0.05) higher during the HITA than during the other two kinds of HIT. The greater contribution of aerobic metabolism to the energy yield during the high-intensity exercise bouts with active recovery was due to faster VO2 kinetics (p< 0.01) and a higher VO2peak during the exercise bouts preceded by active recovery (p < 0.05). In contrast, the anaerobic energy yield (oxygen deficit and peak blood lactate concentrations) was similar in all HITs. Therefore, this study shows that active recovery facilitates performance by increasing aerobic contribution to the whole energy yield turnover during high-intensity intermittent exercise.

Submarine emissions on the seafloor : from cold seeps to hot vents

[EN] The main types of submarine geological emissions are classified as cold seeps (hydrocarbons and brines) and hot vents. These processes result in the emission of geological fluids: brine, gases (mainly hydro-carbons), sediments and rocks. Submarine emissions are associated with an intensive geological, geo-chemical, thermal and biological activity (Judd and Hovland, 2007), and constitute a key process in the dynamics of the global cycles of the planet Earth.

Herbivory intensity over two coexisting macrophytes on soft bottoms of Gran Canaria Island

Máster en Oceanografía

Statistical properties of Radio Halos and the re-acceleration model

Galaxy clusters occupy a special position in the cosmic hierarchy as they are the largest bound structures in the Universe. There is now general agreement on a hierarchical picture for the formation of cosmic structures, in which galaxy clusters are supposed to form by accretion of matter and merging between smaller units. During merger events, shocks are driven by the gravity of the dark matter in the diffuse barionic component, which is heated up to the observed temperature. Radio and hard-X ray observations have discovered non-thermal components mixed with the thermal Intra Cluster Medium (ICM) and this is of great importance as it calls for a “revision” of the physics of the ICM. The bulk of present information comes from the radio observations which discovered an increasing number of Mpcsized emissions from the ICM, Radio Halos (at the cluster center) and Radio Relics (at the cluster periphery). These sources are due to synchrotron emission from ultra relativistic electrons diffusing through µG turbulent magnetic fields. Radio Halos are the most spectacular evidence of non-thermal components in the ICM and understanding the origin and evolution of these sources represents one of the most challenging goal of the theory of the ICM. Cluster mergers are the most energetic events in the Universe and a fraction of the energy dissipated during these mergers could be channelled into the amplification of the magnetic fields and into the acceleration of high energy particles via shocks and turbulence driven by these mergers. Present observations of Radio Halos (and possibly of hard X-rays) can be best interpreted in terms of the reacceleration scenario in which MHD turbulence injected during these cluster mergers re-accelerates high energy particles in the ICM. The physics involved in this scenario is very complex and model details are difficult to test, however this model clearly predicts some simple properties of Radio Halos (and resulting IC emission in the hard X-ray band) which are almost independent of the details of the adopted physics. In particular in the re-acceleration scenario MHD turbulence is injected and dissipated during cluster mergers and thus Radio Halos (and also the resulting hard X-ray IC emission) should be transient phenomena (with a typical lifetime <» 1 Gyr) associated with dynamically disturbed clusters. The physics of the re-acceleration scenario should produce an unavoidable cut-off in the spectrum of the re-accelerated electrons, which is due to the balance between turbulent acceleration and radiative losses. The energy at which this cut-off occurs, and thus the maximum frequency at which synchrotron radiation is produced, depends essentially on the efficiency of the acceleration mechanism so that observations at high frequencies are expected to catch only the most efficient phenomena while, in principle, low frequency radio surveys may found these phenomena much common in the Universe. These basic properties should leave an important imprint in the statistical properties of Radio Halos (and of non-thermal phenomena in general) which, however, have not been addressed yet by present modellings. The main focus of this PhD thesis is to calculate, for the first time, the expected statistics of Radio Halos in the context of the re-acceleration scenario. In particular, we shall address the following main questions: • Is it possible to model “self-consistently” the evolution of these sources together with that of the parent clusters? • How the occurrence of Radio Halos is expected to change with cluster mass and to evolve with redshift? How the efficiency to catch Radio Halos in galaxy clusters changes with the observing radio frequency? • How many Radio Halos are expected to form in the Universe? At which redshift is expected the bulk of these sources? • Is it possible to reproduce in the re-acceleration scenario the observed occurrence and number of Radio Halos in the Universe and the observed correlations between thermal and non-thermal properties of galaxy clusters? • Is it possible to constrain the magnetic field intensity and profile in galaxy clusters and the energetic of turbulence in the ICM from the comparison between model expectations and observations? Several astrophysical ingredients are necessary to model the evolution and statistical properties of Radio Halos in the context of re-acceleration model and to address the points given above. For these reason we deserve some space in this PhD thesis to review the important aspects of the physics of the ICM which are of interest to catch our goals. In Chapt. 1 we discuss the physics of galaxy clusters, and in particular, the clusters formation process; in Chapt. 2 we review the main observational properties of non-thermal components in the ICM; and in Chapt. 3 we focus on the physics of magnetic field and of particle acceleration in galaxy clusters. As a relevant application, the theory of Alfv´enic particle acceleration is applied in Chapt. 4 where we report the most important results from calculations we have done in the framework of the re-acceleration scenario. In this Chapter we show that a fraction of the energy of fluid turbulence driven in the ICM by the cluster mergers can be channelled into the injection of Alfv´en waves at small scales and that these waves can efficiently re-accelerate particles and trigger Radio Halos and hard X-ray emission. The main part of this PhD work, the calculation of the statistical properties of Radio Halos and non-thermal phenomena as expected in the context of the re-acceleration model and their comparison with observations, is presented in Chapts.5, 6, 7 and 8. In Chapt.5 we present a first approach to semi-analytical calculations of statistical properties of giant Radio Halos. The main goal of this Chapter is to model cluster formation, the injection of turbulence in the ICM and the resulting particle acceleration process. We adopt the semi–analytic extended Press & Schechter (PS) theory to follow the formation of a large synthetic population of galaxy clusters and assume that during a merger a fraction of the PdV work done by the infalling subclusters in passing through the most massive one is injected in the form of magnetosonic waves. Then the processes of stochastic acceleration of the relativistic electrons by these waves and the properties of the ensuing synchrotron (Radio Halos) and inverse Compton (IC, hard X-ray) emission of merging clusters are computed under the assumption of a constant rms average magnetic field strength in emitting volume. The main finding of these calculations is that giant Radio Halos are naturally expected only in the more massive clusters, and that the expected fraction of clusters with Radio Halos is consistent with the observed one. In Chapt. 6 we extend the previous calculations by including a scaling of the magnetic field strength with cluster mass. The inclusion of this scaling allows us to derive the expected correlations between the synchrotron radio power of Radio Halos and the X-ray properties (T, LX) and mass of the hosting clusters. For the first time, we show that these correlations, calculated in the context of the re-acceleration model, are consistent with the observed ones for typical µG strengths of the average B intensity in massive clusters. The calculations presented in this Chapter allow us to derive the evolution of the probability to form Radio Halos as a function of the cluster mass and redshift. The most relevant finding presented in this Chapter is that the luminosity functions of giant Radio Halos at 1.4 GHz are expected to peak around a radio power » 1024 W/Hz and to flatten (or cut-off) at lower radio powers because of the decrease of the electron re-acceleration efficiency in smaller galaxy clusters. In Chapt. 6 we also derive the expected number counts of Radio Halos and compare them with available observations: we claim that » 100 Radio Halos in the Universe can be observed at 1.4 GHz with deep surveys, while more than 1000 Radio Halos are expected to be discovered in the next future by LOFAR at 150 MHz. This is the first (and so far unique) model expectation for the number counts of Radio Halos at lower frequency and allows to design future radio surveys. Based on the results of Chapt. 6, in Chapt.7 we present a work in progress on a “revision” of the occurrence of Radio Halos. We combine past results from the NVSS radio survey (z » 0.05 − 0.2) with our ongoing GMRT Radio Halos Pointed Observations of 50 X-ray luminous galaxy clusters (at z » 0.2−0.4) and discuss the possibility to test our model expectations with the number counts of Radio Halos at z » 0.05 − 0.4. The most relevant limitation in the calculations presented in Chapt. 5 and 6 is the assumption of an “averaged” size of Radio Halos independently of their radio luminosity and of the mass of the parent clusters. This assumption cannot be released in the context of the PS formalism used to describe the formation process of clusters, while a more detailed analysis of the physics of cluster mergers and of the injection process of turbulence in the ICM would require an approach based on numerical (possible MHD) simulations of a very large volume of the Universe which is however well beyond the aim of this PhD thesis. On the other hand, in Chapt.8 we report our discovery of novel correlations between the size (RH) of Radio Halos and their radio power and between RH and the cluster mass within the Radio Halo region, MH. In particular this last “geometrical” MH − RH correlation allows us to “observationally” overcome the limitation of the “average” size of Radio Halos. Thus in this Chapter, by making use of this “geometrical” correlation and of a simplified form of the re-acceleration model based on the results of Chapt. 5 and 6 we are able to discuss expected correlations between the synchrotron power and the thermal cluster quantities relative to the radio emitting region. This is a new powerful tool of investigation and we show that all the observed correlations (PR − RH, PR − MH, PR − T, PR − LX, . . . ) now become well understood in the context of the re-acceleration model. In addition, we find that observationally the size of Radio Halos scales non-linearly with the virial radius of the parent cluster, and this immediately means that the fraction of the cluster volume which is radio emitting increases with cluster mass and thus that the non-thermal component in clusters is not self-similar.

Improved ground motion intensity measures for reliability-based demand analysis of structures

In Performance-Based Earthquake Engineering (PBEE), evaluating the seismic performance (or seismic risk) of a structure at a designed site has gained major attention, especially in the past decade. One of the objectives in PBEE is to quantify the seismic reliability of a structure (due to the future random earthquakes) at a site. For that purpose, Probabilistic Seismic Demand Analysis (PSDA) is utilized as a tool to estimate the Mean Annual Frequency (MAF) of exceeding a specified value of a structural Engineering Demand Parameter (EDP). This dissertation focuses mainly on applying an average of a certain number of spectral acceleration ordinates in a certain interval of periods, Sa,avg (T1,…,Tn), as scalar ground motion Intensity Measure (IM) when assessing the seismic performance of inelastic structures. Since the interval of periods where computing Sa,avg is related to the more or less influence of higher vibration modes on the inelastic response, it is appropriate to speak about improved IMs. The results using these improved IMs are compared with a conventional elastic-based scalar IMs (e.g., pseudo spectral acceleration, Sa ( T(¹)), or peak ground acceleration, PGA) and the advanced inelastic-based scalar IM (i.e., inelastic spectral displacement, Sdi). The advantages of applying improved IMs are: (i ) "computability" of the seismic hazard according to traditional Probabilistic Seismic Hazard Analysis (PSHA), because ground motion prediction models are already available for Sa (Ti), and hence it is possibile to employ existing models to assess hazard in terms of Sa,avg, and (ii ) "efficiency" or smaller variability of structural response, which was minimized to assess the optimal range to compute Sa,avg. More work is needed to assess also "sufficiency" and "scaling robustness" desirable properties, which are disregarded in this dissertation. However, for ordinary records (i.e., with no pulse like effects), using the improved IMs is found to be more accurate than using the elastic- and inelastic-based IMs. For structural demands that are dominated by the first mode of vibration, using Sa,avg can be negligible relative to the conventionally-used Sa (T(¹)) and the advanced Sdi. For structural demands with sign.cant higher-mode contribution, an improved scalar IM that incorporates higher modes needs to be utilized. In order to fully understand the influence of the IM on the seismis risk, a simplified closed-form expression for the probability of exceeding a limit state capacity was chosen as a reliability measure under seismic excitations and implemented for Reinforced Concrete (RC) frame structures. This closed-form expression is partuclarly useful for seismic assessment and design of structures, taking into account the uncertainty in the generic variables, structural "demand" and "capacity" as well as the uncertainty in seismic excitations. The assumed framework employs nonlinear Incremental Dynamic Analysis (IDA) procedures in order to estimate variability in the response of the structure (demand) to seismic excitations, conditioned to IM. The estimation of the seismic risk using the simplified closed-form expression is affected by IM, because the final seismic risk is not constant, but with the same order of magnitude. Possible reasons concern the non-linear model assumed, or the insufficiency of the selected IM. Since it is impossibile to state what is the "real" probability of exceeding a limit state looking the total risk, the only way is represented by the optimization of the desirable properties of an IM.

Adaptive finite element simulation of stack pollutant emissions over complex terrains

[EN]A three-dimensional finite element model for the pollutant dispersion is presented. In these environmental processes over a complex terrain, a mesh generator capable of adapting itself to the topographic characteristics is essential. The first stage of the model consists on the construction of an adaptive tetrahedral mesh of a rectangular region bounded in its lower part by the terrain and in its upper part by a horizontal plane. Once the mesh is constructed, an adaptive local refinement of tetrahedra is used in order to capture the plume rise. Wind measurements are used to compute an interpolated wind field, that is modified by using a mass-consistent model and perturbing its vertical component to introduce the plume rise effect...