Stochastic structure and individual-tree growth models

The majority of past and current individual-tree growth modelling methodologies have failed to characterise and incorporate structured stochastic components. Rather, they have relied on deterministic predictions or have added an unstructured random component to predictions. In particular, spatial stochastic structure has been neglected, despite being present in most applications of individual-tree growth models. Spatial stochastic structure (also called spatial dependence or spatial autocorrelation) eventuates when spatial influences such as competition and micro-site effects are not fully captured in models. Temporal stochastic structure (also called temporal dependence or temporal autocorrelation) eventuates when a sequence of measurements is taken on an individual-tree over time, and variables explaining temporal variation in these measurements are not included in the model. Nested stochastic structure eventuates when measurements are combined across sampling units and differences among the sampling units are not fully captured in the model. This review examines spatial, temporal, and nested stochastic structure and instances where each has been characterised in the forest biometry and statistical literature. Methodologies for incorporating stochastic structure in growth model estimation and prediction are described. Benefits from incorporation of stochastic structure include valid statistical inference, improved estimation efficiency, and more realistic and theoretically sound predictions. It is proposed in this review that individual-tree modelling methodologies need to characterise and include structured stochasticity. Possibilities for future research are discussed. (C) 2001 Elsevier Science B.V. All rights reserved.

Influence of patient age and implantation technique on the probability of re-replacement of the homograft aortic valve

Background and aim of the study: Results of valve re-replacement (reoperation) in 898 patients undergoing aortic valve replacement with cryopreserved homograft valves between 1975 and 1998 are reported. The study aim was to provide estimates of unconditional probability of valve reoperation and cumulative incidence function (actual risk) of reoperation. Methods: Valves were implanted by subcoronary insertion (n = 500), inclusion cylinder (n = 46), and aortic root replacement (n = 352). Probability of reoperation was estimated by adopting a mixture model framework within which estimates were adjusted for two risk factors: patient age at initial replacement, and implantation technique. Results: For a patient aged 50 years, the probability of reoperation in his/her lifetime was estimated as 44% and 56% for non-root and root replacement techniques, respectively. For a patient aged 70 years, estimated probability of reoperation was 16% and 25%, respectively. Given that a reoperation is required, patients with non-root replacement have a higher hazard rate than those with root replacement (hazards ratio = 1.4), indicating that non-root replacement patients tend to undergo reoperation earlier before death than root replacement patients. Conclusion: Younger patient age and root versus non-root replacement are risk factors for reoperation. Valve durability is much less in younger patients, while root replacement patients appear more likely to live longer and hence are more likely to require reoperation.

Computational immunology: The coming of age

The explosive growth in biotechnology combined with major advancesin information technology has the potential to radically transformimmunology in the postgenomics era. Not only do we now have readyaccess to vast quantities of existing data, but new data with relevanceto immunology are being accumulated at an exponential rate. Resourcesfor computational immunology include biological databases and methodsfor data extraction, comparison, analysis and interpretation. Publiclyaccessible biological databases of relevance to immunologists numberin the hundreds and are growing daily. The ability to efficientlyextract and analyse information from these databases is vital forefficient immunology research. Most importantly, a new generationof computational immunology tools enables modelling of peptide transportby the transporter associated with antigen processing (TAP), modellingof antibody binding sites, identification of allergenic motifs andmodelling of T-cell receptor serial triggering.

Halo formation in directional solidification

A new model of halo formation in directional solidification is presented. The model describes halo formation in terms of competitive growth between the halo phase and coupled eutectic in liquid with a nominal composition that follows the primary phase liquidus extension with decreasing temperature. The model distinguishes between the effects of constitutional, capillarity and (where applicable) kinetic undercooling and avoids a number of theoretical inconsistencies associated with previous models. The critical growth rate for halo formation in directionally solidified hypereutectic Al-Si alloys is calculated using the model in conjunction with models of primary phase and coupled eutectic growth from the literature. The calculated result agrees reasonably well with the experimental result of Yilmaz and Elliott (Met. Sci. 18 (1984) 362), given the use of a relatively simple isolated dendrite tip model to calculate the growth undercooling of the halo tip. (C) 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

General properties of cosmological models with an isotropic singularity

Much of the published work regarding the Isotropic Singularity is performed under the assumption that the matter source for the cosmological model is a barotropic perfect fluid, or even a perfect fluid with a gamma-law equation of state. There are, however, some general properties of cosmological models which admit an Isotropic Singularity, irrespective of the matter source. In particular, we show that the Isotropic Singularity is a point-like singularity and that vacuum space-times cannot admit an Isotropic Singularity. The relationships between the Isotropic Singularity, and the energy conditions, and the Hubble parameter is explored. A review of work by the authors, regarding the Isotropic Singularity, is presented.

Kinetics of photoacclimation in corals

Traditional models describing the relationship between photosynthesis (P) and irradiance (I) do not account for photoacclimation to short-term variation in irradiance. Here we develop and test a model that predicts the rate of photosynthesis under fluctuating irradiances at the scale of days to weeks. Using oxygen respirometry, we measured the rates of change in the P-I model parameters P-max (maximum rate of gross photosynthesis) and I-k (sub-saturation irradiance) of the photo-symbiotic coral Turbinaria mesenterina (Lamarck) following large and small increases and decreases in growth irradiance. We analyse the behaviour of the dynamic P-I model in turbid-water conditions using a dataset of 3-month continuous irradiance as the input variable. In response to upward or downward changes in experimental growth irradiance, I-k values decreased or increased exponentially, reaching new and stable levels within 5-10 days. I-k responded 4 times stronger than P-max to changes in growth irradiance. The kinetics of I-k did not show hysteresis, and changed in similar ways when irradiance was increased or decreased in small or large amounts. This suggests that mechanisms associated with photo-protection during increases in irradiance, and the maximisation of photosynthetic efficiency during decreases in irradiance, are equally potent. On the scale of months, the dynamic P-I model did not predict higher rates of photosynthesis than the static P-I model, but buffered the variation in photosynthesis during periods of reduced irradiance. Fourier analysis indicated that the kinetics of I-k closely matches the main periodicities in daily irradiance (1-2 weeks). The recorded kinetics of photoacclimation in the Turbinaria-zooxanthella symbiosis is comparable to that of free-living phytoplankton and faster than that of higher plants.

Algebraic Bethe ansatz for integrable one-dimensional extended Hubbard models with open boundary conditions

Nine classes of integrable open boundary conditions, further extending the one-dimensional U-q (gl (212)) extended Hubbard model, have been constructed previously by means of the boundary Z(2)-graded quantum inverse scattering method. The boundary systems are now solved by using the algebraic Bethe ansatz method, and the Bethe ansatz equations are obtained for all nine cases.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Measuring the decoherence rate in a semiconductor charge qubit

We describe a method by which the decoherence time of a solid-state qubit may be measured. The qubit is coded in the orbital degree of freedom of a single electron bound to a pair of donor impurities in a semiconductor host. The qubit is manipulated by adiabatically varying an external electric field. We show that by measuring the total probability of a successful qubit rotation as a function of the control field parameters, the decoherence rate may be determined. We estimate various system parameters, including the decoherence rates due to electromagnetic fluctuations and acoustic phonons. We find that, for reasonable physical parameters, the experiment is possible with existing technology. In particular, the use of adiabatic control fields implies that the experiment can be performed with control electronics with a time resolution of tens of nanoseconds.

Real exchange rate and elasticity of labour supply in a balance-of-payments-constrained macrodynamics

A macrodynamic model is proposed in which the real exchange rate and the elasticity of labour supply interact defining different trajectories of growth and income distribution in a developing economy. Growth depends on imports of capital goods which are paid with exports (there are no capital flows) and hence is constrained by equilibrium in current account. The role of the elasticity of labour supply is to prevent the real exchange rate from appreciating as the economy grows, thereby sustaining international competitiveness. The model allows for endogenous technological change and considers the impact of migration from the subsistence to the modern sector on the cumulative (Kaldor-Verdoorn) process of learning.

Optimal reserve management and sovereign debt

Most models currently used to determine optimal foreign reserve holdings take the level of international debt as given. However, given the sovereign`s willingness-to-pay incentive problems, reserve accumulation may reduce sustainable debt levels. In addition, assuming constant debt levels does not allow addressing one of the puzzles behind using reserves as a means to avoid the negative effects of crisis: why do not sovereign countries reduce their sovereign debt instead? To study the joint decision of holding sovereign debt and reserves, we construct a stochastic dynamic equilibrium model calibrated to a sample of emerging markets. We obtain that the reserve accumulation does not play a quantitatively important role in this model. In fact, we find the optimal policy is not to hold reserves at all. This finding is robust to considering interest rate shocks, sudden stops, contingent reserves and reserve dependent output costs. (c) 2008 Elsevier B.V. All rights reserved.

A new graphic format to facilitate the understanding of technological innovation models: the seesaw of competitiveness

Many models exist in the literature to explain the success of technological innovation. However, no studies have been made regarding graphic formats representing the technological innovation models and their impact, or on the understanding of these models by non-specialists in technology management. Thus, the main objective of this paper is to propose a new graphic configuration to represent the technological innovation management. Based on the literature, the innovation model is presented in the traditional format. Next, the same model is designed in the graphic format - named `the see-saw of competitiveness` - showing the interfaces among the identified factors. The two graphic formats were compared by a group of graduate students in terms of the ease in understanding the conceptual model of innovation. The statistical analysis shows that the seesaw of competitiveness is preferred.

Revisiting the interest rate puzzle

This article makes a connection between Lucas` (1978) asset pricing model and the macroeconomic dynamics for some selected countries. Both the relative risk aversion and the impatience for postponing consumption by synthesizing the investor behaviour can help to understand some key macroeconomic issues across countries, such as the savings decision and the real interest rate. I find that the government consumption makes worse the so-called `equity premium-interest rate puzzle`. The first root of the quadratic function for explaining the real interest rate can produce this puzzle, but not the second root. Thus, Mehra and Prescott (1985) identified only one possible solution.

Ecosystem process models at multiple scales for mapping tropical forest productivity

Quantifying mass and energy exchanges within tropical forests is essential for understanding their role in the global carbon budget and how they will respond to perturbations in climate. This study reviews ecosystem process models designed to predict the growth and productivity of temperate and tropical forest ecosystems. Temperate forest models were included because of the minimal number of tropical forest models. The review provides a multiscale assessment enabling potential users to select a model suited to the scale and type of information they require in tropical forests. Process models are reviewed in relation to their input and output parameters, minimum spatial and temporal units of operation, maximum spatial extent and time period of application for each organization level of modelling. Organizational levels included leaf-tree, plot-stand, regional and ecosystem levels, with model complexity decreasing as the time-step and spatial extent of model operation increases. All ecosystem models are simplified versions of reality and are typically aspatial. Remotely sensed data sets and derived products may be used to initialize, drive and validate ecosystem process models. At the simplest level, remotely sensed data are used to delimit location, extent and changes over time of vegetation communities. At a more advanced level, remotely sensed data products have been used to estimate key structural and biophysical properties associated with ecosystem processes in tropical and temperate forests. Combining ecological models and image data enables the development of carbon accounting systems that will contribute to understanding greenhouse gas budgets at biome and global scales.