Exact Solution of a Drop-push Model for Percolation

Motivated by a computer science algorithm known as `linear probing with hashing' we study a new type of percolation model whose basic features include a sequential `dropping' of particles on a substrate followed by their transport via a `pushing' mechanism. Our exact solution in one dimension shows that, unlike the ordinary random percolation model, the drop-push model has nontrivial spatial correlations generated by the dynamics itself. The critical exponents in the drop-push model are also different from that of the ordinary percolation. The relevance of our results to computer science is pointed out.Comment: 4 pages revtex, 2 eps figure

Boreal forest CO2 exchange and evapotranspiration predicted by nine ecosystem process models: Intermodel comparisons and relationships to field measurements

Nine ecosystem process models were used to predict CO2 and water vapor exchanges by a 150-year-old black spruce forest in central Canada during 1994–1996 to evaluate and improve the models. Three models had hourly time steps, five had daily time steps, and one had monthly time steps. Model input included site ecosystem characteristics and meteorology. Model predictions were compared to eddy covariance (EC) measurements of whole-ecosystem CO2exchange and evapotranspiration, to chamber measurements of nighttime moss-surface CO2release, and to ground-based estimates of annual gross primary production, net primary production, net ecosystem production (NEP), plant respiration, and decomposition. Model-model differences were apparent for all variables. Model-measurement agreement was good in some cases but poor in others. Modeled annual NEP ranged from −11 g C m−2 (weak CO2source) to 85 g C m−2 (moderate CO2 sink). The models generally predicted greater annual CO2sink activity than measured by EC, a discrepancy consistent with the fact that model parameterizations represented the more productive fraction of the EC tower “footprint.” At hourly to monthly timescales, predictions bracketed EC measurements so median predictions were similar to measurements, but there were quantitatively important model-measurement discrepancies found for all models at subannual timescales. For these models and input data, hourly time steps (and greater complexity) compared to daily time steps tended to improve model-measurement agreement for daily scale CO2 exchange and evapotranspiration (as judged by root-mean-squared error). Model time step and complexity played only small roles in monthly to annual predictions

Quantum curves for Hitchin fibrations and the Eynard-Orantin theory

We generalize the topological recursion of Eynard-Orantin (2007) to the family of spectral curves of Hitchin fibrations. A spectral curve in the topological recursion, which is defined to be a complex plane curve, is replaced with a generic curve in the cotangent bundle $T^*C$ of an arbitrary smooth base curve $C$. We then prove that these spectral curves are quantizable, using the new formalism. More precisely, we construct the canonical generators of the formal $\hbar$-deformation family of $D$-modules over an arbitrary projective algebraic curve $C$ of genus greater than $1$, from the geometry of a prescribed family of smooth Hitchin spectral curves associated with the $SL(2,\mathbb{C})$-character variety of the fundamental group $\pi_1(C)$. We show that the semi-classical limit through the WKB approximation of these $\hbar$-deformed $D$-modules recovers the initial family of Hitchin spectral curves.Comment: 34 page

Ecosystem carbon 7 dioxide fluxes after disturbance in forests of North America

Disturbances are important for renewal of North American forests. Here we summarize more than 180 site years of eddy covariance measurements of carbon dioxide flux made at forest chronosequences in North America. The disturbances included stand-replacing fire (Alaska, Arizona, Manitoba, and Saskatchewan) and harvest (British Columbia, Florida, New Brunswick, Oregon, Quebec, Saskatchewan, and Wisconsin) events, insect infestations (gypsy moth, forest tent caterpillar, and mountain pine beetle), Hurricane Wilma, and silvicultural thinning (Arizona, California, and New Brunswick). Net ecosystem production (NEP) showed a carbon loss from all ecosystems following a stand-replacing disturbance, becoming a carbon sink by 20 years for all ecosystems and by 10 years for most. Maximum carbon losses following disturbance (g C m−2y−1) ranged from 1270 in Florida to 200 in boreal ecosystems. Similarly, for forests less than 100 years old, maximum uptake (g C m−2y−1) was 1180 in Florida mangroves and 210 in boreal ecosystems. More temperate forests had intermediate fluxes. Boreal ecosystems were relatively time invariant after 20 years, whereas western ecosystems tended to increase in carbon gain over time. This was driven mostly by gross photosynthetic production (GPP) because total ecosystem respiration (ER) and heterotrophic respiration were relatively invariant with age. GPP/ER was as low as 0.2 immediately following stand-replacing disturbance reaching a constant value of 1.2 after 20 years. NEP following insect defoliations and silvicultural thinning showed lesser changes than stand-replacing events, with decreases in the year of disturbance followed by rapid recovery. NEP decreased in a mangrove ecosystem following Hurricane Wilma because of a decrease in GPP and an increase in ER

Enumeration of simple random walks and tridiagonal matrices

We present some old and new results in the enumeration of random walks in one dimension, mostly developed in works of enumerative combinatorics. The relation between the trace of the $n$-th power of a tridiagonal matrix and the enumeration of weighted paths of $n$ steps allows an easier combinatorial enumeration of the paths. It also seems promising for the theory of tridiagonal random matrices .Comment: several ref.and comments added, misprints correcte

Evanescence in Coined Quantum Walks

In this paper we complete the analysis begun by two of the authors in a previous work on the discrete quantum walk on the line [J. Phys. A 36:8775-8795 (2003) quant-ph/0303105 ]. We obtain uniformly convergent asymptotics for the "exponential decay'' regions at the leading edges of the main peaks in the Schr{\"o}dinger (or wave-mechanics) picture. This calculation required us to generalise the method of stationary phase and we describe this extension in some detail, including self-contained proofs of all the technical lemmas required. We also rigorously establish the exact Feynman equivalence between the path-integral and wave-mechanics representations for this system using some techniques from the theory of special functions. Taken together with the previous work, we can now prove every theorem by both routes.Comment: 32 pages AMS LaTeX, 5 figures in .eps format. Rewritten in response to referee comments, including some additional references. v3: typos fixed in equations (131), (133) and (134). v5: published versio

Chebyshev type lattice path weight polynomials by a constant term method

We prove a constant term theorem which is useful for finding weight polynomials for Ballot/Motzkin paths in a strip with a fixed number of arbitrary `decorated' weights as well as an arbitrary `background' weight. Our CT theorem, like Viennot's lattice path theorem from which it is derived primarily by a change of variable lemma, is expressed in terms of orthogonal polynomials which in our applications of interest often turn out to be non-classical. Hence we also present an efficient method for finding explicit closed form polynomial expressions for these non-classical orthogonal polynomials. Our method for finding the closed form polynomial expressions relies on simple combinatorial manipulations of Viennot's diagrammatic representation for orthogonal polynomials. In the course of the paper we also provide a new proof of Viennot's original orthogonal polynomial lattice path theorem. The new proof is of interest because it uses diagonalization of the transfer matrix, but gets around difficulties that have arisen in past attempts to use this approach. In particular we show how to sum over a set of implicitly defined zeros of a given orthogonal polynomial, either by using properties of residues or by using partial fractions. We conclude by applying the method to two lattice path problems important in the study of polymer physics as models of steric stabilization and sensitized flocculation.Comment: 27 pages, 14 figure

The Impact of Boreal Forest Fire on Climate Warming

We report measurements and analysis of a boreal forest fire, integrating the effects of greenhouse gases, aerosols, black carbon deposition on snow and sea ice, and postfire changes in surface albedo. The net effect of all agents was to increase radiative forcing during the first year (34 ± 31 Watts per square meter of burned area), but to decrease radiative forcing when averaged over an 80-year fire cycle (–2.3 ± 2.2 Watts per square meter) because multidecadal increases in surface albedo had a larger impact than fire-emitted greenhouse gases. This result implies that future increases in boreal fire may not accelerate climate warming

Looking deeper into the soil : biophysical controls and seasonal lags of soil CO2 production and efflux

Author Posting. © Ecological Society of America, 2010. This article is posted here by permission of Ecological Society of America for personal use, not for redistribution. The definitive version was published in Ecological Applications 20 (2010): 1569–1582, doi:10.1890/09-0693.1.We seek to understand how biophysical factors such as soil temperature (Ts), soil moisture (θ), and gross primary production (GPP) influence CO2 fluxes across terrestrial ecosystems. Recent advancements in automated measurements and remote-sensing approaches have provided time series in which lags and relationships among variables can be explored. The purpose of this study is to present new applications of continuous measurements of soil CO2 efflux (F0) and soil CO2 concentrations measurements. Here we explore how variation in Ts, θ, and GPP (derived from NASA's moderate-resolution imaging spectroradiometer [MODIS]) influence F0 and soil CO2 production (Ps). We focused on seasonal variation and used continuous measurements at a daily timescale across four vegetation types at 13 study sites to quantify: (1) differences in seasonal lags between soil CO2 fluxes and Ts, θ, and GPP and (2) interactions and relationships between CO2 fluxes with Ts, θ, and GPP. Mean annual Ts did not explain annual F0 and Ps among vegetation types, but GPP explained 73% and 30% of the variation, respectively. We found evidence that lags between soil CO2 fluxes and Ts or GPP provide insights into the role of plant phenology and information relevant about possible timing of controls of autotrophic and heterotrophic processes. The influences of biophysical factors that regulate daily F0 and Ps are different among vegetation types, but GPP is a dominant variable for explaining soil CO2 fluxes. The emergence of long-term automated soil CO2 flux measurement networks provides a unique opportunity for extended investigations into F0 and Ps processes in the near future.Data collection was possible thanks to NASA, the NSF Center for Embedded Networked Sensing (CCR-0120778), DOE (DE-FG02-03ER63638), CONACyT, UCMEXUS, NSF (EF-0410408), NSF-LTER, KAKENHI (12878089 and 13480150), the Academy of Finland (213093), the Austrian Science Fund (FWF, P18756-B16), the Kearney Foundation, the Canadian Foundation for Climate and Atmospheric Sciences (CFCAS), and the Natural Science and Engineering Research Council of Canada (NSERC). R. Vargas was supported by grant DEB-0639235 during the preparation of this manuscript