UV transit observations of EUV-heated expanded thermospheres of Earth-like exoplanets around M-stars: testing atmosphere evolution scenarios

The detection and investigation of EUV heated, extended and non-hydrostatic upper atmospheres around terrestrial exoplanets would provide important insights into the interaction of the host stars plasma environment as well as the evolution of Earth-type planets their atmospheres and possible magnetic environments. We discuss different scenarios where one can expect that Earth-like planets should experience non-hydrostatic upper atmosphere conditions so that dynamically outward flowing neutral atoms can interact with the stellar plasma flow so that huge hydrogen coronae and energetic neutral atoms (ENA) can be produced via charge exchange. By observing the size of the extended upper atmospheres and related ENA-clouds and by determining the velocities of the surrounding hydrogen atoms, conclusions can be drawn in respect to the origin of these features. Due to the large number of M-type stars in our neighbourhood and their long periods of strong and moderate stellar activity in comparison to G-stars, we expect that M-type stars represent the most promising candidates for the detection of hydrogen ENA-clouds and the subsequent study of the interaction between the host star and the planets?? upper atmosphere. We show that the low mass of M-type stars also makes them preferable targets to observe extended hydrogen clouds around terrestrial exoplanets with a mass as low as one Earth mass. Transit follow-up observations in the UV-range of terrestrial exoplanets around M-type stars with space observatories such as the World Space Observatory-UV (WSO-UV) would provide a unique opportunity to shed more light on the early evolution of Earth-like planets, including those of our own Solar System.     

Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2009

Every three years the IAU Working Group on Cartographic Coordinates and Rotational Elements revises tables giving the directions of the poles of rotation and the prime meridians of the planets, satellites, minor planets, and comets. This report takes into account the IAU Working Group for Planetary System Nomenclature (WGPSN) and the IAU Committee on Small Body Nomenclature (CSBN) definition of dwarf planets, introduces improved values for the pole and rotation rate of Mercury, returns the rotation rate of Jupiter to a previous value, introduces improved values for the rotation of five satellites of Saturn, and adds the equatorial radius of the Sun for comparison. It also adds or updates size and shape information for the Earth, Mars?? satellites Deimos and Phobos, the four Galilean satellites of Jupiter, and 22 satellites of Saturn. Pole, rotation, and size information has been added for the asteroids (21) Lutetia, (511) Davida, and (2867) ?teins. Pole and rotation information has been added for (2) Pallas and (21) Lutetia. Pole and rotation and mean radius information has been added for (1) Ceres. Pole information has been updated for (4) Vesta. The high precision realization for the pole and rotation rate of the Moon is updated. Alternative orientation models for Mars, Jupiter, and Saturn are noted. The Working Group also reaffirms that once an observable feature at a defined longitude is chosen, a longitude definition origin should not change except under unusual circumstances. It is also noted that alternative coordinate systems may exist for various (e.g. dynamical) purposes, but specific cartographic coordinate system information continues to be recommended for each body. The Working Group elaborates on its purpose, and also announces its plans to occasionally provide limited updates to its recommendations via its website, in order to address community needs for some updates more often than every 3 years. Brief recommendations are also made to the general planetary community regarding the need for controlled products, and improved or consensus rotation models for Mars, Jupiter, and Saturn.     

Sensitivity of ENSO characteristics to a new interactive flux correction scheme in a coupled GCM

A fast coupled global climate model (CGCM) is used to study the sensitivity of El Ni?o Southern Oscillation (ENSO) characteristics to a new interactive flux correction scheme. With no flux correction applied our CGCM reveals typical bias in the background state: for instance, the cold tongue in the tropical east Pacific becomes too cold, thus degrading atmospheric sensitivity to variations of sea surface temperature (SST). Sufficient atmospheric sensitivity is essential to ENSO. Our adjustment scheme aims to sustain atmospheric sensitivity by counteracting the SST drift in the model. With reduced bias in the forcing of the atmosphere, the CGCM displays ENSO-type variability that otherwise is absent. The adjustment approach employs a one-way anomaly coupling from the ocean to the atmosphere: heat fluxes seen by the ocean are based on full SST, while heat fluxes seen by the atmosphere are based on anomalies of SST. The latter requires knowledge of the model??s climatological SST field, which is accumulated interactively in the spin-up phase (??training??). Applying the flux correction already during the training period (by utilizing the evolving SST climatology) is necessary for efficiently reducing the bias. The combination of corrected fluxes seen by the atmosphere and uncorrected fluxes seen by the ocean implies a restoring mechanism that counteracts the bias and allows for long stable integrations in our CGCM. A suite of sensitivity runs with varying training periods is utilized to study the effect of different levels of bias in the background state on important ENSO properties. Increased duration of training amplifies the coupled sensitivity in our model and leads to stronger amplitudes and longer periods of the Nino3.4 index, increased emphasis of warm events that is reflected in enhanced skewness, and more pronounced teleconnections in the Pacific. Furthermore, with longer training durations we observe a mode switch of ENSO in our model that closely resembles the observed mode switch related to the mid-1970s ??climate shift??.     

Compositional Bayesian indicator estimation

Indicator kriging is widely used for mapping spatial binary variables and for estimating the global and local spatial distributions of variables in geosciences. For continuous random variables, indicator kriging gives an estimate of the cumulative distribution function, for a given threshold, which is then the estimate of a probability. Like any other kriging procedure, indicator kriging provides an estimation variance that, although not often used in applications, should be taken into account as it assesses the uncertainty of the estimate. An alternative approach to indicator estimation is proposed in this paper. In this alternative approach the complete probability density function of the indicator estimate is evaluated. The procedure is described in a Bayesian framework, using a multivariate Gaussian likelihood and an a priori distribution which are both combined according to Bayes theorem in order to obtain a posterior distribution for the indicator estimate. From this posterior distribution, point estimates, interval estimates and uncertainty measures can be obtained. Among the point estimates, the median of the posterior distribution is the maximum entropy estimate because there is a fifty-fifty chance of the unknown value of the estimate being larger or smaller than the median; that is, there is maximum uncertainty in the choice between two alternatives. Thus in some sense, the latter is an indicator estimator, alternative to the kriging estimator, that includes its own uncertainty. On the other hand, the mode of the posterior distribution estimator, assuming a uniform prior, is coincidental with the simple kriging estimator. Additionally, because the indicator estimate can be considered as a two-part composition which domain of definition is the simplex, the method is extended to compositional Bayesian indicator estimation. Bayesian indicator estimation and compositional Bayesian indicator estimation are illustrated with an environmental case study in which the probability of the content of a geochemical element in soil being over a particular threshold is of interest. The computer codes and its user guides are public domain and freely available.     

Angular momentum exchange during secular migration of two-planet systems

We investigate the secular dynamics of two-planet coplanar systems evolving under mutual gravitational interactions and dissipative forces. We consider two mechanisms responsible for the planetary migration: star-planet (or planet-satellite) tidal interactions and interactions of a planet with a gaseous disc. We show that each migration mechanism is characterized by a specific law of orbital angular momentum exchange. Calculating stationary solutions of the conservative secular problem and taking into account the orbital angular momentum leakage, we trace the evolutionary routes followed by the planet pairs during the migration process. This procedure allows us to recover the dynamical history of two-planet systems and constrain parameters of the involved physical processes.     

Influence of mid-latitude circulation on upper Indus basin precipitation: the explicit role of irrigation

Since much of the flow of the Indus River originates in the Himalayas, Karakoram and Hindu Kush Mountains, an understanding of weather characteristics leading to precipitation over the region is essential for water resources management. This study examines the influence of upper level mid-latitude circulation on the summer precipitation over upper Indus basin (UIB). Using reanalysis data, a geopotential height index (GH) is defined at 200 hPa over central Asia, which has a significant correlation with the precipitation over UIB. GH has also shown significant correlation with the heat low (over Iran and Afghanistan and adjoining Pakistan), easterly shear of zonal winds (associated with central Asian high) and evapotranspiration (over UIB). It is argued that the geopotential height index has the potential to serve as a precursor for the precipitation over UIB. In order to assess the influence of irrigation on precipitation over UIB, a simplified irrigation scheme has been developed and applied to the regional climate model REMO. It has been shown that both versions of REMO (with and without irrigation) show significant correlations of GH with easterly wind shear and heat low. However contrary to reanalysis and the REMO version with irrigation, the REMO version without irrigation does not show any correlation between GH index and evapotranspiration as well as between geopotential height and precipitation over UIB, which is further confirmed by the quantitative analysis of extreme precipitation events over UIB. It is concluded that although atmospheric moisture over coastal Arabian sea region, triggered by wind shear and advected northward due to heat low, also contribute to the UIB precipitation. However for the availability of necessary moisture for precipitation over UIB, the major role is played by the evapotranspiration of water from irrigation. From the results it may also be inferred that the representation of irrigated water in climate models is unavoidable for studying the impact of global warming over the region.     

Zircon formation during fluid circulation in eclogites (Monviso, Western Alps): implications for Zr and Hf budget in subduction zones

The zircons from an eclogite and an enclosed eclogite-facies vein from the Monviso ophiolite (Western Alps) display contrasting chemical and morphologic features and document different stages of the evolution of the ophiolite. The zircons from the eclogite show a typical magmatic zoning and are enriched in heavy rare earth elements (HREEs) over middle rare earth elements (MREEs) and have an accentuated negative Eu anomaly, which indicates that the grains co-crystallised with plagioclase. These magmatic zircons document the formation of oceanic crust at 163 ± 2 Ma. In contrast, zircons from the vein contain inclusions of garnet, omphacite, and rutile, which indicate that they crystallised under eclogite-facies conditions. The vein zircons have Th/U ratios < 0.09, lack Eu anomalies, and are only weakly enriched in HREE with respect to MREE. These features are consistent with a garnet-bearing, plagioclase-free, i.e., eclogite-facies paragenesis. Vein zircons yield an age of 45 ± 1 Ma, which is evidence for Eocene subduction-zone metamorphism of the Monviso ophiolite.In the vein, the apparent coexistence of zircon, omphacite, and garnet permits the determination of a set of trace element distribution coefficients among these minerals at high pressure. This set of partitioning can demonstrate chemical equilibrium among these phases in rocks that show less clear evidence of textural equilibrium. In addition, zircon age can now be linked to sensors of metamorphic pressure-temperature conditions. The presence of zircon and rutile in the vein is another example of high field strength element (HFSE) mobility over short distances in aqueous fluids at eclogite-facies conditions. However, the concentrations of Zr and Hf in the aqueous fluid are estimated to be at least a factor of 10 less than primitive mantle values.Mass balance calculations demonstrate that zircon hosts > 95% of the bulk Zr, 90% of Hf, and ∼25% of U in the vein. Zircon is a residual phase in subducted basalts and sediments up to temperatures of at least 800 to 900 °C. Therefore, residual zircon in subducted crust, together with rutile, control the HFSE in liberated subduction zone fluids/melts and might be partly responsible for negative Zr and Hf anomalies in subduction zone magmas.     

Bayesian Spatial Modelling of Gamma Ray Count Data

Gamma ray logging is a method routinely employed by geophysicists and environmental engineers in site geology evaluations. Modelling of gamma ray data from individual boreholes assists in the local identification of major lithological changes; modelling these data from a network of boreholes assists with lithological mapping and spatial stratigraphic correlation. In this paper we employ Bayesian spatial partition models to analyse gamma ray data spatially. In particular, a spatial partition is defined via a Voronoi tessellation and the mean intensity is assumed constant in each cell of the partition. The number of vertices generating the tessellation as well as the locations of vertices are assumed unknown, and uncertainty about these quantities is described via a hierarchical prior distribution. We describe the advantages of the spatial partition modelling approach in the context of smoothing gamma ray count data and describe an implementation that may be extended to the fitting of a more general model than a constant mean within each cell of the partition. As an illustration of the methodology we consider a data set collected from a network of eight boreholes, which is part of a geophysical study to assist in mapping the lithology of a site. Gamma ray logs are linked with geological information from cores and the spatial analysis of log data assists with predicting the lithology at unsampled locations.     

Dynamics and Variability of Terra Nova Bay Polynya

Abstract. We present a process study on the dynamics and variability of the Terra Nova Bay polynya in the western sector of the Ross Sea. The air-sea heat exchange is known to be particularly large in polynya during the winter, when differences between air and sea temperatures are large. We apply a 1-D model (Pease, 1987; Van Woert, 1999a, 1999b), which is modified in the latent heat parameterisation in order to account for time-dependent relative humidity and cloud coverage. Furthermore, the Ice Collection Depth is correlated linearly with a variable wind speed. The model is forced with two different meteorological data sets: the operational analysis of the European Center for Medium Range Weather Forecasts atmospheric data set and the meteorological parameters measured by an Automatic Weather Station located on the coast of Terra Nova Bay. The results are compared in terms of polynya extension, ice, and High Salinity Shelf Water production. According to the two different wind velocities, the results obtained from the different data sets clearly differ. Qualitatively, however, the results are in good agreement.