Modelling dispersal behaviour on a fractal landscape

We use a spatially explicit population model to explore the population consequences of different habitat selection mechanisms on landscapes with fractal variation in habitat quality. We consider dispersal strategies ranging from random walks to perfect habitat selectors for two species of arboreal marsupial, the greater glider (Petauroides volans) and the mountain brushtail possum (Trichosurus caninus). In this model increasing habitat selection means individuals obtain higher quality territories, but experience increased mortality during dispersal. The net effect is that population sizes are smaller when individuals actively select habitat. We find positive relationships between habitat quality and population size can occur when individuals do not use information about the entire landscape when habitat quality is spatially autocorrelated. We also find that individual behaviour can mitigate the negative effects of spatial variation on population average survival and fecundity. (C) 1998 Elsevier Science Ltd. All rights reserved.

Modelling the large deformations in stratified media - the Cosserat continuum approach

Methods employing continuum approximation in describing the deformation of layered materials possess a clear advantage over explicit models, However, the conventional implicit models based on the theory of anisotropic continua suffers from certain difficulties associated with interface slip and internal instabilities. These difficulties can be remedied by considering the bending stiffness of the layers. This implies the introduction of moment (couple) stresses and internal rotations, which leads to a Cosserat-type theory. In the present model, the behaviour of the layered material is assumed to be linearly elastic; the interfaces are assumed to be elastic perfectly plastic. Conditions of slip or no slip at the interfaces are detected by a Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformation analysis. The model is incorporated into the finite element program AFENA and validated against analytical solutions of elementary buckling problems in layered medium. A problem associated with buckling of the roof and the floor of a rectangular excavation in jointed rock mass under high horizontal in situ stresses is considered as the main application of the theory. Copyright (C) 1999 John Wiley & Sons, Ltd.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Continuity of care and readmission in two service systems: a comparative Victorian and Groningen case-register study

Objective: We compared service consumption, continuity of care and risk of readmission in a record linkage follow-up study of cohorts of patients with schizophrenia and related disorders in Victoria (Australia) and in Groningen (The Netherlands). These areas are interesting to compare because mental health care is in a different stage of deiustitutionalization. More beds are available in Groningen and more community resources are available in Victoria. Method: The cohorts were followed for 4 years, since discharge from inpatient services using record linkage data available in the psychiatric case-registers in both areas. Survival analysis was used to study continuity of care and risk of readmission. Results: Available indicators showed a higher level of continuity of care in Victoria. While the relative risk of readmission was the same in both areas and not affected by aftercare contact after discharge, the number of days spent in hospital was much higher in the Groningen register area. Conclusion: These findings provide further support for earlier reports that the risk of readmission is predominantly affected by attributes of mental illness. However, the duration of admissions, is strongly affected by service system variables, including the provision of continuity of care.

HIV infection among women admitted to the gynaecology service of a district hospital in South Africa

Our objective was to determine the prevalence of HIV infection and disease-specific HIV prevalence among women admitted to the gynaecology service of a district hospital in South Africa over a 3-month period in 1997. This was done with the goal of developing HIV education and counselling services in this setting. HIV status was determined among 196 (96%) of 205 consecutive admissions; 82 (42%) tested HIV positive. The HIV-infected women were younger than the HIV uninfected women (mean age 27 vs 33 years, P=0.001). The disease-specific HIV prevalence was greater than or equal to 40% among women who had had abortions, pelvic inflammatory disease, or ectopic pregnancy. The length of hospital stay (mean 5.4 days) and mortality (1%) were similar in the 2 groups. Inpatient gynaecology services may be act important setting in Africa, within which to provide HIV education, counselling and care.

Providing evidence-based answers to clinical questions - A pilot information service for general practitioners

Objective: To pilot a clinical information service for general practitioners. Methods: A representative sample of 31 GPs was invited to submit clinical questions to a local academic department of general practice. Their views on the service and the usefulness of the information were obtained by telephone interview. Results: Over one month, nine GPs (29% of the sample, 45% of those stating an interest), submitted 20 enquiries comprising 45 discrete clinical questions. The median time to search for evidence, appraise it and write answers to each enquiry was 2.5 hours (range, 1.0-7.4 hours). The median interval between receipt of questions and dispatch of answers was 3 clays (range, 1-12 days). Conclusions: The GPs found the answers useful in clinical decision making; in four out of 20 cases patient management was altered.

Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.

Identifying regional dust transport pathways: Application of kinematic trajectory modelling to a trans-Tasman case

PI kinematic trajectory model is used to investigate potential pathways of dust transport from Australia to New Zealand. Historically, these have been assumed to follow rather direct west-east trajectories spanning 2 to 3 days, often resulting in red snow events in the Southern Alps of New Zealand. However, results from the present study which examined the route taken by air parcels originating in southern Australia during dust storms on 24 and 25 May 1994, indicate that trans-Tasman dust transport trajectories are more diverse than previously thought, and display considerable variation during single events. These mon divergent pathways tie in more closely with aeolian dust sedimentation patterns identified by ocean coring in the Tasman Sea, and may account for the deposition of Australian dust on sub-Antarctic islands located well south of the Australian continent. Copyright 2000 John Wiley Sons, Ltd.

Preliminary experiences with a single-patient trials service in general practice

Objective: To pilot a single-patient trials (SPTs) service in general practice, designed to improve decision-making about long-term medications for chronic conditions. Design: 12-week within-patient, randomised, double-blind, placebo-controlled, crossover comparison of ibuprofen with paracetamol for osteoarthritis, involving three pairs of two-week treatment periods for each participating patient. Setting and patients: Patients attending an academic general practice with a clinical diagnosis of osteoarthritis, with pain of at least a month's duration severe enough to warrant consideration of long-term non-steroidal anti-inflammatory drug (NSAID) use. Main outcome measures: Pain and stiffness; measures of overall arthritis compared with previous fortnight; preference for NSAID at the end of each two-week treatment period; use of escape analgesia; side effects; and management changes as a result of the SPTs. Results: Eight of 14 patients completed SPTs. One was a clear responder to NSAIDs, five were non-responders, and two were indefinite. Of the five who were using NSAIDs before the SPT, two continued and three ceased using them. Clinically useful information assisted decision-making for all eight participants. Medication management changed for six. Conclusions: Single-patient trials can be successfully implemented in general practice and might be a valuable method for GPs to identify patients who respond to medication for chronic stable conditions such as osteoarthritis, in which individual response to medication is variable.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Robust Mixture Modelling Using the t Distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.