Stellar Winds and Dust Avalanches in the AU Mic Debris Disk

Stellar Winds and Dust Avalanches in the AU Mic Debris Disk

Authors: 

Chiang et al 

Abstract:
We explain the fast-moving, ripple-like features in the edge-on debris disk orbiting the young M dwarf AU Mic. The bright features are clouds of submicron dust repelled by the host star's wind. The clouds are produced by avalanches: radial outflows of dust that gain exponentially more mass as they shatter background disk particles in collisional chain reactions. The avalanches are triggered from a region a few au across—the "avalanche zone"—located on AU Mic's primary "birth" ring at a true distance of ~35 au from the star but at a projected distance more than a factor of 10 smaller: the avalanche zone sits directly along the line of sight to the star, on the side of the ring nearest Earth, launching clouds that disk rotation sends wholly to the southeast, as observed. The avalanche zone marks where the primary ring intersects a secondary ring of debris left by the catastrophic disruption of a progenitor up to Varuna in size, less than tens of thousands of years ago. Only where the rings intersect are particle collisions sufficiently violent to spawn the submicron dust needed to seed the avalanches. We show that this picture works quantitatively, reproducing the masses, sizes, and velocities of the observed escaping clouds. The Lorentz force exerted by the wind's magnetic field, whose polarity reverses periodically according to the stellar magnetic cycle, promises to explain the observed vertical undulations. The timescale between avalanches, about 10 yr, might be set by time variability of the wind mass loss rate or, more speculatively, by some self-regulating limit cycle.

Exoplanets as probes of the winds of host stars: the case of the M dwarf GJ 436

Exoplanets as probes of the winds of host stars: the case of the M dwarf GJ 436

Authors:

Vidotto et al

Abstract:
Winds of cool dwarfs are difficult to observe, with only a few M dwarfs presenting observationally-derived mass-loss rates (Mdot), which span several orders of magnitude. Close-in exoplanets are conveniently positioned in the inner regions of stellar winds and can, thus, be used to probe the otherwise-unobservable local properties of their host-stars' winds. Here, we use local stellar wind characteristics observationally-derived in the studies of atmospheric evaporation of the warm-neptune GJ436 b to derive the global characteristics of the wind of its M-dwarf host. Using an isothermal wind model, we constrain the stellar wind temperature to be in the range [0.36,0.43] MK, with Mdot=[0.5,2.5] x 10^{-15} Msyn/yr. By computing the pressure balance between the stellar wind and the interstellar medium, we derive the size of the astrophere of GJ436 to be around 25 au, significantly more compact than the heliosphere. We demonstrate in this paper that transmission spectroscopy, coupled to planetary atmospheric evaporation and stellar wind models, can be a useful tool for constraining the large-scale wind structure of planet-hosting stars. Extending our approach to future planetary systems discoveries will open new perspectives for the combined characterisation of planetary exospheres and winds of cool dwarf stars.

A wind-driving disc model for the mm-wavelength polarization structure of HL Tau

A wind-driving disc model for the mm-wavelength polarization structure of HL Tau

Authors:

Matsakos et al

Abstract:

The recent advent of spatially resolved mm- and cm-wavelength polarimetry in protostellar accretion discs could help clarify the role of magnetic fields in the angular momentum transport in these systems. The best case to date is that of HL~Tau, where the inability to produce a good fit to the 1.25-mm data with a combination of vertical and azimuthal magnetic field components was interpreted as implying that centrifugally driven winds (CDWs) are probably not a significant transport mechanism on the ∼102au scale probed by the observations. Using synthetic polarization maps of heuristic single-field-component discs and of a post-processed simulation of a wind-driving disc, we demonstrate that a much better fit to the data can be obtained if the radial field component, a hallmark of the CDW mechanism, dominates in the polarized emission region. A similar inference was previously made in modelling the far-infrared polarization map of the pc-scale dust ring in the Galactic centre. To reconcile this interpretation with theoretical models of protostellar discs, which indicate that the wind is launched from a comparatively high elevation above the mid-plane, we propose that most of the polarized emission originates -- with a high ($\ga 10$\%) intrinsic degree of polarization -- in small ($\la 0.1\,$mm) grains that remain suspended above the mid-plane, and that the bulk of the mm-wavelength emission is produced -- with low intrinsic polarization -- by larger grains that have settled to the mid-plane.

The Space Weather of Proxima Centauri b

The Space Weather of Proxima Centauri b

Authors:

Garraffo et al

Abstract:

A planet orbiting in the "habitable zone" of our closest neighboring star, Proxima Centauri, has recently been discovered, and the next natural question is whether or not Proxima b is "habitable". Stellar winds are likely a source of atmospheric erosion that could be particularly severe in the case of M dwarf habitable zone planets that reside close to their parent star. Here we study the stellar wind conditions that Proxima b experiences over its orbit. We construct 3-D MHD models of the wind and magnetic field around Proxima Centauri using a surface magnetic field map for a star of the same spectral type and scaled to match the observed ~ 600 G surface magnetic field strength of Proxima. We examine the wind conditions and dynamic pressure over different plausible orbits that sample the constrained parameters of the orbit of Proxima b. For all the parameter space explored, the planet is subject to stellar wind pressures of more than 2000 times those experienced by Earth from the solar wind. During an orbit, Proxima b is also subject to pressure changes of 1 to 3 orders of magnitude on timescales of a day. Its magnetopause standoff distance consequently undergoes sudden and periodic changes by a factor of 2 to 5. Proxima b will traverse the interplanetary current sheet twice each orbit, and likely crosses into regions of subsonic wind quite frequently. These effects should be taken into account in any physically realistic assessment or prediction of its atmospheric reservoir, characteristics and loss.

Tracing Slow Winds from T Tauri Stars via Low Velocity Forbidden Line Emission

Tracing Slow Winds from T Tauri Stars via Low Velocity Forbidden Line Emission

Authors:

Simon et al

Abstract:

Using Keck/HIRES spectra {\Delta}v ~ 7 km/s, we analyze forbidden lines of [O I] 6300 {\AA}, [O I] 5577 {\AA} and [S II] 6731 {\AA} from 33 T Tauri stars covering a range of disk evolutionary stages. After removing a high velocity component (HVC) associated with microjets, we study the properties of the low velocity component (LVC). The LVC can be attributed to slow disk winds that could be magnetically (MHD) or thermally (photoevaporative) driven. Both of these winds play an important role in the evolution and dispersal of protoplanetary material.

LVC emission is seen in all 30 stars with detected [O I] but only in 2 out of eight with detected [S II] , so our analysis is largely based on the properties of the [O I] LVC. The LVC itself is resolved into broad (BC) and narrow (NC) kinematic components. Both components are found over a wide range of accretion rates and their luminosity is correlated with the accretion luminosity, but the NC is proportionately stronger than the BC in transition disks.

The FWHM of both the BC and NC correlates with disk inclination, consistent with Keplerian broadening from radii of 0.05 to 0.5 AU and 0.5 to 5 AU, respectively. The velocity centroids of the BC suggest formation in an MHD disk wind, with the largest blueshifts found in sources with closer to face-on orientations. The velocity centroids of the NC however, show no dependence on disk inclination. The origin of this component is less clear and the evidence for photoevaporation is not conclusive.

A wind-driving disk model for the mm-wavelength polarization structure of HL Tau

A wind-driving disc model for the mm-wavelength polarization structure of HL Tau

Authors:

Matsakos et al

Abstract:

The recent advent of spatially resolved mm- and cm-wavelength polarimetry in protostellar accretion discs could help clarify the role of magnetic fields in the angular momentum transport in these systems. The best case to date is that of HL Tau, where the inability to produce a good fit to the 1.25-mm data with a combination of vertical and azimuthal magnetic field components was interpreted as implying that centrifugally driven winds (CDWs) are probably not a significant transport mechanism on the ∼102 au scale probed by the observations. Using synthetic polarization maps of heuristic single-field-component discs and of a post-processed simulation of a wind-driving disc, we demonstrate that a much better fit to the data can be obtained if the radial field component, a hallmark of the CDW mechanism, dominates in the polarized emission region. A similar inference was previously made in modelling the far-infrared polarization map of the pc-scale dust ring in the Galactic centre. To reconcile this interpretation with theoretical models of protostellar discs, which indicate that the wind is launched from a comparatively high elevation above the mid-plane, we propose that most of the polarized emission originates – with a high (≳ 10%) intrinsic degree of polarization – in small (≲ 0.1 mm) grains that remain suspended above the mid-plane, and that the bulk of the mm-wavelength emission is produced – with low intrinsic polarization – by larger grains that have settled to the mid-plane.

The Effect of Stellar Evolution on Exoplanetary Atmospheres

Stellar magnetism, winds and their effects on planetary environments

Author:

Vidditto

Abstract:

Here, I review some recent works on magnetism of cool, main-sequence stars, their winds and potential impact on surrounding exoplanets. The winds of these stars are very tenuous and persist during their lifetime. Although carrying just a small fraction of the stellar mass, these magnetic winds carry away angular momentum, thus regulating the rotation of the star. Since cool stars are likely to be surrounded by planets, understanding the host star winds and magnetism is a key step towards characterisation of exoplanetary environments. As rotation and activity are intimately related, the spin down of stars leads to a decrease in stellar activity with age. As a consequence, as stars age, a decrease in high-energy (X-ray, extreme ultraviolet) irradiation is observed, which can a ect the evaporation of exoplanetary atmospheres and, thus, also altering exoplanetary evolution.

The Auroras of hot Jupiters

Stellar wind–magnetosphere interaction at exoplanets: computations of auroral radio powers

Authors:

Nichols et al

Abstract:

We present calculations of the auroral radio powers expected from exoplanets with magnetospheres driven by an Earth-like magnetospheric interaction with the solar wind. Specifically, we compute the twin cell-vortical ionospheric flows, currents, and resulting radio powers resulting from a Dungey cycle process driven by dayside and nightside magnetic reconnection, as a function of planetary orbital distance and magnetic field strength. We include saturation of the magnetospheric convection, as observed at the terrestrial magnetosphere, and we present power-law approximations for the convection potentials, radio powers and spectral flux densities. We specifically consider a solar-age system and a young (1 Gyr) system. We show that the radio power increases with magnetic field strength for magnetospheres with saturated convection potential, and broadly decreases with increasing orbital distance. We show that the magnetospheric convection at hot Jupiters will be saturated, and thus unable to dissipate the full available incident Poynting flux, such that the magnetic Radiometric Bode's Law (RBL) presents a substantial overestimation of the radio powers for hot Jupiters. Our radio powers for hot Jupiters are ∼5–1300 TW for hot Jupiters with field strengths of 0.1–10 BJ orbiting a Sun-like star, while we find that competing effects yield essentially identical powers for hot Jupiters orbiting a young Sun-like star. However, in particular, for planets with weaker magnetic fields, our powers are higher at larger orbital distances than given by the RBL, and there are many configurations of planet that are expected to be detectable using SKA.

Interaction of the Evaporating Hot Neptune GJ 436b & its Host Star's Stellar Wind

An evaporating planet in the wind: stellar wind interactions with the radiatively braked exosphere of GJ436 b

Authors:

Bourrier et al

Abstract

The warm Neptune GJ436b was observed with HST/STIS at three different epochs in the stellar Ly-alpha line, showing deep, repeated transits caused by a giant exosphere of neutral hydrogen. The low radiation pressure from the M-dwarf host star was shown to play a major role in the dynamics of the escaping gas. Yet by itself it cannot explain the time-variable spectral features detected in each transit. Here we investigate the combined role of radiative braking and stellar wind interactions using numerical simulations with the EVaporating Exoplanet code (EVE) and we derive atmospheric and stellar properties through the direct comparison of simulated and observed spectra. Our simulations match the last two epochs well. The observed sharp early ingresses come from the abrasion of the planetary coma by the stellar wind. Spectra observed during the transit can be produced by a dual exosphere of planetary neutrals (escaped from the upper atmosphere of the planet) and neutralized protons (created by charge-exchange with the stellar wind). We find similar properties at both epochs for the planetary escape rate (2.5x108 g/s), the stellar photoionization rate (2x10−5 /s), the stellar wind bulk velocity (85 km/s), and its kinetic dispersion velocity (10 km/s). We find high velocities for the escaping gas (50-60 km/s) that may indicate MHD waves that dissipate in the upper atmosphere and drive the planetary outflow. In the last epoch the high density of the stellar wind (3x103 /cm3) led to the formation of an exospheric tail mainly composed of neutralized protons. The observations of GJ436 b allow for the first time to clearly separate the contributions of radiation pressure and stellar wind and to probe the regions of the exosphere shaped by each mechanism.

The Stellar Winds of HL Tauri

The Winds from HL Tau

Authors:

Klaassen et al

Abstract:

Outflowing motions, whether a wind launched from the disk, a jet launched from the protostar, or the entrained molecular outflow, appear to be an ubiquitous feature of star formation. These outwards motions have a number of root causes, and how they manifest is intricately linked to their environment as well as the process of star formation itself.

Using the ALMA Science Verification data of HL Tau, we investigate the high velocity molecular gas being removed from the system as a result of the star formation process. We aim to place these motions in context with the optically detected jet, and the disk. With these high resolution (∼1") ALMA observations of CO (J=1-0), we quantify the outwards motions of the molecular gas. We find evidence for a bipolar outwards flow, with an opening angle, as measured in the red-shifted lobe, starting off at 90∘, and narrowing to 60∘ further from the disk, likely because of magnetic collimation. Its outwards velocity, corrected for inclination angle is of order 2.4 km s−1.

Hot planetary winds near a star: dynamics, wind-wind interactions, and observational signatures

Hot planetary winds near a star: dynamics, wind-wind interactions, and observational signatures

Authors:

Carroll-Nellenback et al

Abstract:

Signatures of "evaporative" winds from exo-planets on short (hot) orbits around their host star have been observed in a number of systems. In this paper we present global AMR simulations that track the launching of the winds, their expansion through the circumstellar environment, and their interaction with a stellar wind. We focus on purely hydrodynamic flows including the anisotropy of the wind launching and explore the orbital/fluid dynamics of the resulting flows in detail. In particular we find that a combination of the tidal and Coriolis forces strongly distorts the planetary "Parker" wind creating "up-orbit" and "down-orbit" streams. We characterize the flows in terms of their orbital elements which change depending on their launch position on the planet. We find that the anisotropy in the atmospheric temperature leads to significant backflow on to the planet. The planetary wind interacts strongly with the stellar wind creating instabilities that cause eventual deposition of planetary gas onto the star. We present synthetic observations of both transit and absorption line-structure for our simulations. For our initial conditions, we find that the orbiting wind material produces absorption signatures at significant distances from the planet and substantial orbit to orbit variability. Ly-{\alpha} absorption shows red and blueshifted features out to 70 km/s. Finally, using semi-analytic models we constrain the effect of radiation pressure, given the approximation of uniform stellar absorption.

Temporal variability of the wind from the star τ Boötis

Temporal variability of the wind from the star τ Boötis

Authors:

Nicholson et al

Abstract:

We present new wind models for {\tau} Bo\"otis ({\tau} Boo), a hot-Jupiter-host-star whose observable magnetic cycles makes it a uniquely useful target for our goal of monitoring the temporal variability of stellar winds and their exoplanetary impacts. Using spectropolarimetric observations from May 2009 to January 2015, the most extensive information of this type yet available, to reconstruct the stellar magnetic field, we produce multiple 3D magnetohydrodynamic stellar wind models. Our results show that characteristic changes in the large-scale magnetic field as the star undergoes magnetic cycles produce changes in the wind properties, both globally and locally at the position of the orbiting planet. Whilst the mass loss rate of the star varies by only a minimal amount (∼ 4 percent), the rates of angular momentum loss and associated spin-down timescales are seen to vary widely (up to ∼ 140 percent), findings consistent with and extending previous research. In addition, we find that temporal variation in the global wind is governed mainly by changes in total magnetic flux rather than changes in wind plasma properties. The magnetic pressure varies with time and location and dominates the stellar wind pressure at the planetary orbit. By assuming a Jovian planetary magnetic field for {\tau} Boo b, we nevertheless conclude that the planetary magnetosphere can remain stable in size for all observed stellar cycle epochs, despite significant changes in the stellar field and the resulting local space weather environment.

How the Magnetic Field and Stellar Wind of Kappa Ceti Would Effect the Habitability of Terretrial Worlds

Magnetic field and wind of Kappa Ceti: towards the planetary habitability of the young Sun when life arose on Earth

Authors:

do Nascimento Jr. et al

Abstract:

We report magnetic field measurements for Kappa1~Cet, a proxy of the young Sun when life arose on Earth. We carry out an analysis of the magnetic properties determined from spectropolarimetric observations and reconstruct its large-scale surface magnetic field to derive the magnetic environment, stellar winds and particle flux permeating the interplanetary medium around Kappa1~Cet. Our results show a closer magnetosphere and mass-loss rate of Mdot = 9.7 x 10^{-13} Msol/yr, i.e., a factor 50 times larger than the current solar wind mass-loss rate, resulting in a larger interaction via space weather disturbances between the stellar wind and a hypothetical young-Earth analogue, potentially affecting the planet's habitability. Interaction of the wind from the young Sun with the planetary ancient magnetic field may have affected the young Earth and its life conditions

How dusty are photoevaporative winds?

How dusty are photoevaporative winds?

Authors:

Hutchinson et al

Abstract:

Not very. We perform dusty smoothed particle hydrodrodynamic (SPH) calculations of photoevaporation in protoplanetary discs. We use unequal-mass particles to resolve more than five orders of magnitude in disc/outflow density and a one-fluid formulation to efficiently simulate an equivalent magnitude range in drag stopping time. We find that only micron sized dust grains and smaller can be entrained in photoevaporative winds. This makes it difficult to explain the dust holes seen in transition discs using photoevaporation and implies that only small grains can be transported to the outer disc by this mechanism. A pileup of micron sized dust grains can occur in the upper atmosphere at critical radii in the disc as grains decouple from the low-density wind. Entrainment is a strong function of location in the disc, resulting in a size sorting of grains in the outflow---the largest grain being carried out between 10--20AU. The peak dust density for each grain occurs at the inner edge of its entrainment region.

Axisymmetric Simulations of Hot Jupiter-Stellar Wind Hydrodynamic Interaction

Axisymmetric Simulations of Hot Jupiter-Stellar Wind Hydrodynamic Interaction

Authors:

Christie et al

Abstract:

Gas giant exoplanets orbiting at close distances to the parent star are subjected to large radiation and stellar wind fluxes. In this paper, hydrodynamic simulations of the planetary upper atmosphere and its interaction with the stellar wind are carried out to understand the possible flow regimes and how they affect the Lyman-alpha transmission spectrum. Following Tremblin and Chiang, charge exchange reactions are included to explore the role of energetic atoms as compared to thermal particles. In order to understand the role of the tail as compared to the leading edge of the planetary gas, the simulations were carried out under axisymmetry, and photoionization and stellar wind electron impact ionization reactions were included to limit the extent of the neutrals away from the planet. By varying the planetary gas temperature, two regimes are found. At high temperature, a supersonic planetary wind is found, which is turned around by the stellar wind and forms a tail behind the planet. At lower temperatures, the planetary wind is shut off when the stellar wind penetrates inside where the sonic point would have been. In this regime mass is lost by viscous interaction at the boundary between planetary and stellar wind gases. Absorption by cold hydrogen atoms is large near the planetary surface, and decreases away from the planet as expected. The hot hydrogen absorption is in an annulus and typically dominated by the tail, at large impact parameter, rather than by the thin leading edge of the mixing layer near the substellar point.

Simulating the Environment Around Planet-Hosting Stars HD 1237, HD 22049 and HD 147513

Simulating the Environment Around Planet-Hosting Stars - I. Coronal Structure

Authors:

Alvarado-Gómez et al

Abstract:

We present the results of a detailed numerical simulation of the circumstellar environment around three exoplanet-hosting stars. A state-of-the-art global magnetohydrodynamic (MHD) model is considered, including Alfv\'en wave dissipation as a self-consistent coronal heating mechanism. This paper contains the description of the numerical set-up, evaluation procedure, and the simulated coronal structure of each system (HD 1237, HD 22049 and HD 147513). The simulations are driven by surface magnetic field maps, recovered with the observational technique of Zeeman Doppler Imaging (ZDI). A detailed comparison of the simulations is performed, where two different implementations of this mapping routine are used to generate the surface field distributions. Quantitative and qualitative descriptions of the coronae of these systems are presented, including synthetic high-energy emission maps in the Extreme Ultra-Violet (EUV) and Soft X-rays (SXR) ranges. Using the simulation results, we are able to recover similar trends as in previous observational studies, including the relation between the magnetic flux and the coronal X-ray emission. Furthermore, for HD 1237 we estimate the rotational modulation of the high-energy emission due to the various coronal features developed in the simulation. We obtain variations, during a single stellar rotation cycle, up to 15\% for the EUV and SXR ranges. The results presented here will be used, in a follow-up paper, to self-consistently simulate the stellar winds and inner astrospheres of these systems.

Stellar Wind/Magnetosphere Interactions in Hot Jupiters

Stellar Wind -- Magnetosphere Interactions in Hot Jupiters

Authors:

Buzasi et al

Abstract:

One potential star-planet interaction mechanism for hot Jupiters involves planetary heating via currents set up by interactions between the stellar wind and planetary magnetosphere. Early modeling results indicate that such currents, which are analogous to the terrestrial global electric circuit (GEC), have the potential to provide sufficient heating to account for the additional radius inflation seen in some hot Jupiters. Here we present a more detailed model of this phenomenon, exploring the scale of the effect, the circumstances under which it is likely to be significant, implications for the planetary magnetospheric structure, and observational signatures.

Do hot Jupiters Inflate Because of Solar Wind/Magnetic Field Interaction Driven Heating?

Extended Heat Deposition in Hot Jupiters: Application to Ohmic Heating

Authors:

Ginzburg et al

Abstract:

Many giant exoplanets in close orbits have observed radii which exceed theoretical predictions. One suggested explanation for this discrepancy is heat deposited deep inside the atmospheres of these "hot Jupiters". Here, we study extended power sources which distribute heat from the photosphere to the deep interior of the planet. Our analytical treatment is a generalization of a previous analysis of localized "point sources". We model the deposition profile as a power law in the optical depth and find that planetary cooling and contraction halt when the internal luminosity (i.e. cooling rate) of the planet drops below the heat deposited in the planet's convective region. A slowdown in the evolutionary cooling prior to equilibrium is possible only for sources which do not extend to the planet's center. We estimate the Ohmic dissipation resulting from the interaction between the atmospheric winds and the planet's magnetic field, and apply our analytical model to Ohmically heated planets. Our model can account for the observed radii of many inflated planets which have equilibrium temperatures ≈1500 K−2500 K, and are inflated to a radius ≈1.5RJ. However, some extremely inflated planets remain unexplained by our model. We also argue that Ohmically inflated planets have already reached their equilibrium phase, and no longer contract. Following Wu & Lithwick (2013) who argued that Ohmic heating could only suspend and not reverse contraction, we calculate the time it takes Ohmic heating to re-inflate a cold planet to its equilibrium configuration. We find that while it is possible to re-inflate a cold planet, the re-inflation timescales are longer by a factor of ≈30 than the cooling time.

Planet Wind/Stellar Wind Bow Shocks Influence Exoplanetary Evaporation

Planetary Evaporation and the Dynamics of Planet Wind/Stellar Wind Bow Shocks

Authors:

Frank et al

Abstract:

We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone &; Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a strength that depends on the escape parameter. When a stellar outflow is included, we see unstable bow waves forming through the outflow's interaction with the planetary wind.

Atmospheric Escape by Magnetically Driven Wind From Hot Jupiters

Atmospheric Escape by Magnetically Driven Wind from Gaseous Planets II --Effects of Magnetic Diffusion--

Authors:

Tanaka et al

Abstract:

We investigate roles of Alfvenic waves in the weakly-ionized atmosphere of hot Jupiters by carrying out non-ideal magnetohydrodynamic (MHD) simulations with Ohmic diffusion in one-dimensional magnetic flux tubes. Turbulence at the surface excites Alfven waves and they propagate upward to drive hot (~ 10^4 K) outflows. The magnetic diffusion plays an important role in the dissipation of the Alfvenic waves in the weakly ionized atmosphere of hot Jupiters. The mass-loss rate of the spontaneously driven planetary wind is considerably reduced, in comparison with that obtained from ideal MHD simulations because the Alfvenic waves are severely damped at low altitudes in the atmosphere, whereas the wave heating is still important in the heating of the upper atmosphere. Dependence on the surface temperature, planetary radius, and velocity dispersion at the surface is also investigated. We find an inversion phenomenon of the transmitted wave energy flux; the energy flux carried by Alfven waves in the upper atmosphere has a nonmonotonic correlation with the input energy flux from the surface in a certain range of the surface temperature because the resistivity is determined by the global physical properties of the atmosphere in a complicated manner. We also point out that the heating and mass loss are expected only in limited zones if the open magnetic field is confined in the limited regions.