Swapping space for time and unfair tests of ecological models

Testing ecological models for management is an increasingly important part of the maturation of ecology as an applied science. Consequently, we need to work at applying fair tests of models with adequate data. We demonstrate that a recent test of a discrete time, stochastic model was biased towards falsifying the predictions. If the model was a perfect description of reality, the test falsified the predictions 84% of the time. We introduce an alternative testing procedure for stochastic models, and show that it falsifies the predictions only 5% of the time when the model is a perfect description of reality. The example is used as a point of departure to discuss some of the philosophical aspects of model testing.

PMCA1 mRNA expression in rat aortic myocytes: A real-time RT-PCR study

The plasmalemmal Ca2+ adenosine triphosphatase (PMCA) is a key regulator of Ca2+ efflux in vascular smooth muscle. In these studies are developed a realtime reverse transcriptase-polymerase chain reaction (real-time RT-PCR) assay for assessing PMCA1 mRNA levels in rat primary cultured aortic myocytes. This assay detected fetal bovine serum-induced increases in PMCA1 mRNA (relative to 18S rRNA) 4, 8, and 24 h after stimulation. Early fetal bovine serum-induced increases in PMCA1 mRNA were insensitive to the Ca2+ channel blockers nifedipine, flunarizine, and SKF-96365. These studies demonstrate the feasibility of real-time RT-PCR to assess mRNA levels of PMCA1 and illustrate dynamic regulation of this Ca2+ pump isoform in rat primary cultured aortic myocytes, (C) 2000 Academic Press.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

Heat flow for Yang-Mills-Higgs fields, part I

The Yang-Mills-Higgs field generalizes the Yang-Mills field. The authors establish the local existence and uniqueness of the weak solution to the heat flow for the Yang-Mills-Higgs field in a vector bundle over a compact Riemannian 4-manifold, and show that the weak solution is gauge-equivalent to a smooth solution and there are at most finite singularities at the maximum existing time.

Occupation, hours worked, and leisure-time physical activity

Background. International research indicates that blue-collar employees typically exhibit lower rates of leisure-time physical activity. While lack of time and work demands are commonly reported barriers to activity, the extent to which time-at-work mediates the relationship between occupation and leisure-time physical activity is unclear. This study investigated the association between occupation, time spent in paid employment, and participation in leisure-time physical activity. Methods. This was a secondary analysis of cross-sectional data from the 1995 Australian Health Survey, focusing on employed persons ages 18-64 years (n = 24,454), Occupation was coded as per the Australian Standard Classification of Occupations and collapsed into three categories (professional, white-collar, blue-collar). Hours worked was categorized into eight levels, ranging from 1-14 to more than 50 h per week. Participation in leisure-time physical activity was categorized as either insufficient or sufficient for health, consistent with recommended levels of energy expenditure (1600 METS-min/fortnight). The relationship between occupation, hours worked, and leisure-time physical activity was examined using logistic regression. Analyses were conducted separately for male and female, and the results are presented as a series of models that successively adjust for a range of potential covariates: age, living arrangement, smoking status, body mass index, and self-reported health. Results. Individuals in blue-collar occupations were approximately 50% more likely to be classified as insufficiently active. This occupational variability in leisure-time physical activity was not explained by hours worked. There was a suggested relationship between hours worked and leisure-time physical activity; however, this differed between men and women, and was difficult to interpret. Conclusions. Occupational variability in leisure-time physical activity cannot be explained by hours worked. Therefore, reports that work constitutes a barrier to participation should be explored further. Identification of the factors contributing to occupational variability in leisure-time physical activity will add to our understanding of why population subgroups differ in their health risk profiles, and assist in the development of health promotion strategies to reduce rates of sedentariness and health inequalities. (C) 2000 American Health Foundation and Academic Press.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.

Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

Residence time distributions of gas flowing through rotating drum bioreactors

Residence time distribution studies of gas through a rotating drum bioreactor for solid-state fermentation were performed using carbon monoxide as a tracer gas. The exit concentration as a function of time differed considerably from profiles expected for plug flow, plug flow with axial dispersion, and continuous stirred tank reactor (CSTR) models. The data were then fitted by least-squares analysis to mathematical models describing a central plug flow region surrounded by either one dead region (a three-parameter model) or two dead regions (a five-parameter model). Model parameters were the dispersion coefficient in the central plug flow region, the volumes of the dead regions, and the exchange rates between the different regions. The superficial velocity of the gas through the reactor has a large effect on parameter values. Increased superficial velocity tends to decrease dead region volumes, interregion transfer rates, and axial dispersion. The significant deviation from CSTR, plug flow, and plug flow with axial dispersion of the residence time distribution of gas within small-scale reactors can lead to underestimation of the calculation of mass and heat transfer coefficients and hence has implications for reactor design and scaleup. (C) 2001 John Wiley & Sons, Inc.

Autonomic activity during human sleep as a function of time and sleep stage

While there is a developing understanding of the influence of sleep on cardiovascular autonomic activity in humans, there remain unresolved issues. In particular, the effect of time within the sleep period, independent of sleep stage, has not been investigated. Further, the influence of sleep on central sympathetic nervous system (SNS) activity is uncertain because results using the major method applicable to humans, the low frequency (LF) component of heart rate Variability (HRV), have been contradictory, and because the method itself is open to criticism. Sleep and cardiac activity were measured in 14 young healthy subjects on three nights. Data was analysed in 2-min epochs. All epochs meeting specified criteria were identified, beginning 2 h before, until 7 h after, sleep onset. Epoch values were allocated to 30-min bins and during sleep were also classified into stage 2, slow wave sleep (SWS) and rapid eye movement (REM) sleep. The measures of cardiac activity were heart irate (HR), blood pressure (BP), high frequency (HF) and LF components of HRV and pre-ejection period (PEP). During non-rapid eye movement (NREM) sleep autonomic balance shifted from sympathetic to parasympathetic dominance, although this appeared to be more because of a shift in parasympathetic nervous system (PNS) activity. Autonomic balance during REM was in general similar to wakefulness. For BP and the HF and LF components the change occurred abruptly at sleep onset and was then constant over time within each stage of sleep, indicating that any change in autonomic balance over the sleep period is a consequence of the changing distribution of sleep stages. Two variables, HR and PEP, did show time effects reflecting a circadian influence over HR and perhaps time asleep affecting PEP. While both the LF component and PEP showed changes consistent with reduced sympathetic tone during sleep, their pattern of change over time differed.

A time-harmonic inverse methodology for the design of RF coils in MRI

An inverse methodology is described to assist in the design of radio-frequency (RF) coils for magnetic resonance imaging (MRI) applications. The time-harmonic electromagnetic Green's functions are used to calculate current on the coil and shield cylinders that will generate a specified internal magnetic field. Stream function techniques and the method of moments are then used to implement this theoretical current density into an RF coil. A novel asymmetric coil operating for a 4.5 T MRI machine was designed and constructed using this methodology and the results are presented.

Mathematical models in percutaneous absorption

A number of mathematical models have been used to describe percutaneous absorption kinetics. In general, most of these models have used either diffusion-based or compartmental equations. The object of any mathematical model is to a) be able to represent the processes associated with absorption accurately, b) be able to describe/summarize experimental data with parametric equations or moments, and c) predict kinetics under varying conditions. However, in describing the processes involved, some developed models often suffer from being of too complex a form to be practically useful. In this chapter, we attempt to approach the issue of mathematical modeling in percutaneous absorption from four perspectives. These are to a) describe simple practical models, b) provide an overview of the more complex models, c) summarize some of the more important/useful models used to date, and d) examine sonic practical applications of the models. The range of processes involved in percutaneous absorption and considered in developing the mathematical models in this chapter is shown in Fig. 1. We initially address in vitro skin diffusion models and consider a) constant donor concentration and receptor conditions, b) the corresponding flux, donor, skin, and receptor amount-time profiles for solutions, and c) amount- and flux-time profiles when the donor phase is removed. More complex issues, such as finite-volume donor phase, finite-volume receptor phase, the presence of an efflux. rate constant at the membrane-receptor interphase, and two-layer diffusion, are then considered. We then look at specific models and issues concerned with a) release from topical products, b) use of compartmental models as alternatives to diffusion models, c) concentration-dependent absorption, d) modeling of skin metabolism, e) role of solute-skin-vehicle interactions, f) effects of vehicle loss, a) shunt transport, and h) in vivo diffusion, compartmental, physiological, and deconvolution models. We conclude by examining topics such as a) deep tissue penetration, b) pharmacodynamics, c) iontophoresis, d) sonophoresis, and e) pitfalls in modeling.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.