Feasibility of agent-based modelling of articular cartilage including a conceptual representation of its structure

Articular cartilage is a complex structure with an architecture in which fluid-swollen proteoglycans constrained within a 3D network of collagen fibrils. Because of the complexity of the cartilage structure, the relationship between its mechanical behaviours at the macroscale level and its components at the micro-scale level are not completely understood. The research objective in this thesis is to create a new model of articular cartilage that can be used to simulate and obtain insight into the micro-macro-interaction and mechanisms underlying its mechanical responses during physiological function. The new model of articular cartilage has two characteristics, namely: i) not use fibre-reinforced composite material idealization ii) Provide a framework for that it does probing the micro mechanism of the fluid-solid interaction underlying the deformation of articular cartilage using simple rules of repartition instead of constitutive / physical laws and intuitive curve-fitting. Even though there are various microstructural and mechanical behaviours that can be studied, the scope of this thesis is limited to osmotic pressure formation and distribution and their influence on cartilage fluid diffusion and percolation, which in turn governs the deformation of the compression-loaded tissue. The study can be divided into two stages. In the first stage, the distributions and concentrations of proteoglycans, collagen and water were investigated using histological protocols. Based on this, the structure of cartilage was conceptualised as microscopic osmotic units that consist of these constituents that were distributed according to histological results. These units were repeated three-dimensionally to form the structural model of articular cartilage. In the second stage, cellular automata were incorporated into the resulting matrix (lattice) to simulate the osmotic pressure of the fluid and the movement of water within and out of the matrix; following the osmotic pressure gradient in accordance with the chosen rule of repartition of the pressure. The outcome of this study is the new model of articular cartilage that can be used to simulate and study the micromechanical behaviours of cartilage under different conditions of health and loading. These behaviours are illuminated at the microscale level using the socalled neighbourhood rules developed in the thesis in accordance with the typical requirements of cellular automata modelling. Using these rules and relevant Boundary Conditions to simulate pressure distribution and related fluid motion produced significant results that provided the following insight into the relationships between osmotic pressure gradient and associated fluid micromovement, and the deformation of the matrix. For example, it could be concluded that: 1. It is possible to model articular cartilage with the agent-based model of cellular automata and the Margolus neighbourhood rule. 2. The concept of 3D inter connected osmotic units is a viable structural model for the extracellular matrix of articular cartilage. 3. Different rules of osmotic pressure advection lead to different patterns of deformation in the cartilage matrix, enabling an insight into how this micromechanism influences macromechanical deformation. 4. When features such as transition coefficient were changed, permeability (representing change) is altered due to the change in concentrations of collagen, proteoglycans (i.e. degenerative conditions), the deformation process is impacted. 5. The boundary conditions also influence the relationship between osmotic pressure gradient and fluid movement at the micro-scale level. The outcomes are important to cartilage research since we can use these to study the microscale damage in the cartilage matrix. From this, we are able to monitor related diseases and their progression leading to potential insight into drug-cartilage interaction for treatment. This innovative model is an incremental progress on attempts at creating further computational modelling approaches to cartilage research and other fluid-saturated tissues and material systems.

A strong quantitative trait locus for wing length on chromosome 2 in a wild population of great reed warblers

Wing length is a key character for essential behaviours related to bird flight such as migration and foraging. In the present study, we initiate the search for the genes underlying wing length in birds by studying a long-distance migrant, the great reed warbler (Acrocephalus arundinaceus). In this species wing length is an evolutionary interesting trait with pronounced latitudinal gradient and sex-specific selection regimes in local populations. We performed a quantitative trait locus (QTL) scan for wing length in great reed warblers using phenotypic, genotypic, pedigree and linkage map data from our long-term study population in Sweden. We applied the linkage analysis mapping method implemented in GRIDQTL (a new web-based software) and detected a genome-wide significant QTL for wing length on chromosome 2, to our knowledge, the first detected QTL in wild birds. The QTL extended over 25 cM and accounted for a substantial part (37%) of the phenotypic variance of the trait. A genome scan for tarsus length (a bodysize-related trait) did not show any signal, implying that the wing-length QTL on chromosome 2 was not associated with body size. Our results provide a first important step into understanding the genetic architecture of avian wing length, and give opportunities to study the evolutionary dynamics of wing length at the locus level. This journal is© 2010 The Royal Society.

Evaluating the role of deep granitic rocks in generating anomalous temperatures in SW-Queensland

Across central Australia and south-west Queensland, a large (~800,000km2) subsurface temperature anomaly occurs (Figure 1). Temperatures are interpreted to be greater than 235°C at 5km depth, ca. 85°C higher than the average geothermal gradient for the upper continental crust (Chopra & Holgate, 2005; Holgate & Gerner, 2011). This anomaly has driven the development of Engineered Geothermal Systems (EGS) at Innamincka, where high temperatures have been related to the radiogenic heat production of High Heat Producing Granites (HHPG) at depth, below thermally insulative sedimentary cover (Chopra & Holgate, 2005; Draper & D’Arcy, 2006; Meixner & Holgate, 2009). To evaluate the role of granitic rocks at depth in generating the broader temperature anomaly in SW-Queensland, we sampled 25 granitic rocks from basement intervals of petroleum drill cores below thermal insulative cover along two transects (WNW–ESE and NNE–SSW — Figure 1) and performed a multidisciplinary study involving petrography, whole-rock chemistry, zircon dating and thermal conductivity measurements.

Exponential integrators and a dual-scale model for wood drying

For the timber industry, the ability to simulate the drying of wood is invaluable for manufacturing high quality wood products. Mathematically, however, modelling the drying of a wet porous material, such as wood, is a diffcult task due to its heterogeneous and anisotropic nature, and the complex geometry of the underlying pore structure. The well{ developed macroscopic modelling approach involves writing down classical conservation equations at a length scale where physical quantities (e.g., porosity) can be interpreted as averaged values over a small volume (typically containing hundreds or thousands of pores). This averaging procedure produces balance equations that resemble those of a continuum with the exception that effective coeffcients appear in their deffnitions. Exponential integrators are numerical schemes for initial value problems involving a system of ordinary differential equations. These methods differ from popular Newton{Krylov implicit methods (i.e., those based on the backward differentiation formulae (BDF)) in that they do not require the solution of a system of nonlinear equations at each time step but rather they require computation of matrix{vector products involving the exponential of the Jacobian matrix. Although originally appearing in the 1960s, exponential integrators have recently experienced a resurgence in interest due to a greater undertaking of research in Krylov subspace methods for matrix function approximation. One of the simplest examples of an exponential integrator is the exponential Euler method (EEM), which requires, at each time step, approximation of φ(A)b, where φ(z) = (ez - 1)/z, A E Rnxn and b E Rn. For drying in porous media, the most comprehensive macroscopic formulation is TransPore [Perre and Turner, Chem. Eng. J., 86: 117-131, 2002], which features three coupled, nonlinear partial differential equations. The focus of the first part of this thesis is the use of the exponential Euler method (EEM) for performing the time integration of the macroscopic set of equations featured in TransPore. In particular, a new variable{ stepsize algorithm for EEM is presented within a Krylov subspace framework, which allows control of the error during the integration process. The performance of the new algorithm highlights the great potential of exponential integrators not only for drying applications but across all disciplines of transport phenomena. For example, when applied to well{ known benchmark problems involving single{phase liquid ow in heterogeneous soils, the proposed algorithm requires half the number of function evaluations than that required for an equivalent (sophisticated) Newton{Krylov BDF implementation. Furthermore for all drying configurations tested, the new algorithm always produces, in less computational time, a solution of higher accuracy than the existing backward Euler module featured in TransPore. Some new results relating to Krylov subspace approximation of '(A)b are also developed in this thesis. Most notably, an alternative derivation of the approximation error estimate of Hochbruck, Lubich and Selhofer [SIAM J. Sci. Comput., 19(5): 1552{1574, 1998] is provided, which reveals why it performs well in the error control procedure. Two of the main drawbacks of the macroscopic approach outlined above include the effective coefficients must be supplied to the model, and it fails for some drying configurations, where typical dual{scale mechanisms occur. In the second part of this thesis, a new dual{scale approach for simulating wood drying is proposed that couples the porous medium (macroscale) with the underlying pore structure (microscale). The proposed model is applied to the convective drying of softwood at low temperatures and is valid in the so{called hygroscopic range, where hygroscopically held liquid water is present in the solid phase and water exits only as vapour in the pores. Coupling between scales is achieved by imposing the macroscopic gradient on the microscopic field using suitably defined periodic boundary conditions, which allows the macroscopic ux to be defined as an average of the microscopic ux over the unit cell. This formulation provides a first step for moving from the macroscopic formulation featured in TransPore to a comprehensive dual{scale formulation capable of addressing any drying configuration. Simulation results reported for a sample of spruce highlight the potential and flexibility of the new dual{scale approach. In particular, for a given unit cell configuration it is not necessary to supply the effective coefficients prior to each simulation.

Formation of single-walled carbon nanotube via the interaction of graphene nanoribbons : ab initio density functional calculations

The interaction of bare graphene nanoribbons (GNRs) was investigated by ab initio density functional theory calculations with both the local density approximation (LDA) and the generalized gradient approximation (GGA). Remarkably, two bare 8-GNRs with zigzag-shaped edges are predicted to form an (8, 8) armchair single-wall carbon nanotube (SWCNT) without any obvious activation barrier. The formation of a (10, 0) zigzag SWCNT from two bare 10-GNRs with armchair-shaped edges has activation barriers of 0.23 and 0.61 eV for using the LDA and the revised PBE exchange correlation functional, respectively, Our results suggest a possible route to control the growth of specific types SWCNT via the interaction of GNRs.

The Glenn A. Fry award lecture 2011 : peripheral optics of the human eye

There has been a low level of interest in peripheral aberrations and corresponding image quality for over 200 years. Most work has been concerned with the second-order aberrations of defocus and astigmatism that can be corrected with conventional lenses. Studies have found high levels of aberration, often amounting to several dioptres, even in eyes with only small central defocus and astigmatism. My investigations have contributed to understanding shape changes in the eye with increases in myopia, changes in eye optics with ageing, and how surgical interventions intended to correct central refractive errors have unintended effects on peripheral optics. My research group has measured peripheral second- and higher-order aberrations over a 42° horizontal × 32° vertical diameter visual field. There is substantial variation in individual aberrations with age and pathology. While the higher-order aberrations in the periphery are usually small compared with second-order aberrations, they can be substantial and change considerably after refractive surgery. The thrust of my research in the next few years is to understand more about the peripheral aberrations of the human eye, to measure visual performance in the periphery and determine whether this can be improved by adaptive optics correction, to use measurements of peripheral aberrations to learn more about the optics of the eye and in particular the gradient index structure of the lens, and to investigate ways of increasing the size of the field of good retinal image quality.

Gradients of refractive index in the crystalline lens and transient changes in refraction among patients with diabetes

Transient hyperopic refractive shifts occur on a timescale of weeks in some patients after initiation of therapy for hyperglycemia, and are usually followed by recovery to the original refraction. Possible lenticular origin of these changes is considered in terms of a paraxial gradient index model. Assuming that the lens thickness and curvatures remain unchanged, as observed in practice, it appears possible to account for initial hyperopic refractive shifts of up to a few diopters by reduction in refractive index near the lens center and alteration in the rate of change between center and surface, so that most of the index change occurs closer to the lens surface. Restoration of the original refraction depends on further change in the refractive index distribution with more gradual changes in refractive index from the lens center to its surface. Modeling limitations are discussed.

Plasma cytokine changes in relation to exercise intensity and muscle damage

The purpose of this study was to compare the effects of exercise intensity and exercise-induced muscle damage on changes in anti-inflammatory cytokines and other inflammatory mediators. Nine well-trained male runners completed three different exercise trials on separate occasions: (1) level treadmill running at 60% VO2max (moderate-intensity trial) for 60 min; (2) level treadmill running at 85% VO2max (high-intensity trial) for 60 min; (3) downhill treadmill running (-10% gradient) at 60% VO2max (downhill running trial) for 45 min. Blood was sampled before, immediately after and 1 h after exercise. Plasma was analyzed for interleukin-1 receptor antagonist (IL-1ra), IL-4, IL-5, IL-10, IL-12p40, IL-13, monocyte chemotactic protein-1 (MCP-1), prostaglandin E(2), leukotriene B(4) and heat shock protein 70 (HSP70). The plasma concentrations of IL-1ra, IL-12p40, MCP-1 and HSP70 increased significantly (P<0.05) after all three trials. Plasma prostaglandin E(2) concentration increased significantly after the downhill running and high-intensity trials, while plasma IL-10 concentration increased significantly only after the high-intensity trial. IL-4 and leukotriene B(4) did not increase significantly after exercise. Plasma IL-1ra and IL-10 concentrations were significantly higher (P<0.05) after the high-intensity trial than after both the moderate-intensity and downhill running trials. Therefore, following exercise up to 1 h duration, exercise intensity appears to have a greater effect on anti-inflammatory cytokine production than exercise-induced muscle damage

Exercise-induced muscle damage, plasma cytokines and markers of neutrophil activation

Introduction: Unaccustomed eccentric exercise often results in muscle damage and neutrophil activation. We examined changes in plasma cytokines stress hormones, creatine kinase activity and myoglobin concentration, neutrophil surface receptor expression, degranulation, and the capacity of neutrophils to generate reactive oxygen species in response to in vitro stimulation after downhill running. Methods: Ten well-trained male runners ran downhill on a treadmill at a gradient of -10% for 45 min at 60% V̇O2max. Blood was sampled immediately before (PRE) and after (POST), 1 h (1 h POST), and 24 h (24 h POST) after exercise. Results: At POST, there were significant increases (P < 0.01) in neutrophil count (32%), plasma interleukin (IL)-6 concentration (460%), myoglobin (Mb) concentration (1100%), and creatine kinase (CK) activity (40%). At 1 h POST, there were further increases above preexercise values for neutrophil count (85%), plasma Mb levels (1800%), and CK activity (56%), and plasma IL-6 concentration remained above preexercise values (410%) (P < 0.01). At 24 h POST, neutrophil counts and plasma IL-6 levels had returned to baseline, whereas plasma Mb concentration (100%) and CK activity (420%) were elevated above preexercise values (P < 0.01). There were no significant changes in neutrophil receptor expression, degranulation and respiratory burst activity, and plasma IL-8 and granulocyte-colony stimulating factor concentrations at any time after exercise. Neutrophil count correlated with plasma Mb concentration at POST (r = 0.64, P < 0.05), and with plasma CK activity at POST (r = 0.83, P < 0.01) and 1 h POST (r = 0.78, P < 0.01). Conclusion: Neutrophil activation remains unchanged after downhill running in well-trained runners, despite increases in plasma markers of muscle damage.

Metabolic and ecological study of environmental pentose utilizing bacteria (E-PUB)

Lignocellulosic materials, such as sugar cane bagasse, a waste product of the sugarcane processing industry, agricultural residues and herbaceous crops, may serve as an abundant and comparatively cheap feedstock for largescale industrial fermentation, resulting in the production of marketable end-products. However, the complex structure of lignocellulosic materials, the presence of various hexose and pentose sugars in the hemicellulose component, and the presence of various compounds that inhibit the organisms selected for the fermentation process, all constitute barriers that add to the production costs and make full scale industrial production economically less feasible. The work presented in this thesis was conducted in order to screen microorganisms for ability to utilize pentose sugars derived from the sugar mill industrial waste. A large number of individual bacterial strains were investigated from hemi-cellulose rich material collected at the Proserpine and Maryborough sugar mills, notably soil samples from the mill sites. The research conducted to isolation of six pentose-capable Gram-positive organisms from the actinomycetes group by using pentose as a sole carbon source in the cultivation process. The isolates were identified as Corynebacterium glutamicum, Actinomyces odontolyticus, Nocardia elegans, and Propionibacterium freudenreichii all of which were isolated from the hemicellulose-enriched soil. Pentose degrading microbes are very rare in the environment, so this was a significant discovery. Previous research indicated that microbes could degrade pentose after genetic modification but the microbes discovered in this research were able to naturally utilize pentose. Six isolates, identified as four different genera, were investigated for their ability to utilize single sugars as substrates (glucose, xylose, arabinose or ribose), and also dual sugars as substrates (a hexose plus a pentose). The results demonstrated that C. glutamicum, A. odontolyticus, N. elegans, and P. freudenreichii were pentose-capable (able to grow using xylose or other pentose sugar), and also showed diauxie growth characteristics during the dual-sugar (glucose, in combination with xylose, arabinose or ribose) carbon source tests. In addition, it was shown that the isolates displayed very small differences in growth rates when grown on dual sugars as compared to single sugars, whether pentose or hexose in nature. The anabolic characteristics of C. glutamicum, A. odontolyticus, N. elegans and P. freudenreichii were subsequently investigated by qualitative analysis of their end-products, using high performance liquid chromatography (HPLC). All of the organisms produced arginine and cysteine after utilization of the pentose substrates alone. In addition, P. freudenreichii produced alanine and glycine. The end-product profile arising from culture with dual carbon sources was also tested. Interestingly, this time the product was different. All of them produced the amino acid glycine, when grown on a combination substrate-mix of glucose with xylose, and also glucose with arabinose. Only N. elegans was able to break down ribose, either singly or in combination with glucose, and the end-product of metabolism of the glucose plus ribose substrate combination was glutamic acid. The ecological analysis of microbial abundance in sugar mill waste was performed using denaturing gradient gel electrophoresis (DGGE) and also the metagenomic microarray PhyloChip method. Eleven solid samples and seven liquid samples were investigated. A very complex bacterial ecosystem was demonstrated in the seven liquid samples after testing with the PhyloChip method. It was also shown that bagasse leachate was the most different, compared to all of the other samples, by virtue of its richness in variety of taxa and the complexity of its bacterial community. The bacterial community in solid samples from Proserpine, Mackay and Maryborough sugar mills showed huge diversity. The information found from 16S rDNA sequencing results was that the bacterial genera Brevibacillus, Rhodospirillaceae, Bacillus, Vibrio and Pseudomonas were present in greatest abundance. In addition, Corynebacterium was also found in the soil samples. The metagenomic studies of the sugar mill samples demonstrate two important outcomes: firstly that the bagasse leachate, as potentially the most pentose-rich sample tested, had the most complex and diverse bacterial community; and secondly that the pentose-capable isolates that were initially discovered at the beginning of this study, were not amongst the most abundant taxonomic groups discovered in the sugar mill samples, and in fact were, as suspected, very rare. As a bioprospecting exercise, therefore, the study has discovered organisms that are naturally present, but in very small numbers, in the appropriate natural environment. This has implications for the industrial application of E-PUB, in that a seeding process using a starter culture will be necessary for industrial purposes, rather than simply assuming that natural fermentation might occur.

Synthesis and characterisation of titanium sol-gels in varied gravity environments

Sol-gel synthesis in varied gravity is only a relatively new topic in the literature and further investigation is required to explore its full potential as a method to synthesise novel materials. Although trialled for systems such as silica, the specific application of varied gravity synthesis to other sol-gel systems such as titanium has not previously been undertaken. Current literature methods for the synthesis of sol-gel material in reduced gravity could not be applied to titanium sol-gel processing, thus a new strategy had to be developed in this study. To successfully conduct experiments in varied gravity a refined titanium sol-gel chemical precursor had to be developed which allowed the single solution precursor to remain un-reactive at temperatures up to 50oC and only begin to react when exposed to a pressure decrease from a vacuum. Due to the new nature of this precursor, a thorough characterisation of the reaction precursors was subsequently undertaken with the use of techniques such as Nuclear Magnetic Resonance, Infra-red and UV-Vis spectroscopy in order to achieve sufficient understanding of precursor chemistry and kinetic stability. This understanding was then used to propose gelation reaction mechanisms under varied gravity conditions. Two unique reactor systems were designed and built with the specific purpose to allow the effects of varied gravity (high, normal, reduced) during synthesis of titanium sol-gels to be studied. The first system was a centrifuge capable of providing high gravity environments of up to 70 g’s for extended periods, whilst applying a 100 mbar vacuum and a temperature of 40-50oC to the reaction chambers. The second system to be used in the QUT Microgravity Drop Tower Facility was also required to provide the same thermal and vacuum conditions used in the centrifuge, but had to operate autonomously during free fall. Through the use of post synthesis characterisation techniques such as Raman Spectroscopy, X-Ray diffraction (XRD) and N2 adsorption, it was found that increased gravity levels during synthesis, had the greatest effect on the final products. Samples produced in reduced and normal gravity appeared to form amorphous gels containing very small particles with moderate surface areas. Whereas crystalline anatase (TiO2), was found to form in samples synthesised above 5 g with significant increases in crystallinity, particle size and surface area observed when samples were produced at gravity levels up to 70 g. It is proposed that for samples produced in higher gravity, an increased concentration gradient of water is forms at the bottom of the reacting film due to forced convection. The particles formed in higher gravity diffuse downward towards this excess of water, which favours the condensation reaction of remaining sol gel precursors with the particles promoting increased particle growth. Due to the removal of downward convection in reduced gravity, particle growth due to condensation reaction processes are physically hindered hydrolysis reactions favoured instead. Another significant finding from this work was that anatase could be produced at relatively low temperatures of 40-50oC instead of the conventional method of calcination above 450oC solely through sol-gel synthesis at higher gravity levels. It is hoped that the outcomes of this research will lead to an increased understanding of the effects of gravity on chemical synthesis of titanium sol-gel, potentially leading to the development of improved products suitable for diverse applications such as semiconductor or catalyst materials as well as significantly reducing production and energy costs through manufacturing these materials at significantly lower temperatures.

Improved ultra-performance liquid chromatography with electrospray ionization quadrupole-time-of-flight high-definition mass spectrometry method for the rapid analysis of the chemical constituents of a typical medical formula: Liuwei Dihuang Wan

Liuwei Dihuang Wan (LWD), a classic Chinese medicinal formulae, has been used to improve or restore declined functions related to aging and geriatric diseases, such as impaired mobility, vision, hearing, cognition and memory. It has attracted increasingly much attention as one of the most popular and valuable herbal medicines. However, the systematic analysis of the chemical constituents of LDW is difficult and thus has not been well established. In this paper, a rapid, sensitive and reliable ultra-performance liquid chromatography with electrospray ionization quadrupole time-of-flight high-definition mass spectrometry (UPLC-ESI-Q-TOF-MS) method with automated MetaboLynx analysis in positive and negative ion mode was established to characterize the chemical constituents of LDW. The analysis was performed on a Waters UPLCTM HSS T3 using a gradient elution system. MS/MS fragmentation behavior was proposed for aiding the structural identification of the components. Under the optimized conditions, a total of 50 peaks were tentatively characterized by comparing the retention time and MS data. It is concluded that a rapid and robust platform based on UPLC-ESI-Q-TOF-MS has been successfully developed for globally identifying multiple-constituents of traditional Chinese medicine prescriptions. This is the first report on systematic analysis of the chemical constituents of LDW. This article is protected by copyright. All rights reserved.

Identification and simultaneous analysis of harmane, harmol, isovitexin and vitexin in Passiflora incarnata extracts with a novel HLPC method

A high performance liquid chromatographic method for the simultaneous analysis of two flavonoids (iso-vitexin and vitexin), and three indole alkaloids (harmane, harmine, and harmol) was developed. This method was then utilised to quantitate levels of these five constituents in methanolic extracts of Australian Passiflora incarnata. HPLC analysis was performed using a Waters™ Novapak C18 (150 × 4 mm, 4 μm) column, with a gradient solvent system of methanol-water-acetic acid. Detection was achieved by PDA UV (254 nm) and fluorescence (excitation 254 nm, emission 414 nm), utilising the external standard method to obtain quantification.

Histogram of weighted local directions for gait recognition

In this paper, we explore the effectiveness of patch-based gradient feature extraction methods when applied to appearance-based gait recognition. Extending existing popular feature extraction methods such as HOG and LDP, we propose a novel technique which we term the Histogram of Weighted Local Directions (HWLD). These 3 methods are applied to gait recognition using the GEI feature, with classification performed using SRC. Evaluations on the CASIA and OULP datasets show significant improvements using these patch-based methods over existing implementations, with the proposed method achieving the highest recognition rate for the respective datasets. In addition, the HWLD can easily be extended to 3D, which we demonstrate using the GEV feature on the DGD dataset, observing improvements in performance.

Simultaneous determination of aldrin, dieldrin, endrin, heptachlor, and p,p'-DDT in medicinal plant extracts using a novel high performance liquid chromatography method

A high performance liquid chromatographic method for the simultaneous determination of five organochlorine pesticides (aldrin, p,p’-DDT, dieldrin, endrin, and heptachlor) was developed. The method was used to determine the levels of these pesticides in medicinal plant samples. Analysis was carried out using a Merck LiChrospher 100 RP C18 (5 μm) column with a gradient solvent system of acetonitrile-water and PDA UV detection (224 nm). Quantification was carried out by the external standard method. The limit of detection for the utilized method was below the local legal limits (ANZFA) for similar plant materials for all 5 pesticides excepting endrin. Medicinal plant extracts were further analyzed by conventional GC-ECD and GC-NPD means using SPE and GPC cleanup as required.