Bow shocks and related density enhancements produced by the winds of massive stars moving through the interstellar medium provide important information regarding the motions of the stars, the properties of their stellar winds, and the characteristics of the local medium. Since bow-shock nebulae are aspherical structures, light scattering within them produces a net polarization signal even if the region is spatially unresolved. Scattering opacity arising from free electrons and dust leads to a distribution of polarized intensity across the bow-shock structure. That polarization encodes information about the shape, composition, opacity, density, and ionization state of the material within the structure. In this paper, we use the Monte Carlo radiative transfer code SLIP to investigate the polarization created when photons scatter in a bow-shock-shaped region of enhanced density surrounding a stellar source. We present results for electron scattering, and investigate the polarization behaviour as a function of optical depth, temperature, and source of photons for two different cases: pure scattering and scattering with absorption. In both regimes, we consider resolved and unresolved cases. We discuss the implications of these results as well as their possible use along with observational data to constrain the properties of observed bow-shock systems. In different situations and under certain assumptions, our simulations can constrain viewing angle, optical depth and temperature of the

Non-thermal radiation has been predicted within bow shocks around runaway stars by recent theoretical works. We present X-ray observations toward the runaway stars ζ Oph by Chandra and Suzaku and of BD+43°3654 by XMM-Newton to search for the presence of non-thermal X-ray emission. We found no evidence of non-thermal emission spatially coincident with the bow shocks; nonetheless, diffuse emission was detected in the vicinity of ζ Oph. After a careful analysis of its spectral characteristics, we conclude that this emission has a thermal nature with a plasma temperature of T ≈ 2 × 106 K. The cometary shape of this emission seems to be in line with recent predictions of radiation-hydrodynamic models of runaway stars. The case of BD+43°3654 is puzzling, as non-thermal emission has been reported in a previous work for this source.

Observations with powerful X-ray telescopes, such as XMM-Newton and Chandra, significantly advance our understanding of massive stars. Nearly all early-type stars are X-ray sources. Studies of their X-ray emission provide important diagnostics of stellar winds. High-resolution X-ray spectra of O-type stars are well explained when stellar wind clumping is taking into account, providing further support to a modern picture of stellar winds as non-stationary, inhomogeneous outflows. X-ray variability is detected from such winds, on time scales likely associated with stellar rotation. High-resolution X-ray spectroscopy indicates that the winds of late O-type stars are predominantly in a hot phase. Consequently, X-rays provide the best observational window to study these winds. X-ray spectroscopy of evolved, Wolf-Rayet type, stars allows to probe their powerful metal enhanced winds, while the mechanisms responsible for the X-ray emission of these stars are not yet understood.

Recombination line profile shapes are derived for ionized spherical stellar winds at radio wavelengths. It is assumed that the wind is optically thick owing to free-free opacity. Emission lines of arbitrary optical depth are obtained assuming that the free-free photosphere forms in the outer, constant expansion portion of the wind. Previous works have derived analytic results for isothermal winds when the line and continuum source functions are equal. Here, semi-analytic results are derived for unequal source functions to reveal that line shapes can be asymmetric about line center. A parameter study is presented and applications discussed.

We present new, contemporaneous Hubble Space Telescope STIS and XMM-Newton observations of the O7 III(n)((f)) star ξ Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of ξ Per. The variable wind of this star is known to have a 2.086-day periodicity. We use a simple, heuristic spot model that fits the low-velocity (near-surface) IUE wind line variability very well, to demonstrate that the low-velocity absorption in the new STIS spectra of N iv λ1718 and Si iv λ1402 vary with the same 2.086-day period. It is remarkable that the period and amplitude of the STIS data agree with those of the IUE spectra obtained 22 yr earlier. We also show that the time variability of the new XMM-Newton fluxes is also consistent with the 2.086-day period. Thus, our new, multiwavelength coordinated observations demonstrate that the mechanism that causes the UV wind line variability is also responsible for a significant fraction of the X-rays in single O stars. The sequence of events for the multiwavelength light-curve minima is Si iv λ1402, N iv λ1718, and X-ray flux, each separated by a phase of about 0.06 relative to the 2.086-day period. Analysis of the X-ray fluxes shows that they become softer as they weaken. This is contrary to expectations if the variability is caused by periodic excess absorption. Furthermore, the high-resolution X-ray spectra suggest that the individual emission lines at maximum are more strongly blueshifted. If we interpret the low-velocity wind line light curves in terms of our model, it implies that there are two bright regions, i.e., regions with less absorption, separated by 180°, on the surface of the star. We note that the presence and persistence of two spots separated by 180° suggest that a weak dipole magnetic field is responsible for the variability of the UV wind line absorption and X-ray flux in ξ Per.

Les surcotes de tempête sont souvent sous-estimées dans les modèles hydrodynamiques, ainsi que les grandes vagues dans les modèles de vagues. Les causes possibles sont une sous-estimation des vents dans les modèles atmosphériques et/ou une formulation incorrecte de la tension de vent. Les objectifs de cette thèse sont (1) d’estimer les biais par vents forts dans les modèles atmosphériques (2) de développer une nouvelle paramétrisation du coefficient de traı̂née permettant de réduire ce biais (3) d’étudier l’impact des vagues sur la tension de vent. La méthode consiste à étudier la réponse de l’atmosphère et de l’océan à la tension de vent. Dans une première partie, nous utilisons le modèle couplé vagues-atmosphère d’ECMWF. Nous montrons que les vents forts sont sous-estimés, avec un biais de l’ordre de -7 m/s à 30 m/s. Des écarts significatifs existent aussi entre les observations, les bouées et les vents issus de ASCAT-KNMI étant généralement inférieurs à ceux des plateformes et des autres données satellites utilisées dans cette étude (AMSR2, ASCAT-RSS, WindSat, SMOS et JASON-2). La nouvelle paramétrisation développée permet d’obtenir des vents plus forts qu’avec celle d’ECMWF par défaut. Dans une deuxième partie (réponse de l’océan), nous utilisons le modèle global océanique TUGO du LEGOS forcé par le modèle couplé vagues-atmosphère d’ECMWF. Nous montrons qu’une paramétrisation de la tension de vent dépendant des vagues plutôt que du vent est plus appropriée quand l’état de mer est jeune. Elle conduit à des surcotes plus proches des observations (marégraphes et traces altimétriques de JASON-2). L’impact des vagues sur la surcote est significatif, et peut atteindre 20 cm. / Storm surges may be underestimated in hydrodynamic models, as well as large wave heights in wave models. This could come from an underestimation of strong winds in atmospheric models and/or an inappropriate wind stress formulation. The objectives of the present work are (1) to estimate how strong are the biases for high winds in atmospheric models (2) to develop a new drag parameterization that could reduce this bias (3) to investigate the impact of the waves on the wind stress. The method consists of studying the response of the atmosphere and the ocean to the wind stress.In a first part, we use the coupled wave-atmosphere model from ECMWF. We show that strong winds may be underestimated, as much as -7 m/s at 30 m/s.Significant differences also exist between observations, with buoys and ASCAT-KNMI generally showing lower wind speeds than the platforms and other remote-sensing data used in this study(AMSR2, ASCAT-RSS, WindSat, SMOS and JASON-2).The newly empirically adjusted Charnock parameterization leads to higher winds compared to the default ECMWF parameterization. In a second part, we use the global ocean model TUGO fromLEGOS forced with ECMWF coupled wave-atmopshere model. We show that a wave-dependent rather than wind-dependent stress formulation is more appropriate, when the sea state is young and the sea rougher. It yields to simulated surges closer to observations (i.e. tide gauges and JASON-2 altimeter tracks). The wave impact on the surges is significant, and may reach 20 cm.

Thesis: Ph. D., Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2019 / Cataloged from PDF version of thesis. / Includes bibliographical references (pages 113-121). / Submesoscale ocean dynamics and instabilities, with characteristic scales 0.1-10 kin, can play a critical role in setting the ocean's surface boundary layer thickness and associated density stratification. Submesoscale instabilities contribute to lateral stirring and tracer dispersal. These dynamics are investigated in the Bay of Bengal, motivated by the upper ocean's potentially coupled interactions with Monsoon winds and convection. The region's excess precipitation and runoff generates strong salinity gradients that typically set density fronts and stratification in the upper 50 m. Since we cannot synoptically measure currents containing fast-evolving and oscillating components across the submesoscale range, we instead analyze passive tracer distributions (spice = density-compensated temperature (T) and salinity (S) anomalies), identifying signatures of flows and testing dynamical theories. / The analysis is based on over 9000 vertical profiles of T and S measured along ~4800 km of ship tracks in the Bay of Bengal during ASIRI and MISO-BOB expeditions in 2013, 2015, and 2018. Observations in the surface mixed layer reveal ~1 km scale-selective correlation of surface T and S, with compensation reducing cross-front density gradients by ~50%. Using a process study ocean model, we show this is caused by submesoscale instabilities slumping fronts, plus surface cooling over the resultant enhanced salinity stratification, potentially thwarting the forward cascade of energy. In the stratified interior, we present a spectral analysis of horizontal spice variance statistics from wavenumber k ~0.01 cpkm to ~1 cpkm. At scales <10 km, stratified layers that are closer to the surface exhibit redder passive tracer spectra (power spectra k⁻³, gradient spectra k⁻¹) than predicted by quasi-geostrophic or frontogenetic theories. / Complimentary observations reveal spice patterns with multiple, parallel, ~10 m thin layers, crossing isopycnals with O(10⁻⁴) slopes, coherent over at least 30-80 kin, with coincident layers of stratification anomalies. Comparison with shear measurements, and a numerical process study, suggest that both submesoscale sheared eddies, and thin near-inertial waves, form such layers. Fast formation timescales and large aspect ratios suggest they enhance horizontal mixing by shear dispersion, reducing variance at ~1-10 km scales. / by Gualtiero Victor Rudi Spiro Jaeger. / Ph. D. / Ph.D. Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution)

A time dependent model of the shallow, tropical convective boundary layer is developed and tested. To simplify the treatment of thermodynamic processes and return to first principles of physics, conserved quantities of entropy and total water density are used as primary model variables. In addition, a new shallow cumulus parameterization scheme is developed and is based on the use of a time dependent cloud kinetic energy equation combining local concepts of cloud processes with the use of a special buoyancy length scale. Two model simulations are performed in an attempt to assess the model's performance and the effectiveness of the parameterization scheme. Results indicate that the model does a reasonable job in both representing the equilibrium structure of a shallow convective boundary layer and in generating a realistic boundary layer structure from an initial state consisting of a shallow moist layer with dry air aloft. The cumulus parameterization scheme appears to adequately represent the transport of thermodynamic quantities associated with convective activity and the use of conserved variables provides an effective way of representing the boundary layer structure and treating the mixing processes associated with cloud processes. This work illustrates the usefulness of generalized conserved variables, particularly entropy and total water density, and indicates that the general approach of using a time dependent cloud kinetic energy equation may be effective for representing thermodynamic processes in the tropical boundary layer. / Graduation date: 1992

Stellar winds play an important role for the evolution of massive stars and their cosmic environment. Multiple lines of evidence, coming from spectroscopy, polarimetry, variability, stellar ejecta, and hydrodynamic modeling, suggest that stellar winds are non-stationary and inhomogeneous. This is referred to as 'wind clumping'. The urgent need to understand this phenomenon is boosted by its far-reaching implications. Most importantly, all techniques to derive empirical mass-loss rates are more or less corrupted by wind clumping. Consequently, mass-loss rates are extremely uncertain. Within their range of uncertainty, completely different scenarios for the evolution of massive stars are obtained. Settling these questions for Galactic OB, LBV and Wolf-Rayet stars is prerequisite to understanding stellar clusters and galaxies, or predicting the properties of first-generation stars. In order to develop a consistent picture and understanding of clumped stellar winds, an international workshop on 'Clumping in Hot Star Winds' was held in Potsdam, Germany, from 18. - 22. June 2007. About 60 participants, comprising almost all leading experts in the field, gathered for one week of extensive exchange and discussion. The Scientific Organizing Committee (SOC) included John Brown (Glasgow), Joseph Cassinelli (Madison), Paul Crowther (Sheffield), Alex Fullerton (Baltimore), Wolf-Rainer Hamann (Potsdam, chair), Anthony Moffat (Montreal), Stan Owocki (Newark), and Joachim Puls (Munich). These proceedings contain the invited and contributed talks presented at the workshop, and document the extensive discussions.

The evolutionary paths of the massive O and B type stars are often defined by angular momentum transformations that involve circumstellar gas disks. This circumstellar gas is revealed in several kinds of observations, and here I describe a series of investigations of the hydrogen line emission from such disk using detailed studies of five massive binaries and a survey of 128 Be stars. By examining three sets of spectra of the active mass-transfer binary system RY Scuti, I determined masses of 7.1 +/- 1.2 M_sun for the bright supergiant and 30.0 +/- 2.1 M_sun for the massive companion that is hidden by an accretion torus. I also present a cartoon model of the complex mass flows in the system. Using optical spectroscopy and X-ray flux data, I investigated the mass transfer processes in four massive X-ray binaries (a massive B star with mass flowing onto a compact, neutron star companion). The B-supergiant system LS I +65 010 transfers mass via stellar winds. I find the X-ray flux modulates with the orbital period. In the other three X-ray binary systems (LS I +61 303, HDE 245770, and X Persei), an outflowing circumstellar disk is responsible for the mass transfer, and in all three systems, the disk appears to be truncated by gravitational interactions with the compact companion. The disk in the microquasar system LS I +61 303 is limited in radius by the periastron separation and an increase in both H-alpha equivalent width and X-ray flux following periastron may be due to a density wave in the disk induced by tidal forces. Observations of HDE 245770 document what appears to be the regeneration of a circumstellar disk. The disk of X Persei appears to have grown to near record proportions and the X-ray flux has dramatically increased. Tidal interaction may generate a spiral density wave in the disk and cause an increase in H-alpha equivalent width and mass transfer to the compact companion. During the course of the analysis of the X-ray binaries, I developed numerical models for estimating the size of the Be star disks using just the H-alpha equivalent width. Finally, I present the results of a three year spectroscopic survey of both the H-alpha and H-gamma regions of 128 Be stars. I find that the median fractional variation in the equivalent width of the disk emission lines is 15% over a two year period. I also find that two-thirds of the sample displays evidence of Fe II emission or absorption resulting from surrounding circumstellar material. Many candidates for non-radial pulsation and binary systems are also found. Spectra and notes for all of the sample stars are presented in an appendix.