The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.

Anthropogenic and Natural Organohalogen Compounds in Blubber of Dolphins and Dugongs (Dugong dugon) from Northeastern Australia

A range of organohalogen compounds (10 polychlorinated biphenyl [PCB] congeners, DDT and metabolites, chlordane-related compounds, the potential natural organochlorine compound Q1, toxaphene, hexachlorobenzene, hexachlorocyclohexanes, dieldrin, and several yet unidentified brominated compounds) were detected in the blubber of four bottlenose dolphins (Tursiops truncatus), one common dolphin (Delphinus delphis), and seven dugongs (Dugong dugon), as well as in adipose tissue of a green turtle (Chelonia mydas) and a python (Morelia spilota) from northeast Queensland (Australia). The green turtle and dugongs accumulated lower organohalogen levels than the dolphins. Lower levels in dugongs were expected because this species is exclusively herbivorous. Highest PCB and DDT levels recorded in dugongs were 209 and 173 mug/kg lipids, respectively. Levels of the nonanthropogenic heptachlorinated compound Q1 (highest level in dugongs was 160 mug/kg lipids) were estimated using the ECD response factor of trans-nonachlor. Highest organohalogen levels were found in blubber of dolphins for sumDDT (575-52,500 mug/kg) and PCBs (600-25,500 mug/kg lipids). Furthermore, Q1 was a major organohalogen detected in all samples analyzed, ranging from 450 -9,100 mug/kg lipids. The highest concentration of Q1 determined in this study represents the highest concentration reported to date in an environmental sample. Levels of chlordane-related compounds were also high (280-7,700 mug/kg, mainly derived from trans-nonachlor), but concentrations of hexachlorobenzene, hexachlorocyclohexanes, dieldrin, and toxaphene were relatively low and contributed little to the overall organohalogen contamination. Furthermore, a series of three major (BC-1, BC-2, and BC-3) and six minor (BC-4 through BC-9) unknown brominated compounds were observable by extracting m/z 79 and m/z 81 from the GC/ECNI-MS full scan run. Structural proposals were made for the two major recalcitrant compounds (referred to as BC-1 and BC-2). BC-2 appears to be a tetrabromo-methoxy-diphenylether (512 u) and BC-1 has 14 u (corresponding with an additional CH2 group) more relative to BC-1. In general the organohalogen pattern observed in blubber of dolphins was different compared to similar samples from other locations in the world, which is apparent from the fact that the four major abundant signals in the GC/ECD chromatogram. of D. delphis originated from the four unknown compounds Q1, BC-1, BC-2, and BC-3.

Pattern-based phylogenetic distance estimation and tree reconstruction

We have developed an alignment-free method that calculates phylogenetic distances using a maximum-likelihood approach for a model of sequence change on patterns that are discovered in unaligned sequences. To evaluate the phylogenetic accuracy of our method, and to conduct a comprehensive comparison of existing alignment-free methods (freely available as Python package decaf+py at http://www.bioinformatics.org.au), we have created a data set of reference trees covering a wide range of phylogenetic distances. Amino acid sequences were evolved along the trees and input to the tested methods; from their calculated distances we infered trees whose topologies we compared to the reference trees. We find our pattern-based method statistically superior to all other tested alignment-free methods. We also demonstrate the general advantage of alignment-free methods over an approach based on automated alignments when sequences violate the assumption of collinearity. Similarly, we compare methods on empirical data from an existing alignment benchmark set that we used to derive reference distances and trees. Our pattern-based approach yields distances that show a linear relationship to reference distances over a substantially longer range than other alignment-free methods. The pattern-based approach outperforms alignment-free methods and its phylogenetic accuracy is statistically indistinguishable from alignment-based distances.

Towards a Modeling Environment for Earth Systems Simulations

The developments of models in Earth Sciences, e.g. for earthquake prediction and for the simulation of mantel convection, are fare from being finalized. Therefore there is a need for a modelling environment that allows scientist to implement and test new models in an easy but flexible way. After been verified, the models should be easy to apply within its scope, typically by setting input parameters through a GUI or web services. It should be possible to link certain parameters to external data sources, such as databases and other simulation codes. Moreover, as typically large-scale meshes have to be used to achieve appropriate resolutions, the computational efficiency of the underlying numerical methods is important. Conceptional this leads to a software system with three major layers: the application layer, the mathematical layer, and the numerical algorithm layer. The latter is implemented as a C/C++ library to solve a basic, computational intensive linear problem, such as a linear partial differential equation. The mathematical layer allows the model developer to define his model and to implement high level solution algorithms (e.g. Newton-Raphson scheme, Crank-Nicholson scheme) or choose these algorithms form an algorithm library. The kernels of the model are generic, typically linear, solvers provided through the numerical algorithm layer. Finally, to provide an easy-to-use application environment, a web interface is (semi-automatically) built to edit the XML input file for the modelling code. In the talk, we will discuss the advantages and disadvantages of this concept in more details. We will also present the modelling environment escript which is a prototype implementation toward such a software system in Python (see www.python.org). Key components of escript are the Data class and the PDE class. Objects of the Data class allow generating, holding, accessing, and manipulating data, in such a way that the actual, in the particular context best, representation is transparent to the user. They are also the key to establish connections with external data sources. PDE class objects are describing (linear) partial differential equation objects to be solved by a numerical library. The current implementation of escript has been linked to the finite element code Finley to solve general linear partial differential equations. We will give a few simple examples which will illustrate the usage escript. Moreover, we show the usage of escript together with Finley for the modelling of interacting fault systems and for the simulation of mantel convection.