A fast Bayesian approach to discrete object detection in astronomical data sets – PowellSnakes I

A new fast Bayesian approach is introduced for the detection of discrete objects immersed in a diffuse background. This new method, called PowellSnakes, speeds up traditional Bayesian techniques by (i) replacing the standard form of the likelihood for the parameters characterizing the discrete objects by an alternative exact form that is much quicker to evaluate; (ii) using a simultaneous multiple minimization code based on Powell's direction set algorithm to locate rapidly the local maxima in the posterior and (iii) deciding whether each located posterior peak corresponds to a real object by performing a Bayesian model selection using an approximate evidence value based on a local Gaussian approximation to the peak. The construction of this Gaussian approximation also provides the covariance matrix of the uncertainties in the derived parameter values for the object in question. This new approach provides a speed up in performance by a factor of '100' as compared to existing Bayesian source extraction methods that use Monte Carlo Markov chain to explore the parameter space, such as that presented by Hobson & McLachlan. The method can be implemented in either real or Fourier space. In the case of objects embedded in a homogeneous random field, working in Fourier space provides a further speed up that takes advantage of the fact that the correlation matrix of the background is circulant. We illustrate the capabilities of the method by applying to some simplified toy models. Furthermore, PowellSnakes has the advantage of consistently defining the threshold for acceptance/rejection based on priors which cannot be said of the frequentist methods. We present here the first implementation of this technique (version I). Further improvements to this implementation are currently under investigation and will be published shortly. The application of the method to realistic simulated Planck observations will be presented in a forthcoming publication.     

Surface magnetic fields on two accreting T Tauri stars: CV Cha and CR Cha

We have produced brightness and magnetic field maps of the surfaces of CV Cha and CR Cha: two actively accreting G- and K-type T Tauri stars in the Chamaeleon I star-forming cloud with ages of 3–5 Myr. Our magnetic field maps show evidence for strong, complex multipolar fields similar to those obtained for young rapidly rotating main-sequence stars. Brightness maps indicate the presence of dark polar caps and low-latitude spots – these brightness maps are very similar to those obtained for other pre-main-sequence and rapidly rotating main-sequence stars.
Only two other classical T Tauri stars have been studied using similar techniques so far: V2129 Oph and BP Tau. CV Cha and CR Cha show magnetic field patterns that are significantly more complex than those recovered for BP Tau, a fully convective T Tauri star.
We discuss possible reasons for this difference and suggest that the complexity of the stellar magnetic field is related to the convection zone; with more complex fields being found in T Tauri stars with radiative cores (V2129 Oph, CV Cha and CR Cha). However, it is clearly necessary to conduct magnetic field studies of T Tauri star systems, exploring a wide range of stellar parameters in order to establish how they affect magnetic field generation, and thus how these magnetic fields are likely to affect the evolution of T Tauri star systems as they approach the main sequence.     

The role of massive AGB stars in the early solar system composition

Abstract— We demonstrate that a massive asymptotic giant branch (AGB) star is a good candidate as the main source of short‐lived radionuclides in the early solar system. Recent identification of massive (4–8 M⊙) AGB stars in the galaxy, which are both lithium‐ and rubidium‐rich, demonstrates that these stars experience proton captures at the base of the convective envelope (hot bottom burning), together with high‐neutron density nucleosynthesis with ²²Ne as a neutron source in the He shell and efficient dredge‐up of the processed material. A model of a 6.5 M⊙ star of solar metallicity can simultaneously match the abundances of ²⁶Al, ⁴¹Ca, ⁶⁰Fe, and ¹⁰⁷Pd inferred to have been present in the solar nebula by using a dilution factor of 1 part of AGB material per 300 parts of original solar nebula material, and taking into account a time interval between injection of the short‐lived nuclides and consolidation of the first meteorites equal to 0.53 Myr. Such a polluting source does not overproduce ⁵³Mn, as supernova models do, and only marginally affects isotopic ratios of stable elements. It is usually argued that it is unlikely that the short‐lived radionuclides in the early solar system came from an AGB star because these stars are rarely found in star forming regions, however, we think that further interdisciplinary studies are needed to address the fundamental problem of the birth of our solar system.     

C-, Sr-isotope and Hg chemostratigraphy of Neoproterozoic cap carbonates of the Sergipano Belt,Northeastern Brazil

Two cap carbonates overlying glaciogenic diamictites crop out extensively in the eastern Vaza Barris Domain of the Sergipano Belt, northeastern Brazil. They are represented by carbonates of the Jacoca Formation, resting on top of diamictite of the Ribeiropolis Formation, and by the Olhos D’Agua Formation (carbonates, organic-rich towards the top), which overlies diamictite of the Palestina Formation. These two sequences were deformed and metamorphosed at sub-greenschist facies-conditions during the Brasiliano cycle (650–600 Ma).     

Tsunamis on Mars: Earth analogues of projected Martian sediment

Bolide impacts on Mars, within the proposed ocean boundaries (“contacts 1 and 2”) in the northern lowlands, would certainly have generated ultra high energy waves similar to tsunamis on Earth. Impacts into putative Noachian and Hesperian seas of variable areal extents and depths would have experienced high-energy inundations (transgressions), which would have left an imprint in the stack of deposits adjacent to the proposed shorelines. On Earth, the principal influencing factors for tsunami-wave energy are the character of shoreline topography and coastal water depth, which control wave compression and shoreline friction. Shorelines with narrow embayments and steep offshore gradients produce wave compression and increased collision of grains within the carried load contrasted with linear shorelines and shallow offshore gradients that dissipate energy. Steep offshore gradients produce concentrated major wave friction with the bed engendering high kinetic energy in the wave during emplacement of tsunami-generated sediment, which differs from shallow offshore beds that produce lower frictional effects over a wider area and drawdown of wave energy. Thus, overprinting of transported quartz grains on Earth is greatest where wave energy is highest, attenuated down to minor or nil overprinting where wave energy is less. Such grain overprinting in the form of energy-induced microtextures would also be observed in other grain types such as olivine and plagioclase, as such mineralogies are expected to dominate the Martian landscape based on orbital and local field (lander and rover) perspectives. Kinetic energy variation in tsunamis is controlled more by the square of velocity than mass, the resulting collisional effects of which produce swarms of v-shaped percussion microfeatures on quartz and other silicate mineral surfaces when velocity and compression are highest. This work indicates that a valid test for the ocean hypothesis is targeting “coastal” areas adjacent to narrow embayments where offshore depths are known to be highest, as possible tsunami-emplaced sediments, especially those that have been protected from atmospheric conditions through relatively rapid burial, may reveal a high frequency of percussion cracks, features of which appear to be unique to such terrestrial environments.     

Daily scale modelling of aquifer–river connectivity in the urban alluvial aquifer in Langreo, Spain

Evaluating aquifer–river interactions is naturally complex, particularly within urban settings. This is largely due to the difficulties involved in quantifying most elements of the water balance. The ability of numerical models to deal with several dynamic variables simultaneously makes them valuable tools to address this kind of problem. An applied, modeling-based approach to investigate the spatial and temporal variations of aquifer–river connectivity within a shallow urban aquifer is presented. Model development is based on comprehensive field campaigns in Langreo, Spain. Two calibration runs (for summer and winter conditions) were carried out in order to evaluate the spatial distribution of recharge rates. The model suggests that baseflows are largely negligible in comparison with total streamflows. This is mostly attributed to the abrupt nature of the catchment, which prevents the existence of sufficiently large alluvial systems to a great extent. Modelling results also show that aquifer–river connectivity at the study site is constrained by urban pumping as well as by seasonal fluctuations.     

The Sierra Ballena Shear Zone in the southernmost Dom Feliciano Belt (Uruguay): evolution, kinematics, and deformation conditions

The Sierra Ballena Shear Zone (SBSZ) is part of a high-strain transcurrent system that divides the Neoproterozoic Dom Feliciano Belt of South America into two different domains. The basement on both sides of the SBSZ shows a deformation stage preceding that of the transcurrent deformation recognized as a high temperature mylonitic foliation associated with migmatization. Grain boundary migration and fluid-assisted grain boundary diffusion enhanced by partial melting were the main deformation mechanisms associated with this foliation. Age estimate of this episode is >658 Ma. The second stage corresponds to the start of transpressional deformation and the nucleation and development of the SBSZ. During this stage, pure shear dominates the deformation, and is characterized by the development of conjugate dextral and sinistral shear zones and the emplacement of syntectonic granites. This event dates to 658–600 Ma based on the age of these intrusions. The third stage was a second transpressional event at about 586 to <560 Ma that was associated with the emplacement of porphyry dikes and granites that show evidence of flattening. Deformation in the SBSZ took place, during the late stages, under regional low-grade conditions, as indicated by the metamorphic paragenesis in the supracrustals of the country rocks. Granitic mylonites show plastic deformation of quartz and brittle behavior of feldspar. A transition from magmatic to solid-state microstructures is also frequently observed in syntectonic granites. Mylonitic porphyries and quartz mylonites resulted from the deformation of alkaline porphyries and quartz veins emplaced in the shear zone. Quartz veins reflect the release of silica associated with the breakdown of feldspar to white mica during the evolution of the granitic mylonites to phyllonites, which resulted in shear zone weakening. Quartz microstructures characteristic of the transition between regime 2 and regime 3, grain boundary migration and incipient recrystallization in feldspar indicate deformation under lower amphibolite to upper greenschist conditions (550–400°C). On the other hand, the mylonitic porphyries display evidence of feldspar recrystallization suggesting magmatic or high-T solid-state deformation during cooling of the dikes.     

Geo-environmental mapping using physiographic analysis: constraints on the evaluation of land instability and groundwater pollution hazards in the Metropolitan District of Campinas,Brazil

Geo-environmental terrain assessments and territorial zoning are useful tools for the formulation and implementation of environmental management instruments (including policy-making, planning, and enforcement of statutory regulations). They usually involve a set of procedures and techniques for delimitation, characterisation and classification of terrain units. However, terrain assessments and zoning exercises are often costly and time-consuming, particularly when encompassing large areas, which in many cases prevent local agencies in developing countries from properly benefiting from such assessments. In the present paper, a low-cost technique based on the analysis of texture of satellite imagery was used for delimitation of terrain units. The delimited units were further analysed in two test areas situated in Southeast Brazil to provide estimates of land instability and the vulnerability of groundwater to pollution hazards. The implementation incorporated procedures for inferring the influences and potential implications of tectonic fractures and other discontinuities on ground behaviour and local groundwater flow. Terrain attributes such as degree of fracturing, bedrock lithology and weathered materials were explored as indicators of ground properties. The paper also discusses constraints on- and limitations of- the approaches taken.     

High-pressure thermal aureoles around two Neoproterozoic synorogenic magmatic epidote-bearing granitoids,Northeastern Brazil

Unusual high-pressure inner thermal aureoles are described from the Minador and Angico Torto epidote-bearing tonalitic plutons that emplaced into greenschist-facies metasedimentary rocks of the Neoproterozoic Cachoeirinha-Salgueiro belt, northeastern Brazil. The foliated pelitic hornfelses display the mineral assemblage garnet, kyanite, staurolite, muscovite, biotite, plagioclase ± quartz. Rare fibrolite is only found very close to the contacts. Hornfelses display steep mineral lineations and steeply-dipping foliations concordant with magmatic contacts. Leucocratic veinlets containing quartz, oligoclase, garnet, kyanite, staurolite, rutile and ilmenite suggest that limited melting conditions were reached very close to magmatic contacts (T ? 650 °C, P around 8 kbar). These high-pressure hornfelses form a few meters thick, rigid envelopes around the two plutons. Contrary to known examples of kyanite-bearing hornfelses that recorded high-temperature decompression, the nearly isobaric cooling down to ca. 450 °C is constrained by 3.20–3.30 Si contents of retrogressive phengites from both inner hornfelses and ductilely-deformed tonalite at the pluton margins. Isograds and bathograds are, therefore, apparently telescoped due to HP/LT shearing, possibly caused by subsequent differential vertical movements affecting these two solidified plutons. The unusual depth of emplacement of these syn-kinematic calc-alkaline plutons is explained by a tentative geodynamic model involving a pre-620 Ma-subduction setting.ResumenLas aureolas internas que rodean dos plutones tonalíticos emplazados dentro de rocas cajas en facies esquistos verdes del Cinturón-plegado Cachoeirinha-Salgueiro al noreste de Brasil, contienen hornfelses pelíticos foliados con granate, kyanita, estaurolita, muscovita, biotita, plagioclasa ± cuarzo. Fibrolita es rara ó es encontrada solamente cerca de las zonas de contacto. Los hornfelses desarrollaron foliaciones concordantes con buzamiento fuerte cerca de los contactos magmáticos y muestran lineaciones minerales casi verticales. Venillas leucocraticas contienen cuarzo, oligoclasa, granate, kyanita, estaurolita, rutilo e ilmenita sugiriendo que el límite en las condiciones de fusión fueron cercanas a los contactos magmáticos (T ? 650 °C, P entorno 8 kbar). Estos hornfelses de alta presión forman una cobertura rígida de poco metros entorno a los dos plutones. Contrariamente a lo conocido en diferentes ejemplos de hornfelses con kyanita con registro de descompresión, la isóbara de enfriamiento para cerca de 450 °C está controlada por el contenido de Si 3.20–3.30 de fengitas retrógradas, en los hornfelses internos y en las márgenes de los plutones donde las tonalitas están deformadas ductilmente. Isógradas y batógradas son, por lo tanto, aparentemente telescopadas debido a cizalla de alta presión/baja temperatura, posiblemente causado por subsecuentes movimientos verticales diferenciales afectando estos dos plutones solidificados. El emplazamiento a una profundidad inusual de estos plutones calci-alcalinos syn-cinemáticos es explicado por un modelo geodinámico tentativo que involucra un escenario de subducción pre-620 Ma.     

A special perturbation method in orbital dynamics

The special perturbation method considered in this paper combines simplicity of computer implementation, speed and precision, and can propagate the orbit of any material particle. The paper describes the evolution of some orbital elements based in Euler parameters, which are constants in the unperturbed problem, but which evolve in the time scale imposed by the perturbation. The variation of parameters technique is used to develop expressions for the derivatives of seven elements for the general case, which includes any type of perturbation. These basic differential equations are slightly modified by introducing one additional equation for the time, reaching a total order of eight. The method was developed in the Grupo de Dinámica de Tethers (GDT) of the UPM, as a tool for dynamic simulations of tethers. However, it can be used in any other field and with any kind of orbit and perturbation. It is free of singularities related to small inclination and/or eccentricity. The use of Euler parameters makes it robust. The perturbation forces are handled in a very simple way: the method requires their components in the orbital frame or in an inertial frame. A comparison with other schemes is performed in the paper to show the good performance of the method.