Wavelets And Adaptive Grids For The Discontinuous Galerkin Method

In this paper, space adaptivity is introduced to control the error in the numerical solution of hyperbolic systems of conservation laws. The reference numerical scheme is a new version of the discontinuous Galerkin method, which uses an implicit diffusive term in the direction of the streamlines, for stability purposes. The decision whether to refine or to unrefine the grid in a certain location is taken according to the magnitude of wavelet coefficients, which are indicators of local smoothness of the numerical solution. Numerical solutions of the nonlinear Euler equations illustrate the efficiency of the method. © Springer 2005.

Classificação de sistemas meteorológicos e comparação da precipitação estimada pelo radar e medida pela rede telemétrica na bacia hidrográfica do alto Tietê

Foram analisadas características da precipitação estimada a partir de 145.194 campos de refletividade, de um total de 827 dias entre 1998 e 2003, obtidos do Radar Meteorológico de São Paulo (RSP). Os eventos foram classificados de acordo com intensidades de precipitação; em Convectivos (EC) e Estratiformes (EE). Quanto à morfologia, cinco tipos de sistemas foram identificados; Convecção Isolada (CI), Brisa Marítima (BM), Linhas de Instabilidade (LI), Bandas Dispersas (BD) e Frentes Frias (FF). Eventos convectivos dominam na primavera e verão e estratiformes no outono e inverno. A CI e a BM tiveram maiores picos de atuação entre outubro e março enquanto as FF de abril a setembro. BD atuam durante todo o ano e as LI só não foram observadas nos meses de junho e julho. Uma comparação pontual entre a precipitação medida pela telemetria e estimada com o radar foi realizada e, mostrou haver, na maioria dos casos, um viés positivo do RSP, para acumulações de 10, 30 e 60 minutos. Com o objetivo de integrar as estimativas de precipitação do radar com as medidas da rede telemétrica, por meio de uma análise objetiva estatística, foram obtidas dos campos de precipitação do radar as estruturas das correlações espaciais em função da distância para acumulações de chuva de 15, 30, 60 e 120 minutos para os cinco tipos de sistemas precipitantes que foram caracterizados. As curvas das correlações espaciais médias de todos os eventos de precipitação de cada sistema foram ajustadas por funções polinomiais de sexta ordem. Os resultados indicam diferenças significativas na estrutura espacial das correlações entre os sistemas precipitantes.

Propagacão das ondas de rossby nos invernos de máxima freqüência de ocorrência de geadas na pampa úmida

Neste trabalho estudou-se a influência dos padrões de onda extratropicais, que favorecem o desenvolvimento de eventos extremos frios no sudeste Sul-Americano, e em particular na região conhecida como Pampa Úmida. O aquecimento anômalo observado na região do oceano Pacífico tropical ocidental a nordeste da Austrália, durante os invernos de máxima freqüência de ocorrência de Geadas Generalizadas (GG) no centro-leste da Argentina, (região conhecida como Pampa Úmida - PU), atua como disparador de ondas de Rossby, as quais se propagam até o continente, favorecendo assim a ocorrência daqueles eventos. O padrão de propagação obtido nas simulações numéricas com um modelo baroclínico global, mostra o predomínio de um número de onda 3. Adicionalmente, foram analisadas as correlações do vento meridional em altos e baixos níveis observados para os eventos de GG, selecionados dentro dos invernos de máxima freqüência de ocorrência desses eventos. O vento meridional global em 250hPa apresenta regiões com correlação estatisticamente significativa com o vento meridional médio na PU. A configuração obtida no caso do vento meridional global em 250hPa, correlacionado com o vento meridional na PU, pode estar associada ao padrão de propagação das ondas simuladas numericamente a partir da forçante tropical. Igualmente importantes e significativos são os valores de correlação do vento sul nos baixos níveis, em particular para toda região da PU. O padrão de ondas simulado está bem representado pelas significativas correlações entre o vento meridional hemisférico em altos níveis e a temperatura no dia de evento de GG.

Convecção linearmente organizada na área de Petrolina, semi-árido do Nordeste do Brasil: aspectos em meso e grande escala

É apresentado um estudo sobre sistemas convectivos linearmente organizados e observados por um radar meteorológico banda-C na região semi-árida do Nordeste do Brasil. São analisados três dias (27 a 29) de março de 1985, com ênfase na investigação do papel desempenhado por fatores locais e de grande escala no desenvolvimento dos sistemas. No cenário de grande escala, a área de cobertura do radar foi influenciada por um cavado de ar superior austral no dia 27 e por um vórtice ciclônico de altos níveis no dia 29. A convergência de umidade próxima à superfície favoreceu a atividade convectiva nos dias 27 e 29, enquanto que divergência de umidade próxima à superfície inibiu a atividade convectiva no dia 28. No cenário de mesoescala, foi observado que o aquecimento diurno é um fator importante para a formação de células convectivas, somando-se a ele o papel determinante da orografia na localização dos ecos. De maneira geral, as imagens de radar mostram os sistemas convectivos linearmente organizados em áreas elevadas e núcleos convectivos intensos envolvidos por uma área de precipitação estratiforme. Os resultados indicam que convergência do fluxo de umidade em grande escala e aquecimento radiativo, são fatores determinantes na evolução e desenvolvimento dos ecos na área de estudo.

Contribuição ao estudo do magmatismo Paleo e Mesoproterozóico do SW do Craton Amazônico através da aplicação de isótopos estáveis de O, H e S

Os isótopos estáveis de O, H e S foram utilizados para investigar a origem das rochas magmáticas nos Terrenos Jauru e Pontes e Lacerda do SW do Craton Amazônico, estado de Mato Grosso, Brasil. No Terreno Jauru as rochas granitóides do Greenstone belt Alto Jauru e da Suíte Cachoeirinha apresentam valores de δ18O entre +9,0‰ e +6,3‰ que indicam derivação a partir de magmas juvenis. Na Suíte Intrusiva Rio Branco valores de δ18O para rochas básicas estão entre +5,4‰ e +5,8‰ e para rochas félsicas entre +8,7‰ e +9,0‰; rochas intermediárias apresentam valores entre +7,3‰ e +8,3‰. Os valores mais baixos de δ18O, obtidos nas rochas básicas, são compatíveis com derivação mantélica, porém as rochas félsicas apresentam valores de δ18O compatíveis com origem crustais. Análises de isótopos estáveis de H (rocha total) forneceram valores de δD entre - 83‰ e -92‰, diferente das assinaturas de rochas metamórficas e de águas meteóricas. Resultados em sulfetos para isótopos estáveis de S em rochas básicas e intermediárias desta suíte apresentam valores de δ34S coerentes com uma fonte mantélica (entre + 0,7‰ e +3,8‰), enquanto os valores de δ34S (entre +5,2‰ e +6,1‰) obtidos nas rochas félsicas sugerem participação crustal na sua gênese. Na Suíte Santa Helena (Terreno Pontes e Lacerda) os resultados obtidos para δ18O se agrupam entre +4,4‰ e +8,9‰ indicando uma origem mantélica. O presente estudo confirma a importância da aplicação de isótopos estáveis para a compreensão de processos magmáticos e evolução crustal.

Anomalous thermodynamic properties in ferropericlase throughout its spin crossover transition

The thermodynamics properties of ferropericlase (Mg(1-x)Fe(x)O where x=0.1875) (Fp) throughout its spin crossover were investigated by first principles. Fp was treated as an ideal solid solution of pure high-spin and low-spin states. The Gibbs free energies of the pure states were addressed using the LDA+U method. A vibrational virtual-crystal model was developed to address the vibrational properties of the pure spin cases and used in conjunction with quasiharmonic theory to compute their vibrational free energies. The thermodynamics properties of Fp display significant anomalies that should be typical of spin crossover systems in general. In Fp, in particular, they are fundamental for understanding the state of earth's interior, where the pressure and temperature conditions of the crossover are realized.

Quasiharmonic elastic constants corrected for deviatoric thermal stresses

The quasiharmonic approximation (QHA), in its simplest form also called the statically constrained (SC) QHA, has been shown to be a straightforward method to compute thermoelastic properties of crystals. Recently we showed that for noncubic solids SC-QHA calculations develop deviatoric thermal stresses at high temperatures. Relaxation of these stresses leads to a series of corrections to the free energy that may be taken to any desired order, up to self-consistency. Here we show how to correct the elastic constants obtained using the SC-QHA. We exemplify the procedure by correcting to first order the elastic constants of MgSiO(3) perovskite and MgSiO(3) postperovskite, the major phases of the Earth's lower mantle. We show that this first-order correction is quite satisfactory for obtaining the aggregated elastic averages of these minerals and their velocities in the lower mantle. This type of correction is also shown to be applicable to experimental measurements of elastic constants in situations where deviatoric stresses can develop, such as in diamond-anvil cells.

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.

Turbulent magnetic pumping in a Babcock-Leighton solar dynamo model

Context. The turbulent pumping effect corresponds to the transport of magnetic flux due to the presence of density and turbulence gradients in convectively unstable layers. In the induction equation it appears as an advective term and for this reason it is expected to be important in the solar and stellar dynamo processes. Aims. We explore the effects of turbulent pumping in a flux-dominated Babcock-Leighton solar dynamo model with a solar-like rotation law. Methods. As a first step, only vertical pumping has been considered through the inclusion of a radial diamagnetic term in the induction equation. In the second step, a latitudinal pumping term was included and then, a near-surface shear was included. Results. The results reveal the importance of the pumping mechanism in solving current limitations in mean field dynamo modeling, such as the storage of the magnetic flux and the latitudinal distribution of the sunspots. If a meridional flow is assumed to be present only in the upper part of the convective zone, it is the full turbulent pumping that regulates both the period of the solar cycle and the latitudinal distribution of the sunspot activity. In models that consider shear near the surface, a second shell of toroidal field is generated above r = 0.95 R(circle dot) at all latitudes. If the full pumping is also included, the polar toroidal fields are efficiently advected inwards, and the toroidal magnetic activity survives only at the observed latitudes near the equator. With regard to the parity of the magnetic field, only models that combine turbulent pumping with near-surface shear always converge to the dipolar parity. Conclusions. This result suggests that, under the Babcock-Leighton approach, the equartorward motion of the observed magnetic activity is governed by the latitudinal pumping of the toroidal magnetic field rather than by a large scale coherent meridional flow. Our results support the idea that the parity problem is related to the quadrupolar imprint of the meridional flow on the poloidal component of the magnetic field and the turbulent pumping positively contributes to wash out this imprint.

The viscosity parameter's dependence on the profile of the disk scale height

It is shown that, for accretion disks, the height scale is a constant whenever hydrostatic equilibrium and the subsonic turbulence regime hold in the disk. In order to have a variable height scale, processes are needed that contribute an extra term to the continuity equation. This contribution makes the viscosity parameter much greater in the outer region and much smaller in the inner region. Under these circumstances, turbulence is the presumable source of viscosity in the disk.

Beryllium abundances along the evolutionary sequence of the open cluster IC 4651-A new test for hydrodynamical stellar models

Context. Previous analyses of lithium abundances in main sequence and red giant stars have revealed the action of mixing mechanisms other than convection in stellar interiors. Beryllium abundances in stars with Li abundance determinations can offer valuable complementary information on the nature of these mechanisms. Aims. Our aim is to derive Be abundances along the whole evolutionary sequence of an open cluster. We focus on the well-studied open cluster IC 4651. These Be abundances are used with previously determined Li abundances, in the same sample stars, to investigate the mixing mechanisms in a range of stellar masses and evolutionary stages. Methods. Atmospheric parameters were adopted from a previous abundance analysis by the same authors. New Be abundances have been determined from high-resolution, high signal-to-noise UVES spectra using spectrum synthesis and model atmospheres. The careful synthetic modeling of the Be lines region is used to calculate reliable abundances in rapidly rotating stars. The observed behavior of Be and Li is compared to theoretical predictions from stellar models including rotation-induced mixing, internal gravity waves, atomic diffusion, and thermohaline mixing. Results. Beryllium is detected in all the main sequence and turn-off sample stars, both slow- and fast-rotating stars, including the Li-dip stars, but is not detected in the red giants. Confirming previous results, we find that the Li dip is also a Be dip, although the depletion of Be is more modest than for Li in the corresponding effective temperature range. For post-main-sequence stars, the Be dilution starts earlier within the Hertzsprung gap than expected from classical predictions, as does the Li dilution. A clear dispersion in the Be abundances is also observed. Theoretical stellar models including the hydrodynamical transport processes mentioned above are able to reproduce all the observed features well. These results show a good theoretical understanding of the Li and Be behavior along the color-magnitude diagram of this intermediate-age cluster for stars more massive than 1.2 M(circle dot).

The close circumstellar environment of Betelgeuse II. Diffraction-limited spectro-imaging from 7.76 to 19.50 mu m with VLT/VISIR

Context. Mass-loss occurring in red supergiants (RSGs) is a major contributor to the enrichment of the interstellar medium in dust and molecules. The physical mechanism of this mass loss is however relatively poorly known. Betelgeuse is the nearest RSG, and as such a prime object for high angular resolution observations of its surface (by interferometry) and close circumstellar environment. Aims. The goal of our program is to understand how the material expelled from Betelgeuse is transported from its surface to the interstellar medium, and how it evolves chemically in this process. Methods. We obtained diffraction-limited images of Betelgeuse and a calibrator (Aldebaran) in six filters in the N band (7.76 to 12.81 mu m) and two filters in the Q band (17.65 and 19.50 mu m), using the VLT/VISIR instrument. Results. Our images show a bright, extended and complex circumstellar envelope at all wavelengths. It is particularly prominent longwards of approximate to 9-10 mu m, pointing at the presence of O-rich dust, such as silicates or alumina. A partial circular shell is observed between 0.5 and 1.0 '' from the star, and could correspond to the inner radius of the dust envelope. Several knots and filamentary structures are identified in the nebula. One of the knots, located at a distance of 0.9 '' west of the star, is particularly bright and compact. Conclusions. The circumstellar envelope around Betelgeuse extends at least up to several tens of stellar radii. Its relatively high degree of clumpiness indicates an inhomogeneous spatial distribution of the material lost by the star. Its extension corresponds to an important intermediate scale, where most of the dust is probably formed, between the hot and compact gaseous envelope observed previously in the near infrared and the interstellar medium.

Reduction of chaotic particle transport driven by drift waves in sheared flows

Investigations of chaotic particle transport by drift waves propagating in the edge plasma of tokamaks with poloidal zonal flow are described. For large aspect ratio tokamaks, the influence of radial electric field profiles on convective cells and transport barriers, created by the nonlinear interaction between the poloidal flow and resonant waves, is investigated. For equilibria with edge shear flow, particle transport is seen to be reduced when the electric field shear is reversed. The transport reduction is attributed to the robust invariant tori that occur in nontwist Hamiltonian systems. This mechanism is proposed as an explanation for the transport reduction in Tokamak Chauffage Alfven Bresilien [R. M. O. Galvao , Plasma Phys. Controlled Fusion 43, 1181 (2001)] for discharges with a biased electrode at the plasma edge.

The impact of the subtropical South Atlantic SST on South American precipitation

The Community Climate Model (CCM3) from the National Center for Atmospheric Research (NCAR) is used to investigate the effect of the South Atlantic sea surface temperature (SST) anomalies on interannual to decadal variability of South American precipitation. Two ensembles composed of multidecadal simulations forced with monthly SST data from the Hadley Centre for the period 1949 to 2001 are analysed. A statistical treatment based on signal-to-noise ratio and Empirical Orthogonal Functions (EOF) is applied to the ensembles in order to reduce the internal variability among the integrations. The ensemble treatment shows a spatial and temporal dependence of reproducibility. High degree of reproducibility is found in the tropics while the extratropics is apparently less reproducible. Austral autumn (MAM) and spring (SON) precipitation appears to be more reproducible over the South America-South Atlantic region than the summer (DJF) and winter (JJA) rainfall. While the Inter-tropical Convergence Zone (ITCZ) region is dominated by external variance, the South Atlantic Convergence Zone (SACZ) over South America is predominantly determined by internal variance, which makes it a difficult phenomenon to predict. Alternatively, the SACZ over western South Atlantic appears to be more sensitive to the subtropical SST anomalies than over the continent. An attempt is made to separate the atmospheric response forced by the South Atlantic SST anomalies from that associated with the El Nino - Southern Oscillation (ENSO). Results show that both the South Atlantic and Pacific SSTs modulate the intensity and position of the SACZ during DJF. Particularly, the subtropical South Atlantic SSTs are more important than ENSO in determining the position of the SACZ over the southeast Brazilian coast during DJF. On the other hand, the ENSO signal seems to influence the intensity of the SACZ not only in DJF but especially its oceanic branch during MAM. Both local and remote influences, however, are confounded by the large internal variance in the region. During MAM and JJA, the South Atlantic SST anomalies affect the magnitude and the meridional displacement of the ITCZ. In JJA, the ENSO has relatively little influence on the interannual variability of the simulated rainfall. During SON, however, the ENSO seems to counteract the effect of the subtropical South Atlantic SST variations on convection over South America.

Agronomic evaluation of calcined crandallite from three Brazilian phosphate deposits

Samples from the weathering mantle containing crandallite of three Brazilian phosphate deposits, Tapira, Catalo, and Juqui, were characterized, calcined, and agronomically evaluated. The calcination process increased total phosphorus (P) and neutral ammonium citrate soluble (NAC) P contents of all samples. The NAC solubility of original Tapira, Catalo, and Juqui was about 5% of total P, whereas for calcined samples it was 54, 16, and 53%, respectively. In a greenhouse study, rates of P were applied at 0, 10, 20, 40, 80, and 120mg P kg(-1) from the calcined materials and MCP (monocalcium phosphate) to an Ultisol cropped with upland and flooded rice for 65 days. The results showed that the calcined P samples increased dry-matter yield and P uptake with increasing rates of P applied for both crops. Tapira and Juqui were more effective for flooded than for upland rice. The calculated values of relative agronomic effectiveness of Tapira, Catalo, and Juqui with respect to MCP were 57, 48, and 53% in dry-matter yield for upland rice and 64, 50, and 69% for flooded rice, respectively.