Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

Polarization effects in molecular dynamics simulations of glass-formers Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8

Thermodynamics, equilibrium structure, and dynamics of glass-forming liquids Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8, have been investigated by molecular dynamics (MD) simulations. A polarizable model was considered for H(2)O and NO(3)- on the basis of previous fluctuating charge models for pure water and the molten salt 2Ca(NO(3))(2)center dot 3KNO(3). Similar thermodynamic properties have been obtained with nonpolarizable and polarizable models. The glass transition temperature, T(g), estimated from MD simulations was dependent on polarization, in particular the dependence of T(g) with electrolyte concentration. Significant polarization effects on equilibrium structure were observed in cation-cation, cation-anion, and water-water structures. Polarization increases the diffusion coefficient of H(2)O, but does not change significantly the diffusion coefficients of ions. Viscosity decreases upon inclusion of polarization, but the conductivity calculated with the polarizable model is smaller than the nonpolarizable model because polarization enhances anion-cation interactions.

Fragility and glassy dynamics of 2Ca(NO(3))(2)center dot 3KNO(3) under pressure: Molecular dynamics simulations

Molecular dynamics simulations of the glass-forming liquid 2Ca(NO(3))(2)center dot 3KNO(3) (CKN) were performed from high temperature liquid states down to low temperature glassy states at six different pressures from 10(-4) to 5.0 GPa. The temperature dependence of the structural relaxation time indicates that the fragility of liquid CKN changes with pressure. In line with recent proposal [Scopigno , Science 302, 849 (2003)], the change on liquid fragility is followed by a proportional change of the nonergodicity factor of the corresponding glass at low temperature. (c) 2008 American Institute of Physics.

Nitrogen dynamics during ecosystem development in tropical forest restoration

We considered whether ecological restoration using high diversity of native tree species serves to restore nitrogen dynamics in the Brazilian Atlantic Forest. We measured delta(15)N and N content in green foliage and soil; vegetation N:P ratio; and soil N mineralization in a preserved natural forest and restored forests of ages 21 and 52 years. Green foliage delta(15)N values, N content, N:P ratio, inorganic N and net mineralization and nitrification rates were all higher, the older the forest. Our findings indicate that the recuperation of N cycling has not been achieved yet in the restored forests even after 52 years, but show that they are following a trajectory of development that is characterized by their N cycling intensity becoming similar to a natural mature forest of the same original forest formation. This study demonstrated that some young restored forests are more limited by N compared to mature natural forests. We document that the recuperation of N cycling in tropical forests can be achieved through ecological restoration actions. (C) 2011 Elsevier B.V. All rights reserved.

Conversion to soy on the Amazonian agricultural frontier increases streamflow without affecting stormflow dynamics

Large-scale soy agriculture in the southern Brazilian Amazon now rivals deforestation for pasture as the region`s predominant form of land use change. Such landscape-level change can have substantial consequences for local and regional hydrology, but these effects remain relatively unstudied in this ecologically and economically important region. We examined how the conversion to soy agriculture influences water balances and stormflows using stream discharge (water yields) and the timing of discharge (stream hydrographs) in small (2.5-13.5 km2) forested and soy headwater watersheds in the Upper Xingu Watershed in the state of Mato Grosso, Brazil. We monitored water yield for 1 year in three forested and four soy watersheds. Mean daily water yields were approximately four times higher in soy than forested watersheds, and soy watersheds showed greater seasonal variability in discharge. The contribution of stormflows to annual streamflow in all streams was low (< 13% of annual streamflow), and the contribution of stormflow to streamflow did not differ between land uses. If the increases in water yield observed in this study are typical, landscape-scale conversion to soy substantially alters water-balance, potentially altering the regional hydrology over large areas of the southern Amazon.

Historical carbon emissions and uptake from the agricultural frontier of the Brazilian Amazon

Tropical ecosystems play a large and complex role in the global carbon cycle. Clearing of natural ecosystems for agriculture leads to large pulses of CO(2) to the atmosphere from terrestrial biomass. Concurrently, the remaining intact ecosystems, especially tropical forests, may be sequestering a large amount of carbon from the atmosphere in response to global environmental changes including climate changes and an increase in atmospheric CO(2). Here we use an approach that integrates census-based historical land use reconstructions, remote-sensing-based contemporary land use change analyses, and simulation modeling of terrestrial biogeochemistry to estimate the net carbon balance over the period 1901-2006 for the state of Mato Grosso, Brazil, which is one of the most rapidly changing agricultural frontiers in the world. By the end of this period, we estimate that of the state`s 925 225 km(2), 221 092 km(2) have been converted to pastures and 89 533 km(2) have been converted to croplands, with forest-to-pasture conversions being the dominant land use trajectory but with recent transitions to croplands increasing rapidly in the last decade. These conversions have led to a cumulative release of 4.8 Pg C to the atmosphere, with similar to 80% from forest clearing and 20% from the clearing of cerrado. Over the same period, we estimate that the residual undisturbed ecosystems accumulated 0.3 Pg C in response to CO2 fertilization. Therefore, the net emissions of carbon from Mato Grosso over this period were 4.5 Pg C. Net carbon emissions from Mato Grosso since 2000 averaged 146 Tg C/yr, on the order of Brazil`s fossil fuel emissions during this period. These emissions were associated with the expansion of croplands to grow soybeans. While alternative management regimes in croplands, including tillage, fertilization, and cropping patterns promote carbon storage in ecosystems, they remain a small portion of the net carbon balance for the region. This detailed accounting of a region`s carbon balance is the type of foundation analysis needed by the new United Nations Collaborative Programmme for Reducing Emissions from Deforestation and Forest Degradation (REDD).

Greenhouse gas emissions from alternative futures of deforestation and agricultural management in the southern Amazon

The Brazilian Amazon is one of the most rapidly developing agricultural areas in the world and represents a potentially large future source of greenhouse gases from land clearing and subsequent agricultural management. In an integrated approach, we estimate the greenhouse gas dynamics of natural ecosystems and agricultural ecosystems after clearing in the context of a future climate. We examine scenarios of deforestation and postclearing land use to estimate the future (2006-2050) impacts on carbon dioxide (CO(2)), methane (CH(4)), and nitrous oxide (N(2)O) emissions from the agricultural frontier state of Mato Grosso, using a process-based biogeochemistry model, the Terrestrial Ecosystems Model (TEM). We estimate a net emission of greenhouse gases from Mato Grosso, ranging from 2.8 to 15.9 Pg CO(2)-equivalents (CO(2)-e) from 2006 to 2050. Deforestation is the largest source of greenhouse gas emissions over this period, but land uses following clearing account for a substantial portion (24-49%) of the net greenhouse gas budget. Due to land-cover and land-use change, there is a small foregone carbon sequestration of 0.2-0.4 Pg CO(2)-e by natural forests and cerrado between 2006 and 2050. Both deforestation and future land-use management play important roles in the net greenhouse gas emissions of this frontier, suggesting that both should be considered in emissions policies. We find that avoided deforestation remains the best strategy for minimizing future greenhouse gas emissions from Mato Grosso.

Stable isotopes of bulk organic matter to trace carbon and nitrogen dynamics in an estuarine ecosystem in Babitonga Bay (Santa Catarina, Brazil)

The biogeochemical processes affecting the transport and cycling of terrestrial organic carbon in coastal and transition areas are still not fully understood One means of distinguishing between the sources of organic materials contributing to particulate organic matter (POM) in Babitonga Bay waters and sediments is by the direct measurement of delta(13)C of dissolved inorganic carbon (DIC) and delta(13)C and delta(15)N in the organic constituents. An isotopic survey was taken from samples collected in the Bay in late spring of 2004. The results indicate that the delta(13)C and delta(15)N compositions of OM varied from -21.7 parts per thousand to -26 2 parts per thousand. and from + 9 2 parts per thousand. to -0 1 parts per thousand, respectively. delta(13)C from DIC ranges from +0.04 parts per thousand to -12.7 parts per thousand The difference in the isotope compositions enables the determination of three distinct end-members terrestrial, marine and urban Moreover, the evaluation of source contribution to the particulate organic matter (POM) in the Bay, enables assessment of the anthropogenic impact. Comparing the depleted values of delta(13)C(DIC) and delta(13)C(POC) it is possible to further understand the carbon dynamic within Babitonga Bay (C) 2010 Elsevier BV All rights reserved

Simulation of sugarcane residue decomposition and aboveground growth

Due to the worldwide increase in demand for biofuels, the area cultivated with sugarcane is expected to increase. For environmental and economic reasons, an increasing proportion of the areas are being harvested without burning, leaving the residues on the soil surface. This periodical input of residues affects soil physical, chemical and biological properties, as well as plant growth and nutrition. Modeling can be a useful tool in the study of the complex interactions between the climate, residue quality, and the biological factors controlling plant growth and residue decomposition. The approach taken in this work was to parameterize the CENTURY model for the sugarcane crop, to simulate the temporal dynamics of aboveground phytomass and litter decomposition, and to validate the model through field experiment data. When studying aboveground growth, burned and unburned harvest systems were compared, as well as the effect of mineral fertilizer and organic residue applications. The simulations were performed with data from experiments with different durations, from 12 months to 60 years, in Goiana, TimbaA(0)ba and Pradpolis, Brazil; Harwood, Mackay and Tully, Australia; and Mount Edgecombe, South Africa. The differentiation of two pools in the litter, with different decomposition rates, was found to be a relevant factor in the simulations made. Originally, the model had a basically unlimited layer of mulch directly available for decomposition, 5,000 g m(-2). Through a parameter optimization process, the thickness of the mulch layer closer to the soil, more vulnerable to decomposition, was set as 110 g m(-2). By changing the layer of mulch at any given time available for decomposition, the sugarcane residues decomposition simulations where close to measured values (R (2) = 0.93), contributing to making the CENTURY model a tool for the study of sugarcane litter decomposition patterns. The CENTURY model accurately simulated aboveground carbon stalk values (R (2) = 0.76), considering burned and unburned harvest systems, plots with and without nitrogen fertilizer and organic amendment applications, in different climates and soil conditions.

Simulation of Soil Carbon Dynamics under Sugarcane with the CENTURY Model

Currently there is a trend for the expansion of the area cropped with sugarcane (Saccharum officinarum L.), driven by an increase in the world demand for biofuels, due to economical, environmental, and geopolitical issues. Although sugarcane is traditionally harvested by burning dried leaves and tops, the unburned, mechanized harvest has been progressively adopted. The use of process based models is useful in understanding the effects of plant litter in soil C dynamics. The objective of this work was to use the CENTURY model in evaluating the effect of sugarcane residue management in the temporal dynamics of soil C. The approach taken in this work was to parameterize the CENTURY model for the sugarcane crop, to simulate the temporal dynamics of soil C, validating the model through field experiment data, and finally to make predictions in the long term regarding soil C. The main focus of this work was the comparison of soil C stocks between the burned and unburned litter management systems, but the effect of mineral fertilizer and organic residue applications were also evaluated. The simulations were performed with data from experiments with different durations, from 1 to 60 yr, in Goiana and Timbauba, Pernambuco, and Pradopolis, Sao Paulo, all in Brazil; and Mount Edgecombe, Kwazulu-Natal, South Africa. It was possible to simulate the temporal dynamics of soil C (R(2) = 0.89). The predictions made with the model revealed that there is, in the long term, a trend for higher soil C stocks with the unburned management. This increase is conditioned by factors such as climate, soil texture, time of adoption of the unburned system, and N fertilizer management.

DISTRIBUTION OF CAPYBARAS IN AN AGROECOSYSTEM, SOUTHEASTERN BRAZIL, BASED ON ECOLOGICAL NICHE MODELING

Southeastern Brazil has seen dramatic landscape modifications in recent decades, due to expansion of agriculture and urban areas; these changes have influenced the distribution and abundance of vertebrates. We developed predictive models of ecological and spatial distributions of capybaras (Hydrochoerus hydrochaeris) using ecological niche modeling. Most Occurrences of capybaras were in flat areas with water bodies Surrounded by sugarcane and pasture. More than 75% of the Piracicaba River basin was estimated as potentially habitable by capybara. The models had low omission error (2.3-3.4%), but higher commission error (91.0-98.5%); these ""model failures"" seem to be more related to local habitat characteristics than to spatial ones. The potential distribution of capybaras in the basin is associated with anthropogenic habitats, particularly with intensive land use for agriculture.

Dynamics of carbon, biomass, and structure in two Amazonian forests

Amazon forests are potentially globally significant sources or sinks for atmospheric carbon dioxide. In this study, we characterize the spatial trends in carbon storage and fluxes in both live and dead biomass (necromass) in two Amazonian forests, the Biological Dynamic of Forest Fragments Project (BDFFP), near Manaus, Amazonas, and the Tapajos National Forest (TNF) near Santarem, Para. We assessed coarse woody debris (CWD) stocks, tree growth, mortality, and recruitment in ground-based plots distributed across the terra firme forest at both sites. Carbon dynamics were similar within each site, but differed significantly between the sites. The BDFFP and the TNF held comparable live biomass (167 +/- 7.6 MgC.ha(-1) versus 149 +/- 6.0 MgC.ha(-1), respectively), but stocks of CWD were 2.5 times larger at TNF (16.2 +/- 1.5 MgC.ha(-1) at BDFFP, versus 40.1 +/- 3.9 MgC.ha(-1) at TNF). A model of current forest dynamics suggests that the BDFFP was close to carbon balance, and its size class structure approximated a steady state. The TNF, by contrast, showed rapid carbon accrual to live biomass (3.24 +/- 0.22 MgC.ha(-1).a(-1) in TNF, 2.59 +/- 0.16 MgC.ha(-1).a(-1) in BDFFP), which was more than offset by losses from large stocks of CWD, as well as ongoing shifts of biomass among size classes. This pattern in the TNF suggests recovery from a significant disturbance. The net loss of carbon from the TNF will likely last 10 - 15 years after the initial disturbance (controlled by the rate of decay of coarse woody debris), followed by uptake of carbon as the forest size class structure and composition continue to shift. The frequency and longevity of forests showing such disequilibruim dynamics within the larger matrix of the Amazon remains an essential question to understanding Amazonian carbon balance.

Multifractal regime transition in a modified minority game model

The search for more realistic modeling of financial time series reveals several stylized facts of real markets. In this work we focus on the multifractal properties found in price and index signals. Although the usual minority game (MG) models do not exhibit multifractality, we study here one of its variants that does. We show that the nonsynchronous MG models in the nonergodic phase is multifractal and in this sense, together with other stylized facts, constitute a better modeling tool. Using the structure function (SF) approach we detected the stationary and the scaling range of the time series generated by the MG model and, from the linear (non-linear) behavior of the SF we identified the fractal (multifractal) regimes. Finally, using the wavelet transform modulus maxima (WTMM) technique we obtained its multifractal spectrum width for different dynamical regimes. (C) 2009 Elsevier Ltd. All rights reserved.

An opinion dynamics model for the diffusion of innovations

We study the dynamics of the adoption of new products by agents with continuous opinions and discrete actions (CODA). The model is such that the refusal in adopting a new idea or product is increasingly weighted by neighbor agents as evidence against the product. Under these rules, we study the distribution of adoption times and the final proportion of adopters in the population. We compare the cases where initial adopters are clustered to the case where they are randomly scattered around the social network and investigate small world effects on the final proportion of adopters. The model predicts a fat tailed distribution for late adopters which is verified by empirical data. (C) 2009 Elsevier B.V. All rights reserved.

Continuous opinions and discrete actions in opinion dynamics problems

A model where agents show discrete behavior regarding their actions, but have continuous opinions that are updated by interacting with other agents is presented. This new updating rule is applied to both the voter and Sznajd models for interaction between neighbors, and its consequences are discussed. The appearance of extremists is naturally observed and it seems to be a characteristic of this model.