Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Groundwater Variability Across Temporal and Spatial Scales in the Central and Northeastern U.S.

Depth-to-water measurements from 181 monitoring wells in unconfined or semi-confined aquifers in nine regions of the central and northeastern U.S. were analyzed. Groundwater storage exhibited strong seasonal variations in all regions, with peaks in spring and lows in autumn, and its interannual variability was nearly unbounded, such that the impacts of droughts, floods, and excessive pumping could persist for many years. We found that the spatial variability of groundwater storage anomalies (deviations from the long term mean) increases as a power function of extent scale (square root of area). That relationship, which is linear on a log-log graph, is common to other hydrological variables but had never before been shown with groundwater data. We describe how the derived power function can be used to determine the number of wells needed to estimate regional mean groundwater storage anomalies with a desired level of accuracy, or to assess uncertainty in regional mean estimates from a set number of observations. We found that the spatial variability of groundwater storage anomalies within a region often increases with the absolute value of the regional mean anomaly, the opposite of the relationship between soil moisture spatial variability and mean. Recharge (drainage from the lowest model soil layer) simulated by the Variable Infiltration Capacity (VIC) model was compatible with observed monthly groundwater storage anomalies and month-to-month changes in groundwater storage.

water cycle↗

Spatial and Temporal Variability of Carbon Dioxide Using Structure Functions in Urban Areas: Insights for Future Active Remote CO2 Sensors

High resolution in-situ CO2 measurements were recorded onboard the NASA P-3B during the DISCOVER-AQ (Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality) Field Campaigns during July 2011 over Washington DC/Baltimore, MD; January-February 2013 over the San Joaquin Valley, CA; September 2013 over Houston, TX; and July-August 2014 over Denver, CO. Each of these campaigns have approximately two hundred vertical soundings of CO2 within the lower troposphere (surface to about 5 kilometers) at 6-8 different sites in each of the urban areas. In this study, we used structure function analysis, which is a useful way to quantify spatial and temporal variability, by displaying differences with average observations, to evaluate the variability of CO2 in the 0-2 kilometers range (representative of the planetary boundary layer). These results can then be used to provide guidance in the development of science requirements for the future ASCENDS (Active Sensing of CO2 Emissions over Nights, Days, and Seasons) mission to measure near-surface CO2 variability in different urban areas. We also compare the observed in-situ CO2 variability with the variability of the CO2 column-averaged optical depths in the 0-1 kilometer and 0-3.5 kilometers altitude ranges in the four geographically different urban areas, using vertical weighting functions for potential future ASCENDS lidar CO2 sensors operating in the 1.57 and 2.05 millimeter measurement regions. In addition to determining the natural variability of CO2 near the surface and in the column, radiocarbon method using continuous CO2 and CO measurements are used to examine the variation of emission quantification between anthropogenic and biogenic sources in the DC/Maryland urban site.

Choi, Yonghoon↗

A42F-03: Small-Scale Variability in Tropical Tropopause Layer Humidity

Recent advances in statistical parameterizations of cirrus cloud processes for use in global models are highlighting the need for information about small-scale fluctuations in upper tropospheric humidity and the physical processes that control the humidity variability. To address these issues, we have analyzed high-resolution airborne water vapor measurements obtained in the Airborne Tropical TRopopause EXperiment (ATTREX) over the tropical Pacific between 14 and 20 km. Using accurate and precise 1-Hz water vapor measurements along approximately-level aircraft flight legs, we calculate structure functions spanning horizontal scales ranging from about 0.2 to 50 km, and we compare the water vapor variability in the lower (about 14 km) and upper (16-19 km) Tropical Tropopause Layer (TTL). We also compare the magnitudes and scales of variability inside TTL cirrus versus in clear-sky regions. The measurements show that in the upper TTL, water vapor concentration variance is stronger inside cirrus than in clear-sky regions. Using simulations of TTL cirrus formation, we show that small variability in clear-sky humidity is amplified by the strong sensitivity of ice nucleation rate to supersaturation, which results in highly-structured clouds that subsequently drive variability in the water vapor field. In the lower TTL, humidity variability is correlated with recent detrainment from deep convection. The structure functions indicate approximately power-law scaling with spectral slopes ranging from about minus 5 divided by 3, to minus 2.

water vapor↗

Significant Spatial Variability in Radar-Derived West Antarctic Accumulation Linked to Surface Winds and Topography

Across the Antarctic Ice Sheet, accumulation heavily influences firn compaction and surface height changes. Therefore, accumulation varies over short distances (<25 km), complicating the derivation of ice sheet mass changes from altimetry and reducing how accurately field measurements can be spatially extrapolated. However, current atmospheric reanalyses have grid spacings (>25 km) that are too coarse to resolve this variability. To address this limitation, we construct a fine-scale accumulation product from airborne snow radar observations by superimposing along-track fluctuations in accumulation onto an atmospheric reanalysis product. Our resulting airborne product reflects large-scale (>25 km) orographic precipitation patterns while providing robust and unprecedented insight into Antarctic accumulation variability on subgrid scales. On these smaller scales, we find significant, regionally dependent accumulation variability (𝜎(sub relative) > 40%). This variability in accumulation is correlated with variability in topographic surface slope in the wind direction (p < 0.01), confirming that subgrid-scale accumulation variability is driven by snow redistribution by wind.

Dattler, Marissa E.↗

Publishing Variables Archived at GES DISC to Earth System Grid Federation (ESGF)

We present a straightforward and low-cost approach to publish variables archived at NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) to the Earth System Grid Federation (ESGF). An ESGF publication requires a single standard-name variable aggregated over time to facilitate data inter-comparison. It also contains significant metadata to enable searching in ESGF. We look up standard names on high demand in ESGF search history, and using OPeNDAP and NcML technologies we aggregate the corresponding variables available in the GES DISC archive with augmented metadata required by CMIP6 and obs4MIPs Data Specification version 2.1. At this writing 10 variables from a standard product of the Atmospheric Infrared Sounder along with the Tech Notes are published in ESGF by NASA Center for Climate Simulation (NCCS). Users can view, analyze, and subset remotely, and download these aggregated variables via links in any ESGF node after searching. We plan to work on and publish more variables and data from different NASA missions and experiments in our archive.

Fan Fang↗

Radiography in High Mass X-Ray Binaries: Micro-Structure of the Stellar Wind Through Variability of the Column Density

Context. In high mass X-ray binaries, an accreting compact object orbits a high mass star, which loses mass through a dense and in homogeneous wind. Aims. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass-loss rates for high mass stars. Methods. We developed a simplified representation of the stellar wind where all the matter is gathered in spherical “clumps” that are radially advected away from the star. This model enables us to explore the connections between the stochastic properties of the wind and the variability of the column density for a comprehensive set of parameters related to the orbit and to the wind micro-structure, such as the size of the clumps and their individual mass. In particular, we focus on the evolution with the orbital phase of the standard deviation of the column density and of the characteristic duration of enhanced absorption episodes. Using the porosity length, we derive analytical predictions and compare them to the standard deviations and coherence time scales that were obtained. Results. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a single clump in front of the compact object. We thus provide a procedure to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. Conclusions. The coherence time scale grants direct access to the size of the clumps, while their mass can be deduced separately from the amplitude of the variability. We further show how monitoring the variability at superior conjunctions can probe the onset of the clump-forming region above the stellar photosphere. If the high column density variations in some high mass X-ray binaries are due to unaccreted clumps which are passing by the line-of-sight, this would require high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are marginally compatible with the ones derived from radiative-hydrodynamics simulations. Alternatively, the following components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism that has yet to be identified or a dense environment in the immediate vicinity of the accretor, such as an accretion disk, an outflow, or a spherical shell surrounding the magnetosphere of the accreting neutron star.

stars: mass loss↗

A search for variability in exoplanet analogues and low-gravity brown dwarfs

We report the results of a J-band survey for photometric variability in a sample of young, low-gravity objects using the New Technology Telescope (NTT) and the United Kingdom Infrared Telescope (UKIRT). Surface gravity is a key parameter in the atmospheric properties of brown dwarfs and this is the first large survey that aims to test the gravity dependence of variability properties. We do a full analysis of the spectral signatures of youth and assess the group membership probability of each target using membership tools from the literature. This results in a 30 object sample of young low-gravity brown dwarfs. Since we are lacking in objects with spectral types later than L9, we focus our statistical analysis on the L0–L8.5 objects. We find that the variability occurrence rate of L0–L8.5 low-gravity brown dwarfs in this survey is 30(+16,−8) per cent⁠. We reanalyse the results of Radigan (2014) and find that the field dwarfs with spectral types L0–L8.5 have a variability occurrence rate of 11(+13,−4) per cent⁠. We determine a probability of 98 per cent that the samples are drawn from different distributions. This is the first quantitative indication that the low-gravity objects are more likely to be variable than the field dwarf population. Furthermore, we present follow-up J(S) and K(S) observations of the young, planetary-mass variable object PSO 318.5–22 over three consecutive nights. We find no evidence of phase shifts between the J(S) and K(S) bands and find higher J(S) amplitudes. We use the J(S) light curves to measure a rotational period of 8.45 ± 0.05 h for PSO 318.5–22.

Johanna M. Vos↗

The Stochastic X-ray Variability of the Accreting Millisecond Pulsar IGR J17062-6143

This work presents an investigation of the stochastic X-ray variability from the 164 Hz accreting millisecond pulsar IGR J17062–6143, based on regular observations collected with the Neutron Star Interior Composition Explorer between 2017 July and 2020 August. Over this period, the power density spectrum showed a stable morphology, with broad ∼ 25% rms band-limited noise below 16 Hz. Quasiperiodic oscillations (QPOs) were occasionally observed, with the most notably detections including a low-frequency QPO centered at 2.7 Hz and a sharp QPO centered at 115 Hz that may be a 2:3 resonance with the spin frequency. Further, the energy dependence of the band-limited noise is studied through a spectroscopic analysis of the complex covariance in two frequency intervals. It is found that the power law continuum is the primary driver for the observed variability, although the thermal (blackbody) emission also appears to be intrinsically variable in area (5% rms) and temperature (1% rms). Notably, the 1 keV emission feature seen in all X-ray spectra of IGR J17062–6143 varies with the same amplitude as the power law, but systematically lags behind that continuum emission. These results appear consistent with a scenario in which a time variable Compton-scattering corona is the primary source for the observed stochastic variability, with the variability observed in the emission feature and at the lowest photon energies being due to the disk reflection of the power law continuum.

Petrus Bult↗

Atmospheric Variability and Measurement Uncertainty: Pitfalls in Averaging in situ Data

Satellite measurements and atmospheric models, two essential components of the integrated global observing system, provide crucial tools for monitoring and predicting regional and global foci, spanning numerous Earth Science fields. Unfortunately, models can lack the spatial and temporal resolution needed to resolve finer scale structure. Satellite measurements also have similar tempo-spatial restriction issues, but additionally can only measure certain species and can have biases that must be evaluated. Ground measurement networks are critical components of this system, by providing both independent inputs of species needed by models (including those that satellites do not provide) and assisting with investigation of biases in satellite products. Similarly, aircraft measurements play a vital role in closing the gaps between satellite measurements, model products, and ground monitoring networks by providing high accuracy, high-resolution data on local to regional spatial scales. Therefore, both ground-based and airborne observations are widely used to assess model predictions and satellite observations. One of the great challenges in using both ground and aircraft data in this fashion is matching the data temporally and spatially to the model/satellite data. A meaningful comparison with model or satellite requires a solid assessment of the variability of the in-situ measurements, which include both the instrument uncertainty and the statistical uncertainty due to atmospheric variability. While instrument uncertainty is generally more straightforwardly characterized, it can be challenging to accurately capture this variability uncertainty as it often presents in a non-Gaussian manner (e.g. emission plumes, frontal passages). We will present results examining spatial and temporal variability over a selection of scales relevant to satellite measurements and models of several in situ measurement species spanning both airborne and ground measurements. The extent of non-Gaussian variability will be quantified, and we will discuss additional statistical parameters that help assess the fitness of gaussian variability assumption when temporally or spatially averaging.

satellite validation↗

Near-infrared Variability in the 2MASS Calibration Fields: A Search for Planetary Transit Candidates

The Two Micron All Sky Survey (2MASS) photometric calibration observations cover approximately 6 square degrees on the sky in 35 'calibration fields,' each sampled in nominal photometric conditions between 562 and 3692 times during the 4 years of the 2MASS mission. We compile a catalog of variables from the calibration observations to search for M dwarfs transited by extrasolar planets. We present our methods for measuring periodic and nonperiodic flux variability. From 7554 sources with apparent K(sub s) magnitudes between 5.6 and 16.1, we identify 247 variables, including extragalactic variables and 23 periodic variables. We have discovered three M dwarf eclipsing systems, including two candidates for transiting extrasolar planets.

statistical methods↗

On the Lack of Stratospheric Dynamical Variability in Low-top Versions of the CMIP5 Models

We describe the main differences in simulations of stratospheric climate and variability by models within the fifth Coupled Model Intercomparison Project (CMIP5) that have a model top above the stratopause and relatively fine stratospheric vertical resolution (high-top), and those that have a model top below the stratopause (low-top). Although the simulation of mean stratospheric climate by the two model ensembles is similar, the low-top model ensemble has very weak stratospheric variability on daily and interannual time scales. The frequency of major sudden stratospheric warming events is strongly underestimated by the low-top models with less than half the frequency of events observed in the reanalysis data and high-top models. The lack of stratospheric variability in the low-top models affects their stratosphere-troposphere coupling, resulting in short-lived anomalies in the Northern Annular Mode, which do not produce long-lasting tropospheric impacts, as seen in observations. The lack of stratospheric variability, however, does not appear to have any impact on the ability of the low-top models to reproduce past stratospheric temperature trends. We find little improvement in the simulation of decadal variability for the high-top models compared to the low-top, which is likely related to the fact that neither ensemble produces a realistic dynamical response to volcanic eruptions.

troposphere↗

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data. Part II: Using Maximum Covariance Analysis to Effectively Compare Spatiotemporal Variability of Satellite and AERONET Measured Aerosol Optical Depth

Moderate Resolution Imaging SpectroRadiometer (MODIS) and Multi-angle Imaging Spectroradiomater (MISR) provide regular aerosol observations with global coverage. It is essential to examine the coherency between space- and ground-measured aerosol parameters in representing aerosol spatial and temporal variability, especially in the climate forcing and model validation context. In this paper, we introduce Maximum Covariance Analysis (MCA), also known as Singular Value Decomposition analysis as an effective way to compare correlated aerosol spatial and temporal patterns between satellite measurements and AERONET data. This technique not only successfully extracts the variability of major aerosol regimes but also allows the simultaneous examination of the aerosol variability both spatially and temporally. More importantly, it well accommodates the sparsely distributed AERONET data, for which other spectral decomposition methods, such as Principal Component Analysis, do not yield satisfactory results. The comparison shows overall good agreement between MODIS/MISR and AERONET AOD variability. The correlations between the first three modes of MCA results for both MODIS/AERONET and MISR/ AERONET are above 0.8 for the full data set and above 0.75 for the AOD anomaly data. The correlations between MODIS and MISR modes are also quite high (greater than 0.9). We also examine the extent of spatial agreement between satellite and AERONET AOD data at the selected stations. Some sites with disagreements in the MCA results, such as Kanpur, also have low spatial coherency. This should be associated partly with high AOD spatial variability and partly with uncertainties in satellite retrievals due to the seasonally varying aerosol types and surface properties.

aerosols↗

Relative Contributions of Mean-State Shifts and ENSO-Driven Variability to Precipitation Changes in a Warming Climate

El Niño-Southern Oscillation (ENSO) is an important driver of regional hydroclimate variability through far-reaching teleconnections. This study uses simulations performed with coupled general circulation models (CGCMs) to investigate how regional precipitation in the twenty-first century may be affected by changes in both ENSO-driven precipitation variability and slowly evolving mean rainfall. First, a dominant, time-invariant pattern of canonical ENSO variability (cENSO) is identified in observed SST data. Next, the fidelity with which 33 state-of-the-art CGCMs represent the spatial structure and temporal variability of this pattern (as well as its associated precipitation responses) is evaluated in simulations of twentieth-century climate change. Possible changes in both the temporal variability of this pattern and its associated precipitation teleconnections are investigated in twenty-first-century climate projections. Models with better representation of the observed structure of the cENSO pattern produce winter rainfall teleconnection patterns that are in better accord with twentieth-century observations and more stationary during the twenty-first century. Finally, the model-predicted twenty-first-century rainfall response to cENSO is decomposed into the sum of three terms: 1) the twenty-first-century change in the mean state of precipitation, 2) the historical precipitation response to the cENSO pattern, and 3) a future enhancement in the rainfall response to cENSO, which amplifies rainfall extremes. By examining the three terms jointly, this conceptual framework allows the identification of regions likely to experience future rainfall anomalies that are without precedent in the current climate.

Simulation↗

Sensitivity of Water Scarcity Events to ENSO-Driven Climate Variability at the Global Scale

Globally, freshwater shortage is one of the most dangerous risks for society. Changing hydro-climatic and socioeconomic conditions have aggravated water scarcity over the past decades. A wide range of studies show that water scarcity will intensify in the future, as a result of both increased consumptive water use and, in some regions, climate change. Although it is well-known that El Niño- Southern Oscillation (ENSO) affects patterns of precipitation and drought at global and regional scales, little attention has yet been paid to the impacts of climate variability on water scarcity conditions, despite its importance for adaptation planning. Therefore, we present the first global-scale sensitivity assessment of water scarcity to ENSO, the most dominant signal of climate variability. We show that over the time period 1961-2010, both water availability and water scarcity conditions are significantly correlated with ENSO-driven climate variability over a large proportion of the global land area (> 28.1 %); an area inhabited by more than 31.4% of the global population. We also found, however, that climate variability alone is often not enough to trigger the actual incidence of water scarcity events. The sensitivity of a region to water scarcity events, expressed in terms of land area or population exposed, is determined by both hydro-climatic and socioeconomic conditions. Currently, the population actually impacted by water scarcity events consists of 39.6% (CTA: consumption-to-availability ratio) and 41.1% (WCI: water crowding index) of the global population, whilst only 11.4% (CTA) and 15.9% (WCI) of the global population is at the same time living in areas sensitive to ENSO-driven climate variability. These results are contrasted, however, by differences in growth rates found under changing socioeconomic conditions, which are relatively high in regions exposed to water scarcity events. Given the correlations found between ENSO and water availability and scarcity conditions, and the relative developments of water scarcity impacts under changing socioeconomic conditions, we suggest that there is potential for ENSO-based adaptation and risk reduction that could be facilitated by more research on this emerging topic.

risk↗

The Role of Low-Level, Terrain-Induced Jets in Rainfall Variability in Tigris Euphrates Headwaters

Rainfall variability in the Tigris Euphrates headwaters is a result of interaction between topography and meteorological features at a range of spatial scales. Here, the Weather Research and Forecasting (WRF) Model, driven by the NCEP-DOE AMIP-II reanalysis (R-2), has been implemented to better understand these interactions. Simulations were performed over a domain covering most of the Middle East. The extended simulation period (1983 - 2013) enables us to study seasonality, interannual variability, spatial variability, and extreme events of rainfall. Results showed that the annual cycle of precipitation produced by WRF agrees much more closely with observations than does R-2. This was particularly evident during the transition months of April and October, which were further examined to study the underlying physical mechanisms. In both months, WRF improves representation of interannual variability relative to R-2, with a substantially larger benefit in April. This improvement results primarily from WRFs ability to resolve two low-level, terrain-induced flows in the region that are either absent or weak in R-2: one parallel to the western edge of the Zagros Mountains, and one along the east Turkish highlands. The first shows a complete reversal in its direction during wet and dry days, when flowing southeasterly it transports moisture from the Persian Gulf to the region, and when flowing northwesterly it blocks moisture and transports it away from the region. The second is more directly related to synoptic-scale systems and carries moist, warm air from the Mediterranean and Red Seas toward the region. The combined contribution of these flows explains about 50 of interannual variability in both WRF and observations for April and October precipitation.

rainfall variability↗

Detection of Vegetation Trends in Highly Variable Environments After Grazing Exclusion in Mongolia

Aims: Environmental variability (e.g. in precipitation) has a large effect on vegetation dynamics, and this often makes it difficult to assess the recovery of vegetation after a disturbance. In this study, we assessed vegetation recovery trends in response to grazing exclusion while considering for the annual environmental variability. Location: Two regions with highly variable precipitation: a steppe near Mandalgobi, in Mongolia's Central Gobi province, and a desert steppe near Bulgan, in Mongolia's South Gobi province. Methods: Changes in vegetation were observed along grazing gradients at the above two sites, and vegetation thresholds were identified. We established reference plots in pre-threshold areas along the grazing gradients. We removed the impact of livestock grazing at various locations along the grazing gradients by establishing exclosures, and investigated vegetation from 2005 to 2013. We developed a smoothed hierarchical model within a Bayesian framework, and examined the effect of grazing exclusion on vegetation, focusing especially on the extent of grass cover recovery. In addition, we compared soil nutrient conditions in the reference plots and inside and outside each exclosure along the grazing gradients. Results: Temporal trends in the cover of perennial grass in each plot inside and outside of the exclosures largely coincided, irrespective of grazing intensity, and exclosure had no effect (Bulgan) or a negative effect (Mandalgobi) on vegetation recovery. Soil nutrient content was not significantly affected by exclosure at a given distance from the grazing source, but decreased significantly with decreasing grazing intensity. Thus, recovery of the land from a post-threshold state may not be apparent even after 9 yr of grazing exclusion in environments with highly variable precipitation. Conclusion: The effect of exclosure duration on perennial grass cover was limited, even after controlling for environmental variability. Once a vegetation threshold has been crossed, merely removing livestock from the landscape may not be sufficient for that area to recover.

Bayesian statistics;Drought;Environmental variabil↗

Satellite Observations and Modeling to Understand the Lower Mekong River Basin Streamflow Variability

In this work, we have used the Soil & Water Assessment Tool (SWAT) to examine streamflow variability of the Lower Mekong River Basin (LMRB) associated with changes in the Upper Mekong River Basin (UMRB) inflows. Two hypothetical experiments were formulated and evaluated for the LMRB, where we conducted runoff simulations with multiple inflow changes that include upstream runoff yield increase and decrease scenarios. Streamflow variability of the LMRB was quantified by two streamflow metrics that explain flow variability and predictability, and high flow disturbance. The model experiments were performed for the Lower Mekong River Basin with identical climate, soil, and other watershed characteristics data. Remote sensing precipitation (Tropical Rainfall Measurement Mission, TRMM, and Global Precipitation Measurement mission, GPM), meteorological data as well as spatial data that include a digital elevation model, newly developed soil information (Harmonized World Soil Database, HWSD), and land use and land cover were processed as input to the LMRB model simulations. Observed daily streamflow data along the Lower Mekong River from Chiang Sean, Thailand to Kratie, Cambodia were used for calibration and validation. Our work results suggest that the Lower Mekong River streamflow is highly variable and has a low predictability (Colwell index of about 32%). We found that releasing more water from upstream Mekong during rainfall months by 30% would result in a reduction in the Lower Mekong streamflow predictability by about 21%. This reduction in predictability is mainly attributed to a decrease in the Contingency index. Our work shows that the ability to predict floods/droughts at the Lower Mekong River would be reduced if there is any anticipated change (i.e., increase/decrease) from UMRB releases. Our results also show that releasing more flows from the upstream Mekong would also affect flood duration and the frequency of flood occurrences downstream. The results of this work thus help to quantify the sensitivity of streamflow variability at the Lower Mekong River Basin to upstream anthropogenic changes.

Streamflow Predictability↗

Luminous Blue Variable Outbursts from the Variations of Helium Opacity

Luminous blue variables are evolved massive stars that exhibit poorly understood variability on timescales from months to years with significant mass loss rates. They typically reach an effective temperature of ≈ 9000K during outbursts. These variations and large mass loss rates are believed to be caused by the radiation force on the cooler, more opaque outer layers that can balance or even exceed the force of gravity, although the exact mechanisms for the outburst and variability are unknown and cannot be calculated in one dimension. Here we report three dimensional simulations of radiation-dominated massive stars, which show that helium opacity plays an important role for the outburst and setting the observed effective temperature ≈ 9000K. It also likely triggers the episodic mass loss at rates of 10(exp -7) – 10(exp -5)M(sub ⨀)/yr. The helium opacity peak only shows up in our three-dimensional simulations but not in one dimensional calculation because the envelope structures need to be determined self-consistently with convection. The simulations not only reproduce similar observed long time scale variabilities, but also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days and reaches larger amplitudes for even cooler stars. These short time scale variabilities should be observed with future high cadence observations.

Helium opacities↗