Modelling the effects of chorus species composition and caller density on acoustic masking interference in multispecies choruses of crickets and katydids

Natural multispecies acoustic choruses such as the dusk chorus of a tropical rain forest consist of simultaneously signalling individuals of different species whose calls travel through a common shared medium before reaching their `intended' receivers. This causes masking interference between signals and impedes signal detection, recognition and localization. The levels of acoustic overlap depend on a number of factors, including call structure, intensity, habitat-dependent signal attenuation and receiver tuning. In addition, acoustic overlaps should also depend on caller density and the species composition of choruses, including relative and absolute abundance of the different calling species. In this study, we used simulations to examine the effects of chorus species relative abundance and caller density on the levels of effective heterospecific acoustic overlap in multispecies choruses composed of the calls of five species of crickets and katydids that share the understorey of a rain forest in southern India. We found that on average species-even choruses resulted in higher levels of effective heterospecific acoustic overlap than choruses with strong dominance structures. This effect was found consistently across dominance levels ranging from 0.4 to 0.8 for larger choruses of forty individuals. For smaller choruses of twenty individuals, the effect was seen consistently for dominance levels of 0.6 and 0.8 but not 0.4. Effective acoustic overlap (EAO) increased with caller density but the manner and extent of increase depended both on the species' call structure and the acoustic context provided by the composition scenario. The Phaloria sp. experienced very low levels of EAO and was highly buffered to changes in acoustic context whereas other species experienced high FAO across contexts or were poorly buffered. These differences were not simply predictable from call structures. These simulation-based findings may have important implications for acoustic biodiversity monitoring and for the study of acoustic masking interference in natural environments. (C) 2013 Elsevier B.V. All rights reserved.

Syntheses and structural investigation of some alkali metal ion-mediated (LVO2-)-O-V (L2- = tridentate ONO ligands) species: DNA binding, photoinduced DNA cleavage and cytotoxic activities

Eight alkali metal ion-mediated dioxidovanadium(V), {(VO2L1-6)-O-V} A(H2O)n]proportional to, complexes for A = Li+, Na+, K+ and Cs+, containing tridentate aroylhydrazonate ligands coordinating via ONO donor atoms, are described. All the synthesised ligands and the metal complexes were successfully characterised by elemental analysis, IR, UV-Vis and NMR spectroscopy. X-ray crystallographic investigation of 3, 5-7 shows the presence of distorted NO4 coordination geometries for LVO2- in each case, and varying mu-oxido and/ or mu-aqua bridging with interesting variations correlated with the size of the alkali metal ions: with small Li+, no bridging-O is found but four ion aggregates are found with Na+, chains for K+ and finally, layers for Cs+. Two (5) or three-dimensional (3, 6 and 7) architectures are consolidated by hydrogen bonding. The dioxidovanadium(V) complexes were found to exhibit DNA binding activity due to their interaction with CT-DNA by the groove binding mode, with binding constants ranging from 10(3) to 10(4) M-1. Complexes 1-8 were also tested for DNA nuclease activity against pUC19 plasmid DNA which showed that 6 and 7 had the best DNA binding and photonuclease activity; these results support their good protein binding and cleavage activity with binding constants ranging from 104 to 105 M-1. Finally, the in vitro antiproliferative activity of all complexes was assayed against the HeLa cell line. Some of the complexes (2, 5, 6 and 7) show considerable activity compared to commonly used chemotherapeutic drugs. The variation in cytotoxicity of the complexes is influenced by the various functional groups attached to the aroylhydrazone derivative.

Temporal variability of forest communities: empirical estimates of population change in 4000 tree species

Long-term surveys of entire communities of species are needed to measure fluctuations in natural populations and elucidate the mechanisms driving population dynamics and community assembly. We analysed changes in abundance of over 4000 tree species in 12 forests across the world over periods of 6-28years. Abundance fluctuations in all forests are large and consistent with population dynamics models in which temporal environmental variance plays a central role. At some sites we identify clear environmental drivers, such as fire and drought, that could underlie these patterns, but at other sites there is a need for further research to identify drivers. In addition, cross-site comparisons showed that abundance fluctuations were smaller at species-rich sites, consistent with the idea that stable environmental conditions promote higher diversity. Much community ecology theory emphasises demographic variance and niche stabilisation; we encourage the development of theory in which temporal environmental variance plays a central role.

Cryptic species and Miocene diversification of Palaearctic naked-toed geckos (Squamata: Gekkonidae) in the Indian dry zone

We sampled Palaearctic naked-toed geckos from across their range in India and used two mitochondrial and two nuclear genes to reconstruct relationships within a global phylogeny. Published sequences of Peninsular Indian Hemidactylus allow us to contrast these two groups in dating analyses - providing insights into the history of the Indian dry zone. Palaearctic naked-toed geckos first moved onto the Indian Plate in the Oligocene, with higher-level diversification probably linked to collision of the Indian and Eurasian plates, and subsequent dispersal into-India and diversification with increasing Miocene aridity. An independent gekkonid radiation with species in the dry zone, Hemidactylus diversified during the same period in Peninsular India. Our results demonstrate that dry zone taxa across India may date back to at least the Miocene, with a potential historical climatic barrier between the Indus and Peninsular Indian Divisions. `Cyrtopodion' aravallense is revealed to be a complex with seven genetically and environmentally divergent lineages that began diversifying in the late Miocene, congruent with increased aridity in north-western India. This discovery of cryptic diversity in the Indian dry zone represents the first terrestrial vertebrate radiation from north-western central India and highlights how little we understand of the regions' biodiversity, emphasizing the need for systematic geographic sampling and multiline evidence to reveal true patterns of diversity. The ancestor of `Cyrtopodion' aravallense came into the nascent Indian dry zone in the Miocene and has since diversified, potentially in the absence of any sympatric scansorial rupicolous geckos. Cyrtopodion scabrum represents a unique case of commensalism and shows phylogeographic structure in its presumed native range. The taxonomic implications of our study include a number of undescribed species, recognition of `Cyrtopodion' as a distinct lineage and the non-monophyly of Altiphylax.

Halogen Bonds in Crystal Engineering: Like Hydrogen Bonds yet Different

CONSPECTUS: The halogen bond is an attractive interaction in which an electrophilic halogen atom approaches a negatively polarized species. Short halogen atom contacts in crystals have been known for around 50 years. Such contacts are found in two varieties: type I, which is symmetrical, and type II, which is bent. Both are influenced by geometric and chemical considerations. Our research group has been using halogen atom interactions as design elements in crystal engineering, for nearly 30 years. These interactions include halogen center dot center dot center dot halogen interactions (X center dot center dot center dot X) and halogen center dot center dot center dot heteroatom interactions (X center dot center dot center dot B). Many X center dot center dot center dot X and almost all X center dot center dot center dot B contacts can be classified as halogen bonds. In this Account, we illustrate examples of crystal engineering where one can build up from previous knowledge with a focus that is provided by the modern definition of the halogen bond. We also comment on the similarities and differences between halogen bonds and hydrogen bonds. These interactions are similar because the protagonist atoms halogen and hydrogen are both electrophilic in nature. The interactions are distinctive because the size of a halogen atom is of consequence when compared with the atomic sizes of, for example, C, N, and O, unlike that of a hydrogen atom. Conclusions may be drawn pertaining to the nature of X center dot center dot center dot X interactions from the Cambridge Structural Database (CSD). There is a clear geometric and chemical distinction between type I and type II, with only type II being halogen bonds. Cl/Br isostructurality is explained based on a geometric model. In parallel, experimental studies on 3,4-dichlorophenol and its congeners shed light on the nature of halogen center dot center dot center dot halogen interactions and reveal the chemical difference between Cl and Br. Variable temperature studies also show differences between type I and type II contacts. In terms of crystal design, halogen bonds offer a unique opportunity in the strength, atom size and interaction gradation; this may be used in the design of ternary cocrystals. Structural modularity in which an entire crystal structure is defined as a combination of modules is rationalized on the basis of the intermediate strength of a halogen bond. The specific directionality of the halogen bond makes it a good tool to achieve orthogonality in molecular crystals. Mechanical properties can be tuned systematically by varying these orthogonally oriented halogen center dot center dot center dot halogen interactions. In a further development, halogen bonds are shown to play a systematic role in organization of LSAMs (long range synthon aufbau module), which are bigger structural units containing multiple synthons. With a formal definition in place, this may be the right time to look at differences between halogen bonds and hydrogen bonds and exploit them in more subtle ways in crystal engineering.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

Vermicompost reduces the invasiveness of the earthworm species Dichogaster bolaui compared to the use of compost in a degraded tropical soil in Northern Vietnam

In Northern Vietnam, organic fertilization programmes are being tested to restore soil fertility and reduce soil erosion. However, the amendment of organic matter in soil is also associated with the development of the invasive earthworm species Dichogaster bolaui. The objective of this study was to investigate the influence of organic matter amendment quality (compost vs. vermicompost) on D. bolaui. Our study confirmed D. bolaui development in organic patches in the field. However, we also observed that the flat-backed millipede Asiomorpha coarctata proliferated in these organic patches. Native to Asia, this millipede species is also considered as invasive in America. Both D. bolaui and A. coarctata more rapidly colonized compost than vermicompost patches. A laboratory experiment confirmed this trend and showed the limited development of D. bolaui in vermicompost. This is probably because of the decreased palatability of this substrate to soil fauna. In conclusion, any restoration practice that aims to increase the organic stocks in soils degraded by erosion should consider the quality of the organic amendment. In Northern Vietnam, vermicompost may be the preferred substrate for restoring soils while limiting the spread of D. bolaui. (C) 2014 Elsevier Masson SAS. All rights reserved.

Interferon-gamma induced cell death: Regulation and contributions of nitric oxide, cJun N-terminal kinase, reactive oxygen species and peroxynitrite

Interferon-gamma (Ifn gamma), a known immunomodulatory cytokine, regulates cell proliferation and survival. In this study, the mechanisms leading to the selective susceptibility of some tumor cells to Ifn gamma were deciphered. Seven different mouse tumor cell lines tested demonstrated upregulation of MHC class I to variable extents with Ifn gamma; however, only the cell lines, H6 hepatoma and L929 fibrosarcoma, that produce higher amounts of nitric oxide (NO) and reactive oxygen species (ROS) are sensitive to Ifn gamma-induced cell death. NO inhibitors greatly reduce Ifn gamma-induced ROS; however, ROS inhibitors did not affect the levels of Ifn gamma-induced NO, demonstrating that NO regulates ROS. Consequently, NO inhibitors are more effective, compared to ROS inhibitors, in reducing Ifn gamma-induced cell death. Further analysis revealed that Ifn gamma induces peroxynitrite and 3-nitrotyrosine amounts and a peroxynitrite scavenger, FeTPPS, reduces cell death. Ifn gamma treatment induces the phosphorylation of c-jun N-terminal kinase (Jnk) in H6 and L929 but not CT26, a colon carcinoma cell line, which is resistant to Ifn gamma-mediated death. Jnk activation downstream to NO leads to induction of ROS, peroxynitrite and cell death in response to Ifn gamma. Importantly, three cell lines tested, i.e. CT26, EL4 and Neuro2a, that are resistant to cell death with Ifn gamma alone become sensitive to the combination of Ifn gamma and NO donor or ROS inducer in a peroxynitrite-dependent manner. Overall, this study delineates the key roles of NO as the initiator and Jnk, ROS, and peroxynitrite as the effectors during Ifn gamma-mediated cell death. The implications of these findings in the Ifn gamma-mediated treatment of malignancies are discussed. (C) 2014 Elsevier B.V. All rights reserved.

Interferon-gamma induced cell death: Regulation and contributions of nitric oxide, cJun N-terminal kinase, reactive oxygen species and peroxynitrite

Interferon-gamma (Ifn gamma), a known immunomodulatory cytokine, regulates cell proliferation and survival. In this study, the mechanisms leading to the selective susceptibility of some tumor cells to Ifn gamma were deciphered. Seven different mouse tumor cell lines tested demonstrated upregulation of MHC class I to variable extents with Ifn gamma; however, only the cell lines, H6 hepatoma and L929 fibrosarcoma, that produce higher amounts of nitric oxide (NO) and reactive oxygen species (ROS) are sensitive to Ifn gamma-induced cell death. NO inhibitors greatly reduce Ifn gamma-induced ROS; however, ROS inhibitors did not affect the levels of Ifn gamma-induced NO, demonstrating that NO regulates ROS. Consequently, NO inhibitors are more effective, compared to ROS inhibitors, in reducing Ifn gamma-induced cell death. Further analysis revealed that Ifn gamma induces peroxynitrite and 3-nitrotyrosine amounts and a peroxynitrite scavenger, FeTPPS, reduces cell death. Ifn gamma treatment induces the phosphorylation of c-jun N-terminal kinase (Jnk) in H6 and L929 but not CT26, a colon carcinoma cell line, which is resistant to Ifn gamma-mediated death. Jnk activation downstream to NO leads to induction of ROS, peroxynitrite and cell death in response to Ifn gamma. Importantly, three cell lines tested, i.e. CT26, EL4 and Neuro2a, that are resistant to cell death with Ifn gamma alone become sensitive to the combination of Ifn gamma and NO donor or ROS inducer in a peroxynitrite-dependent manner. Overall, this study delineates the key roles of NO as the initiator and Jnk, ROS, and peroxynitrite as the effectors during Ifn gamma-mediated cell death. The implications of these findings in the Ifn gamma-mediated treatment of malignancies are discussed. (C) 2014 Elsevier B.V. All rights reserved.

A Comparative Assessment of Hand Preference in Captive Red Howler Monkeys, Alouatta seniculus and Yellow-Breasted Capuchin Monkeys, Sapajus xanthosternos

There are two major theories that attempt to explain hand preference in non-human primates-the `task complexity' theory and the `postural origins' theory. In the present study, we proposed a third hypothesis to explain the evolutionary origin of hand preference in non-human primates, stating that it could have evolved owing to structural and functional adaptations to feeding, which we refer to as the `niche structure' hypothesis. We attempted to explore this hypothesis by comparing hand preference across species that differ in the feeding ecology and niche structure: red howler monkeys, Alouatta seniculus and yellow-breasted capuchin monkeys, Sapajus xanthosternos. The red howler monkeys used the mouth to obtain food more frequently than the yellow-breasted capuchin monkeys. The red howler monkeys almost never reached for food presented on the opposite side of a wire mesh or inside a portable container, whereas the yellow-breasted capuchin monkeys reached for food presented in all four spatial arrangements (scattered, on the opposite side of a wire mesh, inside a suspended container, and inside a portable container). In contrast to the red howler monkeys that almost never acquired bipedal and clinging posture, the yellow-breasted capuchin monkeys acquired all five body postures (sitting, bipedal, tripedal, clinging, and hanging). Although there was no difference between the proportion of the red howler monkeys and the yellow-breasted capuchin monkeys that preferentially used one hand, the yellow-breasted capuchin monkeys exhibited an overall weaker hand preference than the red howler monkeys. Differences in hand preference diminished with the increasing complexity of the reaching-for-food tasks, i.e., the relatively more complex tasks were perceived as equally complex by both the red howler monkeys and the yellow-breasted capuchin monkeys. These findings suggest that species-specific differences in feeding ecology and niche structure can influence the perception of the complexity of the task and, consequently, hand preference.

Importance of intraspecifically gregarious species in a tropical bird community

In both single- and mixed-species social groups, certain participants are known to play important roles in providing benefits. Identifying these participants is critical for understanding group dynamics, but is often difficult with large roving social groups in the wild. Here, we develop a new approach to characterize roles in social groups and apply it to mixed-species bird flocks (flocks hereafter) in an Indian tropical evergreen forest. Two types of species, namely intraspecifically gregarious and sallying species, are thought to play important roles in flocks because studies have shown they attract other flock participants. However, it is unclear why these types are attractive and whether they are essential for flock formation. We address these questions by focusing on the composition of the subset of flocks containing only two species each. In two-species flocks, it is reasonable to assume that at least one species obtains some kind of benefit. Therefore, only those species combinations that result in benefit to at least one species should occur as two-species flocks. Using data from 540 flocks overall, of which 158 were two-species flocks, we find that intraspecifically gregarious species are disproportionately represented in two-species flocks and always lead flocks when present, and that flocks containing them are joined significantly more by other species. Our results suggest that intraspecifically gregarious species are likely to be the primary benefit providers in flocks and are important for tropical flock formation. Our study also provides a new approach to understanding importance in other mixed-species and single-species social groups.

Species and acoustic diversity of bats in a palaeotropical wet evergreen forest in southern India

The Western Ghats of India is among the top 25 biodiversity hotspots in the world. About 43% of the reported 117 bat species in India are found in this region, but few quantitative studies of bat echolocation calls and diversity have been carried out here thus far. A quantitative study of bat diversity was therefore conducted using standard techniques, including mist-netting, acoustical and roost surveys in the wet evergreen forests of Kudremukh National Park in the Western Ghats of Karnataka. A total of 106 bats were caught over 108 sampling nights, representing 17 species, 3 belonging to Megachiroptera and 14 to Microchiroptera. Acoustical and roost surveys added three more species, two from Microchiroptera and one from Megachiroptera. Of these 20 species, 4 belonged to the family Pteropodidae, 10 to Vespertilionidae, 3 to Rhinolophidae, 2 to Megadermatidae and 1 to Hipposideridae. We recorded the echolocation calls of 13 of the 16 microchiropteran species, of which the calls of 4 species (Pipistrellus coromandra, Pipistrellus affinis, Pipistrellus ceylonicus and Harpiocephalus harpia) have been recorded for the first time. Discriminant function analyses of the calls of 11 species provided 91.7% correct classification of individuals to their respective species, indicating that the echolocation calls could be used successfully for non-invasive acoustic surveys and monitoring of bat species in the future.

A phylogeny of the only ground-dwelling radiation of Cyrtodactylus (Squamata, Gekkonidae): diversification of Geckoella across peninsular India and Sri Lanka

The subgenus Geckoella, the only ground-dwelling radiation within Cyrtodactylus, closely overlaps in distribution with brookii group Hemidactylus in peninsular India and Sri Lanka. Both groups have Oligocene origins, the latter with over thrice as many described species. The striking difference in species richness led us to believe that Geckoella diversity is underestimated, and we sampled for Geckoella across peninsular India. A multi-locus phylogeny reveals Geckoella diversity is hugely underestimated, with at least seven undescribed species, doubling previously known richness. Strikingly, the new species correspond to cryptic lineages within described Indian species (complexes); a number of these endemic lineages from the hills of peninsular India outside the Western Ghats, highlighting the undocumented diversity of the Indian dry zone. The Geckoella phylogeny demonstrates deep splits between the Indian species and Sri Lankan G. triedrus, and between Indian dry and wet zone clades, dating back to the late Oligocene. Geckoella and brookii group Hemidactylus show contrasting diversification patterns. Geckoella shows signals of niche conservatism and appears to have retained its ancestral forest habitat. The late Miocene burst in speciation in Geckoella may be linked to the expansion of rain forests during the mid-Miocene climatic optimum and subsequent fragmentation with increasing late Miocene aridification. (C) 2014 Elsevier Inc. All rights reserved.

Lineage delimitation and description of nine new species of bush frogs (Anura: Raorchestes, Rhacophoridae) from the Western Ghats Escarpment

Bush frogs of the genus Raorchestes are distributed mainly in the Western Ghats Escarpment of Peninsular India. The inventory of species in this genus is incomplete and there is ambiguity in the systematic status of species recognized by morphological criteria. To address the dual problem of taxon sampling and systematic uncertainty in bush frogs, we used a large-scale spatial sampling design, explicitly incorporating the geographic and ecological heterogeneity of the Western Ghats. We then used a hierarchical multi-criteria approach by combining mitochondrial phylogeny, genetic distance, geographic range, morphology and advertisement call to delimit bush frog lineages. Our analyses revealed the existence of a large number of new lineages with varying levels of genetic divergence. Here, we provide diagnoses and descriptions for nine lineages that exhibit divergence across multiple axes. The discovery of new lineages that exhibit high divergence across wide ranges of elevation and across the major massifs highlights the large gaps in historical sampling. These discoveries underscore the significance of addressing inadequate knowledge of species distribution, namely the ``Wallacean shortfall'', in addressing the problem of taxon sampling and unknown diversity in tropical hotspots. A biogeographically informed sampling and analytical approach was critical in detecting and delineating lineages in a consistent manner across the genus. Through increased taxon sampling, we were also able to discern a number of well-supported sub-clades that were either unresolved or absent in earlier phylogenetic reconstructions and identify a number of shallow divergent lineages which require further examination for assessment of their taxonomic status.

And yet it shrinks: A novel method for correcting bias in forest tree growth estimates caused by water-induced fluctuations

Accuracy in tree woody growth estimates is important to global carbon budget estimation and climate-change science. Tree growth in permanent sampling plots (PSPs) is commonly estimated by measuring stem diameter changes, but this method is susceptible to bias resulting from water-induced reversible stem shrinkage. In the absence of bias correction, temporal variability in growth is likely to be overestimated and incorrectly attributed to fluctuations in resource availability, especially in forests with high seasonal and inter-annual variability in water. We propose and test a novel approach for estimating and correcting this bias at the community level. In a 50-ha PSP from a seasonally dry tropical forest in southern India, where tape measurements have been taken every four years from 1988 to 2012, for nine trees we estimated bias due to reversible stem shrinkage as the difference between woody growth measured using tree rings and that estimated from tape. We tested if the bias estimated from these trees could be used as a proxy to correct bias in tape-based growth estimates at the PSP scale. We observed significant shrinkage-related bias in the growth estimates of the nine trees in some censuses. This bias was strongly linearly related to tape-based growth estimates at the level of the PSP, and could be used as a proxy. After bias was corrected, the temporal variance in growth rates of the PSP decreased, while the effect of exceptionally dry or wet periods was retained, indicating that at least a part of the temporal variability arose from reversible shrinkage-related bias. We also suggest that the efficacy of the bias correction could be improved by measuring the proxy on trees that belong to different size classes and census timing, but not necessarily to different species. Our approach allows for reanalysis - and possible reinterpretation of temporal trends in tree growth, above ground biomass change, or carbon fluxes in forests, and their relationships with resource availability in the context of climate change. (C) 2014 Elsevier B.V. All rights reserved.