Results of the pollen analysis of the Dyanushka peat sediments K7/P2

A 415cm thick permafrost peat section from the Verkhoyansk Mountains was radiocarbon-dated and studied using palaeobotanical and sedimentological approaches. Accumulation of organic-rich sediment commenced in a former oxbow lake, detached from a Dyanushka River meander during the Younger Dryas stadial, at ~12.5 kyr BP. Pollen data indicate that larch trees, shrub alder and dwarf birch were abundant in the vegetation at that time. Local presence of larch during the Younger Dryas is documented by well-preserved and radiocarbon-dated needles and cones. The early Holocene pollen assemblages reveal high percentages of Artemisia pollen, suggesting the presence of steppe-like communities around the site, possibly in response to a relatively warm and dry climate ~11.4-11.2 kyr BP. Both pollen and plant macrofossil data demonstrate that larch woods were common in the river valley. Remains of charcoal and pollen of Epilobium indicate fire events and mark a hiatus ~11.0-8.7 kyr BP. Changes in peat properties, C31/C27 alkane ratios and radiocarbon dates suggest that two other hiatuses occurred ~8.2-6.9 and ~6.7-0.6 kyr BP. Prior to 0.6 kyr BP, a major fire destroyed the mire surface. The upper 60 cm of the studied section is composed of aeolian sands modified in the uppermost part by the modern soil formation. For the first time, local growth of larch during the Younger Dryas has been verified in the western foreland of the Verkhoyansk Mountains (~170km south of the Arctic Circle), thus increasing our understanding of the quick reforestation of northern Eurasia by the early Holocene.

Pollen records from Bol'shoy Lyakhovsky Island, Siberia

Cryolithological, ground ice and fossil bioindicator (pollen, diatoms, plant macrofossils, rhizopods, insects, mammal bones) records from Bol'shoy Lyakhovsky Island permafrost sequences (73°20'N, 141°30'E) document the environmental history in the region for the past c. 115 kyr. Vegetation similar to modern subarctic tundra communities prevailed during the Eemian/Early Weichselian transition with a climate warmer than the present. Sparse tundra-like vegetation and harsher climate conditions were predominant during the Early Weichselian. The Middle Weichselian deposits contain peat and peaty soil horizons with bioindicators documenting climate amelioration. Although dwarf willows grew in more protected places, tundra and steppe vegetation prevailed. Climate conditions became colder and drier c. 30 kyr BP. No sediments dated between c. 28.5 and 12.05 14C kyr BP were found, which may reflect active erosion during that time. Herb and shrubby vegetation were predominant 11.6-11.3 14C kyr BP. Summer temperatures were c. 4 °C higher than today. Typical arctic environments prevailed around 10.5 14C kyr BP. Shrub alder and dwarf birch tundra were predominant between c. 9 and 7.6 kyr BP. Reconstructed summer temperatures were at least 4 °C higher than present. However, insect remains reflect that steppe-like habitats existed until c. 8 kyr BP. After 7.6 kyr BP, shrubs gradually disappeared and the vegetation cover became similar to that of modern tundra. Pollen and beetles indicate a severe arctic environment c. 3.7 kyr BP. However, Betula nana, absent on the island today, was still present. Together with our previous study on Bol'shoy Lyakhovsky Island covering the period between about 200 and 115 kyr, a comprehensive terrestrial palaeoenvironmental data set from this area in western Beringia is now available for the past two glacial-interglacial cycles.

Type Ia supernovae from violent mergers of carbon-oxygen white dwarfs: polarization signatures

The violent merger of two carbon-oxygen white dwarfs has been proposed as a viable progenitor for some Type Ia supernovae. However, it has been argued that the strong ejecta asymmetries produced by this model might be inconsistent with the low degree of polarization typically observed in Type Ia supernova explosions. Here, we test this claim by carrying out a spectropolarimetric analysis for the model proposed by Pakmor et al. for an explosion triggered during the merger of a 1.1 and 0.9 M⊙ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarization signal varies significantly with viewing angle. We find that polarization levels for observers in the equatorial plane are modest (≲1 per cent) and show clear evidence for a dominant axis, as a consequence of the ejecta symmetry about the orbital plane. In contrast, orientations out of the plane are associated with higher degrees of polarization and departures from a dominant axis. While the particular model studied here gives a good match to highly polarized events such as SN 2004dt, it has difficulties in reproducing the low polarization levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarization features that are not observed in the majority of currently available spectropolarimetric data of Type Ia supernovae. Future studies will map out the parameter space of the merger scenario to investigate if alternative models can provide better agreement with observations.

Calcium-rich gap transients: Tidal detonations of white dwarfs?

We hypothesize that at least some of the recently discovered class of calcium-rich gap transients are tidal detonation events of white dwarfs (WDs) by black holes (BHs) or possibly neutron stars. We show that the properties of the calcium-rich gap transients agree well with the predictions of the tidal detonation model. Under the predictions of this model, we use a follow-up X-ray observation of one of these transients, SN 2012hn, to place weak upper limits on the detonator mass of this system that include all intermediate-mass BHs (IMBHs). As these transients are preferentially in the stellar haloes of galaxies, we discuss the possibility that these transients are tidal detonations of WDs caused by random flyby encounters with IMBHs in dwarf galaxies or globular clusters. This possibility has been already suggested in the literature but without connection to the calcium-rich gap transients. In order for the random flyby cross-section to be high enough, these events would have to be occurring inside these dense stellar associations. However, there is a lack of evidence for IMBHs in these systems, and recent observations have ruled out all but the very faintest dwarf galaxies and globular clusters for a few of these transients. Another possibility is that these are tidal detonations caused by three-body interactions, where a WD is perturbed towards the detonator in isolated multiple star systems. We highlight a number of ways this could occur, even in lower mass systems with stellar-mass BHs or neutron stars. Finally, we outline several new observational tests of this scenario, which are feasible with current instrumentation.

The crowded magnetosphere of the post-common-envelope binary QS Virginis

We present high-speed photometry and high-resolution spectroscopy of the eclipsing post-common-envelope binary QS Virginis (QS Vir). Our Ultraviolet and Visual Echelle Spectrograph (UVES) spectra span multiple orbits over more than a year and reveal the presence of several large prominences passing in front of both the M star and its white dwarf companion, allowing us to triangulate their positions. Despite showing small variations on a time-scale of days, they persist for more than a year and may last decades. One large prominence extends almost three stellar radii from the M star. Roche tomography reveals that the M star is heavily spotted and that these spots are long-lived and in relatively fixed locations, preferentially found on the hemisphere facing the white dwarf. We also determine precise binary and physical parameters for the system. We find that the 14 220 ± 350 K white dwarf is relatively massive, 0.782 ± 0.013 M⊙, and has a radius of 0.010 68 ± 0.000 07 R⊙, consistent with evolutionary models. The tidally distorted M star has a mass of 0.382 ± 0.006 M⊙ and a radius of 0.381 ± 0.003 R⊙, also consistent with evolutionary models. We find that the magnesium absorption line from the white dwarf is broader than expected. This could be due to rotation (implying a spin period of only ˜700 s), or due to a weak (˜100 kG) magnetic field, we favour the latter interpretation. Since the M star's radius is still within its Roche lobe and there is no evidence that it is overinflated, we conclude that QS Vir is most likely a pre-cataclysmic binary just about to become semidetached.

The O stars in the VLT-FLAMES Tarantula Survey

The VLT-FLAMES Tarantula Survey (VFTS) has secured mid-resolution spectra of over 300 O-type stars in the 30 Doradus region of the Large Magellanic Cloud. A homogeneous analysis of such a large sample requires automated techniques, an approach that will also be needed for the upcoming analysis of the Gaia surveys of the Northern and Southern Hemisphere supplementing the Gaia measurements. We point out the importance of Gaia for the study of O stars, summarize the O star science case of VFTS and present a test of the automated modeling technique using synthetically generated data. This method employs a genetic algorithm based optimization technique in combination with fastwind model atmospheres. The method is found to be robust and able to recover the main photospheric parameters accurately. Precise wind parameters can be obtained as well, however, as expected, for dwarf stars the rate of acceleration of the ow is poorly constrained.

Searching for swept-up hydrogen and helium in the late-time spectra of 11 nearby Type Ia supernovae

The direct detection of a stellar system that explodes as a Type Ia supernova (SN Ia) has not yet been successful. Various indirect methods have been used to investigate SN Ia progenitor systems but none have produced conclusive results. A prediction of single-degenerate models is that H- (or He-) rich material from the envelope of the companion star should be swept up by the SN ejecta in the explosion. Seven SNe Ia have been analysed to date looking for signs of H-rich material in their late-time spectra and none were detected. We present results from new late-time spectra of 11 SNe Ia obtained at the Very Large Telescope using XShooter and FORS2. We present the tentative detection of Hα emission for SN 2013ct, corresponding to ∼0.007 M_⊙ of stripped/ablated companion star material (under the assumptions of the spectral modelling). This mass is significantly lower than expected for single-degenerate scenarios, suggesting that >0.1 M_⊙ of H-rich is present but not observed. We do not detect Hα emission in the other 10 SNe Ia. This brings the total sample of normal SNe Ia with non-detections (<0.001–0.058 M_⊙) of H-rich material to 17 events. The simplest explanation for these non-detections is that these objects did not result from the explosion of a CO white dwarf accreting matter from a H-rich companion star via Roche lobe overflow or symbiotic channels. However, further spectral modelling is needed to confirm this. We also find no evidence of He-emission features, but models with He-rich companion stars are not available to place mass limits.

Predicting polarization signatures for double-detonation and delayed-detonation models of Type Ia supernovae

Calculations of synthetic spectropolarimetry are one means to test multidimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M⊙ white dwarf binary system is too asymmetric to explain the low polarization levels commonly observed in normal Type Ia supernovae. Here, we present polarization simulations for two alternative scenarios: the sub-Chandrasekhar mass double-detonation and the Chandrasekhar mass delayed-detonation model. Specifically, we study a 2D double-detonation model and a 3D delayed-detonation model, and calculate polarization spectra for multiple observer orientations in both cases. We find modest polarization levels (<1 per cent) for both explosion models. Polarization in the continuum peaks at ∼0.1–0.3 per cent and decreases after maximum light, in excellent agreement with spectropolarimetric data of normal Type Ia supernovae. Higher degrees of polarization are found across individual spectral lines. In particular, the synthetic Si II λ6355 profiles are polarized at levels that match remarkably well the values observed in normal Type Ia supernovae, while the low degrees of polarization predicted across the O I λ7774 region are consistent with the non-detection of this feature in current data. We conclude that our models can reproduce many of the characteristics of both flux and polarization spectra for well-studied Type Ia supernovae, such as SN 2001el and SN 2012fr. However, the two models considered here cannot account for the unusually high level of polarization observed in extreme cases such as SN 2004dt.

Three-dimensional simulations of gravitationally confined detonations compared to observations of SN 1991T

The gravitationally confined detonation (GCD) model has been proposed as a possible explosion mechanism for Type Ia supernovae in the single-degenerate evolution channel. It starts with ignition of a deflagration in a single off-centre bubble in a near-Chandrasekhar-mass white dwarf. Driven by buoyancy, the deflagration flame rises in a narrow cone towards the surface. For the most part, the main component of the flow of the expanding ashes remains radial, but upon reaching the outer, low-pressure layers of the white dwarf, an additional lateral component develops. This causes the deflagration ashes to converge again at the opposite side, where the compression heats fuel and a detonation may be launched. We first performed five three-dimensional hydrodynamic simulations of the deflagration phase in 1.4 M⊙ carbon/oxygen white dwarfs at intermediate-resolution (2563computational zones). We confirm that the closer the initial deflagration is ignited to the centre, the slower the buoyant rise and the longer the deflagration ashes takes to break out and close in on the opposite pole to collide. To test the GCD explosion model, we then performed a high-resolution (5123 computational zones) simulation for a model with an ignition spot offset near the upper limit of what is still justifiable, 200 km. This high-resolution simulation met our deliberately optimistic detonation criteria, and we initiated a detonation. The detonation burned through the white dwarf and led to its complete disruption. For this model, we determined detailed nucleosynthetic yields by post-processing 106 tracer particles with a 384 nuclide reaction network, and we present multi-band light curves and time-dependent optical spectra. We find that our synthetic observables show a prominent viewing-angle sensitivity in ultraviolet and blue wavelength bands, which contradicts observed SNe Ia. The strong dependence on the viewing angle is caused by the asymmetric distribution of the deflagration ashes in the outer ejecta layers. Finally, we compared our model to SN 1991T. The overall flux level of the model is slightly too low, and the model predicts pre-maximum light spectral features due to Ca, S, and Si that are too strong. Furthermore, the model chemical abundance stratification qualitatively disagrees with recent abundance tomography results in two key areas: our model lacks low-velocity stable Fe and instead has copious amounts of high-velocity 56Ni and stable Fe. We therefore do not find good agreement of the model with SN 1991T.

Pulling Back the Veil: The Characterization and Habitability of Enshrouded Worlds

Thesis (Ph.D.)--University of Washington, 2016-08

Resolved Stars' Insights into Galaxy Physics

Thesis (Ph.D.)--University of Washington, 2016-08

Traitement de données observationnelles et développement d’un code de simulation pour SpIOMM et SITELLE

L’arrivée du spectromètre imageur à transformée de Fourier SITELLE au télescope Canada-France-Hawaï souligne la nécessité d’un calculateur de temps d’exposition permettant aux utilisateurs de l’instrument de planifier leurs observations et leurs demandes de temps de télescope. Une grande partie de mon projet est ainsi le développement d’un code de simulation capable de reproduire les résultats de SITELLE et de son prédecesseur SpIOMM, installé à l’Observatoire du Mont-Mégantic. La précision des simulations est confirmée par une comparaison avec des données SpIOMM et les premières observations de SITELLE. La seconde partie de mon projet consiste en une analyse spectrale de données observationelles. Prenant avantage du grand champ de vue de SpIOMM, les caractéristiques du gaz ionisé (vitesse radiale et intensité) sont étudiées pour l’ensemble de la paire de galaxies en interaction Arp 72. La courbe de rotation dans le visible ainsi que le gradient de métallicité de NGC 5996, la galaxie principale d’Arp 72, sont obtenues ici pour la première fois. La galaxie spirale NGC 7320 est également étudiée à partir d’observations faites à la fois avec SpIOMM et SITELLE.

Inhomogeneous cloud coverage through the Coulomb explosion of dust in substellar atmospheres

Context. Recent observations of brown dwarf spectroscopic variability in the infrared infer the presence of patchy cloud cover. Aims. This paper proposes a mechanism for producing inhomogeneous cloud coverage due to the depletion of cloud particles through the Coulomb explosion of dust in atmospheric plasma regions. Charged dust grains Coulomb-explode when the electrostatic stress of the grain exceeds its mechanical tensile stress, which results in grains below a critical radius a < a Coul crit being broken up. Methods. This work outlines the criteria required for the Coulomb explosion of dust clouds in substellar atmospheres, the effect on the dust particle size distribution function, and the resulting radiative properties of the atmospheric regions. Results. Our results show that for an atmospheric plasma region with an electron temperature of Te = 10 eV (≈105 K), the critical grain radius varies from 10−7 to 10−4 cm, depending on the grains’ tensile strength. Higher critical radii up to 10−3 cm are attainable for higher electron temperatures. We find that the process produces a bimodal particle size distribution composed of stable nanoscale seed particles and dust particles with a ≥ a Coul crit , with the intervening particle sizes defining a region devoid of dust. As a result, the dust population is depleted, and the clouds become optically thin in the wavelength range 0.1–10 μm, with a characteristic peak that shifts to higher wavelengths as more sub-micrometer particles are destroyed. Conclusions. In an atmosphere populated with a distribution of plasma volumes, this will yield regions of contrasting radiative properties, thereby giving a source of inhomogeneous cloud coverage. The results presented here may also be relevant for dust in supernova remnants and protoplanetary disks.

Ionisation and discharge in cloud-forming atmospheres of brown dwarfs and extrasolar planets

Brown dwarfs and giant gas extrasolar planets have cold atmospheres with rich chemical compositions from which mineral cloud particles form. Their properties, like particle sizes and material composition, vary with height, and the mineral cloud particles are charged due to triboelectric processes in such dynamic atmospheres. The dynamics of the atmospheric gas is driven by the irradiating host star and/or by the rotation of the objects that changes during its lifetime. Thermal gas ionisation in these ultra-cool but dense atmospheres allows electrostatic interactions and magnetic coupling of a substantial atmosphere volume. Combined with a strong magnetic field , a chromosphere and aurorae might form as suggested by radio and x-ray observations of brown dwarfs. Non-equilibrium processes like cosmic ray ionisation and discharge processes in clouds will increase the local pool of free electrons in the gas. Cosmic rays and lighting discharges also alter the composition of the local atmospheric gas such that tracer molecules might be identified. Cosmic rays affect the atmosphere through air showers in a certain volume which was modelled with a 3D Monte Carlo radiative transfer code to be able to visualise their spacial extent. Given a certain degree of thermal ionisation of the atmospheric gas, we suggest that electron attachment to charge mineral cloud particles is too inefficient to cause an electrostatic disruption of the cloud particles. Cloud particles will therefore not be destroyed by Coulomb explosion for the local temperature in the collisional dominated brown dwarf and giant gas planet atmospheres. However, the cloud particles are destroyed electrostatically in regions with strong gas ionisation. The potential size of such cloud holes would, however, be too small and might occur too far inside the cloud to mimic the effect of, e.g. magnetic field induced star spots.

Virus diseases of chickpeas and pulse crops in Australia

Accurate identification of viruses is critical for resistance breeding and for development of management strategies. To this end, we are developing PCR diagnostics for the luteoviruses / poleroviruses that commonly affect chickpea and pulse crops in Australia. This is helping to overcome the shortfalls in virus identifications that often result from cross reactions of viruses to some antibodies. We compared these PCR tests with antibody based Tissue blot immune-assay (TBIA) in virus surveys of chickpea and pulse crops from eastern Australia. We used a multiplex PCR for Beet western yellows virus (BWYV), Bean leaf roll virus (BLRV), Phasey bean virus (PhBV – a new polerovirus species) and Soybean dwarf virus (SbDV) to investigate the importance of each virus and their host range from different locations. Important alternative hosts included Malva parviflora which was commonly found to be infected with BWYV from many locations and Medicago polymorpha was a host for BLRV, PhBV and SbDV. Using the virus species-specific PCR, 49 virus affected plants (mostly crop plants) from surveys in 2013 were screened, revealing the following infections; 38 SbDV, 5 PhBV, 3 BWYV, 2 BLRV and 1 mixed SbDV/BWYV. From the 45 samples that were not BWYV by PCR, 33 were false-positives in the BWYV TBIA. This demonstrates the BWYV antibody used was not useful for identifying BWYV and PCR indicated that SbDV was the dominant virus from the samples tested from the 2013 season. Preliminary results from the 2014 season indicate a significant change, with SbDV being only a minor component of the total virus population. Further work to clarify the Australian luteovirus complex through molecular techniques is in progress.