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.

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.

Brucite Microbialites in Living Coral Skeletons: Indicators of Extreme Microenvironments in Shallow-Marine Settings

Brucite [Mg(OH)2] microbialites occur in vacated interseptal spaces of living scleractinian coral colonies (Acropora, Pocillopora, Porites) from subtidal and intertidal settings in the Great Barrier Reef, Australia, and subtidal Montastraea from the Florida Keys, United States. Brucite encrusts microbial filaments of endobionts (i.e., fungi, green algae, cyanobacteria) growing under organic biofilms; the brucite distribution is patchy both within interseptal spaces and within coralla. Although brucite is undersaturated in seawater, its precipitation was apparently induced in the corals by lowered pCO2 and increased pH within microenvironments protected by microbial biofilms. The occurrence of brucite in shallow-marine settings highlights the importance of microenvironments in the formation and early diagenesis of marine carbonates. Significantly, the brucite precipitates discovered in microenvironments in these corals show that early diagenetic products do not necessarily reflect ambient seawater chemistry. Errors in environmental interpretation may arise where unidentified precipitates occur in microenvironments in skeletal carbonates that are subsequently utilized as geochemical seawater proxies.

Institutional problems of the emerging scientific community in Iran

This article attempts to explore the concept of scientific community at the macro-national level in the context of Iran. Institutionalisation of science and its professional growth has been constrained by several factors. The article first conceptualises the notion of science community as found in the literature in the context of Iran, and attempts to map through some indicators. The main focus, however, lies in mapping some institutional problems through empirical research. This was undertaken in 2002–04 in order to analyse the structure of the scientific community in Iran in the ‘exact sciences’ (mathematics, physics, chemistry, biology and earth sciences). The empirical work was done in two complementary perspectives: through a questionnaire and statistical analysis of it, and through semistructured interviews with the researchers. There are number of problems confronting scientists in Iran. Facilities provided by institutions is one of the major problems of research. Another is the tenuous cooperation among scientists. This is reported by most of the researchers, who deplore the lack of cooperation among their group. Relationships are mostly with the Ph.D. students and only marginally with colleagues. Our research shows that the more brilliant the scientists, the more frustrated they are from scientific institutions in Iran. Medium-range researchers seem to be much happier about the scientific institution to which they belong than the brighter scholars. The scientific institutions in Iran seem to be built for the needs of the former rather than the latter. These institutions seem not to play a positive role in the case of the best scientists. On the whole, many ingredients of the scientific community, at least at its inception, are present among Iranian scientists: the strong desire for scientific achievement in spite of personal, institutional and economic problems.