Comparison of different objective functions for parameterization of simple respiration models

The eddy covariance measurements of carbon dioxide fluxes collected around the world offer a rich source for detailed data analysis. Simple, aggregated models are attractive tools for gap filling, budget calculation, and upscaling in space and time. Key in the application of these models is their parameterization and a robust estimate of the uncertainty and reliability of their predictions. In this study we compared the use of ordinary least squares (OLS) and weighted absolute deviations (WAD, which is the objective function yielding maximum likelihood parameter estimates with a double exponential error distribution) as objective functions within the annual parameterization of two respiration models: the Q10 model and the Lloyd and Taylor model. We introduce a new parameterization method based on two nonparametric tests in which model deviation (Wilcoxon test) and residual trend analyses (Spearman test) are combined. A data set of 9 years of flux measurements was used for this study. The analysis showed that the choice of the objective function is crucial, resulting in differences in the estimated annual respiration budget of up to 40%. The objective function should be tested thoroughly to determine whether it is appropriate for the application for which the model will be used. If simple models are used to estimate a respiration budget, a trend test is essential to achieve unbiased estimates over the year. The analyses also showed that the parameters of the Lloyd and Taylor model are highly correlated and difficult to determine precisely, thereby limiting the physiological interpretability of the parameter

Characteristics of Magnetoplasmas Semiannual Status Report No. 12, May 1 - Oct. 31, 1965

Magnetoplasma characteristics - anomalous diffusion across magnetic field, heat conduction in plasma, cesium plasma generator, and electron velocity distribution function in magnetoplasma

Influence of the Soret effect on convection of binary fluids

Convection in horizontal layers of binary fluids heated from below and in particular the influence of the Soret effect on the bifurcation properties of extended stationary and traveling patterns that occur for negative Soret coupling is investigated theoretically. The fixed points corresponding to these two convection structures are determined for realistic boundary conditions with a many mode Galerkin scheme for temperature and concentration and an accurate one mode truncation of the velocity field. This solution procedure yields the stable and unstable solutions for all stationary and traveling patterns so that complete phase diagrams for the different convection types in typical binary liquid mixtures can easily be computed. Also the transition from weakly to strongly nonlinear states can be analyzed in detail. An investigation of the concentration current and of the relevance of its constituents shows the way for a simplification of the mode representation of temperature and concentration field as well as for an analytically manageable few mode description.Comment: 30 pages, 12 figure

Renewal of Collaborative Research: Economically Viable Forest Harvesting Practices That Increase Carbon Sequestration

Forests provide wildlife habitat, water and air purification, climate moderation, and timber and nontimber products. Concern about climate change has put forests in the limelight as sinks of atmospheric carbon. The C stored in the global vegetation, mostly in forests, is nearly equivalent to the amount present in atmospheric CO{sub 2}. Both voluntary and government-mandated carbon trading markets are being developed and debated, some of which include C sequestration resulting from forest management as a possible tradeable commodity. However, uncertainties regarding sources of variation in sequestration rates, validation, and leakage remain significant challenges for devising strategies to include forest management in C markets. Hence, the need for scientifically-based information on C sequestration by forest management has never been greater. The consequences of forest management on the US carbon budget are large, because about two-thirds of the {approx}300 million hectare US forest resource is classified as 'commercial forest.' In most C accounting budgets, forest harvesting is usually considered to cause a net release of C from the terrestrial biosphere to the atmosphere. However, forest management practices could be designed to meet the multiple goals of providing wood and paper products, creating economic returns from natural resources, while sequestering C from the atmosphere. The shelterwood harvest strategy, which removes about 30% of the basal area of the overstory trees in each of three successive harvests spread out over thirty years as part of a stand rotation of 60-100 years, may improve net C sequestration compared to clear-cutting because: (1) the average C stored on the land surface over a rotation increases, (2) harvesting only overstory trees means that a larger fraction of the harvested logs can be used for long-lived sawtimber products, compared to more pulp resulting from clearcutting, (3) the shelterwood cut encourages growth of subcanopy trees by opening up the forest canopy to increasing light penetration. Decomposition of onsite harvest slash and of wastes created during timber processing releases CO{sub 2} to the atmosphere, thus offsetting some of the C sequestered in vegetation. Decomposition of soil C and dead roots may also be temporarily stimulated by increased light penetration and warming of the forest floor. Quantification of these processes and their net effect is needed. We began studying C sequestration in a planned shelterwood harvest at the Howland Forest in central Maine in 2000. The harvest took place in 2002 by the International Paper Corporation, who assisted us to track the fates of harvest products (Scott et al., 2004, Environmental Management 33: S9-S22). Here we present the results of intensive on-site studies of the decay of harvest slash, soil respiration, growth of the remaining trees, and net ecosystem exchange (NEE) of CO{sub 2} during the first six years following the harvest. These results are combined with calculations of C in persisting off-site harvest products to estimate the net C consequences to date of this commercial shelterwood harvest operation. Tower-based eddy covariance is an ideal method for this study, as it integrates all C fluxes in and out of the forest over a large 'footprint' area and can reveal how the net C flux, as well as gross primary productivity and respiration, change following harvest. Because the size of this experiment precludes large-scale replication, we are use a paired-airshed approach, similar to classic large-scale paired watershed experiments. Measurements of biomass and C fluxes in control and treatment stands were compared during a pre-treatment calibration period, and then divergence from pre-treatment relationships between the two sites measured after the harvest treatment. Forests store carbon (C) as they accumulate biomass. Many forests are also commercial sources of timber and wood fiber. In most C accounting budgets, forest harvesting is usually considered to cause a net release of C from the terrestrial biosphere to the atmosphere. However, it might also be possible for commercial use of forests to contribute to terrestrial sequestration of C. The objective of the our research project is to determine whether shelterwood cutting regimes now being adopted in the commercial forests of Maine and other areas of the country can achieve these multiple goals

Bistability of Slow and Fast Traveling Waves in Fluid Mixtures

The appearence of a new type of fast nonlinear traveling wave states in binary fluid convection with increasing Soret effect is elucidated and the parameter range of their bistability with the common slower ones is evaluated numerically. The bifurcation behavior and the significantly different spatiotemporal properties of the different wave states - e.g. frequency, flow structure, and concentration distribution - are determined and related to each other and to a convenient measure of their nonlinearity. This allows to derive a limit for the applicability of small amplitude expansions. Additionally an universal scaling behavior of frequencies and mixing properties is found. PACS: 47.20.-k, 47.10.+g, 47.20.KyComment: 4 pages including 5 Postscript figure

Structure and oxidation kinetics of the Si(100)-SiO2 interface

We present first-principles calculations of the structural and electronic properties of Si(001)-SiO2 interfaces. We first arrive at reasonable structures for the c-Si/a-SiO2 interface via a Monte-Carlo simulated annealing applied to an empirical interatomic potential, and then relax these structures using first-principles calculations within the framework of density-functional theory. We find a transition region at the interface, having a thickness on the order of 20\AA, in which there is some oxygen deficiency and a corresponding presence of sub-oxide Si species (mostly Si^+2 and Si^+3). Distributions of bond lengths and bond angles, and the nature of the electronic states at the interface, are investigated and discussed. The behavior of atomic oxygen in a-SiO2 is also investigated. The peroxyl linkage configuration is found to be lower in energy than interstitial or threefold configurations. Based on these results, we suggest a possible mechanism for oxygen diffusion in a-SiO2 that may be relevant to the oxidation process.Comment: 7 pages, two-column style with 6 postscript figures embedded. Uses REVTEX and epsf macros. Also available at http://www.physics.rutgers.edu/~dhv/preprints/index.html#ng_sio

Canopy nitrogen, carbon assimilation, and albedo in temperate and boreal forests: Functional relations and potential climate feedbacks

The availability of nitrogen represents a key constraint on carbon cycling in terrestrial ecosystems, and it is largely in this capacity that the role of N in the Earth\u27s climate system has been considered. Despite this, few studies have included continuous variation in plant N status as a driver of broad-scale carbon cycle analyses. This is partly because of uncertainties in how leaf-level physiological relationships scale to whole ecosystems and because methods for regional to continental detection of plant N concentrations have yet to be developed. Here, we show that ecosystem CO2 uptake capacity in temperate and boreal forests scales directly with whole-canopy N concentrations, mirroring a leaf-level trend that has been observed for woody plants worldwide. We further show that both CO2 uptake capacity and canopy N concentration are strongly and positively correlated with shortwave surface albedo. These results suggest that N plays an additional, and overlooked, role in the climate system via its influence on vegetation reflectivity and shortwave surface energy exchange. We also demonstrate that much of the spatial variation in canopy N can be detected by using broad-band satellite sensors, offering a means through which these findings can be applied toward improved application of coupled carbon cycle–climate models

Polarization Dependent Excitation and High Harmonic Generation from Intense Mid-IR Laser Pulses in ZnO

The generation of high order harmonics from femtosecond mid-IR laser pulses in ZnO has shown great potential to reveal new insight into the ultrafast electron dynamics on a few femtosecond timescale. In this work we report on the experimental investigation of photoluminescence and high-order harmonic generation (HHG) in a ZnO single crystal and polycrystalline thin film irradiated with intense femtosecond mid-IR laser pulses. The ellipticity dependence of the HHG process is experimentally studied up to the 17th harmonic order for various driving laser wavelengths in the spectral range 3–4 µm. Interband Zener tunneling is found to exhibit a significant excitation efficiency drop for circularly polarized strong-field pump pulses. For higher harmonics with energies larger than the bandgap, the measured ellipticity dependence can be quantitatively described by numerical simulations based on the density matrix equations. The ellipticity dependence of the below and above ZnO band gap harmonics as a function of the laser wavelength provides an efficient method for distinguishing the dominant HHG mechanism for different harmonic orders

Universal Crossover between Efros-Shklovskii and Mott Variable-Range-Hopping Regimes

A universal scaling function, describing the crossover between the Mott and the Efros-Shklovskii hopping regimes, is derived, using the percolation picture of transport in strongly localized systems. This function is agrees very well with experimental data. Quantitative comparison with experiment allows for the possible determination of the role played by polarons in the transport.Comment: 7 pages + 1 figure, Revte