Beyond the Fragmentation Threshold Hypothesis: Regime Shifts in Biodiversity Across Fragmented Landscapes

Ecological systems are vulnerable to irreversible change when key system properties are pushed over thresholds, resulting in the loss of resilience and the precipitation of a regime shift. Perhaps the most important of such properties in human-modified landscapes is the total amount of remnant native vegetation. In a seminal study Andren proposed the existence of a fragmentation threshold in the total amount of remnant vegetation, below which landscape-scale connectivity is eroded and local species richness and abundance become dependent on patch size. Despite the fact that species patch-area effects have been a mainstay of conservation science there has yet to be a robust empirical evaluation of this hypothesis. Here we present and test a new conceptual model describing the mechanisms and consequences of biodiversity change in fragmented landscapes, identifying the fragmentation threshold as a first step in a positive feedback mechanism that has the capacity to impair ecological resilience, and drive a regime shift in biodiversity. The model considers that local extinction risk is defined by patch size, and immigration rates by landscape vegetation cover, and that the recovery from local species losses depends upon the landscape species pool. Using a unique dataset on the distribution of non-volant small mammals across replicate landscapes in the Atlantic forest of Brazil, we found strong evidence for our model predictions - that patch-area effects are evident only at intermediate levels of total forest cover, where landscape diversity is still high and opportunities for enhancing biodiversity through local management are greatest. Furthermore, high levels of forest loss can push native biota through an extinction filter, and result in the abrupt, landscape-wide loss of forest-specialist taxa, ecological resilience and management effectiveness. The proposed model links hitherto distinct theoretical approaches within a single framework, providing a powerful tool for analysing the potential effectiveness of management interventions.

Transport of North African dust from the Bod,l, depression to the Amazon Basin: a case study

Through long-range transport of dust, the North-African desert supplies essential minerals to the Amazon rain forest. Since North African dust reaches South America mostly during the Northern Hemisphere winter, the dust sources active during winter are the main contributors to the forest. Given that the Bod,l, depression area in southwestern Chad is the main winter dust source, a close link is expected between the Bod,l, emission patterns and volumes and the mineral supply flux to the Amazon. Until now, the particular link between the Bod,l, and the Amazon forest was based on sparse satellite measurements and modeling studies. In this study, we combine a detailed analysis of space-borne and ground data with reanalysis model data and surface measurements taken in the central Amazon during the Amazonian Aerosol Characterization Experiment (AMAZE-08) in order to explore the validity and the nature of the proposed link between the Bod,l, depression and the Amazon forest. This case study follows the dust events of 11-16 and 18-27 February 2008, from the emission in the Bod,l, over West Africa (most likely with contribution from other dust sources in the region) the crossing of the Atlantic Ocean, to the observed effects above the Amazon canopy about 10 days after the emission. The dust was lifted by surface winds stronger than 14 m s(-1), usually starting early in the morning. The lofted dust, mixed with biomass burning aerosols over Nigeria, was transported over the Atlantic Ocean, and arrived over the South American continent. The top of the aerosol layer reached above 3 km, and the bottom merged with the boundary layer. The arrival of the dusty air parcel over the Amazon forest increased the average concentration of aerosol crustal elements by an order of magnitude.

Cloud condensation nuclei in pristine tropical rainforest air of Amazonia: size-resolved measurements and modeling of atmospheric aerosol composition and CCN activity

Atmospheric aerosol particles serving as cloud condensation nuclei (CCN) are key elements of the hydrological cycle and climate. We have measured and characterized CCN at water vapor supersaturations in the range of S=0.10-0.82% in pristine tropical rainforest air during the AMAZE-08 campaign in central Amazonia. The effective hygroscopicity parameters describing the influence of chemical composition on the CCN activity of aerosol particles varied in the range of kappa approximate to 0.1-0.4 (0.16+/-0.06 arithmetic mean and standard deviation). The overall median value of kappa approximate to 0.15 was by a factor of two lower than the values typically observed for continental aerosols in other regions of the world. Aitken mode particles were less hygroscopic than accumulation mode particles (kappa approximate to 0.1 at D approximate to 50 nm; kappa approximate to 0.2 at D approximate to 200 nm), which is in agreement with earlier hygroscopicity tandem differential mobility analyzer (H-TDMA) studies. The CCN measurement results are consistent with aerosol mass spectrometry (AMS) data, showing that the organic mass fraction (f(org)) was on average as high as similar to 90% in the Aitken mode (D <= 100 nm) and decreased with increasing particle diameter in the accumulation mode (similar to 80% at D approximate to 200 nm). The kappa values exhibited a negative linear correlation with f(org) (R(2)=0.81), and extrapolation yielded the following effective hygroscopicity parameters for organic and inorganic particle components: kappa(org)approximate to 0.1 which can be regarded as the effective hygroscopicity of biogenic secondary organic aerosol (SOA) and kappa(inorg)approximate to 0.6 which is characteristic for ammonium sulfate and related salts. Both the size dependence and the temporal variability of effective particle hygroscopicity could be parameterized as a function of AMS-based organic and inorganic mass fractions (kappa(p)=kappa(org) x f(org)+kappa(inorg) x f(inorg)). The CCN number concentrations predicted with kappa(p) were in fair agreement with the measurement results (similar to 20% average deviation). The median CCN number concentrations at S=0.1-0.82% ranged from N(CCN,0.10)approximate to 35 cm(-3) to N(CCN,0.82)approximate to 160 cm(-3), the median concentration of aerosol particles larger than 30 nm was N(CN,30)approximate to 200 cm(-3), and the corresponding integral CCN efficiencies were in the range of N(CCN,0.10/NCN,30)approximate to 0.1 to N(CCN,0.82/NCN,30)approximate to 0.8. Although the number concentrations and hygroscopicity parameters were much lower in pristine rainforest air, the integral CCN efficiencies observed were similar to those in highly polluted megacity air. Moreover, model calculations of N(CCN,S) assuming an approximate global average value of kappa approximate to 0.3 for continental aerosols led to systematic overpredictions, but the average deviations exceeded similar to 50% only at low water vapor supersaturation (0.1%) and low particle number concentrations (<= 100 cm(-3)). Model calculations assuming aconstant aerosol size distribution led to higher average deviations at all investigated levels of supersaturation: similar to 60% for the campaign average distribution and similar to 1600% for a generic remote continental size distribution. These findings confirm earlier studies suggesting that aerosol particle number and size are the major predictors for the variability of the CCN concentration in continental boundary layer air, followed by particle composition and hygroscopicity as relatively minor modulators. Depending on the required and applicable level of detail, the information and parameterizations presented in this paper should enable efficient description of the CCN properties of pristine tropical rainforest aerosols of Amazonia in detailed process models as well as in large-scale atmospheric and climate models.

Melting and freezing of spherical bismuth nanoparticles confined in a homogeneous sodium borate glass

The melting temperature and the crystallization temperature of Bi nanoclusters confined in a sodium borate glass were experimentally determined as functions of the cluster radius. The results indicate that, on cooling, liquid Bi nanodroplets exhibit a strong undercooling effect for a wide range of radii. The difference between the melting temperature and the freezing temperature decreases for decreasing radius and vanishes for Bi nanoparticles with a critical radius R = 1.9 nm. The magnitude of the variation in density across the melting and freezing transitions for Bi nanoparticles with R = 2 nm is 40% smaller than for bulk Bi. These experimental results support a basic core-shell model for the structure of Bi nanocrystals consisting of a central crystalline volume surrounded by a structurally disordered shell. The volume fraction of the crystalline core decreases for decreasing nanoparticle radius and vanishes for R = 1.9 nm. Thus, on cooling, the liquid nanodroplets with R < 1.9 nm preserve, across the liquid-to-solid transformation, their homogeneous and disordered structure without crystalline core.

Remote sensing the vertical profile of cloud droplet effective radius, thermodynamic phase, and temperature

Cloud-aerosol interaction is a key issue in the climate system, affecting the water cycle, the weather, and the total energy balance including the spatial and temporal distribution of latent heat release. Information on the vertical distribution of cloud droplet microphysics and thermodynamic phase as a function of temperature or height, can be correlated with details of the aerosol field to provide insight on how these particles are affecting cloud properties and their consequences to cloud lifetime, precipitation, water cycle, and general energy balance. Unfortunately, today's experimental methods still lack the observational tools that can characterize the true evolution of the cloud microphysical, spatial and temporal structure in the cloud droplet scale, and then link these characteristics to environmental factors and properties of the cloud condensation nuclei. Here we propose and demonstrate a new experimental approach (the cloud scanner instrument) that provides the microphysical information missed in current experiments and remote sensing options. Cloud scanner measurements can be performed from aircraft, ground, or satellite by scanning the side of the clouds from the base to the top, providing us with the unique opportunity of obtaining snapshots of the cloud droplet microphysical and thermodynamic states as a function of height and brightness temperature in clouds at several development stages. The brightness temperature profile of the cloud side can be directly associated with the thermodynamic phase of the droplets to provide information on the glaciation temperature as a function of different ambient conditions, aerosol concentration, and type. An aircraft prototype of the cloud scanner was built and flew in a field campaign in Brazil. The CLAIM-3D (3-Dimensional Cloud Aerosol Interaction Mission) satellite concept proposed here combines several techniques to simultaneously measure the vertical profile of cloud microphysics, thermodynamic phase, brightness temperature, and aerosol amount and type in the neighborhood of the clouds. The wide wavelength range, and the use of multi-angle polarization measurements proposed for this mission allow us to estimate the availability and characteristics of aerosol particles acting as cloud condensation nuclei, and their effects on the cloud microphysical structure. These results can provide unprecedented details on the response of cloud droplet microphysics to natural and anthropogenic aerosols in the size scale where the interaction really happens.

Robust relations between CCN and the vertical evolution of cloud drop size distribution in deep convective clouds

In-situ measurements in convective clouds (up to the freezing level) over the Amazon basin show that smoke from deforestation fires prevents clouds from precipitating until they acquire a vertical development of at least 4 km, compared to only 1-2 km in clean clouds. The average cloud depth required for the onset of warm rain increased by similar to 350 m for each additional 100 cloud condensation nuclei per cm(3) at a super-saturation of 0.5% (CCN0.5%). In polluted clouds, the diameter of modal liquid water content grows much slower with cloud depth (at least by a factor of similar to 2), due to the large number of droplets that compete for available water and to the suppressed coalescence processes. Contrary to what other studies have suggested, we did not observe this effect to reach saturation at 3000 or more accumulation mode particles per cm(3). The CCN0.5% concentration was found to be a very good predictor for the cloud depth required for the onset of warm precipitation and other microphysical factors, leaving only a secondary role for the updraft velocities in determining the cloud drop size distributions. The effective radius of the cloud droplets (r(e)) was found to be a quite robust parameter for a given environment and cloud depth, showing only a small effect of partial droplet evaporation from the cloud's mixing with its drier environment. This supports one of the basic assumptions of satellite analysis of cloud microphysical processes: the ability to look at different cloud top heights in the same region and regard their r(e) as if they had been measured inside one well developed cloud. The dependence of r(e) on the adiabatic fraction decreased higher in the clouds, especially for cleaner conditions, and disappeared at r(e)>=similar to 10 mu m. We propose that droplet coalescence, which is at its peak when warm rain is formed in the cloud at r(e)=similar to 10 mu m, continues to be significant during the cloud's mixing with the entrained air, cancelling out the decrease in r(e) due to evaporation.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: A statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types-each one associated, respectively, with the polar-headgroup and the acyl-chain states-which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

Combined Monte Carlo and quantum mechanics study of the solvatochromism of phenol in water. The origin of the blue shift of the lowest pi-pi* transition

A combined and sequential use of Monte Carlo simulations and quantum mechanical calculations is made to analyze the spectral shift of the lowest pi-pi* transition of phenol in water. The solute polarization is included using electrostatic embedded calculations at the MP2/aug-cc-pVDZ level giving a dipole moment of 2.25 D, corresponding to an increase of 76% compared to the calculated gas-phase value. Using statistically uncorrelated configurations sampled from the MC simulation,first-principle size-extensive calculations are performed to obtain the solvatochromic shift. Analysis is then made of the origin of the blue shift. Results both at the optimized geometry and in room-temperature liquid water show that hydrogen bonds of water with phenol promote a red shift when phenol is the proton-donor and a blue shift when phenol is the proton-acceptor. In the case of the optimized clusters the calculated shifts are in very good agreement with results obtained from mass-selected free jet expansion experiments. In the liquid case the contribution of the solute-solvent hydrogen bonds partially cancels and the total shift obtained is dominated by the contribution of the outer solvent water molecules. Our best result, including both inner and outer water molecules, is 570 +/- 35 cm(-1), in very good agreement with the small experimental shift of 460 cm(-1) for the absorption maximum.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Approximation to density functional theory for the calculation of band gaps of semiconductors

The local-density approximation (LDA) together with the half occupation (transitionstate) is notoriously successful in the calculation of atomic ionization potentials. When it comes to extended systems, such as a semiconductor infinite system, it has been very difficult to find a way to half ionize because the hole tends to be infinitely extended (a Bloch wave). The answer to this problem lies in the LDA formalism itself. One proves that the half occupation is equivalent to introducing the hole self-energy (electrostatic and exchange correlation) into the Schrodinger equation. The argument then becomes simple: The eigenvalue minus the self-energy has to be minimized because the atom has a minimal energy. Then one simply proves that the hole is localized, not infinitely extended, because it must have maximal self-energy. Then one also arrives at an equation similar to the self- interaction correction equation, but corrected for the removal of just 1/2 electron. Applied to the calculation of band gaps and effective masses, we use the self- energy calculated in atoms and attain a precision similar to that of GW, but with the great advantage that it requires no more computational effort than standard LDA.

A preliminary study of the incorparation of NPK fertilizer into chitosan nanoparticles

The use of slow release fertilizer has become a new trend to save fertilizer consumption and to minimize environmental pollution. Due to its polymeric cationic, biodegradable, bioabsorbable, and bactericidal characteristics, chitosan (CS) nanoparticle is an interesting material for use in controlled release systems. However, there are no attempts to explore the potential of chitosan nanoparticles as controlled release for NPK fertilizers. In this work chitosan nanoparticles were obtained by polymerizing methacrylic acid for the incorporation of NPK fertilizers. The interaction and stability of chitosan nanoparticle suspensions containing nitrogen (N), phosphorus (P) and potassium (K) were evaluated by FTIR spectroscopy, particle size analysis and zeta-potential. The FTIR results indicated the existence of electrostatic interactions between chitosan nanoparticles and the elements N, P and K. The stability of the CS-PMAA colloidal suspension was higher with the addition of nitrogen and potassium than with the addition of phosphorus, due to the higher anion charge from the calcium phosphate than the anion charges from the potassium chloride and urea. The mean diameter increase of the CS-PMAA nanoparticles in suspension with the addition of different compounds indicated that the elements are being aggregated on the surface of the chitosan nanoparticles.

Nitrogen availability patterns in white-sand vegetations of Central Brazilian Amazon

Addressing spatial variability in nitrogen (N) availability in the Central Brazilian Amazon, we hypothesized that N availability varies among white-sand vegetation types (campina and campinarana) and lowland tropical forests (dense terra-firme forests) in the Central Brazilian Amazon, under the same climate conditions. Accordingly, we measured soil and foliar N concentration and N isotope ratios (delta(15)N) throughout the campina-campinarana transect and compared to published dense terra-firme forest results. There were no differences between white-sand vegetation types in regard to soil N concentration, C:N ratio and delta(15)N across the transect. Both white-sand vegetation types showed very low foliar N concentrations and elevated foliar C:N ratios, and no significant difference between site types was observed. Foliar delta(15)N was depleted, varying from -9.6 to 1.6aEuro degrees in the white-sand vegetations. The legume Aldina heterophylla had the highest average delta(15)N values (-1.5aEuro degrees) as well as the highest foliar N concentration (2.1%) while the non-legume species had more depleted delta(15)N values and the average foliar N concentrations varied from 0.9 to 1.5% among them. Despite the high variation in foliar delta(15)N among plants, a significant and gradual (15)N-enrichment in foliar isotopic signatures throughout the campina-campinarana transect was observed. Individual plants growing in the campinarana were significantly enriched in (15)N compared to those in campina. In the white-sand N-limited ecosystems, the differentiation of N use seems to be a major cause of variations observed in foliar delta(15)N values throughout the campina-campinarana transect.

Understanding the Influences of Spatial Patterns on N Availability Within the Brazilian Amazon Forest

Nitrogen variations at different spatial scales and integrated across functional groups were addressed for lowland tropical forests in the Brazilian Amazon as follows: (1) how does N availability vary across the region over different spatial scales (regional x landscape scale); ( 2) how are these variations in N availability integrated across plant functional groups ( legume 9 non-legume trees). Leaf N, P, and Ca concentrations as well the leaf N isotope ratios (delta(15)N) from a large set of legume and non-legume tree species were measured. Legumes had higher foliar N/Ca ratios than non-legumes, consistent with the high energetic costs in plant growth associated with higher foliar P/Ca ratios found in legumes than in non-legumes. At the regional level, foliar delta(15)N decreased with increasing rainfall. At the landscape level, N availability was higher in the forests on clayey soils on the plateau than in forests on sandier soils. The isotope as well as the non-isotope data relationships here documented, explain to a large extent the variation in delta(15)N signatures across gradients of rainfall and soil. Although at the regional level, the precipitation regime is a major determinant of differences in N availability, at the landscape level, under the same precipitation regime, soil type seems to be a major factor influencing the availability of N in the Brazilian Amazon forest.

Spatio-temporal patterns of throughfall and solute deposition in an open tropical rain forest

The brief interaction of precipitation with a forest canopy can create a high spatial variability of both throughfall and solute deposition. We hypothesized that (i) the variability in natural forest systems is high but depends on system-inherent stability, (ii) the spatial variability of solute deposition shows seasonal dynamics depending on the increase in rainfall frequency, and (iii) spatial patterns persist only in the short-term. The study area in the north-western Brazilian state of Rondonia is subject to a climate with a distinct wet and dry season. We collected rain and throughfall on an event basis during the early wet season (n = 14) and peak of the wet season (n = 14) and analyzed the samples for pH and concentrations of NH4+, Na+, K+, Ca2+ Mg2+,, Cl-, NO3-, SO42- and DOC. The coefficient 3 4 cient of variation for throughfall based on both sampling intervals was 29%, which is at the lower end of values reported from other tropical forest sites, but which is higher than in most temperate forests. Coefficients of variation of solute deposition ranged from 29% to 52%. This heterogeneity of solute deposition is neither particularly high nor particularly tow compared with a range of tropical and temperate forest ecosystems. We observed an increase in solute deposition variability with the progressing wet season, which was explained by a negative correlation between heterogeneity of solute deposition and antecedent dry period. The temporal stability of throughfall. patterns was Low during the early wet season, but gained in stability as the wet season progressed. We suggest that rapid plant growth at the beginning of the rainy season is responsible for the lower stability, whereas less vegetative activity during the later rainy season might favor the higher persistence of ""hot"" and ""cold"" spots of throughfall. quantities. The relatively high stability of throughfall patterns during later stages of the wet season may influence processes at the forest floor and in the soil. Solute deposition patterns showed less clear trends but all patterns displayed a short-term stability only. The weak stability of those patterns is apt to impede the formation of solute deposition -induced biochemical microhabitats in the soil. (C) 2008 Elsevier B.V. All rights reserved.