A case study of the direct radiative effect of biomass burning aerosols on precipitation in the Eastern Amazon

Numerical experiments with the Brazilian additions to the Regional Atmospheric Modeling System were performed with two nested grids (50 and 10 km horizontal resolution, respectively) with and without the effect of biomass burning for 8 different situations for 96 h integrations. Only the direct radiative effect of aerosols is considered. The results were analyzed in large areas encompassing the BR163 road (one of the main areas of deforestation in the Amazon). mainly where most of the burning takes place. The precipitation change due to the direct radiative impact of biomass burning is generally negative (i.e., there is a decrease of precipitation). However, there are a few cases with a positive impact. Two opposite forcing mechanisms were explored: (a) the thermodynamic forcing that is generally negative in the sense that the aerosol tends to stabilize the lower atmosphere and (b) the dynamic impact associated with the low level horizontal pressure gradients produced by the aerosol plumes. In order to understand the non-linear relationship between the two effects, experiments were performed with 4-fold emissions. In these cases, the dynamic effect overcomes the stabilization produced by the radiative forcing and precipitation increase is observed in comparison with the control experiment. This study suggests that. in general, the biomass burning radiative forcing decreases the precipitation. However, very large concentrations of aerosols may lead to an increase of precipitation due to the dynamical forcing associated with the horizontal pressure gradients. (C) 2009 Elsevier B.V. All rights reserved.

Precipitation diurnal cycle and summer climatology assessment over South America: An evaluation of Regional Climate Model version 3 simulations

Regional Climate Model version 3 (RegCM3) simulations of 17 summers (1988-2004) over part of South America south of 5 degrees S were evaluated to identify model systematic errors. Model results were compared to different rainfall data sets (Climate Research Unit (CRU), Climate Prediction Center (CPC), Global Precipitation Climatology Project (GPCP), and National Centers for Environmental Prediction (NCEP) reanalysis), including the five summers mean (1998-2002) precipitation diurnal cycle observed by the Tropical Rainfall Measuring Mission (TRMM)-Precipitation Radar (PR). In spite of regional differences, the RegCM3 simulates the main observed aspects of summer climatology associated with the precipitation (northwest-southeast band of South Atlantic Convergence Zone (SACZ)) and air temperature (warmer air in the central part of the continent and colder in eastern Brazil and the Andes Mountains). At a regional scale, the main RegCM3 failures are the underestimation of the precipitation in the northern branch of the SACZ and some unrealistic intense precipitation around the Andes Mountains. However, the RegCM3 seasonal precipitation is closer to the fine-scale analyses (CPC, CRU, and TRMM-PR) than is the NCEP reanalysis, which presents an incorrect north-south orientation of SACZ and an overestimation of its intensity. The precipitation diurnal cycle observed by TRMM-PR shows pronounced contrasts between Tropics and Extratropics and land and ocean, where most of these features are simulated by RegCM3. The major similarities between the simulation and observation, especially the diurnal cycle phase, are found over the continental tropical and subtropical SACZ regions, which present afternoon maximum (1500-1800 UTC) and morning minimum (0900-1200 UTC). More specifically, over the core of SACZ, the phase and amplitude of the simulated precipitation diurnal cycle are very close to the TRMM-PR observations. Although there are amplitude differences, the RegCM3 simulates the observed nighttime rainfall in the eastern Andes Mountains, over the Atlantic Ocean, and also over northern Argentina. The main simulation deficiencies are found in the Atlantic Ocean and near the Andes Mountains. Over the Atlantic Ocean the convective scheme is not triggered; thus the rainfall arises from the grid-scale scheme and therefore differs from the TRMM-PR. Near the Andes, intense (nighttime and daytime) simulated precipitation could be a response of an incorrect circulation and topographic uplift. Finally, it is important to note that unlike most reported bias of global models, RegCM3 does not trigger the moist convection just after sunrise over the southern part of the Amazon.

Impact of biomass burning aerosols on precipitation in the Amazon: A modeling case study

A study of the potential role of aerosols in modifying clouds and precipitation is presented using a numerical atmospheric model. Measurements of cloud condensation nuclei (CCN) and cloud size distribution properties taken in the southwestern Amazon region during the transition from dry to wet seasons were used as guidelines to define the microphysical parameters for the simulations. Numerical simulations were carried out using the Brazilian Development on Regional Atmospheric Modeling System, and the results presented considerable sensitivity to changes in these parameters. High CCN concentrations, typical of polluted days, were found to result in increases or decreases in total precipitation, depending on the level of pollution used as a reference, showing a complexity that parallels the aerosol-precipitation interaction. Our results show that on the grids evaluated, higher CCN concentrations reduced low-to-moderate rainfall rates and increased high rainfall rates. The principal consequence of the increased pollution was a change from a warm to a cold rain process, which affected the maximum and overall mean accumulated precipitation. Under polluted conditions, cloud cover diminished, allowing greater amounts of solar radiation to reach the surface. Aerosol absorption of radiation in the lower layers of the atmosphere delayed convective evolution but produced higher maximum rainfall rates due to increased instability. In addition, the intensity of the surface sensible heat flux, as well as that of the latent heat flux, was reduced by the lower temperature difference between surface and air, producing greater energy stores at the surface.

Precipitation variability in Sao Paulo State, Brazil

The State of Sao Paulo is the richest in Brazil, responsible for over 30% of the Brazilian gross rate. It has a population of around 30 million and its economy is based on agriculture and industrial products. Any change in climate can have a profound influence on the socio-economics of the State. In order to determine changes in total and extreme rainfall over Sao Paulo State, climate change indices derived from daily precipitation data were calculated using specially designed software. Maps of trends for a subset of 59 rain gauge stations were analysed for the period 1950-1999 and also for a subset of this period, 1990-1999, representing more recent climate. A non-parametric Mann-Kendall test was applied to the time series. Maps of trends for six annual precipitation indices (annual total precipitation (PRCPTOT), very heavy precipitation days (R20mm), events greater than the 95th percentile (R95p), maximum five days precipitation total (RX5day), the length of the largest wet spell (CWD) and the length of the largest dry spell (CDD)) were analysed for the entire period. These exhibited statistically significant trends associated with a wetter climate. A significant increase in PRCPTOT, associated with very heavy precipitation days, were observed at more than 45% of the rain gauge stations. The Mann-Kendall test identified that the positive trend in PRCPTOT is possibly related to the increase in the R95p and R20mm indices. Therefore, the results suggest that there has been a change in precipitation intensity. In contrast, the indices for the more recent shorter time series are significantly different to the longer term indices. The results indicate that intense precipitation is becoming concentrated in a few days and spread over the period when the CDD and R20mm indices show positive trends, while negative ones are seen in the RX5day index. The trends found could be related to many anthropogenic aspects such as biomass burning aerosols and land use.

Decomposition in tropical forests: a pan-tropical study of the effects of litter type, litter placement and mesofaunal exclusion across a precipitation gradient

1. Litter decomposition recycles nutrients and causes large fluxes of carbon dioxide into the atmosphere. It is typically assumed that climate, litter quality and decomposer communities determine litter decay rates, yet few comparative studies have examined their relative contributions in tropical forests. 2. We used a short-term litterbag experiment to quantify the effects of litter quality, placement and mesofaunal exclusion on decomposition in 23 tropical forests in 14 countries. Annual precipitation varied among sites (760-5797 mm). At each site, two standard substrates (Raphia farinifera and Laurus nobilis) were decomposed in fine- and coarse-mesh litterbags both above and below ground for approximately 1 year. 3. Decomposition was rapid, with >95% mass loss within a year at most sites. Litter quality, placement and mesofaunal exclusion all independently affected decomposition, but the magnitude depended upon site. Both the average decomposition rate at each site and the ratio of above- to below-ground decay increased linearly with annual precipitation, explaining 60-65% of among-site variation. Excluding mesofauna had the largest impact on decomposition, reducing decomposition rates by half on average, but the magnitude of decrease was largely independent of climate. This suggests that the decomposer community might play an important role in explaining patterns of decomposition among sites. Which litter type decomposed fastest varied by site, but was not related to climate. 4. Synthesis. A key goal of ecology is to identify general patterns across ecological communities, as well as relevant site-specific details to understand local dynamics. Our pan-tropical study shows that certain aspects of decomposition, including average decomposition rates and the ratio of above- to below-ground decomposition are highly correlated with a simple climatic index: mean annual precipitation. However, we found no relationship between precipitation and effects of mesofaunal exclusion or litter type, suggesting that site-specific details may also be required to understand how these factors affect decomposition at local scales.

Discontinuous local semiflows for Kurzweil equations leading to LaSalle`s invariance principle for differential systems with impulses at variable times

In this paper, we consider an initial value problem for a class of generalized ODEs, also known as Kurzweil equations, and we prove the existence of a local semidynamical system there. Under certain perturbation conditions, we also show that this class of generalized ODEs admits a discontinuous semiflow which we shall refer to as an impulsive semidynamical system. As a consequence, we obtain LaSalle`s invariance principle for such a class of generalized ODEs. Due to the importance of LaSalle`s invariance principle in studying stability of differential systems, we include an application to autonomous ordinary differential systems with impulse action at variable times. (C) 2011 Elsevier Inc. All rights reserved.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Rainforest Aerosols as Biogenic Nuclei of Clouds and Precipitation in the Amazon

The Amazon is one of the few continental regions where atmospheric aerosol particles and their effects on climate are not dominated by anthropogenic sources. During the wet season, the ambient conditions approach those of the pristine pre-industrial era. We show that the fine submicrometer particles accounting for most cloud condensation nuclei are predominantly composed of secondary organic material formed by oxidation of gaseous biogenic precursors. Supermicrometer particles, which are relevant as ice nuclei, consist mostly of primary biological material directly released from rainforest biota. The Amazon Basin appears to be a biogeochemical reactor, in which the biosphere and atmospheric photochemistry produce nuclei for clouds and precipitation sustaining the hydrological cycle. The prevailing regime of aerosol-cloud interactions in this natural environment is distinctly different from polluted regions.

Critical discontinuous phase transition in the threshold contact process

We analyze a threshold contact process on a square lattice in which particles are created on empty sites with at least two neighboring particles and are annihilated spontaneously. We show by means of Monte Carlo simulations that the process undergoes a discontinuous phase transition at a definite value of the annihilation parameter, in accordance with the Gibbs phase rule, and that the discontinuous transition exhibits critical behavior. The simulations were performed by using boundary conditions in which the sites of the border of the lattice are permanently occupied by particles.

Orbitally driven east-west antiphasing of South American precipitation

The variations of tropical precipitation are antiphased between the hemispheres on orbital timescales. This antiphasing arises through the alternating strength of incoming solar radiation in the two hemispheres, which affects monsoon intensity and hence the position of the meridional atmospheric circulation of the Hadley cells(1-4). Here we compare an oxygen isotopic record recovered from a speleothem from northeast Brazil for the past 26,000 years with existing reconstructions of precipitation in tropical South America(5-8). During the Holocene, we identify a similar, but zonally oriented, antiphasing of precipitation within the same hemisphere: northeast Brazil experiences humid conditions during low summer insolation and aridity when summer insolation is high, whereas the rest of southern tropical South America shows opposite characteristics. Simulations with a general circulation model that incorporates isotopic variations support this pattern as well as the link to insolation-driven monsoon activity. Our results suggest that convective heating over tropical South America and associated adjustments in large-scale subsidence over northeast Brazil lead to a remote forcing of the South American monsoon, which determines most of the precipitation changes in the region on orbital timescales.

Computational aspects of harmonic wavelet Galerkin methods and an application to a precipitation front propagation model

This article is dedicated to harmonic wavelet Galerkin methods for the solution of partial differential equations. Several variants of the method are proposed and analyzed, using the Burgers equation as a test model. The computational complexity can be reduced when the localization properties of the wavelets and restricted interactions between different scales are exploited. The resulting variants of the method have computational complexities ranging from O(N(3)) to O(N) (N being the space dimension) per time step. A pseudo-spectral wavelet scheme is also described and compared to the methods based on connection coefficients. The harmonic wavelet Galerkin scheme is applied to a nonlinear model for the propagation of precipitation fronts, with the front locations being exposed in the sizes of the localized wavelet coefficients. (C) 2011 Elsevier Ltd. All rights reserved.

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

We design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters. (C) 2010 Elsevier B.V. All rights reserved.

Accurate velocity reconstruction for Discontinuous Galerkin approximations of two-phase porous media flows

We consider Discontinuous Galerkin approximations of two-phase, immiscible porous media flows in the global pressure/fractional flow formulation with capillary pressure. A sequential approach is used with a backward Euler step for the saturation equation, equal-order interpolation for the pressure and the saturation, and without any limiters. An accurate total velocity field is recovered from the global pressure equation to be used in the saturation equation. Numerical experiments show the advantages of the proposed reconstruction. To cite this article: A. Ern et al., C R. Acad. Sci. Paris, Ser. 1347 (2009). (C) 2009 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.