Long-term analysis of Hubbard Brook stable oxygen isotope ratios of streamwater and precipitation sulfate

In response to decreasing atmospheric emissions of sulfur (S) since the 1970s there has been a concomitant decrease in S deposition to watersheds in the Northeastern U.S. Previous study at the Hubbard Brook Experimental Forest, NH (USA) using chemical and isotopic analyzes ( $ \delta^{34} {\text{S}}_{{{\text{SO}}_{4} }} $ ) combined with modeling has suggested that there is an internal source of S within these watersheds that results in a net loss of S via sulfate in drainage waters. The current study expands these previous investigations by the utilization of δ¹⁸O analyzes of precipitation sulfate and streamwater sulfate. Archived stream and bulk precipitation samples at the Hubbard Brook Experimental Forest from 1968–2004 were analyzed for stable oxygen isotope ratios of sulfate ( $ \delta^{18} {\text{O}}_{{{\text{SO}}_{4} }} $ ). Overall decreasing temporal trends and seasonally low winter values of $ \delta^{18} {\text{O}}_{{{\text{SO}}_{4} }} $ in bulk precipitation are most likely attributed to similar trends in precipitation $ \delta^{18} {\text{O}}_{{{\text{H}}_{2} {\text{O}}}} $ values. Regional climate trends and changes in temperature control precipitation $ \delta^{18} {\text{O}}_{{{\text{H}}_{2} {\text{O}}}} $ values that are reflected in the $ \delta^{18} {\text{O}}_{{{\text{SO}}_{4} }} $ values of precipitation. The significant relationship between ambient temperature and the $ \delta^{18} {\text{O}}_{{{\text{H}}_{2} {\text{O}}}} $ values of precipitation is shown from a nearby site in Ottawa, Ontario (Canada). Although streamwater $ \delta^{18} {\text{O}}_{{{\text{SO}}_{4} }} $ values did not reveal temporal trends, a large difference between precipitation and streamwater $ \delta^{18} {\text{O}}_{{{\text{SO}}_{4} }} $ values suggest the importance of internal cycling of S especially through the large organic S pool and the concomitant effect on the $ \delta^{18} {\text{O}}_{{{\text{SO}}_{4} }} $ values in drainage waters.     

The effects of temperature and exercise training on swimming performance in juvenile qingbo (Spinibarbus sinensis)

To investigate the effects of temperature and exercise training on swimming performance in juvenile qingbo (Spinibarbus sinensis), we measured the following: (1) the resting oxygen consumption rate $ \left( {{\dot{\text{M}}\text{O}}_{{ 2 {\text{rest}}}} } \right) $ , critical swimming speed (U _crit) and active oxygen consumption rate $ \left( {{\dot{\text{M}}\text{O}}_{{ 2 {\text{active}}}} } \right) $ of fish at acclimation temperatures of 10, 15, 20, 25 and 30 °C and (2) the $ \dot{M}{\text{O}}_{{ 2 {\text{rest}}}} $ , U _crit and $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ of both exercise-trained (exhaustive chasing training for 14 days) and control fish at both low and high acclimation temperatures (15 and 25 °C). The relationship between U _crit and temperature (T) approximately followed a bell-shaped curve as temperature increased: U _crit = 8.21/{1 + [(T ? 27.2)/17.0]²} (R ² = 0.915, P < 0.001, N = 40). The optimal temperature for maximal U _crit (8.21 BL s^?1) in juvenile qingbo was 27.2 °C. Both the $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ and the metabolic scope (MS, $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} - \dot{M}{\text{O}}_{{ 2 {\text{rest}}}} $ ) of qingbo increased with temperature from 10 to 25 °C (P < 0.05), but there were no significant differences between fish acclimated to 25 and 30 °C. The relationships between $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ or MS and temperature were described as $ {\dot{\text{M}}\text{O}}_{{ 2 {\text{active}}}} = 1,214.29/\left\{ {1 + \left[ {\left( {T - 28.8} \right)/10.6} \right]^{2} } \right\}\;\left( {R^{2} = 0.911,\;P < 0.001,\;N = 40} \right) $ and MS = 972.67/{1 + [(T ? 28.0)/9.34]²} (R ² = 0.878, P < 0.001, N = 40). The optimal temperatures for $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ and MS in juvenile qingbo were 28.8 and 28.0 °C, respectively. Exercise training resulted in significant increases in both U _crit and $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ at a low temperature (P < 0.05), but training exhibited no significant effect on either U _crit or $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ at a high temperature. These results suggest that exercise training had different effects on swimming performance at different temperatures. These differences may be related to changes in aerobic metabolic capability, arterial oxygen delivery, available dissolved oxygen, imbalances in ion fluxes and stimuli to remodel tissues with changes in temperature.     

Retention of Dissolved Inorganic Nitrogen by Foliage and Twigs of Four Temperate Tree Species

Nitrogen (N) retention by tree canopies is believed to be an important process for tree nutrient uptake, and its quantification is a key issue in determining the impact of atmospheric N deposition on forest ecosystems. Due to dry deposition and retention by other canopy elements, the actual uptake and assimilation by the tree canopy is often obscured in throughfall studies. In this study, ¹⁵N-labeled solutions ( $ ^{15} {\text{NH}}_{4}^{ + } $ and $ ^{15} {\text{NO}}_{3}^{ - } $ ) were used to assess dissolved inorganic N retention by leaves/needles and twigs of European beech, pedunculate oak, silver birch, and Scots pine saplings. The effects of N form, tree species, leaf phenology, and applied $ {\text{NO}}_{3}^{ - } $ to $ {\text{NH}}_{4}^{ + } $ ratio on the N retention were assessed. Retention patterns were mainly determined by foliar uptake, except for Scots pine. In twigs, a small but significant ¹⁵N enrichment was detected for $ {\text{NH}}_{4}^{ + } $ , which was found to be mainly due to physicochemical adsorption to the woody plant surface. The mean $ {{^{15} {\text{NH}}_{4}^{ + } } \mathord{\left/ {\vphantom {{^{15} {\text{NH}}_{4}^{ + } } {^{15} {\text{NO}}_{3}^{ - } }}} \right. \kern-0em} {^{15} {\text{NO}}_{3}^{ - } }} $ retention ratio varied considerably among species and phenological stadia, which indicates that the use of a fixed ratio in the canopy budget model could lead to an over- or underestimation of the total N retention. In addition, throughfall water under each branch was collected and analyzed for $ ^{15} {\text{NH}}_{4}^{ + } $ , $ ^{15} {\text{NO}}_{3}^{ - } $ , and all major ions. Net throughfall of $ ^{15} {\text{NH}}_{4}^{ + } $ was, on average, 20 times higher than the actual retention of $ ^{15} {\text{NH}}_{4}^{ + } $ by the plant material. This difference in $ ^{15} {\text{NH}}_{4}^{ + } $ retention could not be attributed to pools and fluxes measured in this study. The retention of $ ^{15} {\text{NH}}_{4}^{ + } $ was correlated with the net throughfall of K⁺, Mg²⁺, Ca²⁺, and weak acids during leaf development and the fully leafed period, while no significant relationships were found for $ ^{15} {\text{NO}}_{3}^{ - } $ retention. This suggests that the main driving factors for $ {\text{NH}}_{4}^{ + } $ retention might be ion exchange processes during the start and middle of the growing season and passive diffusion at leaf senescence. Actual assimilation or abiotic uptake of N through leaves and twigs was small in this study, for example, 1–5% of the applied dissolved ¹⁵N, indicating that the impact of canopy N retention from wet deposition on forest productivity and carbon sequestration is likely limited.     

The average local ionization energy as a tool for identifying reactive sites on defect-containing model graphene systems

In a continuing effort to further explore the use of the average local ionization energy $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ as a computational tool, we have investigated how well $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ computed on molecular surfaces serves as a predictive tool for identifying the sites of the more reactive electrons in several nonplanar defect-containing model graphene systems, each containing one or more pentagons. They include corannulene (C₂₀H₁₀), two inverse Stone-Thrower-Wales defect-containing structures C₂₆H₁₂ and C₄₂H₁₆, and a nanotube cap model C₂₂H₆, whose end is formed by three fused pentagons. Coronene (C₂₄H₁₂) has been included as a reference planar defect-free graphene model. We have optimized the structures of these systems as well as several monohydrogenated derivatives at the B3PW91/6-31G* level, and have computed their $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ on molecular surfaces corresponding to the 0.001 au, 0.003 au and 0.005 au contours of the electronic density. We find that (1) the convex sides of the interior carbons of the nonplanar models are more reactive than the concave sides, and (2) the magnitudes of the lowest $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ surface minima (the $ {{\overline{\mathrm{I}}}_{{\mathrm{S}\text{,}\min }}} $ ) correlate well with the interaction energies for hydrogenation at these sites. These $ {{\overline{\mathrm{I}}}_{{\mathrm{S}\text{,}\min }}} $ values decrease in magnitude as the nonplanarity of the site increases, consistent with earlier studies. A practical benefit of the use of $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ is that a single calculation suffices to characterize the numerous sites on a large molecular system, such as graphene and defect-containing graphene models.

On the page number of RNA secondary structures with pseudoknots

Let ${\mathcal {S}}$ denote the set of (possibly noncanonical) base pairs {i, j} of an RNA tertiary structure; i.e. ${\{i, j\} \in \mathcal {S}}$ if there is a hydrogen bond between the ith and jth nucleotide. The page number of ${\mathcal {S}}$ , denoted ${\pi(\mathcal {S})}$ , is the minimum number k such that ${\mathcal {S}}$ can be decomposed into a disjoint union of k secondary structures. Here, we show that computing the page number is NP-complete; we describe an exact computation of page number, using constraint programming, and determine the page number of a collection of RNA tertiary structures, for which the topological genus is known. We describe an approximation algorithm from which it follows that ${\omega(\mathcal {S}) \leq \pi(\mathcal {S}) \leq \omega(\mathcal {S}) \cdot \log n}$ , where the clique number of ${\mathcal {S}, \omega(\mathcal {S})}$ , denotes the maximum number of base pairs that pairwise cross each other.     

Release of non-exchangeable {}^{{{\text{15}}}}{\text{NH}}^{{\text{ + }}}_{{\text{4}}} from subgrade, decomposed granite substrates and uptake by non-mycorrhizal and mycorrhizal California native annual grass, Vulpia microstachys

Release rates of recently fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ from non-exchangeable interlayer sites in 2:1 silicate minerals were determined for decomposed granite (DG) saprolites from three locations in California, USA. Recently-fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ release from the DG substrate was quantified by extracting diffused $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ with H-resin, as well as a native, annual grass Vulpia microstachys. The $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ release data varied with via the method of extraction, which included H-resin pre-treatments (Na⁺ or H⁺) and V. microstachys uptake (mycorrhizal inoculated or uninoculated). After 6 weeks (1008 h), more $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ was recovered from fixed interlayer positions by the H-resins as compared to uptake by V. microstachys. The H⁺ treated H-resins recovered more released $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ (≈94 mg ${\text{NH}}^{{\text{ + }}}_{{\text{4}}} - {\text{N}}\;{\text{kg}}^{1} $ or (12%) of total fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ ) in two of the three DG samples as compared to the Na⁺ treated resins, (which recovered ≈70–78 mg ${\text{NH}}^{{\text{ + }}}_{{\text{4}}} - {\text{N}}\;{\text{kg}}^{{{\text{ - 1}}}} $ (or 9–10%) of the total fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ ). The V. microstachys assimilated 8–9% of the total fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ with mycorrhizal inoculum as compared to only 2% without a mycorrhizal inoculum, over the same time period. The fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ release kinetics from the H-resin experiments were most accurately described by first order and power function models, and can be characterized as biphasic using a heterogeneous diffusion model. Uptake of both the ¹⁵N and ambient, unlabelled N from the soils was closely related to plant biomass. There was no significant difference in percent of N per unit of biomass between the control and mycorrhizal treatments. The findings presented here indicate that observed, long-term $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ release rates from DG in studies utilizing resins, may overestimate the levels of fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ made available to plants and microorganisms. Additionally, the study suggested that mycorrhizae facilitate the acquisition and plant uptake of fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ , resulting in markedly increased plant biomass production.     

Bone refilling in cortical basic multicellular units: insights into tetracycline double labelling from a computational model

Bone remodelling is carried out by ‘bone multicellular units’ ( $\text{ BMU }$ s) in which active osteoclasts and active osteoblasts are spatially and temporally coupled. The refilling of new bone by osteoblasts towards the back of the $\text{ BMU }$ occurs at a rate that depends both on the number of osteoblasts and on their secretory activity. In cortical bone, a linear phenomenological relationship between matrix apposition rate and $\text{ BMU }$ cavity radius is found experimentally. How this relationship emerges from the combination of complex, nonlinear regulations of osteoblast number and secretory activity is unknown. Here, we extend our previous mathematical model of cell development within a single cortical $\text{ BMU }$ to investigate how osteoblast number and osteoblast secretory activity vary along the $\text{ BMU }$ ’s closing cone. The mathematical model is based on biochemical coupling between osteoclasts and osteoblasts of various maturity and includes the differentiation of osteoblasts into osteocytes and bone lining cells, as well as the influence of $\text{ BMU }$ cavity shrinkage on osteoblast development and activity. Matrix apposition rates predicted by the model are compared with data from tetracycline double labelling experiments. We find that the linear phenomenological relationship observed in these experiments between matrix apposition rate and $\text{ BMU }$ cavity radius holds for most of the refilling phase simulated by our model, but not near the start and end of refilling. This suggests that at a particular bone site undergoing remodelling, bone formation starts and ends rapidly, supporting the hypothesis that osteoblasts behave synchronously. Our model also suggests that part of the observed cross-sectional variability in tetracycline data may be due to different bone sites being refilled by $\text{ BMU }$ s at different stages of their lifetime. The different stages of a $\text{ BMU }$ ’s lifetime (such as initiation stage, progression stage, and termination stage) depend on whether the cell populations within the $\text{ BMU }$ are still developing or have reached a quasi-steady state whilst travelling through bone. We find that due to their longer lifespan, active osteoblasts reach a quasi-steady distribution more slowly than active osteoclasts. We suggest that this fact may locally enlarge the Haversian canal diameter (due to a local lack of osteoblasts compared to osteoclasts) near the $\text{ BMU }$ ’s point of origin.     

Provenance variation and genetic parameters of Eucalyptus viminalis in Argentina

Genetic parameters for growth, stem straightness, pilodyn penetration, relative bark thickness and survival were estimated in a base-population of five open-pollinated provenance/progeny trials of Eucalyptus viminalis. The trials, located in northern, central and southern Buenos Aires Province, Argentina, comprised 148 open-pollinated families from 13 Australian native provenances and eight local Argentinean seedlots. The Australian native provenances come from a limited range of the natural distribution. Overall survival, based on the latest assessment of each trial, was 62.4%. Single-site analyses showed that statistically significant provenances differences (p?<?0.05) for at least one of the studied traits in three out of the five trials analyzed. The local land race performed inconsistently in this study. The average narrow-sense individual-tree heritability estimate $ \left( {{{\hat{h}}^2}} \right) $ was 0.27 for diameter and 0.17 for total height. Values of $ {\hat{h}^2} $ also increased with age. Pilodyn penetration, assessed at only one site, was more heritable $ \left( {{{\hat{h}}^2} = 0.32} \right) $ than the average of growth traits. Estimated individual-tree heritabilities were moderate to low for stem straightness (average of 0.20) and relative bark thickness (0.16). The estimated additive genetic correlations $ \left( {{{{r}}_{{A}}}} \right) $ between diameter and height were consistently high and positive ( $ {{r}_{^A}} $ average of 0.90). High additive genetic correlations were observed between growth variables and pilodyn penetration ( $ {{r}_{^A}} $ average of 0.58). Relative bark thickness showed a negative correlation with diameter $ \left( {{{{r}}_{^A}} = - 0.39} \right) $ and height $ \left( {{{{r}}_{^A}} = - 0.51} \right) $ . The average estimated additive genetic correlation between sites was high for diameter (0.67). The implications of all these parameter estimates for genetic improvement of E. viminalis in Argentina are discussed.     

Longitudinal distribution of nitrate δ15N and δ18O in two contrasting tropical rivers: implications for instream nitrogen cycling

The longitudinal variations in the nitrogen (δ¹⁵N) and oxygen (δ¹⁸O) isotopic compositions of nitrate (NO₃ ^?), the carbon isotopic composition (δ¹³C) of dissolved inorganic carbon (DIC) and the δ¹³C and δ¹⁵N of particulate organic matter were determined in two Southeast Asian rivers contrasting in the watershed geology and land use to understand internal nitrogen cycling processes. The $ \delta^{15} {\text{N}}_{{{\text{NO}}_{3} }} $ became higher longitudinally in the freshwater reach of both rivers. The $ \delta^{18} {\text{O}}_{{{\text{NO}}_{3} }} $ also increased longitudinally in the river with a relatively steeper longitudinal gradient and a less cultivated watershed, while the $ \delta^{18} {\text{O}}_{{{\text{NO}}_{3} }} $ gradually decreased in the other river. A simple model for the $ \delta^{15} {\text{N}}_{{{\text{NO}}_{3} }} $ and the $ \delta^{18} {\text{O}}_{{{\text{NO}}_{3} }} $ that accounts for simultaneous input and removal of NO₃ ^? suggested that the dynamics of NO₃ ^? in the former river were controlled by the internal production by nitrification and the removal by denitrification, whereas that in the latter river was significantly affected by the anthropogenic NO₃ ^? loading in addition to the denitrification and/or assimilation. In the freshwater-brackish transition zone, heterotrophic activities in the river water were apparently elevated as indicated by minimal dissolved oxygen, minimal δ¹³C_DIC and maximal pCO₂. The δ¹⁵N of suspended particulate nitrogen (PN) varied in parallel to the $ \delta^{15} {\text{N}}_{{{\text{NO}}_{3} }} $ there, suggesting that the biochemical recycling processes (remineralization of PN coupled to nitrification, and assimilation of NO₃ ^?-N back to PN) played dominant roles in the instream nitrogen transformation. In the brackish zone of both rivers, the $ \delta^{15} {\text{N}}_{{{\text{NO}}_{3} }} $ displayed a declining trend while the $ \delta^{18} {\text{O}}_{{{\text{NO}}_{3} }} $ increased sharply. The redox cycling of NO₃ ^?/NO₂ ^? and/or deposition of atmospheric nitrogen oxides may have been the major controlling factor for the estuarine $ \delta^{15} {\text{N}}_{{{\text{NO}}_{3} }} $ and $ \delta^{18} {\text{O}}_{{{\text{NO}}_{3} }} $ , however, the exact mechanism behind the observed trends is currently unresolved.     

Investigation of the possibility of exceeding the Greenwald density limit during ECRH in T-10

In T-10 experiments, attempts were made to significantly exceed the Greenwald limit $\bar n_{Gr} $ during high-power (P _ab=750 kW) electron-cyclotron resonance heating (ECRH) and gas puffing. Formally, the density limit $(\bar n_e )_{\lim } $ exceeding the Greenwald limit $({{(\bar n_e )_{\lim } } \mathord{\left/ {\vphantom {{(\bar n_e )_{\lim } } {\bar n_{Gr} }}} \right. \kern-0em} {\bar n_{Gr} }} = 1.8)$ was achieved for q _L=8.2. However, as q _L decreased, the ratio ${{(\bar n_e )_{\lim } } \mathord{\left/ {\vphantom {{(\bar n_e )_{\lim } } {\bar n_{Gr} }}} \right. \kern-0em} {\bar n_{Gr} }}$ also decreased, approaching unity at q _L≈3. It was suggested that the “current radius” (i.e., the radius of the magnetic surface enclosing the bulk of the plasma current I _p), rather than the limiter radius, was the parameter governing the value of $(\bar n_e )_{\lim } $ . In the ECRH experiments, no substantial degradation of plasma confinement was observed up to $\bar n_e \sim 0.9(\bar n_e )_{\lim } $ regardless of the ratio ${{(\bar n_e )_{\lim } } \mathord{\left/ {\vphantom {{(\bar n_e )_{\lim } } {\bar n_{Gr} }}} \right. \kern-0em} {\bar n_{Gr} }}$ . In different scenarios of the density growth up to $(\bar n_e )_{\lim } $ , two types of disruptions related to the density limit were observed.     

Estimation of viscosity from passage time of liquids flowing through a microchannel array

Recently, a microchannel flow analyzer (MC-FAN) has been used to study the flow properties of blood. However, the correlation between blood passage time measured by use of the MC-FAN and hemorheology has not been clarified. In this study, a simple model is proposed for estimation of liquid viscosity from the passage time t _p of liquids. The t _p data for physiological saline were well represented by the model. According to the model, the viscosity of Newtonian fluids was estimated reasonably well from the t _p data. For blood samples, although the viscosity $ \eta_{\text{mc}} $ estimated from t _p was shown to be smaller than the viscosity $ \eta_{{450{\text{s}}^{ - 1} }} $ measured by use of a rotatory viscometer at a shear rate of 450 s^?1, $ \eta_{\text{mc}} $ was correlated with $ \eta_{{450{\text{s}}^{ - 1} }} $ . An empirical equation for estimation of $ \eta_{{450{\text{s}}^{ - 1} }} $ from $ \eta_{\text{mc}} $ of blood samples is proposed.     

Theoretical study of the reaction of CH2XO (X = F, Cl, Br) radicals with the NO radical

In this paper, we focus on the multiple-channel reactions of CH₂XO (X = F, Cl, Br) radicals with the NO radical by means of direct dynamic methods. All structures of the stationary points were obtained at the MP2/6-311+G(d,p) level and vibrational frequency analysis was also performed at this level of theory. The minimum energy path (MEP) was obtained via the intrinsic reaction coordinate (IRC) theory at the MP2/6-311+G(d,p) level, and higher-level energetic information was refined by the MC-QCISD method. The rate constants for the three hydrogen abstraction reaction channels over the temperature range 200–1,500 K were calculated by the improved canonical variational transition state theory (ICVT) with a correction for small-curvature tunneling (SCT). The rate constants calculated in this manner were in good agreement with the available experimental data, and the three-parameter rate–temperature formulae for the temperature range 200–1,500 K were $ {k_{1{\text{a}} }}(T)=0.32\times {10^{-18 }}{T^{1.83 }}\exp \left( {1748.54/T} \right) $ , $ {k_{2{\text{a}} }}(T)=0.22\times {10^{-19 }}{T^{2.19 }}\exp \left( {1770.19/T} \right) $ , $ {k_{3{\text{a}} }}(T)=0.88\times {10^{-20 }}{T^{2.20 }}\exp \left( {1513.82/T} \right) $ (in units of cm³ molecule^?1?s^?1).     

Coupling-induced population synchronization in an excitatory population of subthreshold Izhikevich neurons

We consider an excitatory population of subthreshold Izhikevich neurons which exhibit noise-induced firings. By varying the coupling strength J, we investigate population synchronization between the noise-induced firings which may be used for efficient cognitive processing such as sensory perception, multisensory binding, selective attention, and memory formation. As J is increased, rich types of population synchronization (e.g., spike, burst, and fast spike synchronization) are found to occur. Transitions between population synchronization and incoherence are well described in terms of an order parameter $\mathcal{O}$ . As a final step, the coupling induces oscillator death (quenching of noise-induced spikings) because each neuron is attracted to a noisy equilibrium state. The oscillator death leads to a transition from firing to non-firing states at the population level, which may be well described in terms of the time-averaged population spike rate $\overline{R}$ . In addition to the statistical-mechanical analysis using $\mathcal{O}$ and $\overline{R}$ , each population and individual state are also characterized by using the techniques of nonlinear dynamics such as the raster plot of neural spikes, the time series of the membrane potential, and the phase portrait. We note that population synchronization of noise-induced firings may lead to emergence of synchronous brain rhythms in a noisy environment, associated with diverse cognitive functions.     

Variability in temporary emigration rates of individually marked female Weddell seals prior to first reproduction

In many species, temporary emigration (TE) is a phenomenon, often indicative of life-history characteristics such as dormancy, skipped reproduction, or partial migration, whereby certain individuals in a population are temporarily unobservable at a particular site. TE may be a flexible condition-dependent strategy that allows individuals to mitigate effects of adverse conditions. Consequently, TE rates ought to be highly variable, but sources of variations are poorly understood for most species. We used data from known-aged female Weddell seals (Leptonychotes weddellii) tagged in Erebus Bay, Antarctica, to investigate sources of variation in TE rates prior to reproduction and to evaluate possible implications for age-specific probability of first reproduction. TE rates were near 1 the year after birth, decreased to an average of 0.15 ( $ \widehat{\text{SE}} $  = 0.01) by age 8, and were similar thereafter. TE rates varied substantially from year-to-year and were lower for seals that attended reproductive colonies the previous year than for seals that did not attend (e.g., $ \overline{{\hat{\psi }_{{i,{\text{age}}\,8}}^{\text{UU}} - \hat{\psi }_{{i,\,{\text{age}}\,8}}^{\text{PU}} }} $  = 0.22). Recruitment rates were marginally greater for seals that did attend than for seals that did not attend colonies the previous year. For Weddell seals specifically, our results suggest that (1) motivation to attend colonies varied temporally, (2) as seals grew older they had increased motivation to attend even before reproductive maturity, and (3) seals appear to follow various attendance strategies. More broadly, our results support the idea of TE as a variable, condition-dependent strategy, and highlight the utility of TE models for providing population and life-history insights for diverse taxa.     

Stochastic population growth in spatially heterogeneous environments

Classical ecological theory predicts that environmental stochasticity increases extinction risk by reducing the average per-capita growth rate of populations. For sedentary populations in a spatially homogeneous yet temporally variable environment, a simple model of population growth is a stochastic differential equation dZ _t = μ Z _t dt + σ Z _t dW _t , t ≥ 0, where the conditional law of Z _t+Δt ? Z _t given Z _t = z has mean and variance approximately z μΔt and z ² σ ²Δt when the time increment Δt is small. The long-term stochastic growth rate ${\lim_{t \to \infty} t^{-1}\log Z_t}$ for such a population equals ${\mu -\frac{\sigma^2}{2}}$ . Most populations, however, experience spatial as well as temporal variability. To understand the interactive effects of environmental stochasticity, spatial heterogeneity, and dispersal on population growth, we study an analogous model ${{\bf X}_t = (X_t^1, \ldots, X_t^n)}$ , t ≥ 0, for the population abundances in n patches: the conditional law of X _t+Δt given X _t = x is such that the conditional mean of ${X_{t+\Delta t}^i - X_t^i}$ is approximately ${[x^i \mu_i + \sum_j (x^j D_{ji} - x^i D_{ij})] \Delta t}$ where μ _i is the per capita growth rate in the ith patch and D _ij is the dispersal rate from the ith patch to the jth patch, and the conditional covariance of ${X_{t+\Delta t}^i - X_t^i}$ and ${X_{t + \Delta t}^j - X_t^j}$ is approximately x ⁱ x ^j σ _ij Δt for some covariance matrix Σ = (σ _ij ). We show for such a spatially extended population that if ${S_t = X_t^1 + \cdots + X_t^n}$ denotes the total population abundance, then Y _t = X _t /S _t , the vector of patch proportions, converges in law to a random vector Y _∞ as ${t \to \infty}$ , and the stochastic growth rate ${\lim_{t \to \infty} t^{-1}\log S_t}$ equals the space-time average per-capita growth rate ${\sum_i \mu_i \mathbb{E}[Y_\infty^i]}$ experienced by the population minus half of the space-time average temporal variation ${\mathbb{E}[\sum_{i,j}\sigma_{ij}Y_\infty^i Y_\infty^j]}$ experienced by the population. Using this characterization of the stochastic growth rate, we derive an explicit expression for the stochastic growth rate for populations living in two patches, determine which choices of the dispersal matrix D produce the maximal stochastic growth rate for a freely dispersing population, derive an analytic approximation of the stochastic growth rate for dispersal limited populations, and use group theoretic techniques to approximate the stochastic growth rate for populations living in multi-scale landscapes (e.g. insects on plants in meadows on islands). Our results provide fundamental insights into “ideal free” movement in the face of uncertainty, the persistence of coupled sink populations, the evolution of dispersal rates, and the single large or several small (SLOSS) debate in conservation biology. For example, our analysis implies that even in the absence of density-dependent feedbacks, ideal-free dispersers occupy multiple patches in spatially heterogeneous environments provided environmental fluctuations are sufficiently strong and sufficiently weakly correlated across space. In contrast, for diffusively dispersing populations living in similar environments, intermediate dispersal rates maximize their stochastic growth rate.     

Growth of Wolinella succinogenes with polysulphide as terminal acceptor of phosphorylative electron transport

Polysulphide was formed according to reaction (1), when tetrathionate was (1) $${\text{S}}_4 {\text{O}}_6^{2 - } + {\text{HS}}^ - \to 2{\text{S}}_2 {\text{O}}_3^{2 - } + {\text{S(O)}} + {\text{H}}^ + $$ added to an anaerobic buffer (pH 8.5) containing excess sulphide. S(O) denotes the zero oxidation state sulphur in the polysulphide mixture S _infn ^sup2- . The addition of formate to the polysulphide solution in the presence of Wolinella succinogenes caused the reduction of polysulphide according to reaction (2). The bacteria grew in a medium containing formate and sulphide, (2) $${\text{HCO}}_2^ - + {\text{S(O)}} + {\text{H}}2{\text{O}} \to {\text{HCO}}_3^ - + {\text{HS}}^ - + {\text{H}}^ + $$ when tetrathionate was continuously added. The cell density increased proportional to reaction (3) which represents the sum of reactions (1) and (3) $${\text{HCO}}_2^ - + {\text{S}}_{\text{4}} {\text{O}}_6^{2 - } + {\text{H}}2{\text{O}} \to {\text{HCO}}_3^ - + 2{\text{S}}_{\text{2}} {\text{O}}_3^{2 - } + 2{\text{H}}^ + $$ (2). The cell yield per mol formate was nearly the same as during growth on formate and elemental sulphur, while the velocity of growth was greater. The specific activities of polysulphide reduction by formate measured with bacteria grown with tetrathionate or with elemental sulphur were consistent with the growth parameters. The results suggest that W. succinogenes grow at the expense of formate oxidation by polysulphide and that polysulphide is an intermediate during growth on formate and elemental sulphur.     

Heritability for resistance to Puccinia psidii Winter rust in Eucalyptus grandis Hill ex Maiden in Southwestern Brazil

The climates of the central and southern regions of São Paulo state in Brazil favor pathogens such as Puccinia psidii Winter, which causes a common and severe disease in Eucalyptus plantations under 2 years old. We studied genetic parameters including genotype by environment interaction (G × E) of resistance to P. psidii rust in Eucalyptus grandis at nine sites in São Paulo State. Open-pollinated progeny from ten ‘provenances’ were established in a randomized complete block design; at individual sites there were from 134 to 160 progenies, from four to eight blocks, and five to six trees per plot. Significant provenance and progeny(provenance) differences were detected, as was G × E involving progeny(provenance). However, the G × E involved little if any rank changes, indicating that selection can be done efficiently at a single site, if the disease level is sufficient. The estimated coefficient of genetic variation among the progeny within provenances $ \left( {{{{\widehat{\mathrm{CV}}}}_{\mathrm{g}}}} \right) $ was high and variable among the sites (ranging from 11 % to 36.7 %), demonstrating different expression of genetic variability among the sites. The estimated heritability at the individual-tree level $ \left( {{{\widehat{h}}^2}} \right) $ and within a plot $ \left( {\widehat{h}_{\mathrm{w}}^2} \right) $ ranged from low to intermediate (ranging from 0.04 to 0.46) and was high at the progeny-mean level $ \left( {\widehat{h}_{\mathrm{f}}^2} \right) $ (ranging from 0.30 to 0.86). Our study shows good prospects of controlling this disease by selection among and within progenies in a single site.     

Synchronicity of movement paths of barren-ground caribou and tundra wolves

Movement patterns of highly mobile animals can reveal life history strategies and ecological relationships. We hypothesized that wolves (Canis lupus) would display similar movement patterns as their prey, barren-ground caribou (Rangifer tarandus groenlandicus), and that movements of the two species would co-vary with season. We tested for interspecific movement dynamics using animal locations from wolves and caribou monitored concurrently from mid-October to June, across the Northwest Territories and Nunavut, Canada. We used a correlated random walk as a null model to test for pattern in movements and the bearing procedure to detect whether movements were consistently directional. There was a statistical difference between the movements of caribou and wolves (F _1,9 = 13.232, P = 0.005), when compared to a correlated random walk, and a significant interaction effect between season and species (F _1,9 = 6.815, P = 0.028). During winter, the movements of caribou were strongly correlated with the 80°–90° ( $\overline{X}$ X ¯ r = 0.859, SE = 0.065) and 270°–280° ( $\overline{X}$ X ¯ r = 0.875, SE = 0.059) bearing classes suggesting an east–west movement gradient. Wolf movements during winter showed large variation in direction, but were generally east to west. Peak mean correlation for caribou movements during spring was distinct at 40°–50° ( $\overline{X}$ X ¯ r = 0.978, SE = 0.006) revealing movement to the north-east calving grounds. During spring, wolf movements correlated with the 80°–90° ( $\overline{X}$ X ¯ r = 0.861, SE = 0.043) and 270°–280° ( $\overline{X}$ X ¯ r = 0.850, SE = 0.064) bearing class. Directionality of movement suggested that during winter, caribou and wolves had a similar distribution at the large spatial scales we tested. During spring migration, however, caribou and wolves employed asynchronous movement strategies. Our findings demonstrate the utility of the correlated random walk and bearing procedure for quantifying the movement patterns of co-occurring species.     

Comparative studies on the photosynthetic responses of three freshwater phytoplankton species to temperature and light regimes

The photosynthetic responses of Microcystis aeruginosa, Scenedesmus obliquus, and Cyclotella meneghiniana to temperature and light regimes were investigated. M. aeruginosa had a higher specific growth rate at 30°C than at 14 and 20°C. Its specific growth rate was the maximum among the three species at 30°C. This suggests that M. aeruginosa could predominate in a water body having high temperature. When exposed to high light, M. aeruginosa showed lower maximal photosystem II (PSII) quantum yield (Φ_M), operational PSII quantum yield ( $ \Phi_{\text{M}}^{\prime } $ ), and active reaction centers per excited cross section (RC/CS_m) than S. obliquus and C. meneghiniana. Moreover, after 2?h low light recovery at 14°C and 20°C, the recovery of Φ_M, $ \Phi_{\text{M}}^{\prime } $ and RC/CS_m in M. aeruginosa were less than the other two species. This indicates that the capacity of high light adaptation of M. aeruginosa is the lowest among the studied species at 14–20°C. When exposed to high light, C. meneghiniana had higher Φ_M and $ \Phi_{\text{M}}^{\prime } $ lost and induced higher nonphotochemical quenching at 14–20°C. This suggests that C. meneghiniana developed a higher resistance to high light under low growth temperatures. M. aeruginosa showed the lowest light compensation point among these three species, which indicates that it could utilize low light more efficiently than the other two species. Cyclic electron flow around PSII may play a role in the photoprotective mechanism of all these three species.