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.

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.

Novel omega-conotoxins from Conus catus discriminate among neuronal calcium channel subtypes

omega -Conotoxins selective for N-type calcium channels are useful in the management of severe pain. In an attempt to expand the therapeutic potential of this class, four new omega -conotoxins (CVIA-D) have been discovered in the venom of the piscivorous cone snail, Conus catus, using assay-guided fractionation and gene cloning. Compared with other omega -conotoxins, CVID has a novel loop 4 sequence and the highest selectivity for N-type over P/Q-type calcium channels in radioligand binding assays. CVIA-D also inhibited contractions of electrically stimulated rat vas deferens. In electrophysiological studies, omega -conotoxins CVID and MVIIA had similar potencies to inhibit current through central (alpha (1B-d)) and peripheral (alpha (1B-b)) splice variants of the rat N-type calcium channels when coexpressed with rat beta (3) in Xenopus oocytes, However, the potency of CVID and MVIIA increased when alpha (1B-d) and alpha (1B-b) were expressed in the absence of rat beta (3), an effect most pronounced for CVID at alpha (1B-d) (up to 540-fold) and least pronounced for MVIIA at alpha (1B-d) (3-fold). The novel selectivity of CVID may have therapeutic implications. H-1 NMR studies reveal that CMD possesses a combination of unique structural features, including two hydrogen bonds that stabilize loop 2 and place loop 2 proximal to loop 4, creating a globular surface that is rigid and well defined.

Quantitative analysis of vascular to cardiac selectivity of L- and T-type voltage-operated calcium channel antagonists in human tissues

1. Classical L-type voltage-operated calcium channel (VOCC) antagonists dilate blood vessels, depress myocardial contractility and slow cardiac conduction. 2. We compared four L-type VOCC antagonists and a novel tetralol derivative, mibefradil, reportedly 10-fold more selective for T- (transient) over L-type VOCC in two in vitro assays of human tissue, namely isolated small arteries from the aortic vasa vasorum in a myograph and right atrial trabeculae muscle under isometric force conditions. 3. In arteries contracted with K+ (62 mmol/L), the relaxation pIC(50) values for the VOCC antagonists felodipine, nifedipine, amlodipine, verapamil and mibefradil were 8.30, 7.78, 6.64, 6.26 and 6.22, respectively. In atrial trabeculae, the pIC(50) values to inhibit the inotropic response to a submaximal concentration of isoprenaline (6 nmol/L) for felodipine, nifedipine, verapamil, amlodipine and mibefradil were 7.21, 6.95, 6.91, 5.94 and 4.61, respectively. 4. Taking the anti-log (pIC(50) vessel - pIC(50) atrium) the vascular relaxation to cardiac depression potency ratios for mibefradil, felodipine, nifedipine, amlodipine and verapamil were 41, 12, 7, 5 and 0.22, respectively. 5. We conclude that, in human tissue assays, perhaps T- over L-type VOCC selectivity confers the most favourable vascular selectivity on mibefradil. Alternatively, splice variants of L-type VOCC in the vasculature (CaV1.2b) may be more sensitive to mibefradil than the splice variants in the heart (CaV1.2a).

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.

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.

Neuronal sodium-channel alpha 1-subunit mutations in generalized epilepsy with febrile seizures plus

Generalized epilepsy with febrile seizures plus (GEFS+) is a familial epilepsy syndrome characterized by the presence of febrile and afebrile seizures. The first gene, GEFS1, was mapped to chromosome 19q and was identified as the sodium-channel beta1-subunit, SCN1B. A second locus on chromosome 2q, GEFS2, was recently identified as the sodium-channel alpha1-subunit, SCN1A. Single-stranded conformation analysis (SSCA) of SCN1A was performed in 53 unrelated index cases to estimate the frequency of mutations in patients with GEFS+. No mutations were found in 17 isolated cases of GEFS+. Three novel SCN1A mutations-D188V, V1353L, and I1656M-were found in 36 familial cases; of the remaining 33 families, 3 had mutations in SCN1B. On the basis of SSCA, the combined frequency of SCN1A and SCN1B mutations in familial cases of GEFS+ was found to be 17%.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Modelling life history strategies with capture-recapture data: Evolutionary demography of the water skink Eulamprus tympanum

Matrix population models, elasticity analysis and loop analysis can potentially provide powerful techniques for the analysis of life histories. Data from a capture-recapture study on a population of southern highland water skinks (Eulamprus tympanum) were used to construct a matrix population model. Errors in elasticities were calculated by using the parametric bootstrap technique. Elasticity and loop analyses were then conducted to identify the life history stages most important to fitness. The same techniques were used to investigate the relative importance of fast versus slow growth, and rapid versus delayed reproduction. Mature water skinks were long-lived, but there was high immature mortality. The most sensitive life history stage was the subadult stage. It is suggested that life history evolution in E. tympanum may be strongly affected by predation, particularly by birds. Because our population declined over the study, slow growth and delayed reproduction were the optimal life history strategies over this period. Although the techniques of evolutionary demography provide a powerful approach for the analysis of life histories, there are formidable logistical obstacles in gathering enough high-quality data for robust estimates of the critical parameters.

Disaggregation of polygons of surficial geology and soil maps using spatial modelling and legacy data

Examples from the Murray-Darling basin in Australia are used to illustrate different methods of disaggregation of reconnaissance-scale maps. One approach for disaggregation revolves around the de-convolution of the soil-landscape paradigm elaborated during a soil survey. The descriptions of soil ma units and block diagrams in a soil survey report detail soil-landscape relationships or soil toposequences that can be used to disaggregate map units into component landscape elements. Toposequences can be visualised on a computer by combining soil maps with digital elevation data. Expert knowledge or statistics can be used to implement the disaggregation. Use of a restructuring element and k-means clustering are illustrated. Another approach to disaggregation uses training areas to develop rules to extrapolate detailed mapping into other, larger areas where detailed mapping is unavailable. A two-level decision tree example is presented. At one level, the decision tree method is used to capture mapping rules from the training area; at another level, it is used to define the domain over which those rules can be extrapolated. (C) 2001 Elsevier Science B.V. All rights reserved.

Modelling socially optimal land allocations for sugar cane growing in North Queensland: a linked mathematical programming and choice modelling study

A modelling framework is developed to determine the joint economic and environmental net benefits of alternative land allocation strategies. Estimates of community preferences for preservation of natural land, derived from a choice modelling study, are used as input to a model of agricultural production in an optimisation framework. The trade-offs between agricultural production and environmental protection are analysed using the sugar industry of the Herbert River district of north Queensland as an example. Spatially-differentiated resource attributes and the opportunity costs of natural land determine the optimal tradeoffs between production and conservation for a range of sugar prices.