The use of genetic algorithms to model protoplanetary discs

The protoplanetary discs of T Tauri and Herbig Ae/Be stars have previously been studied using geometric disc models to fit their spectral energy distribution (SED). The simulations provide a means to reproduce the signatures of various circumstellar structures, which are related to different levels of infrared excess. With the aim of improving our previous model, which assumed a simple flat-disc configuration, we adopt here a reprocessing flared-disc model that assumes hydrostatic, radiative equilibrium. We have developed a method to optimize the parameter estimation based on genetic algorithms (GAs). This paper describes the implementation of the new code, which has been applied to Herbig stars from the Pico dos Dias Survey catalogue, in order to illustrate the quality of the fitting for a variety of SED shapes. The star AB Aur was used as a test of the GA parameter estimation, and demonstrates that the new code reproduces successfully a canonical example of the flared-disc model. The GA method gives a good quality of fit, but the range of input parameters must be chosen with caution, as unrealistic disc parameters can be derived. It is confirmed that the flared-disc model fits the flattened SEDs typical of Herbig stars; however, embedded objects (increasing SED slope) and debris discs (steeply decreasing SED slope) are not well fitted with this configuration. Even considering the limitation of the derived parameters, the automatic process of SED fitting provides an interesting tool for the statistical analysis of the circumstellar luminosity of large samples of young stars.     

Promoting geological heritage through postage stamps

Postage stamps are useful in promoting geological heritage. This concept is explored through an analysis of British and New Zealand stamps, which portray a variety of different types of geological heritage. Coastal, glacial, fluvial and limestone scenery all feature, together with volcanic landforms such as the Giant's Causeway. Other stamps celebrate the work of historically significant figures who have contributed to geology, such as Charles Darwin. The field of palaeontology is represented in stamps, for example, featuring dinosaurs. Others portray economic activities and engineering achievements associated with geology such as oil and gas exploration. Industrial heritage is depicted by stamps featuring, for example, abandoned tin mines, while building stones are represented in stamps of ancient monuments, bridges, cathedrals and domestic architecture. Perhaps subjects for future stamp issues in the UK might draw upon the recently published 100 top geosites list.     

Carbon dioxide and carbon monoxide photoproduction quantum yields in the Delaware Estuary

Photochemical mineralization of dissolved organic matter (DOM) plays an important role in the cycling of carbon in estuarine systems. A key to modeling this process is knowledge of apparent quantum yields (AQYs) for the photochemical products. Here we determined spectral AQYs for carbon dioxide (CO₂) and carbon monoxide (CO), the main products of DOM photomineralization, along the main axis of the Delaware Estuary. Apparent quantum yields for CO₂ photoproduction were determined shipboard using a multi-spectral irradiation system. Carbon monoxide AQYs were determined in stored samples by employing a narrow band spectral irradiation system. A single AQY spectrum described carbon dioxide photochemical production within the estuary whereas CO AQY spectra varied with salinity, suggesting different precursors and mechanisms for the production of these two species. CO₂ AQYs were used along with shipboard measurements of DOM absorbance and solar irradiance to calculate photoproduction rates. Calculated CO₂ photoproduction rates agreed with directly measured rates (2 to 4 μM CO₂ d^? 1) within experimental error, supporting the further development and use of AQYs to calculate regional-scale photochemical fluxes.     

Noble gases in E-chondrites

Noble gas data are reported for 12 E-chondrites. Combined with literature data, they show that K-Ar ages are >4 Æ for 14 out of 18 meteorites, yet U, Th-He ages are often shorter, perhaps due to late, mild reheating. Cosmic-ray exposure ages differ systematically between types 4 and 6, with E4's mostly below 16 Myr and E6's above 30 Myr. This may mean that the E-chondrite parent body contains predominantly a single petrologic type on the (~ 1 km) scale of individual impacts, in contrast to the more thoroughly mixed parent bodies of the ordinary chondrites.The heavy noble gases consist of at least two primordial components: the usual planetary component (${}^{36}\text{Ar}{}^{132}\text{Xe}\text{~ 80}$) and a less fractionated, ‘subsolar’ component ($\text{2700 ≤}{}^{36}\text{Ar}{}^{132}\text{Xe}\text{≤ 3800}$). The latter is found in highest concentration in the E4 chondrite South Oman (³⁶Ar = 760 × 10^?8cc/g, ${}^{36}\text{Ar}{}^{132}\text{Xe}\text{= 2700}$). The isotopic compositions of both components are similar to typical planetary values, indicating that some factor other than mass controlled the noble gas elemental ratios. The heavy Xe isotopes occasionally show some of the lowest ${}^{134}\text{Xe}{}^{132}\text{Xe}$ and ${}^{136}\text{Xe}{}^{132}\text{Xe}$ ratios measured in bulk chondrites, suggestive of nearly fission-free Xe (e.g. ${}^{136}\text{Xe}{}^{132}\text{Xe}\text{= 0.3095 ± 0.0020}$). Amounts of planetary gas in E4 E6 chondrites fall in the range for ordinary chondrites of types 4–6, but, in contrast to the ordinary chondrites. fail to correlate with petrologic type or volatile trace element contents. Another unusual feature of E-chondrites is that primordial Ne is present even in most 4's and 5's (²⁰Ne_p ~ 1 to 7 × 10^?8cc/g). with an isotopic composition consistent with planetary Ne.Analyses of mineral separates show that the planetary gases are concentrated in an HF- and HCl-insoluble mineral similar to phase Q, the poorly characterized, HNO₃-soluble carrier of primordial gases in carbonaceous and ordinary chondrites. The subsolar gases, on the other hand, are located in an HCl- and HNO₃-resistant phase, possibly enstatite or a minor phase included in enstatite. Much of the ¹²⁹Xe_r (50% for E4's, > 70% for E6's) is in HCl-resistant but HF-soluble sites, suggestive of a silicate.A similar subsolar component may be responsible for the high ${}^{36}\text{Ar}{}^{132}\text{Xe}$ ratios of some C3's, unequilibrated ordinary chondrites, and the unique aubrite Shallowater. The planet Venus also has a high $\text{Ar}\text{Kr}$ ratio, well above the planetary range, and hence may have acquired its noble gases from an E-chondrite-like material, similar to South Oman.     

Bicyclic sesquiterpenoids and diterpenoids in Australian crude oils

Bicyclanes previously reported only in heavily biodegraded Texas Gulf Coast crudes have been found to be ubiquitous in Australian crude oils of non-marine origin from four different basins. The compounds are present in oils, thought to be derived from the same or similar sources, that have undergone varying degrees of biodegradation. They are also found to be present in oils of different geological age. In addition a series of tricyclic diterpenoid hydrocarbons was common to four oils from the Gippsland Basin. Four of these compounds had the molecular formula C₂₀H₃₄ and mass spectral fragmentation patterns suggested they were mono-unsaturated diterpenoids. The presence of unsaturated diterpenoids in crude oils appears to be a unique observation. It is proposed that the diterpenoids may be the source of the bicyclanes also observed in these oils.     

High pressure metamorphism in the Scourian of NW Scotland: Evidence from garnet granulites

Ultramafic and mafic granulites from Archaean gneisses in N.W. Scotland (the Scourian) show evidence of two periods of granulite facies mineral growth. The first produced a high pressure clinopyroxene +garnet±plagioclase assemblage at an estimatedP-T of 12–15 kb and 1,000° C. Uplift of the complex caused partial breakdown of the garnet by reaction with clinopyroxene to produce orthopyroxene +plagioclase ±spinel±amphibole symplectites, at an estimatedP-T of 10–14 kb and 800°–900° C. Garnet stability is shown to depend on both whole-rock Fe/Mg ratios and onP-T conditions. The pressures imply crustal thicknesses in the Archaean of least 35–45 km.     

Extinct I129 in C3 chondrites

Eight C3 chondrites were examined by the I¹²⁹Xe¹²⁹ dating method, to see whether their IXe “ages” (better, initial $\text{I}{}^{129}\text{I}{}^{127}$ratios ≡ R₀) correlate with any other properties. The R₀'s range from 1.60 × 10^?4 to 1.09 × 10^?4, corresponding to IXe ages from 2.0 Myr before to 6.7 Myr after Murchison magnetite. Three C3O's (Lancé, Felix, Ornans) have essentially indistinguishable R₀'s of (1.41 ± 0.13) to (1.17 ± 0.10) × 10^?4; the fourth C3O, Warrenton, is undatable owing to homogenization of radiogenic and trapped Xe.Four C3V's show a distinct spread: Vigarano and Grosnaja are highest [R₀ = (1.60 ± 0.07) and (1.57 ± 0.14) × 10^?4], Mokoia is intermediate, and Kaba is lowest [R₀ = (1.38 ± 0.06) and (1.09 ± 0.10) × 10^?4]. Literature values for Allende place it near Kaba. These R₀'s correlate inversely with 4 other properties: I-, Br-, and Cd-content, and olivine composition, both percent mean deviation (PMD) and proportion of iron-poor olivine grains (≤2% fayalite).It is difficult to accept the ～9 Myr spread in R₀ as a true age, reflecting either nebular or parent-body processes. This time span is more than an order of magnitude longer than the lifetime of the solar nebula inferred from astronomical evidence. Nor does the degree of thermal metamorphism, which is slight for C3's anyway, correlate with R₀. A more plausible interpretation is that the variations in R₀ reflect mainly isotopic heterogeneity of iodine. The simplest model that accounts for the correlations with R₀ involves mixing of two iodine components in the solar nebula, associated with gas and grains, respectively. The second, of lower $\text{I}{}^{129}\text{I}{}^{127}$ ratio, predominated at later times and thus became enriched in late-formed meteorites, along with other volatiles such as Cd and Br. The low Fe content and large PMD of olivine may reflect either less metamorphism owing to shallow location in the parent body, or greater reduction of Fe²⁺ during chondrule formation.     

Quantitative P-T paths from zoned minerals: Theory and tectonic applications

An analytical approach to the analysis of zoning profiles in minerals is presented that simultaneously accounts for all of the possible continuous reactions that may be operative in a given assemblage. The method involves deriving a system of simultaneous linear differential equations consisting of a Gibbs-Duhem equation for each phase, a set of linearly independent stoichiometric relations among the chemical potentials of phase components in the assemblage, and a set of equations describing the total differential of the slope of the tangent plane to the Gibbs free energy surface of solid solution phases. The variables are the differentials of T, P, chemical potentials of all phase components, and independent compositional terms of solid solution phases. The required input data are entropies, volumes, the compositions of coexisting phases at a reference P and T, and an expression for the curvature of the Gibbs functions for solid solution phases. Results derived are slopes of isopleths (dP/dT, dX/dT or dX/dP) which can be used to contour P-T diagrams with mineral composition.To interpret mineral zoning, T and P can be expressed as functions of n independent composition parameters, where n is the variance of the mineral assemblage. The total differentials of P and T are differential equations that can be solved by finite difference techniques using the derivatives obtained from the analytical formulation of phase equilibria.Results calculated from Zone I and Zone IV garnets of Tracy et al. (1976) indicate that Zone I garnets grew while T increased (T+72° C) and P decreased sharply (P–3 kb). Zone IV garnets zoned in response to decreasing T (T–17° C) and P (P–1 kb). A P-T path calculated for a zoned garnet from the Greinerschiefer series, western Tauern Window, Austria, also indicates growth during decompression (–3kb) and heating (T+15° C). A P-T path calculated for the Wissahickon schist (Crawford and Mark 1982) indicates growth during cooling and compression (T–25 C, P+2.2 kb). The calculated P-T paths differ according to structural environment and can be used to relate mineral growth to tectonic processes.