Scalar field description for parametric models of dark energy

We investigate theoretical and observational aspects of a time-dependent parameterization for the dark energy equation of state w(z), which is a well behaved function of the redshift z over the entire cosmological evolution, i.e., z is an element of [-1, infinity). By using a theoretical algorithm of constructing the quintes-sence potential directly from the w(z) function, we derive and discuss the general features of the resulting potential for the cases in which dark energy is separately conserved and when it is coupled to dark matter. Since the parameterization here discussed allows us to divide the parametric plane in defined regions associated to distinct classes of dark energy models, we use some of the most recent observations from type Ia supernovae, baryon acoustic oscillation peak and Cosmic Microwave Background shift parameter to check which class is observationally preferred. We show that the largest portion of the confidence contours lies into the region corresponding to a possible crossing of the so-called phantom divide line at some point of the cosmic evolution.

STRUCTURAL CHARACTERIZATION OF CANOPIES OF Eucalyptus spp. USING RADIOMETRIC DATA FROM TM/Landsat 5

Empirical approaches and, more recently, physical approaches, have grounded the establishment of logical connections between radiometric variables derived from remote data and biophysical variables derived from vegetation cover. This study was aimed at evaluating correlations of dendrometric and density data from canopies of Eucalyptus spp., as collected in Capao Bonito forest unit, with radiometric data from imagery acquired by the TM/Landsat-5 sensor on two orbital passages over the study site (dates close to field data collection). Results indicate that stronger correlations were identified between crown dimensions and canopy height with near-infrared spectral band data (rho(s)4), irrespective of the satellite passage date. Estimates of spatial distribution of dendrometric data and canopy density (D) using spectral characterization were consistent with the spatial distribution of tree ages during the study period. Statistical tests were applied to evaluate performance disparities of empirical models depending on which date data were acquired. Results indicated a significant difference between models based on distinct data acquisition dates.

Low-Resolution Molecular Models Reveal the Oligomeric State of the PPAR and the Conformational Organization of Its Domains in Solution

The peroxisome proliferator-activated receptors (PPARs) regulate genes involved in lipid and carbohydrate metabolism, and are targets of drugs approved for human use. Whereas the crystallographic structure of the complex of full length PPAR gamma and RXR alpha is known, structural alterations induced by heterodimer formation and DNA contacts are not well understood. Herein, we report a small-angle X-ray scattering analysis of the oligomeric state of hPPAR gamma alone and in the presence of retinoid X receptor (RXR). The results reveal that, in contrast with other studied nuclear receptors, which predominantly form dimers in solution, hPPAR gamma remains in the monomeric form by itself but forms heterodimers with hRXR alpha. The low-resolution models of hPPAR gamma/RXR alpha complexes predict significant changes in opening angle between heterodimerization partners (LBD) and extended and asymmetric shape of the dimer (LBD-DBD) as compared with X-ray structure of the full-length receptor bound to DNA. These differences between our SAXS models and the high-resolution crystallographic structure might suggest that there are different conformations of functional heterodimer complex in solution. Accordingly, hydrogen/deuterium exchange experiments reveal that the heterodimer binding to DNA promotes more compact and less solvent-accessible conformation of the receptor complex.

New Optimization Algorithms for Structural Reliability Analysis

Solution of structural reliability problems by the First Order method require optimization algorithms to find the smallest distance between a limit state function and the origin of standard Gaussian space. The Hassofer-Lind-Rackwitz-Fiessler (HLRF) algorithm, developed specifically for this purpose, has been shown to be efficient but not robust, as it fails to converge for a significant number of problems. On the other hand, recent developments in general (augmented Lagrangian) optimization techniques have not been tested in aplication to structural reliability problems. In the present article, three new optimization algorithms for structural reliability analysis are presented. One algorithm is based on the HLRF, but uses a new differentiable merit function with Wolfe conditions to select step length in linear search. It is shown in the article that, under certain assumptions, the proposed algorithm generates a sequence that converges to the local minimizer of the problem. Two new augmented Lagrangian methods are also presented, which use quadratic penalties to solve nonlinear problems with equality constraints. Performance and robustness of the new algorithms is compared to the classic augmented Lagrangian method, to HLRF and to the improved HLRF (iHLRF) algorithms, in the solution of 25 benchmark problems from the literature. The new proposed HLRF algorithm is shown to be more robust than HLRF or iHLRF, and as efficient as the iHLRF algorithm. The two augmented Lagrangian methods proposed herein are shown to be more robust and more efficient than the classical augmented Lagrangian method.

Comparison of two experimental models of pulmonary hypertension

Objective: To compare two models of pulmonary hypertension (monocrotaline and monocrotaline+pneumonectomy) regarding hemodynamic severity, structure of pulmonary arteries, inflammatory markers (IL-1 and PDGF), and 45-day survival. Methods: We used 80 Sprague-Dawley rats in two study protocols: structural analysis; and survival analysis. The rats were divided into four groups: control; monocrotaline (M), pneumonectomy (P), and monocrotaline+pneumonectomy (M+P). In the structural analysis protocol, 40 rats (10/group) were catheterized for the determination of hemodynamic variables, followed by euthanasia for the removal of heart and lung tissue. The right ventricle (RV) was dissected from the interventricular septum (IS), and the ratio between RV weight and the weight of the left ventricle (LV) plus IS (RV/LV+IS) was taken as the index of RV hypertrophy. In lung tissues, we performed histological analyses, as well as using ELISA to determine IL-1 and PDGF levels. In the survival protocol, 40 animals (10/group) were followed for 45 days. Results: The M and M+P rats developed pulmonary hypertension, whereas the control and P rats did not. The RV/LV+IS ratio was significantly higher in M+P rats than in M rats, as well as being significantly higher in M and M+P rats than in control and P rats. There were no significant differences between the M and M+P rats regarding the area of the medial layer of the pulmonary arteries; IL-1 and PDGF levels; or survival. Conclusions: On the basis of our results, we cannot conclude that the monocrotaline+pneumonectomy model is superior to the monocrotaline model.

Purification, characterization and structural determination of chitinases produced by Moniliophthora perniciosa

The enzyme chitinase from Moniliophthora perniciosa the causative agent of the witches' broom disease in Theobroma cacao, was partially purified with ammonium sulfate and filtration by Sephacryl S-200 using sodium phosphate as an extraction buffer. Response surface methodology (RSM) was used to determine the optimum pH and temperature conditions. Four different isoenzymes were obtained: ChitMp I, ChitMp II, ChitMp III and ChitMp IV. ChitMp I had an optimum temperature at 44-73ºC and an optimum pH at 7.0-8.4. ChitMp II had an optimum temperature at 45-73ºC and an optimum pH at 7.0-8.4. ChitMp III had an optimum temperature at 54-67ºC and an optimum pH at 7.3-8.8. ChitMp IV had an optimum temperature at 60ºC and an optimum pH at 7.0. For the computational biology, the primary sequence was determined in silico from the database of the Genome/Proteome Project of M. perniciosa, yielding a sequence with 564 bp and 188 amino acids that was used for the three-dimensional design in a comparative modeling methodology. The generated models were submitted to validation using Procheck 3.0 and ANOLEA. The model proposed for the chitinase was subjected to a dynamic analysis over a 1 ns interval, resulting in a model with 91.7% of the residues occupying favorable places on the Ramachandran plot and an RMS of 2.68.

Effect of endogenous hydrogen sulfide inhibition on structural and functional renal disturbances induced by gentamicin

Animal models of gentamicin nephrotoxicity present acute tubular necrosis associated with inflammation, which can contribute to intensify the renal damage. Hydrogen sulfide (H2S) is a signaling molecule involved in inflammation. We evaluated the effect of DL-propargylglycine (PAG), an inhibitor of endogenous H2S formation, on the renal damage induced by gentamicin. Male Wistar rats (N = 8) were injected with 40 mg/kg gentamicin (im) twice a day for 9 days, some of them also received PAG (N = 8, 10 mg·kg-1·day-1, ip). Control rats (N = 6) were treated with saline or PAG only (N = 4). Twenty-four-hour urine samples were collected one day after the end of these treatments, blood samples were collected, the animals were sacrificed, and the kidneys were removed for quantification of H2S formation and histological and immunohistochemical studies. Gentamicin-treated rats presented higher sodium and potassium fractional excretion, increased plasma creatinine [4.06 (3.00; 5.87) mg%] and urea levels, a greater number of macrophages/monocytes, and a higher score for tubular interstitial lesions [3.50 (3.00; 4.00)] in the renal cortex. These changes were associated with increased H2S formation in the kidneys from gentamicin-treated rats (230.60 ± 38.62 µg·mg protein-1·h-1) compared to control (21.12 ± 1.63) and PAG (11.44 ± 3.08). Treatment with PAG reduced this increase (171.60 ± 18.34), the disturbances in plasma creatinine levels [2.20 (1.92; 4.60) mg%], macrophage infiltration, and score for tubular interstitial lesions [2.00 (2.00; 3.00)]. However, PAG did not interfere with the increase in fractional sodium excretion provoked by gentamicin. The protective effect of PAG on gentamicin nephrotoxicity was related, at least in part, to decreased H2S formation.

High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

Abstract Background The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery. Results Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer. Conclusions Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.

Statistical models for experimental model membranes.

Biological membranes are constituted from lipid bilayers and proteins. Investigation of protein-membrane interaction, essential for biological function of cells, must rest upon solid knowledge of lipid bilayer behavior. Thus, extensive studies of an experimental model for membranes, lipid bilayers in water solution, have been undertaken in the last decades. These systems present structural, thermal and electrical properties which depend on temperature, ionic strength or concentration. In this talk, we shall discuss statistical models for lipid bilayers, as well as the relation between their properties and results for properties of lipid dispersions investigated by the laboratories supervised by Teresa Lamy (IF-USP) and Amando Ito (FFCL-USP).

Structural stability of flexible lines in catenary configuration under torsion

Catenary risers can present during installation a very low tension close to seabed, which combined with torsion moment can lead to a structural instability, resulting in a loop. This is undesirable once it is possible that the loop turns into a kink, creating damage. This work presents a numerical methodology to analyze the conditions of loop formation in catenary risers. Stability criteria were applied to finite element models, including geometric nonlinearities and contact constraint due to riser-seabed interaction. The classical Greenhill's formula was used to predict the phenomenon and parametric analysis shows a “universal plot” able to predict instability in catenaries using a simple equation that can be applied for typical risers installation conditions and, generically, for catenary lines under torsion.

Application of 137Cs and 210Pb radionuclides to determine sedimentation rates of recent sediments from admiraty Bay, Antarctica Peninsula

Studies about natural and artificial radionuclides in areas such as the Antarctic are key to understand natural and dynamic processes in marine environments. These studies are important to determine levels of radioactive elements and local sedimentation rates. Five marine sediment cores were collected in different points of Admiralty Bay, in the Antarctic Peninsula. The purpose of this study was to determine 137Cs, 226Ra and 210Pb and sedimentation rates at each site. 137Cs, 210Pb and 226Ra were assayed by gamma-counting through direct measurement of the peak at 661 keV, 47 keV and 609 keV, respectively. Sedimentation rates were obtained by 137Cs and 210Pb (CIC and CRS). The activities for 137Cs ranged from 0.84 to 7.09 Bq kg-1; to 226Ra from 6.77 to 31.07 Bq kg-1 and for 210Pb ranged from 1.10 to 36.90 Bq kg-1. The sedimentation rates obtained by the three models ranged from 0.11±0.01 cm y-1 to 0.46±0.05 cm y-1. The levels of 137Cs registered in this study, as well as in other studies in the Antarctic region indicate that global fallout is the main cause of artificial radionuclides present in this environment, since the Antarctic has not suffered a direct action of human activities that released radioactive elements. The possible grain size variations that occur in the studied points of Admiralty Bay may explain the differences found in the vertical distribution of radionuclides, because of the different values of sedimentation rates and respective dating determined in their profiles

Inhibitors of Trypanosoma brucei trypanothione reductase: comparative molecular field analysis modeling and structural basis for selective inhibition

Background: Sleeping sickness is a major cause of death in Africa. Since no secure treatment is available, the development of novel therapeutic agents is urgent. In this context, the enzyme trypanothione reductase (TR) is a prominent molecular target that has been investigated in drug design for sleeping sickness. Results: In this study, comparative molecular field analysis models were generated for a series of Trypanosoma brucei TR inhibitors. Statistically significant results were obtained and the models were applied to predict the activity of external test sets, with good correlation between predicted and experimental results. We have also investigated the structural requirements for the selective inhibition of the parasite's enzyme over the human glutathione reductase. Conclusion: The quantitative structure-activity relationship models provided valuable information regarding the essential molecular requirements for the inhibitory activity upon the target protein, providing important insights into the design of more potent and selective TR inhibitors.

Stem Cells as a therapy for myocardial infarction in animal models

Advances in stem cell biology have challenged the notion that infarcted myocardium is irreparable. The pluripotent ability of stem cells to differentiate into specialized cell lines began to garner intense interest within cardiology when it was shown in animal models that intramyocardial injection of bone marrow stem cells (MSCs), or the mobilization of bone marrow stem cells with spontaneous homing to myocardium, could improve cardiac function and survival after induced myocardial infarction (MI) [1, 2]. Furthermore, the existence of stem cells in myocardium has been identified in animal heart [3, 4], and intense research is under way in an attempt to clarify their potential clinical application for patients with myocardial infarction. To date, in order to identify the best one, different kinds of stem cells have been studied; these have been derived from embryo or adult tissues (i.e. bone marrow, heart, peripheral blood etc.). Currently, three different biologic therapies for cardiovascular diseases are under investigation: cell therapy, gene therapy and the more recent “tissue-engineering” therapy . During my Ph.D. course, first I focalised my study on the isolation and characterization of Cardiac Stem Cells (CSCs) in wild-type and transgenic mice and for this purpose I attended, for more than one year, the Cardiovascular Research Institute of the New York Medical College, in Valhalla (NY, USA) under the direction of Doctor Piero Anversa. During this period I learnt different Immunohistochemical and Biomolecular techniques, useful for investigating the regenerative potential of stem cells. Then, during the next two years, I studied the new approach of cardiac regenerative medicine based on “tissue-engineering” in order to investigate a new strategy to regenerate the infracted myocardium. Tissue-engineering is a promising approach that makes possible the creation of new functional tissue to replace lost or failing tissue. This new discipline combines isolated functioning cells and biodegradable 3-dimensional (3D) polymeric scaffolds. The scaffold temporarily provides the biomechanical support for the cells until they produce their own extracellular matrix. Because tissue-engineering constructs contain living cells, they may have the potential for growth and cellular self-repair and remodeling. In the present study, I examined whether the tissue-engineering strategy within hyaluron-based scaffolds would result in the formation of alternative cardiac tissue that could replace the scar and improve cardiac function after MI in syngeneic heterotopic rat hearts. Rat hearts were explanted, subjected to left coronary descending artery occlusion, and then grafted into the abdomen (aorta-aorta anastomosis) of receiving syngeneic rat. After 2 weeks, a pouch of 3 mm2 was made in the thickness of the ventricular wall at the level of the post-infarction scar. The hyaluronic scaffold, previously engineered for 3 weeks with rat MSCs, was introduced into the pouch and the myocardial edges sutured with few stitches. Two weeks later we evaluated the cardiac function by M-Mode echocardiography and the myocardial morphology by microscope analysis. We chose bone marrow-derived mensenchymal stem cells (MSCs) because they have shown great signaling and regenerative properties when delivered to heart tissue following a myocardial infarction (MI). However, while the object of cell transplantation is to improve ventricular function, cardiac cell transplantation has had limited success because of poor graft viability and low cell retention, that’s why we decided to combine MSCs with a biopolimeric scaffold. At the end of the experiments we observed that the hyaluronan fibres had not been substantially degraded 2 weeks after heart-transplantation. Most MSCs had migrated to the surrounding infarcted area where they were especially found close to small-sized vessels. Scar tissue was moderated in the engrafted region and the thickness of the corresponding ventricular wall was comparable to that of the non-infarcted remote area. Also, the left ventricular shortening fraction, evaluated by M-Mode echocardiography, was found a little bit increased when compared to that measured just before construct transplantation. Therefore, this study suggests that post-infarction myocardial remodelling can be favourably affected by the grafting of MSCs delivered through a hyaluron-based scaffold

Probabilistic approaches in civil engineering: generation of random fields and structural identification with genetic algorithms

The inherent stochastic character of most of the physical quantities involved in engineering models has led to an always increasing interest for probabilistic analysis. Many approaches to stochastic analysis have been proposed. However, it is widely acknowledged that the only universal method available to solve accurately any kind of stochastic mechanics problem is Monte Carlo Simulation. One of the key parts in the implementation of this technique is the accurate and efficient generation of samples of the random processes and fields involved in the problem at hand. In the present thesis an original method for the simulation of homogeneous, multi-dimensional, multi-variate, non-Gaussian random fields is proposed. The algorithm has proved to be very accurate in matching both the target spectrum and the marginal probability. The computational efficiency and robustness are very good too, even when dealing with strongly non-Gaussian distributions. What is more, the resulting samples posses all the relevant, welldefined and desired properties of “translation fields”, including crossing rates and distributions of extremes. The topic of the second part of the thesis lies in the field of non-destructive parametric structural identification. Its objective is to evaluate the mechanical characteristics of constituent bars in existing truss structures, using static loads and strain measurements. In the cases of missing data and of damages that interest only a small portion of the bar, Genetic Algorithm have proved to be an effective tool to solve the problem.

Retrieval of trace gases vertical profile in the lower atmosphere combining. Differential Optical Absorption Spectroscopy with radiative transfer models

The motivation for the work presented in this thesis is to retrieve profile information for the atmospheric trace constituents nitrogen dioxide (NO2) and ozone (O3) in the lower troposphere from remote sensing measurements. The remote sensing technique used, referred to as Multiple AXis Differential Optical Absorption Spectroscopy (MAX-DOAS), is a recent technique that represents a significant advance on the well-established DOAS, especially for what it concerns the study of tropospheric trace consituents. NO2 is an important trace gas in the lower troposphere due to the fact that it is involved in the production of tropospheric ozone; ozone and nitrogen dioxide are key factors in determining the quality of air with consequences, for example, on human health and the growth of vegetation. To understand the NO2 and ozone chemistry in more detail not only the concentrations at ground but also the acquisition of the vertical distribution is necessary. In fact, the budget of nitrogen oxides and ozone in the atmosphere is determined both by local emissions and non-local chemical and dynamical processes (i.e. diffusion and transport at various scales) that greatly impact on their vertical and temporal distribution: thus a tool to resolve the vertical profile information is really important. Useful measurement techniques for atmospheric trace species should fulfill at least two main requirements. First, they must be sufficiently sensitive to detect the species under consideration at their ambient concentration levels. Second, they must be specific, which means that the results of the measurement of a particular species must be neither positively nor negatively influenced by any other trace species simultaneously present in the probed volume of air. Air monitoring by spectroscopic techniques has proven to be a very useful tool to fulfill these desirable requirements as well as a number of other important properties. During the last decades, many such instruments have been developed which are based on the absorption properties of the constituents in various regions of the electromagnetic spectrum, ranging from the far infrared to the ultraviolet. Among them, Differential Optical Absorption Spectroscopy (DOAS) has played an important role. DOAS is an established remote sensing technique for atmospheric trace gases probing, which identifies and quantifies the trace gases in the atmosphere taking advantage of their molecular absorption structures in the near UV and visible wavelengths of the electromagnetic spectrum (from 0.25 μm to 0.75 μm). Passive DOAS, in particular, can detect the presence of a trace gas in terms of its integrated concentration over the atmospheric path from the sun to the receiver (the so called slant column density). The receiver can be located at ground, as well as on board an aircraft or a satellite platform. Passive DOAS has, therefore, a flexible measurement configuration that allows multiple applications. The ability to properly interpret passive DOAS measurements of atmospheric constituents depends crucially on how well the optical path of light collected by the system is understood. This is because the final product of DOAS is the concentration of a particular species integrated along the path that radiation covers in the atmosphere. This path is not known a priori and can only be evaluated by Radiative Transfer Models (RTMs). These models are used to calculate the so called vertical column density of a given trace gas, which is obtained by dividing the measured slant column density to the so called air mass factor, which is used to quantify the enhancement of the light path length within the absorber layers. In the case of the standard DOAS set-up, in which radiation is collected along the vertical direction (zenith-sky DOAS), calculations of the air mass factor have been made using “simple” single scattering radiative transfer models. This configuration has its highest sensitivity in the stratosphere, in particular during twilight. This is the result of the large enhancement in stratospheric light path at dawn and dusk combined with a relatively short tropospheric path. In order to increase the sensitivity of the instrument towards tropospheric signals, measurements with the telescope pointing the horizon (offaxis DOAS) have to be performed. In this circumstances, the light path in the lower layers can become very long and necessitate the use of radiative transfer models including multiple scattering, the full treatment of atmospheric sphericity and refraction. In this thesis, a recent development in the well-established DOAS technique is described, referred to as Multiple AXis Differential Optical Absorption Spectroscopy (MAX-DOAS). The MAX-DOAS consists in the simultaneous use of several off-axis directions near the horizon: using this configuration, not only the sensitivity to tropospheric trace gases is greatly improved, but vertical profile information can also be retrieved by combining the simultaneous off-axis measurements with sophisticated RTM calculations and inversion techniques. In particular there is a need for a RTM which is capable of dealing with all the processes intervening along the light path, supporting all DOAS geometries used, and treating multiple scattering events with varying phase functions involved. To achieve these multiple goals a statistical approach based on the Monte Carlo technique should be used. A Monte Carlo RTM generates an ensemble of random photon paths between the light source and the detector, and uses these paths to reconstruct a remote sensing measurement. Within the present study, the Monte Carlo radiative transfer model PROMSAR (PROcessing of Multi-Scattered Atmospheric Radiation) has been developed and used to correctly interpret the slant column densities obtained from MAX-DOAS measurements. In order to derive the vertical concentration profile of a trace gas from its slant column measurement, the AMF is only one part in the quantitative retrieval process. One indispensable requirement is a robust approach to invert the measurements and obtain the unknown concentrations, the air mass factors being known. For this purpose, in the present thesis, we have used the Chahine relaxation method. Ground-based Multiple AXis DOAS, combined with appropriate radiative transfer models and inversion techniques, is a promising tool for atmospheric studies in the lower troposphere and boundary layer, including the retrieval of profile information with a good degree of vertical resolution. This thesis has presented an application of this powerful comprehensive tool for the study of a preserved natural Mediterranean area (the Castel Porziano Estate, located 20 km South-West of Rome) where pollution is transported from remote sources. Application of this tool in densely populated or industrial areas is beginning to look particularly fruitful and represents an important subject for future studies.