The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 1: Model description and evaluation

We introduce the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS). CATT-BRAMS is an on-line transport model fully consistent with the simulated atmospheric dynamics. Emission sources from biomass burning and urban-industrial-vehicular activities for trace gases and from biomass burning aerosol particles are obtained from several published datasets and remote sensing information. The tracer and aerosol mass concentration prognostics include the effects of sub-grid scale turbulence in the planetary boundary layer, convective transport by shallow and deep moist convection, wet and dry deposition, and plume rise associated with vegetation fires in addition to the grid scale transport. The radiation parameterization takes into account the interaction between the simulated biomass burning aerosol particles and short and long wave radiation. The atmospheric model BRAMS is based on the Regional Atmospheric Modeling System (RAMS), with several improvements associated with cumulus convection representation, soil moisture initialization and surface scheme tuned for the tropics, among others. In this paper the CATT-BRAMS model is used to simulate carbon monoxide and particulate material (PM(2.5)) surface fluxes and atmospheric transport during the 2002 LBA field campaigns, conducted during the transition from the dry to wet season in the southwest Amazon Basin. Model evaluation is addressed with comparisons between model results and near surface, radiosondes and airborne measurements performed during the field campaign, as well as remote sensing derived products. We show the matching of emissions strengths to observed carbon monoxide in the LBA campaign. A relatively good comparison to the MOPITT data, in spite of the fact that MOPITT a priori assumptions imply several difficulties, is also obtained.

Continental passive microwave-based rainfall estimation algorithm: Application to the Amazon Basin

This paper presents a new statistical algorithm to estimate rainfall over the Amazon Basin region using the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI). The algorithm relies on empirical relationships derived for different raining-type systems between coincident measurements of surface rainfall rate and 85-GHz polarization-corrected brightness temperature as observed by the precipitation radar (PR) and TMI on board the TRMM satellite. The scheme includes rain/no-rain area delineation (screening) and system-type classification routines for rain retrieval. The algorithm is validated against independent measurements of the TRMM-PR and S-band dual-polarization Doppler radar (S-Pol) surface rainfall data for two different periods. Moreover, the performance of this rainfall estimation technique is evaluated against well-known methods, namely, the TRMM-2A12 [ the Goddard profiling algorithm (GPROF)], the Goddard scattering algorithm (GSCAT), and the National Environmental Satellite, Data, and Information Service (NESDIS) algorithms. The proposed algorithm shows a normalized bias of approximately 23% for both PR and S-Pol ground truth datasets and a mean error of 0.244 mm h(-1) ( PR) and -0.157 mm h(-1)(S-Pol). For rain volume estimates using PR as reference, a correlation coefficient of 0.939 and a normalized bias of 0.039 were found. With respect to rainfall distributions and rain area comparisons, the results showed that the formulation proposed is efficient and compatible with the physics and dynamics of the observed systems over the area of interest. The performance of the other algorithms showed that GSCAT presented low normalized bias for rain areas and rain volume [0.346 ( PR) and 0.361 (S-Pol)], and GPROF showed rainfall distribution similar to that of the PR and S-Pol but with a bimodal distribution. Last, the five algorithms were evaluated during the TRMM-Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA) 1999 field campaign to verify the precipitation characteristics observed during the easterly and westerly Amazon wind flow regimes. The proposed algorithm presented a cumulative rainfall distribution similar to the observations during the easterly regime, but it underestimated for the westerly period for rainfall rates above 5 mm h(-1). NESDIS(1) overestimated for both wind regimes but presented the best westerly representation. NESDIS(2), GSCAT, and GPROF underestimated in both regimes, but GPROF was closer to the observations during the easterly flow.

Newtonian perturbations on models with matter creation

Creation of cold dark matter (CCDM) can macroscopically be described by a negative pressure, and, therefore, the mechanism is capable to accelerate the Universe, without the need of an additional dark energy component. In this framework, we discuss the evolution of perturbations by considering a Neo-Newtonian approach where, unlike in the standard Newtonian cosmology, the fluid pressure is taken into account even in the homogeneous and isotropic background equations (Lima, Zanchin, and Brandenberger, MNRAS 291, L1, 1997). The evolution of the density contrast is calculated in the linear approximation and compared to the one predicted by the Lambda CDM model. The difference between the CCDM and Lambda CDM predictions at the perturbative level is quantified by using three different statistical methods, namely: a simple chi(2)-analysis in the relevant space parameter, a Bayesian statistical inference, and, finally, a Kolmogorov-Smirnov test. We find that under certain circumstances, the CCDM scenario analyzed here predicts an overall dynamics (including Hubble flow and matter fluctuation field) which fully recovers that of the traditional cosmic concordance model. Our basic conclusion is that such a reduction of the dark sector provides a viable alternative description to the accelerating Lambda CDM cosmology.

Chemical potential and the nature of dark energy: The case of a phantom field

The influence of a possible nonzero chemical potential mu on the nature of dark energy is investigated by assuming that the dark energy is a relativistic perfect simple fluid obeying the equation of state, p=omega rho (omega < 0, constant). The entropy condition, S >= 0, implies that the possible values of omega are heavily dependent on the magnitude, as well as on the sign of the chemical potential. For mu > 0, the omega parameter must be greater than -1 (vacuum is forbidden) while for mu < 0 not only the vacuum but even a phantomlike behavior (omega <-1) is allowed. In any case, the ratio between the chemical potential and temperature remains constant, that is, mu/T=mu(0)/T(0). Assuming that the dark energy constituents have either a bosonic or fermionic nature, the general form of the spectrum is also proposed. For bosons mu is always negative and the extended Wien's law allows only a dark component with omega <-1/2, which includes vacuum and the phantomlike cases. The same happens in the fermionic branch for mu < 0. However, fermionic particles with mu > 0 are permitted only if -1 dynamics and statistical arguments constrain the equation-of-state parameter to be omega <-1/2, a result surprisingly close to the maximal value required to accelerate a Friedmann-Robertson-Walker-type universe dominated by matter and dark energy (omega less than or similar to-10/21).

NGC 7097: THE ACTIVE GALACTIC NUCLEUS AND ITS MIRROR, REVEALED BY PRINCIPAL COMPONENT ANALYSIS TOMOGRAPHY

Three-dimensional spectroscopy techniques are becoming more and more popular, producing an increasing number of large data cubes. The challenge of extracting information from these cubes requires the development of new techniques for data processing and analysis. We apply the recently developed technique of principal component analysis (PCA) tomography to a data cube from the center of the elliptical galaxy NGC 7097 and show that this technique is effective in decomposing the data into physically interpretable information. We find that the first five principal components of our data are associated with distinct physical characteristics. In particular, we detect a low-ionization nuclear-emitting region (LINER) with a weak broad component in the Balmer lines. Two images of the LINER are present in our data, one seen through a disk of gas and dust, and the other after scattering by free electrons and/or dust particles in the ionization cone. Furthermore, we extract the spectrum of the LINER, decontaminated from stellar and extended nebular emission, using only the technique of PCA tomography. We anticipate that the scattered image has polarized light due to its scattered nature.

The metal content of bulge field stars from FLAMES-GIRAFFE spectra - I. Stellar parameters and iron abundances

Aims. We determine the iron distribution function (IDF) for bulge field stars, in three different fields along the Galactic minor axis and at latitudes b = -4 degrees, b = -6 degrees, and b = -12 degrees. A fourth field including NGC 6553 is also included in the discussion. Methods. About 800 bulge field K giants were observed with the GIRAFFE spectrograph of FLAMES@VLT at spectral resolution R similar to 20 000. Several of them were observed again with UVES at R similar to 45 000 to insure the accuracy of the measurements. The LTE abundance analysis yielded stellar parameters and iron abundances that allowed us to construct an IDF for the bulge that, for the first time, is based on high-resolution spectroscopy for each individual star. Results. The IDF derived here is centered on solar metallicity, and extends from [Fe/H] similar to -1.5 to [Fe/H] similar to + 0.5. The distribution is asymmetric, with a sharper cutoff on the high-metallicity side, and it is narrower than previously measured. A variation in the mean metallicity along the bulge minor axis is clearly between b = -4 degrees and b = -6 degrees ([Fe/H] decreasing similar to by 0.6 dex per kpc). The field at b = -12 degrees. is consistent with the presence of a gradient, but its quantification is complicated by the higher disk/bulge fraction in this field. Conclusions. Our findings support a scenario in which both infall and outflow were important during the bulge formation, and then suggest the presence of a radial gradient, which poses some challenges to the scenario in which the bulge would result solely from the vertical heating of the bar.

Trapping of strangelets in the geomagnetic field

Strangelets arriving from the interstellar medium are an interesting target for experiments searching for evidence of this hypothetical state of hadronic matter. We entertain the possibility of a trapped strangelet population, quite analogous to ordinary nuclei and electron belts. For a population of strangelets to be trapped by the geomagnetic field, these incoming particles would have to fulfill certain conditions, namely, having magnetic rigidities above the geomagnetic cutoff and below a certain threshold for adiabatic motion to hold. We show in this work that, for fully ionized strangelets, there is a narrow window for stable trapping. An estimate of the stationary population is presented and the dominant loss mechanisms discussed. It is shown that the population would be substantially enhanced with respect to the interstellar medium flux (up to 2 orders of magnitude) due to quasistable trapping.

Dynamics and constraints of the massive graviton dark matter flat cosmologies

We discuss the dynamics of the Universe within the framework of the massive graviton cold dark matter scenario (MGCDM) in which gravitons are geometrically treated as massive particles. In this modified gravity theory, the main effect of the gravitons is to alter the density evolution of the cold dark matter component in such a way that the Universe evolves to an accelerating expanding regime, as presently observed. Tight constraints on the main cosmological parameters of the MGCDM model are derived by performing a joint likelihood analysis involving the recent supernovae type Ia data, the cosmic microwave background shift parameter, and the baryonic acoustic oscillations as traced by the Sloan Digital Sky Survey red luminous galaxies. The linear evolution of small density fluctuations is also analyzed in detail. It is found that the growth factor of the MGCDM model is slightly different (similar to 1-4%) from the one provided by the conventional flat Lambda CDM cosmology. The growth rate of clustering predicted by MGCDM and Lambda CDM models are confronted to the observations and the corresponding best fit values of the growth index (gamma) are also determined. By using the expectations of realistic future x-ray and Sunyaev-Zeldovich cluster surveys we derive the dark matter halo mass function and the corresponding redshift distribution of cluster-size halos for the MGCDM model. Finally, we also show that the Hubble flow differences between the MGCDM and the Lambda CDM models provide a halo redshift distribution departing significantly from the those predicted by other dark energy models. These results suggest that the MGCDM model can observationally be distinguished from Lambda CDM and also from a large number of dark energy models recently proposed in the literature.

Kinematics of galaxies in compact groups Studying the B-band Tully-Fischer relation

We obtained new Fabry-Perot data cubes and derived velocity fields, monochromatic, and velocity dispersion maps for 28 galaxies in the Hickson compact groups 37, 40, 47, 49, 54, 56, 68, 79, and 93. We also derived rotation curves for 9 of the studied galaxies, 6 of which are strongly asymmetric. Combining these new data with previously published 2D kinematic maps of compact group galaxies, we investigated the differences between the kinematic and morphological position angles for a sample of 46 galaxies. We find that one third of the unbarred compact group galaxies have position angle misalignments between the stellar and gaseous components. This and the asymmetric rotation curves are clear signatures of kinematic perturbations, probably because of interactions among compact group galaxies. A comparison between the B-band Tully-Fisher relation for compact group galaxies and for the GHASP field-galaxy sample shows that, despite the high fraction of compact group galaxies with asymmetric rotation curves, these lay on the TF relation defined by galaxies in less dense environments, although with more scatter. This agrees with previous results, but now confirmed for a larger sample of 41 galaxies. We confirm the tendency for compact group galaxies at the low-mass end of the Tully-Fisher relation (HCG 49b, 89d, 96c, 96d, and 100c) to have either a magnitude that is too bright for its mass (suggesting brightening by star formation) and/or a low maximum rotational velocity for its luminosity (suggesting tidal stripping). These galaxies are outside the Tully Fisher relation at the 1 sigma level, even when the minimum acceptable values of inclinations are used to compute their maximum velocities. Including such galaxies with nu < 100 km s(-1) in the determination of the zero point and slope of the compact group B-band Tully-Fisher relation would strongly change the fit, making it different from the relation for field galaxies, which has to be kept in mind when studying scaling relations of interacting galaxies, especially at high redshifts.

Insights on the Milky Way bulge formation from the correlations between kinematics and metallicity

Context. Two main scenarios for the formation of the Galactic bulge are invoked, the first one through gravitational collapse or hierarchical merging of subclumps, the second through secular evolution of the Galactic disc. Aims. We aim to constrain the formation of the Galactic bulge through studies of the correlation between kinematics and metallicities in Baade's Window (l = 1 degrees, b = -4 degrees) and two other fields along the bulge minor axis (l = 0 degrees, b = -6 degrees and b = -12 degrees). Methods. We combine the radial velocity and the [Fe/H] measurements obtained with FLAMES/GIRAFFE at the VLT with a spectral resolution of R = 20 000, plus for the Baade's Window field the OGLE-II proper motions, and compare these with published N-body simulations of the Galactic bulge. Results. We confirm the presence of two distinct populations in Baade's Window found in Hill et al. (2010, A&A, submitted): the metal-rich population presents bar-like kinematics while the metal-poor population shows kinematics corresponding to an old spheroid or a thick disc. In this context the metallicity gradient along the bulge minor axis observed by Zoccali et al. (2008, A&A, 486, 177), visible also in the kinematics, can be related to a varying mix of these two populations as one moves away from the Galactic plane, alleviating the apparent contradiction between the kinematic evidence of a bar and the existence of a metallicity gradient. Conclusions. We show evidence that the two main scenarios for the bulge formation co-exist within the Milky Way bulge.

MASS DISTRIBUTION IN HICKSON COMPACT GROUPS OF GALAXIES

This study presents the mass distribution for a sample of 18 late-type galaxies in nine Hickson compact groups. We used Ha rotation curves (RCs) from high-resolution two-dimensional velocity fields of Fabry-Perot observations and the J-band photometry from the Two Micron All Sky Survey, in order to determine the dark halo and the visible matter distributions. The study compares two halo density profiles, an isothermal core-like distribution, and a cuspy one. We also compare their visible and dark matter distributions with those of galaxies belonging to cluster and field galaxies coming from two samples: 40 cluster galaxies of Barnes et al. and 35 field galaxies of Spano et al. The central halo surface density is found to be constant with respect to the total absolute magnitude similar to what is found for the isolated galaxies. This suggests that the halo density is independent of galaxy type and environment. We have found that core-like density profiles better fit the RCs than cuspy-like ones. No major differences have been found between field, cluster, and compact group galaxies with respect to their dark halo density profiles.

Star formation history of Canis Major R1 I. Wide-Field X-ray study of the young stellar population

Aims. The CMa R1 star-forming region contains several compact clusters as well as many young early-B stars. It is associated with a well-known bright rimmed nebula, the nature of which is unclear (fossil HII region or supernova remnant). To help elucidate the nature of the nebula, our goal was to reconstruct the star-formation history of the CMa R1 region, including the previously unknown older, fainter low-mass stellar population, using X-rays. Methods. We analyzed images obtained with the ROSAT satellite, covering similar to 5 sq. deg. Complementary VRI photometry was performed with the Gemini South telescope. Colour-magnitude and colour-colour diagrams were used in conjunction with pre-main sequence evolutionary tracks to derive the masses and ages of the X-ray sources. Results. The ROSAT images show two distinct clusters. One is associated with the known optical clusters near Z CMa, to which similar to 40 members are added. The other, which we name the ""GU CMa"" cluster, is new, and contains similar to 60 members. The ROSAT sources are young stars with masses down to M(star) similar to 0.5 M(circle dot), and ages up to 10 Myr. The mass functions of the two clusters are similar, but the GU CMa cluster is older than the cluster around Z CMa by at least a few Myr. Also, the GU CMa cluster is away from any molecular cloud, implying that star formation must have ceased; on the contrary (as already known), star formation is very active in the Z CMa region.

Planetary transit candidates in Corot-IRa01 field

Context. CoRoT is a pioneering space mission devoted to the analysis of stellar variability and the photometric detection of extrasolar planets. Aims. We present the list of planetary transit candidates detected in the first field observed by CoRoT, IRa01, the initial run toward the Galactic anticenter, which lasted for 60 days. Methods. We analysed 3898 sources in the coloured bands and 5974 in the monochromatic band. Instrumental noise and stellar variability were taken into account using detrending tools before applying various transit search algorithms. Results. Fifty sources were classified as planetary transit candidates and the most reliable 40 detections were declared targets for follow-up ground-based observations. Two of these targets have so far been confirmed as planets, CoRoT-1b and CoRoT-4b, for which a complete characterization and specific studies were performed.

Social Networks Shape the Transmission Dynamics of Hepatitis C Virus

Hepatitis C virus (HCV) infects 170 million people worldwide, and is a major public health problem in Brazil, where over 1% of the population may be infected and where multiple viral genotypes co-circulate. Chronically infected individuals are both the source of transmission to others and are at risk for HCV-related diseases, such as liver cancer and cirrhosis. Before the adoption of anti-HCV control measures in blood banks, this virus was mainly transmitted via blood transfusion. Today, needle sharing among injecting drug users is the most common form of HCV transmission. Of particular importance is that HCV prevalence is growing in non-risk groups. Since there is no vaccine against HCV, it is important to determine the factors that control viral transmission in order to develop more efficient control measures. However, despite the health costs associated with HCV, the factors that determine the spread of virus at the epidemiological scale are often poorly understood. Here, we sequenced partial NS5b gene sequences sampled from blood samples collected from 591 patients in Sao Paulo state, Brazil. We show that different viral genotypes entered Sao Paulo at different times, grew at different rates, and are associated with different age groups and risk behaviors. In particular, subtype 1b is older and grew more slowly than subtypes 1a and 3a, and is associated with multiple age classes. In contrast, subtypes 1a and 3b are associated with younger people infected more recently, possibly with higher rates of sexual transmission. The transmission dynamics of HCV in Sao Paulo therefore vary by subtype and are determined by a combination of age, risk exposure and underlying social network. We conclude that social factors may play a key role in determining the rate and pattern of HCV spread, and should influence future intervention policies.

Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

Background: The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized. Results: We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for similar to 40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance. Conclusion: These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.