Improving wheat simulation capabilities in Australia from a cropping systems perspective - III. The integrated wheat model (I_WHEAT)

Previous work has identified several short-comings in the ability of four spring wheat and one barley model to simulate crop processes and resource utilization. This can have important implications when such models are used within systems models where final soil water and nitrogen conditions of one crop define the starting conditions of the following crop. In an attempt to overcome these limitations and to reconcile a range of modelling approaches, existing model components that worked demonstrably well were combined with new components for aspects where existing capabilities were inadequate. This resulted in the Integrated Wheat Model (I_WHEAT), which was developed as a module of the cropping systems model APSIM. To increase predictive capability of the model, process detail was reduced, where possible, by replacing groups of processes with conservative, biologically meaningful parameters. I_WHEAT does not contain a soil water or soil nitrogen balance. These are present as other modules of APSIM. In I_WHEAT, yield is simulated using a linear increase in harvest index whereby nitrogen or water limitations can lead to early termination of grainfilling and hence cessation of harvest index increase. Dry matter increase is calculated either from the amount of intercepted radiation and radiation conversion efficiency or from the amount of water transpired and transpiration efficiency, depending on the most limiting resource. Leaf area and tiller formation are calculated from thermal time and a cultivar specific phyllochron interval. Nitrogen limitation first reduces leaf area and then affects radiation conversion efficiency as it becomes more severe. Water or nitrogen limitations result in reduced leaf expansion, accelerated leaf senescence or tiller death. This reduces the radiation load on the crop canopy (i.e. demand for water) and can make nitrogen available for translocation to other organs. Sensitive feedbacks between light interception and dry matter accumulation are avoided by having environmental effects acting directly on leaf area development, rather than via biomass production. This makes the model more stable across environments without losing the interactions between the different external influences. When comparing model output with models tested previously using data from a wide range of agro-climatic conditions, yield and biomass predictions were equal to the best of those models, but improvements could be demonstrated for simulating leaf area dynamics in response to water and nitrogen supply, kernel nitrogen content, and total water and nitrogen use. I_WHEAT does not require calibration for any of the environments tested. Further model improvement should concentrate on improving phenology simulations, a more thorough derivation of coefficients to describe leaf area development and a better quantification of some processes related to nitrogen dynamics. (C) 1998 Elsevier Science B.V.

Capillary network models for transport in packed beds: Considerations of pore aspect ratio

The conventional analysis for the estimation of the tortuosity factor for transport in porous media is modified here to account for the effect of pore aspect ratio. Structural models of the porous medium are also constructed for calculating the aspect ratio as a function of porosity. Comparison of the model predictions with the extensive data of Currie (1960) for the effective diffusivity of hydrogen in packed beds shows good agreement with a network model of randomly oriented intersecting pores for porosities upto about 50 percent, which is the region of practical interest. The predictions based on this network model are also found to be in better agreement with the data of Currie than earlier expressions developed for unconsolidated and grainy media.

Two dimensional computational fluid dynamic models for waste stabilisation ponds

Traditional waste stabilisation pond (WSP) models encounter problems predicting pond performance because they cannot account for the influence of pond features, such as inlet structure or pond geometry, on fluid hydrodynamics. In this study, two dimensional (2-D) computational fluid dynamics (CFD) models were compared to experimental residence time distributions (RTD) from literature. In one of the-three geometries simulated, the 2-D CFD model successfully predicted the experimental RTD. However, flow patterns in the other two geometries were not well described due to the difficulty of representing the three dimensional (3-D) experimental inlet in the 2-D CFD model, and the sensitivity of the model results to the assumptions used to characterise the inlet. Neither a velocity similarity nor geometric similarity approach to inlet representation in 2-D gave results correlating with experimental data. However. it was shown that 2-D CFD models were not affected by changes in values of model parameters which are difficult to predict, particularly the turbulent inlet conditions. This work suggests that 2-D CFD models cannot be used a priori to give an adequate description of the hydrodynamic patterns in WSP. (C) 1998 Elsevier Science Ltd. All rights reserved.

Metabolite mean transit times in the liver as predicted by various models of hepatic elimination

Predicted area under curve (AUC), mean transit time (MTT) and normalized variance (CV2) data have been compared for parent compound and generated metabolite following an impulse input into the liver, Models studied were the well-stirred (tank) model, tube model, a distributed tube model, dispersion model (Danckwerts and mixed boundary conditions) and tanks-in-series model. It is well known that discrimination between models for a parent solute is greatest when the parent solute is highly extracted by the liver. With the metabolite, greatest model differences for MTT and CV2 occur when parent solute is poorly extracted. In all cases the predictions of the distributed tube, dispersion, and tasks-in-series models are between the predictions of the rank and tube models. The dispersion model with mixed boundary conditions yields identical predictions to those for the distributed tube model (assuming an inverse gaussian distribution of tube transit times). The dispersion model with Danckwerts boundary conditions and the tanks-in series models give similar predictions to the dispersion (mixed boundary conditions) and the distributed tube. The normalized variance for parent compound is dependent upon hepatocyte permeability only within a distinct range of permeability values. This range is similar for each model but the order of magnitude predicted for normalized variance is model dependent. Only for a one-compartment system is the MIT for generated metabolite equal to the sum of MTTs for the parent compound and preformed metabolite administered as parent.

Inhibition of phenotype modulation, growth, and migration of vascular smooth muscle cells by a guanosine-rich 30-mer phosphorothioate oligodeoxynucleotide

The testing of a 30-mer dG-rich phosphorothioate oligodeoxynucleotide (LG4PS) for effects on the behaviour of vascular smooth muscle cells (VSMC) in vitro and in vivo is described. LG4PS at 0.3 mu M inhibited significantly the phenotype modulation of freshly isolated rabbit VSMC, and cell outgrowth from pig aortic explants was inhibited similar to 80% by 5 mu M LG4PS. The growth of proliferating rabbit and pig VSMC was inhibited similar to 70% by 0.3 mu M and 5 mu M LG4PS, respectively. Though less marked, the antiproliferative effects of LG4PS on human VSMC were comparable to those obtained with heparin. The cytotoxic effects of LG4PS on VSMC in vitro were low. Despite these promising results, adventitial application of 2-200 nmol LG4PS in pluronic gel failed to reduce vascular hyperplasia in balloon-injured rabbit carotid arteries, and the highest dose caused extensive mortality. (C) 1997 Academic Press Limited.

Modelling fungal growth on surfaces

Fungal growth in time and space at the substrate surface was modelled for a simple system mimicking solid-state fermentation, using a polycarbonate Nucleopore membrane laid over a glucose solution. Biomass production depends on both tip density and the diffusion of glucose within the fungal hyphae. The model predicts early increases in both height and concentration, followed by a period in which the biomass profile moves with a constant wavefront. The rate of increase in height increases as tip diffusivity increases or as the Monod saturation constant for glucose decreases.

Modeling ex vivo hematopoiesis using chemical engineering metaphors

Ex vivo hematopoiesis is increasingly used for clinical applications. Models of ex vivo hematopoiesis are required to better understand the complex dynamics and to optimize hematopoietic culture processes. A general mathematical modeling framework is developed which uses traditional chemical engineering metaphors to describe the complex hematopoietic dynamics. Tanks and tubular reactors are used to describe the (pseudo-) stochastic and deterministic elements of hematopoiesis, respectively. Cells at any point in the differentiation process can belong to either an immobilized, inert phase (quiescent cells) or a mobile, active phase (cycling cells). The model describes five processes: (1) flow (differentiation), (2) autocatalytic formation (growth),(3) degradation (death), (4) phase transition from immobilized to mobile phase (quiescent to cycling transition), and (5) phase transition from mobile to immobilized phase (cycling to quiescent transition). The modeling framework is illustrated with an example concerning the effect of TGF-beta 1 on erythropoiesis. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

Anthocyanin synthesis, growth and nutrient uptake in suspension cultures of strawberry cells

Strawberry (Fragaria ananassa cv. Shikinari) cell suspension cultures carried out in shake flasks for 18 d were closely examined for cell growth, anthocyanin synthesis and the development of pigmented cells in relation to the uptake of carbohydrate, extracellular PO4, NO3, NH4, and calcium. Cell viability, extracellular anthocyanin content, pH and electrical conductivity of the broth were also monitored. The specific growth rate of strawberry cells at exponential phase was 0.27 and 0.28 d(-1) based on fresh and dry weight, respectively. Anthocyanin synthesis was observed to increase continuously to a maximum value of 0.86 mg/g fresh cell weight (FCW) at day 6, and was partially growth-associated. Anthocyanin synthesis was linearly related to the increase in pigmented cell ratio, which increased with time and reached a maximum value of ca. 70% at day 6 due to reduction in cell viability and depletion of substrate. Total carbohydrate uptake was closely associated with increase in cell growth, and glucose was utilized in preference to fructose. Nitrate and ammonia were consumed until 9 d of culture, but phosphate was completely absorbed within 4 d. Calcium was assimilated throughout the growth cycle. After 9 d, cell lysis was observed which resulted in the leakage of intracellular substances and a concomitant pH rise. Anthocyanin was never detected in the broth although the broth became darkly pigmented during the lysis period. This suggests that anthocyanin was synthesized only by viable pigmented cells, and degraded rapidly upon cell death and lysis. Based on the results of kinetic analysis, a model was developed by incorporating governing equations for the ratio of pigmented cells into a Bailey and Nicholson's model. This was verified by comparison with the experimental data. The results suggest Bat the model satisfactorily describes the strawberry cell culture process, and may thus be used for process optimization.

Economic incentives for farmers in the Philippine uplands to adopt hedgerow intercropping

Land degradation in the Philippine uplands is severe and widespread. Most upland areas are steep, and intense rainfall on soils disturbed by intensive agriculture can produce high rates of soil loss. This has serious implications for the economic welfare of a growing upland population with few feasible livelihood alternatives. Hedgerow intercropping can greatly reduce soil loss from annual cropping systems and has been considered an appropriate technology for soil conservation research and extension in the Philippine uplands. However; adoption of hedgerow intercropping has been sporadic and transient, rarely continuing once external support has been withdrawn. The objective of this paper is to investigate the economic incentives for farmers in the Philippine uplands to adopt hedgerow intercropping relative to traditional open-field maize farming. Cost-benefit analysis is used to compare the economic viability of hedgerow intercropping, as it has been promoted to upland farmers, with the viability of traditional methods of open-field farming. The APSIM and SCUAF models were used to predict the effect of soil erosion on maize yields from open-field farming and hedgerow intercropping. The results indicate that there have been strong economic incentives for farmers with limited planning horizons to reject hedgerow intercropping because the benefits of sustained yields are not realized rapidly enough to compensate for high establishment costs. Alternative forms of hedgerow intercropping such as natural vegetation and grass strips reduce establishment and maintenance costs and are therefore more economically attractive to farmers than hedgerow intercropping with shrub legumes. The long-term economic viability of hedgerow intercropping depends on the economic setting and the potential for hedgerow intercropping to sustain maize production relative to traditional open-field farming. (C) 1998 Academic Press.

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Enhanced downregulation of the p75 nerve growth factor receptor by cholesteryl and bis-cholesteryl antisense oligonucleotides

The effects of conjugating cholesterol to either or both ends of a phosphorothioate (PS) oligonucleotide were analyzed in terms of cellular uptake and antisense efficacy. The oligo sequence was directed against the p75 nerve growth factor receptor (p75), and was tested in differentiated PC12 cells, which express high levels of this protein. The addition of a single cholesteryl group to the 5'-end significantly increased cellular uptake and improved p75 mRNA downregulation compared with the unmodified PS oligo, However, only a minor degree of downregulation of p75 protein was obtained with 5' cholesteryl oligos, Three different linkers was used to attach the 5' cholesteryl group but were found not to have any impact on efficacy. Addition of a single cholesteryl group to the 3'-end led to greater p75 mRNA downregulation (31%) and p75 protein downregulation (28%) than occurred with the 5' cholesteryl oligos. The biggest improvement in antisense efficacy, both at the mRNA and protein levels, was obtained from the conjugation of cholesterol to both ends of the oligo. One of the bis-cholesteryl oligos was nearly as effective as cycloheximide at decreasing synthesis of p75, The bis-cholesteryl oligos also displayed significant efficacy at 1 mu M, whereas the other oligos required 5 mu M to be effective. The enhanced efficacy of bis-cholesteryl oligos is likely to be due to a combination of enhanced cellular uptake and resistance to both 5' and 3' exonucleases.

Fibroblast growth factor receptor expression in outer hair cells of rat cochlea

FIBROBLAST growth factors (FGFs) are critical for normal development of the organ of Corti, and may also protect hair cells from ototoxic damage. Four different fibroblast growth factors are known, three of which have different splice variants in the extracellular immunoglobin-like (Ig) III FGF-binding domain, giving different patterns of sensitivity to the different FGFs. Analysis of a cDNA library of rat outer hair cells by the polymerase chain reaction, using isoform specific primers, showed expression only of FGF receptor 3, splice variant IIIc. This allows us to predict the pattern of sensitivity to applied FGFs, may be useful in targeting outer hair cells selectively during an FGF-based strategy for cochlear therapy. (C) 1998 Lippincott Williams & Wilkins.