A comprehensive, multi-process box-model approach to glacial-interglacial carbon cycling

The canonical question of which physical, chemical or biological mechanisms were responsible for oceanic uptake of atmospheric CO2 during the last glacial is yet unanswered. Insight from paleo proxies has led to a multitude of hypotheses but none so far have been convincingly supported in three dimensional numerical modelling experiments. The processes that influence the CO2 uptake and export production are inter-related and too complex to solve conceptually while complex numerical models are time consuming and expensive to run which severely limits the combinations of mechanisms that can be explored. Instead, an intermediate inverse box model approach is used here in which the whole parameter space is explored. The glacial circulation and biological production states are derived from these using proxies of glacial export production and the need to draw down CO2 into the ocean. We find that circulation patterns which explain glacial observations include reduced Antarctic Bottom Water formation and high latitude mixing and to a lesser extent reduced equatorial upwelling. The proposed mechanism of CO2 uptake by an increase of eddies in the Southern Ocean, leading to a reduced residual circulation, is not supported. Regarding biological mechanisms, an increase in the nutrient utilization in either the equatorial regions or the northern polar latitudes can reduce atmospheric CO2 and satisfy proxies of glacial export production. Consistent with previous studies, CO2 is drawn down more easily through increased productivity in the Antarctic region than the sub-Antarctic, but that violates observations of lower export production there

A Standard Input Format for Computer Codes Which Solve Stochastic Programs with Recourse and a Library of Utilities to Simplify its Use

We explain our suggestions for standardizing input formats for computer codes which solve stochastic programs with recourse. The main reason to set some conventions is to allow programs implementing different methods of solution to be used interchangeably. The general philosophy behind our design is a) to remain fairly faithful to the de facto standard for the statement of LP problems established by IBM for use with MPSX and subsequently adopted by the authors of MINOS, b) to provide sufficient flexibility so that a variety of problems may be expressed in the standard format, c) to allow problems originally formulated as deterministic LP to be converted to stochastic problems with a minimum of effort, d) to permit new options to be added as the need arises, and e) to provide some routines to facilitate the task of reading files specified in the standard format

Limaria hians (Mollusca : Limacea): a neglected reef-forming keystone species

1. A key component of physical habitat along braided river systems is the exposed riverine sediment within the active zone. The relatively unmanaged, gravel-bed Fiume Tagliamento, Italy, provides the focus for exploring two ecologically important properties of exposed riverine sediments: their within-patch and between-patch variability in calibre. 2. To characterize between-patch variation in exposed riverine sediments, replicate (within-patch) samples were obtained from three geomorphologically distinct locations along 130 km of the river: bar heads along the margin of the low-flow channel, the heads of major bars across the exposed surface of the active zone, and floodplain surfaces. A photographic technique enabled rapid and consistent field sampling of the coarse sediments at bar heads along the low-flow channel margin and on major bars across the dry bed. 3. A downstream decrease in particle size and an increase in within-patch heterogeneity in sediment size were observed within bar head sediments along the margin of the low-flow channel. Comparisons between major bar and low-flow channel samples revealed greatest within-patch variability in individual sediment size indices (D50, A- and B-axes of the larger particles) at headwater sites, greatest between-patch variability in the three measured indices in the central reaches, and lowest between-patch variability at downstream sites. However, there was a distinct increase in the overall heterogeneity in particle size, which was sustained across all patches, in a downstream direction. 4. There was a clear downstream decrease in the size of floodplain sediments in the headwaters, but thereafter there was no distinct downstream trend in any of the calculated particle size indices. 5. The geomorphological controls on the observed patterns and the potential ecological significance of the patterns, particularly for plant establishment, are discussed in relation to the relative relief of the active zone, and the highly variable hydrological and climatic regime along the river

Phase transitions in the steady state behavior of mechanically perturbed spin glasses and ferromagnets

We analyze the steady state regime of systems interpolating between spin glasses and ferromagnets under a tapping dynamics recently introduced by analogy with the dynamics of mechanically perturbed granular media. A crossover from a second order to first order ferromagnetic transition as a function of the spin coupling distribution is found. The flat measure over blocked states introduced by Edwards for granular media is used to explain this scenario. Annealed calculations of the Edwards entropy are shown to qualitatively explain the nature of the phase transitions. A Monte-Carlo construction of the Edwards measure confirms that this explanation is also quantitatively accurate

Linear response of vibrated granular systems to sudden changes in the vibration intensity

The short-term memory effects recently observed in vibration-induced compaction of granular materials are studied. It is shown that they can be explained by means of quite plausible hypothesis about the mesoscopic description of the evolution of the system. The existence of a critical time separating regimes of ``anomalous'' and ``normal'' responses is predicted. A simple model fitting into the general framework is analyzed in the detail. The relationship between this work and previous studies is discussed.Comment: 10 pages, 6 figures; fixed errata, updtated reference

Numerical Study of Competing Spin-Glass and Ferromagnetic Order

Two and three dimensional random Ising models with a Gaussian distribution of couplings with variance $J$ and non-vanishing mean value $J_0$ are studied using the zero-temperature domain-wall renormalization group (DWRG). The DWRG trajectories in the ($J_0,J$) plane after rescaling can be collapsed on two curves: one for $J_0/J > r_c$ and other for $J_0/J < r_c$. In the first case the DWRG flows are toward the ferromagnetic fixed point both in two and three dimensions while in the second case flows are towards a paramagnetic fixed point and spin-glass fixed point in two and three dimensions respectively. No evidence for an extra phase is found.Comment: a bit more data is taken, 5 pages, 4 eps figures included, to appear in PR

Steady State Behavior of Mechanically Perturbed Spin Glasses and Ferromagnets

A zero temperature dynamics of Ising spin glasses and ferromagnets on random graphs of finite connectivity is considered, like granular media these systems have an extensive entropy of metastable states. We consider the problem of what energy a randomly prepared spin system falls to before becoming stuck in a metastable state. We then introduce a tapping mechanism, analogous to that of real experiments on granular media, this tapping, corresponding to flipping simultaneously any spin with probability $p$, leads to stationary regime with a steady state energy $E(p)$. We explicitly solve this problem for the one dimensional ferromagnet and $\pm J$ spin glass and carry out extensive numerical simulations for spin systems of higher connectivity. The link with the density of metastable states at fixed energy and the idea of Edwards that one may construct a thermodynamics with a flat measure over metastable states is discussed. In addition our simulations on the ferromagnetic systems reveal a novel first order transition, whereas the usual thermodynamic transition on these graphs is second order.Comment: 11 pages, 7 figure

Spin glass transition in a magnetic field: a renormalization group study

We study the transition of short range Ising spin glasses in a magnetic field, within a general replica symmetric field theory, which contains three masses and eight cubic couplings, that is defined in terms of the fields representing the replicon, anomalous and longitudinal modes. We discuss the symmetry of the theory in the limit of replica number n to 0, and consider the regular case where the longitudinal and anomalous masses remain degenerate. The spin glass transitions in zero and non-zero field are analyzed in a common framework. The mean field treatment shows the usual results, that is a transition in zero field, where all the modes become critical, and a transition in non-zero field, at the de Almeida-Thouless (AT) line, with only the replicon mode critical. Renormalization group methods are used to study the critical behavior, to order epsilon = 6-d. In the general theory we find a stable fixed-point associated to the spin glass transition in zero field. This fixed-point becomes unstable in the presence of a small magnetic field, and we calculate crossover exponents, which we relate to zero-field critical exponents. In a finite magnetic field, we find no physical stable fixed-point to describe the AT transition, in agreement with previous results of other authors.Comment: 36 pages with 4 tables. To be published in Phys. Rev.

Quantum Kinetic Theory III: Quantum kinetic master equation for strongly condensed trapped systems

We extend quantum kinetic theory to deal with a strongly Bose-condensed atomic vapor in a trap. The method assumes that the majority of the vapor is not condensed, and acts as a bath of heat and atoms for the condensate. The condensate is described by the particle number conserving Bogoliubov method developed by one of the authors. We derive equations which describe the fluctuations of particle number and phase, and the growth of the Bose-Einstein condensate. The equilibrium state of the condensate is a mixture of states with different numbers of particles and quasiparticles. It is not a quantum superposition of states with different numbers of particles---nevertheless, the stationary state exhibits the property of off-diagonal long range order, to the extent that this concept makes sense in a tightly trapped condensate.Comment: 3 figures submitted to Physical Review

Novel order parameter to describe the critical behavior of Ising spin glass models

A novel order parameter $\Phi$ for spin glasses is defined based on topological criteria and with a clear physical interpretation. $\Phi$ is first investigated for well known magnetic systems and then applied to the Edwards-Anderson $\pm J$ model on a square lattice, comparing its properties with the usual $q$ order parameter. Finite size scaling procedures are performed. Results and analyses based on $\Phi$ confirm a zero temperature phase transition and allow to identify the low temperature phase. The advantages of $\Phi$ are brought out and its physical meaning is established.Comment: 13 pages, 4 figures, to appear in Physica