Source author record

Steven R. Cranmer

Steven R. Cranmer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2020arXiv

Alfvén-wave driven magnetic rotator winds from low-mass stars I: rotation dependences of magnetic braking and mass-loss rate

Observations of stellar rotation show that low-mass stars lose angular momentum during the main sequence. We simulate the winds of Sun-like stars with a range of rotation rates, covering the fast and slow magneto-rotator regimes, including the transition between the two. We generalize an Alfvén-wave driven solar wind model that builds on previous works by including the magneto-centrifugal force explicitly. In this model, the surface-averaged open magnetic flux is assumed to scale as $B_\ast f^{\rm open}_\ast \propto {\rm Ro}^{-1.2}$, where $f^{\rm open}_\ast$ and ${\rm Ro}$ are the surface open-flux filling factor and Rossby number, respectively. We find that, 1. the angular momentum loss rate (torque) of the wind is described as $τ_w \approx 2.59 \times 10^{30} {\rm \ erg} \ \left( Ω_\ast / Ω_\odot \right)^{2.82}$, yielding a spin-down law $Ω_\ast \propto t^{-0.55}$. 2. the mass-loss rate saturates at $\dot{M}_w \sim 3.4 \times 10^{-14} M_\odot {\rm \ yr^{-1}}$, due to the strong reflection and dissipation of Alfvén waves in the chromosphere. This indicates that the chromosphere has a strong impact in connecting the stellar surface and stellar wind. Meanwhile, the wind ram pressure scales as $P_w \propto Ω_\ast^{0.57}$, which is able to explain the lower-envelope of the observed stellar winds by Wood et al. 3. the location of the Alfvén radius is shown to scale in a way that is consistent with 1D analytic theory. Additionally, the precise scaling of the Alfvén radius matches previous works which used thermally-driven winds. Our results suggest that the Alfvén-wave driven magnetic rotator wind plays a dominant role in the stellar spin-down during the main-sequence.

preprint2020arXiv

Heating Rates for Protons and Electrons in Polar Coronal Holes: Empirical Constraints from the Ultraviolet Coronagraph Spectrometer

Ultraviolet spectroscopy of the extended solar corona is a powerful tool for measuring the properties of protons, electrons, and heavy ions in the accelerating solar wind. The large coronal holes that expand up from the north and south poles at solar minimum are low-density collisionless regions in which it is possible to detect departures from one-fluid thermal equilibrium. An accurate characterization of these departures is helpful in identifying the kinetic processes ultimately responsible for coronal heating. In this paper, Ultraviolet Coronagraph Spectrometer (UVCS) measurements of the H I Lyman alpha line are analyzed to constrain values for the solar wind speed, electron density, electron temperature, proton temperature (parallel and perpendicular to the magnetic field) and Alfven-wave amplitude. The analysis procedure involves creating a large randomized ensemble of empirical models, simulating their Lyman alpha profiles, and building posterior probability distributions for only the models that agree with the UVCS data. The resulting temperatures do not exhibit a great deal of radial variation between heliocentric distances of 1.4 and 4 solar radii. Typical values for the electron, parallel proton, and perpendicular proton temperatures are 1.2, 1.8, and 1.9 MK, respectively. Resulting values for the "nonthermal" Alfven wave amplitude show evidence for weak dissipation, with a total energy-loss rate that agrees well with an independently derived total heating rate for the protons and electrons. The moderate Alfven-wave amplitudes appear to resolve some tension in the literature between competing claims of both higher (undamped) and lower (heavily damped) values.

preprint2016arXiv

Constraints on Planetesimal Collision Models in Debris Disks

Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3" resolution and $\sim5$ $μ$Jy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, $n(a)\propto a^{-q}$, of $\langle q \rangle = 3.36\pm0.02$, with a possible trend of decreasing $q$ for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution ($q=3.5$) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system.

preprint2016arXiv

Explaining Inverted Temperature Loops in the Quiet Solar Corona with Magnetohydrodynamic Wave Mode Conversion

Coronal loops trace out bipolar, arch-like magnetic fields above the Sun's surface. Recent measurements that combine rotational tomography, extreme ultraviolet imaging, and potential-field extrapolation have shown the existence of large loops with inverted temperature profiles; i.e., loops for which the apex temperature is a local minimum, not a maximum. These "down loops" appear to exist primarily in equatorial quiet regions near solar minimum. We simulate both these and the more prevalent large-scale "up loops" by modeling coronal heating as a time-steady superposition of: (1) dissipation of incompressible Alfven-wave turbulence, and (2) dissipation of compressive waves formed by mode conversion from the initial population of Alfven waves. We found that when a large percentage (> 99%) of the Alfven waves undergo this conversion, heating is greatly concentrated at the footpoints and stable "down loops" are created. In some cases we found loops with three maxima that are also gravitationally stable. Models that agree with the tomographic temperature data exhibit higher gas pressures for "down loops" than for "up loops," which is consistent with observations. These models also show a narrow range of Alfven wave amplitudes: 3 to 6 km/s at the coronal base. This is low in comparison to typical observed amplitudes of 20 to 30 km/s in bright X-ray loops. However, the large-scale loops we model are believed to comprise a weaker diffuse background that fills much of the volume of the corona. By constraining the physics of loops that underlie quiescent streamers, we hope to better understand the formation of the slow solar wind.

preprint2016arXiv

Predictions for Dusty Mass Loss from Asteroids during Close Encounters with Solar Probe Plus

The Solar Probe Plus (SPP) mission will explore the Sun's corona and innermost solar wind starting in 2018. The spacecraft will also come close to a number of Mercury-crossing asteroids with perihelia less than 0.3 AU. At small heliocentric distances, these objects may begin to lose mass, thus becoming "active asteroids" with comet-like comae or tails. This paper assembles a database of 97 known Mercury-crossing asteroids that may be encountered by SPP, and it presents estimates of their time-dependent visible-light fluxes and mass loss rates. Assuming a similar efficiency of sky background subtraction as was achieved by STEREO, we find that approximately 80% of these asteroids are bright enough to be observed by the Wide-field Imager for SPP (WISPR). A model of gas/dust mass loss from these asteroids is developed and calibrated against existing observations. This model is used to estimate the visible-light fluxes and spatial extents of spherical comae. Observable dust clouds occur only when the asteroids approach the Sun closer than 0.2 AU. The model predicts that during the primary SPP mission between 2018 and 2025, there should be 113 discrete events (for 24 unique asteroids) during which the modeled comae have angular sizes resolvable by WISPR. The largest of these correspond to asteroids 3200 Phaethon, 137924, 155140, and 289227, all with angular sizes of roughly 15 to 30 arcminutes. We note that the SPP trajectory may still change, but no matter the details there should still be multiple opportunities for fruitful asteroid observations.

preprint2016arXiv

Statistical study of network jets observed in the solar transition region: A comparison between coronal holes and quiet sun regions

Recent IRIS observations have revealed a prevalence of intermittent small-scale jets with apparent speeds of 80 - 250 km s$^{-1}$, emanating from small-scale bright regions inside network boundaries of coronal holes. We find that these network jets appear not only in coronal holes but also in quiet-sun regions. Using IRIS 1330A (C II) slit-jaw images, we extract several parameters of these network jets, e.g. apparent speed, length, lifetime and increase in foot-point brightness. Using several observations, we find that some properties of the jets are very similar but others are obviously different between the quiet sun and coronal holes. For example, our study shows that the coronal-hole jets appear to be faster and longer than those in the quiet sun. This can be directly attributed to a difference in the magnetic configuration of the two regions with open magnetic field lines rooted in coronal holes and magnetic loops often present in quiet sun. We have also detected compact bright loops, likely transition region loops, mostly in quiet sun. These small loop-like regions are generally devoid of network jets. In spite of different magnetic structures in the coronal hole and quiet sun in the transition region, there appears to be no substantial difference for the increase in foot-point brightness of the jets, which suggests that the generation mechanism of these network jets is likely the same in both regions.

preprint2015arXiv

Driving Solar Spicules and Jets with Magnetohydrodynamic Turbulence: Testing a Persistent Idea

The solar chromosphere contains thin, highly dynamic strands of plasma known as spicules. Recently, it has been suggested that the smallest and fastest (Type II) spicules are identical to intermittent jets observed by the Interface Region Imaging Spectrograph. These jets appear to expand out along open magnetic field lines rooted in unipolar network regions of coronal holes. In this paper we revisit a thirty-year-old idea that spicules may be caused by upward forces associated with Alfven waves. These forces involve the conversion of transverse Alfven waves into compressive acoustic-like waves that steepen into shocks. The repeated buffeting due to upward shock propagation causes nonthermal expansion of the chromosphere and a transient levitation of the transition region. Some older models of wave-driven spicules assumed sinusoidal wave inputs, but the solar atmosphere is highly turbulent and stochastic. Thus, we model this process using the output of a time-dependent simulation of reduced magnetohydrodynamic turbulence. The resulting mode-converted compressive waves are strongly variable in time, with a higher transition region occurring when the amplitudes are large and a lower transition region when the amplitudes are small. In this picture, the transition region bobs up and down by several Mm on timescales less than a minute. These motions produce narrow, intermittent extensions of the chromosphere that have similar properties as the observed jets and Type II spicules.

preprint2015arXiv

Time-Dependent Turbulent Heating of Open Flux Tubes in the Chromosphere, Corona, and Solar Wind

We investigate several key questions of plasma heating in open-field regions of the corona that connect to the solar wind. We present results for a model of Alfven-wave-driven turbulence for three typical open magnetic field structures: a polar coronal hole, an open flux tube neighboring an equatorial streamer, and an open flux tube near a strong-field active region. We compare time-steady, one-dimensional turbulent heating models (Cranmer et al., 2007) against fully time-dependent three-dimensional reduced-magnetohydrodynamics modeling of BRAID (van Ballegooijen et al., 2011). We find that the time-steady results agree well with time-averaged results from BRAID. The time-dependence allows us to investigate the variability of the magnetic fluctuations and of the heating in the corona. The high-frequency tail of the power spectrum of fluctuations forms a power law whose exponent varies with height, and we discuss the possible physical explanation for this behavior. The variability in the heating rate is bursty and nanoflare-like in nature, and we analyze the amount of energy lost via dissipative heating in transient events throughout the simulation. The average energy in these events is 10^21.91 erg, within the "picoflare" range, and many events reach classical "nanoflare" energies. We also estimated the multithermal distribution of temperatures that would result from the heating-rate variability, and found good agreement with observed widths of coronal differential emission measure (DEM) distributions. The results of the modeling presented in this paper provide compelling evidence that turbulent heating in the solar atmosphere by Alfven waves accelerates the solar wind in open flux tubes.

preprint2014arXiv

Ensemble Simulations of Proton Heating in the Solar Wind via Turbulence and Ion Cyclotron Resonance

Protons in the solar corona and heliosphere exhibit anisotropic velocity distributions, violation of magnetic moment conservation, and a general lack of thermal equilibrium with the other particle species. There is no agreement about the identity of the physical processes that energize non-Maxwellian protons in the solar wind, but a traditional favorite has been the dissipation of ion cyclotron resonant Alfven waves. This paper presents kinetic models of how ion cyclotron waves heat protons on their journey from the corona to interplanetary space. It also derives a wide range of new solutions for the relevant dispersion relations, marginal stability boundaries, and nonresonant velocity-space diffusion rates. A phenomenological model containing both cyclotron damping and turbulent cascade is constructed to explain the suppression of proton heating at low alpha-proton differential flow speeds. These effects are implemented in a large-scale model of proton thermal evolution from the corona to 1 AU. A Monte Carlo ensemble of realistic wind speeds, densities, magnetic field strengths, and heating rates produces a filled region of parameter space (in a plane described by the parallel plasma beta and the proton temperature anisotropy ratio) similar to what is measured. The high-beta edges of this filled region are governed by plasma instabilities and strong heating rates. The low-beta edges correspond to weaker proton heating and a range of relative contributions from cyclotron resonance. On balance, the models are consistent with other studies that find only a small fraction of the turbulent power spectrum needs to consist of ion cyclotron waves.

preprint2014arXiv

Suprathermal Electrons in the Solar Corona: Can Nonlocal Transport Explain Heliospheric Charge States?

There have been several ideas proposed to explain how the Sun's corona is heated and how the solar wind is accelerated. Some models assume that open magnetic field lines are heated by Alfven waves driven by photospheric motions and dissipated after undergoing a turbulent cascade. Other models posit that much of the solar wind's mass and energy is injected via magnetic reconnection from closed coronal loops. The latter idea is motivated by observations of reconnecting jets and also by similarities of ion composition between closed loops and the slow wind. Wave/turbulence models have also succeeded in reproducing observed trends in ion composition signatures versus wind speed. However, the absolute values of the charge-state ratios predicted by those models tended to be too low in comparison with observations. This letter refines these predictions by taking better account of weak Coulomb collisions for coronal electrons, whose thermodynamic properties determine the ion charge states in the low corona. A perturbative description of nonlocal electron transport is applied to an existing set of wave/turbulence models. The resulting electron velocity distributions in the low corona exhibit mild suprathermal tails characterized by "kappa" exponents between 10 and 25. These suprathermal electrons are found to be sufficiently energetic to enhance the charge states of oxygen ions, while maintaining the same relative trend with wind speed that was found when the distribution was assumed to be Maxwellian. The updated wave/turbulence models are in excellent agreement with solar wind ion composition measurements.

preprint2014arXiv

Turbulence-Driven Coronal Heating and Improvements to Empirical Forecasting of the Solar Wind

Forecasting models of the solar wind often rely on simple parameterizations of the magnetic field that ignore the effects of the full magnetic field geometry. In this paper, we present the results of two solar wind prediction models that consider the full magnetic field profile and include the effects of Alfvén waves on coronal heating and wind acceleration. The one-dimensional MHD code ZEPHYR self-consistently finds solar wind solutions without the need for empirical heating functions. Another 1D code, introduced in this paper (The Efficient Modified-Parker-Equation-Solving Tool, TEMPEST), can act as a smaller, stand-alone code for use in forecasting pipelines. TEMPEST is written in Python and will become a publicly available library of functions that is easy to adapt and expand. We discuss important relations between the magnetic field profile and properties of the solar wind that can be used to independently validate prediction models. ZEPHYR provides the foundation and calibration for TEMPEST, and ultimately we will use these models to predict observations and explain space weather created by the bulk solar wind. We are able to reproduce with both models the general anticorrelation seen in comparisons of observed wind speed at 1 AU and the flux tube expansion factor. There is significantly less spread than comparing the results of the two models than between ZEPHYR and a traditional flux tube expansion relation. We suggest that the new code, TEMPEST, will become a valuable tool in the forecasting of space weather.

preprint2013arXiv

Connecting the Sun's High-Resolution Magnetic Carpet to the Turbulent Heliosphere

The solar wind is connected to the Sun's atmosphere by flux tubes that are rooted in an ever-changing pattern of positive and negative magnetic polarities on the surface. Observations indicate that the magnetic field is filamentary and intermittent across a wide range of spatial scales. However, we do not know to what extent the complex flux tube topology seen near the Sun survives as the wind expands into interplanetary space. In order to study the possible long-distance connections between the corona and the heliosphere, we developed new models of turbulence-driven solar wind acceleration along empirically constrained field lines. We used a potential-field model of the Quiet Sun to trace field lines into the ecliptic plane with unprecedented spatial resolution at their footpoints. For each flux tube, a one-dimensional model was created with an existing wave/turbulence code that solves equations of mass, momentum, and energy conservation from the photosphere to 4 AU. To take account of stream-stream interactions between flux tubes, we used those models as inner boundary conditions for a time-steady MHD description of radial and longitudinal structure in the ecliptic. Corotating stream interactions smear out much of the smallest-scale variability, making it difficult to see how individual flux tubes on granular or supergranular scales can survive out to 1 AU. However, our models help clarify the level of "background" variability with which waves and turbulent eddies should be expected to interact. Also, the modeled fluctuations in magnetic field magnitude were seen to match measured power spectra quite well.

preprint2013arXiv

Constraining a Model of Turbulent Coronal Heating for AU Microscopii with X-Ray, Radio, and Millimeter Observations

Many low-mass pre-main-sequence stars exhibit strong magnetic activity and coronal X-ray emission. Even after the primordial accretion disk has been cleared out, the star's high-energy radiation continues to affect the formation and evolution of dust, planetesimals, and large planets. Young stars with debris disks are thus ideal environments for studying the earliest stages of non-accretion-driven coronae. In this paper we simulate the corona of AU Mic, a nearby active M dwarf with an edge-on debris disk. We apply a self-consistent model of coronal loop heating that was derived from numerical simulations of solar field-line tangling and magnetohydrodynamic turbulence. We also synthesize the modeled star's X-ray luminosity and thermal radio/millimeter continuum emission. A realistic set of parameter choices for AU Mic produces simulated observations that agree with all existing measurements and upper limits. This coronal model thus represents an alternative explanation for a recently discovered ALMA central emission peak that was suggested to be the result of an inner "asteroid belt" within 3 AU of the star. However, it is also possible that the central 1.3 mm peak is caused by a combination of active coronal emission and a bright inner source of dusty debris. Additional observations of this source's spatial extent and spectral energy distribution at millimeter and radio wavelengths will better constrain the relative contributions of the proposed mechanisms.

preprint2013arXiv

Stellar Granulation as the Source of High-Frequency Flicker in Kepler Light Curves

A large fraction of cool, low-mass stars exhibit brightness fluctuations that arise from a combination of convective granulation, acoustic oscillations, magnetic activity, and stellar rotation. Much of the short-timescale variability takes the form of stochastic noise, whose presence may limit the progress of extrasolar planet detection and characterization. In order to lay the groundwork for extracting useful information from these quasi-random signals, we focus on the origin of the granulation-driven component of the variability. We apply existing theoretical scaling relations to predict the star-integrated variability amplitudes for 508 stars with photometric light curves measured by the Kepler mission. We also derive an empirical correction factor that aims to account for the suppression of convection in F-dwarf stars with magnetic activity and shallow convection zones. So that we can make predictions of specific observational quantities, we performed Monte Carlo simulations of granulation light curves using a Lorentzian power spectrum. These simulations allowed us to reproduce the so-called "flicker floor" (i.e., a lower bound in the relationship between the full light-curve range and power in short-timescale fluctuations) that was found in the Kepler data. The Monte Carlo model also enabled us to convert the modeled fluctuation variance into a flicker amplitude directly comparable with observations. When the magnetic suppression factor described above is applied, the model reproduces the observed correlation between stellar surface gravity and flicker amplitude. Observationally validated models like these provide new and complementary evidence for a possible impact of magnetic activity on the properties of near-surface convection.

preprint2012arXiv

Proton, Electron, and Ion Heating in the Fast Solar Wind from Nonlinear Coupling Between Alfvenic and Fast-Mode Turbulence

In the parts of the solar corona and solar wind that experience the fewest Coulomb collisions, the component proton, electron, and heavy ion populations are not in thermal equilibrium with one another. Observed differences in temperatures, outflow speeds, and velocity distribution anisotropies are useful constraints on proposed explanations for how the plasma is heated and accelerated. This paper presents new predictions of the rates of collisionless heating for each particle species, in which the energy input is assumed to come from magnetohydrodynamic (MHD) turbulence. We first created an empirical description of the radial evolution of Alfven, fast-mode, and slow-mode MHD waves. This model provides the total wave power in each mode as a function of distance along an expanding flux tube in the high-speed solar wind. Next we solved a set of cascade advection-diffusion equations that give the time-steady wavenumber spectra at each distance. An approximate term for nonlinear coupling between the Alfven and fast-mode fluctuations is included. For reasonable choices of the parameters, our model contains enough energy transfer from the fast mode to the Alfven mode to excite the high-frequency ion cyclotron resonance. This resonance is efficient at heating protons and other ions in the direction perpendicular to the background magnetic field, and our model predicts heating rates for these species that agree well with both spectroscopic and in situ measurements. Nonetheless, the high-frequency waves comprise only a small part of the total Alfvenic fluctuation spectrum, which remains highly two-dimensional as is observed in interplanetary space.

preprint2011arXiv

Testing a Predictive Theoretical Model for the Mass Loss Rates of Cool Stars

The basic mechanisms responsible for producing winds from cool, late-type stars are still largely unknown. We take inspiration from recent progress in understanding solar wind acceleration to develop a physically motivated model of the time-steady mass loss rates of cool main-sequence stars and evolved giants. This model follows the energy flux of magnetohydrodynamic turbulence from a subsurface convection zone to its eventual dissipation and escape through open magnetic flux tubes. We show how Alfven waves and turbulence can produce winds in either a hot corona or a cool extended chromosphere, and we specify the conditions that determine whether or not coronal heating occurs. These models do not utilize arbitrary normalization factors, but instead predict the mass loss rate directly from a star's fundamental properties. We take account of stellar magnetic activity by extending standard age-activity-rotation indicators to include the evolution of the filling factor of strong photospheric magnetic fields. We compared the predicted mass loss rates with observed values for 47 stars and found significantly better agreement than was obtained from the popular scaling laws of Reimers, Schroeder, and Cuntz. The algorithm used to compute cool-star mass loss rates is provided as a self-contained and efficient computer code. We anticipate that the results from this kind of model can be incorporated straightforwardly into stellar evolution calculations and population synthesis techniques.

preprint2010arXiv

A Model for the Stray Light Contamination of the UVCS Instrument on SOHO

We present a detailed model of stray-light suppression in the spectrometer channels of the Ultraviolet Coronagraph Spectrometer (UVCS) on the SOHO spacecraft. The control of diffracted and scattered stray light from the bright solar disk is one of the most important tasks of a coronagraph. We compute the fractions of light that diffract past the UVCS external occulter and non-specularly pass into the spectrometer slit. The diffracted component of the stray light depends on the finite aperture of the primary mirror and on its figure. The amount of non-specular scattering depends mainly on the micro-roughness of the mirror. For reasonable choices of these quantities, the modeled stray-light fraction agrees well with measurements of stray light made both in the laboratory and during the UVCS mission. The models were constructed for the bright H I Lyman alpha emission line, but they are applicable to other spectral lines as well.

preprint2010arXiv

Extended Coronal Heating and Solar Wind Acceleration Over the Solar Cycle

This paper reviews our growing understanding of the physics behind coronal heating (in open-field regions) and the acceleration of the solar wind. Many new insights have come from the last solar cycle's worth of observations and theoretical work. Measurements of the plasma properties in the extended corona, where the primary solar wind acceleration occurs, have been key to discriminating between competing theories. We describe how UVCS/SOHO measurements of coronal holes and streamers over the last 14 years have provided clues about the detailed kinetic processes that energize both fast and slow wind regions. We also present a brief survey of current ideas involving the coronal source regions of fast and slow wind streams, and how these change over the solar cycle. These source regions are discussed in the context of recent theoretical models (based on Alfven waves and MHD turbulence) that have begun to successfully predict both the heating and acceleration in fast and slow wind regions with essentially no free parameters. Some new results regarding these models - including a quantitative prediction of the lower density and temperature at 1 AU seen during the present solar minimum in comparison to the prior minimum - are also shown.

preprint2009arXiv

An Efficient Approximation of the Coronal Heating Rate for Use in Global Sun-Heliosphere Simulations

The origins of the hot solar corona and the supersonically expanding solar wind are still the subject of debate. A key obstacle in the way of producing realistic simulations of the Sun-heliosphere system is the lack of a physically motivated way of specifying the coronal heating rate. Recent one-dimensional models have been found to reproduce many observed features of the solar wind by assuming the energy comes from Alfven waves that are partially reflected, then dissipated by magnetohydrodynamic turbulence. However, the nonlocal physics of wave reflection has made it difficult to apply these processes to more sophisticated (three-dimensional) models. This paper presents a set of robust approximations to the solutions of the linear Alfven wave reflection equations. A key ingredient to the turbulent heating rate is the ratio of inward to outward wave power, and the approximations developed here allow this to be written explicitly in terms of local plasma properties at any given location. The coronal heating also depends on the frequency spectrum of Alfven waves in the open-field corona, which has not yet been measured directly. A model-based assumption is used here for the spectrum, but the results of future measurements can be incorporated easily. The resulting expression for the coronal heating rate is self-contained, computationally efficient, and applicable directly to global models of the corona and heliosphere. This paper tests and validates the approximations by comparing the results to exact solutions of the wave transport equations in several cases relevant to the fast and slow solar wind.

preprint2009arXiv

Coronal Holes

Coronal holes are the darkest and least active regions of the Sun, as observed both on the solar disk and above the solar limb. Coronal holes are associated with rapidly expanding open magnetic fields and the acceleration of the high-speed solar wind. This paper reviews measurements of the plasma properties in coronal holes and how these measurements are used to reveal details about the physical processes that heat the solar corona and accelerate the solar wind. It is still unknown to what extent the solar wind is fed by flux tubes that remain open (and are energized by footpoint-driven wave-like fluctuations), and to what extent much of the mass and energy is input intermittently from closed loops into the open-field regions. Evidence for both paradigms is summarized in this paper. Special emphasis is also given to spectroscopic and coronagraphic measurements that allow the highly dynamic non-equilibrium evolution of the plasma to be followed as the asymptotic conditions in interplanetary space are established in the extended corona. For example, the importance of kinetic plasma physics and turbulence in coronal holes has been affirmed by surprising measurements from UVCS that heavy ions are heated to hundreds of times the temperatures of protons and electrons. These observations point to specific kinds of collisionless Alfven wave damping (i.e., ion cyclotron resonance), but complete models do not yet exist. Despite our incomplete knowledge of the complex multi-scale plasma physics, however, much progress has been made toward the goal of understanding the mechanisms responsible for producing the observed properties of coronal holes.

preprint2009arXiv

Empirical Constraints on Proton and Electron Heating in the Fast Solar Wind

We analyze measured proton and electron temperatures in the high-speed solar wind in order to calculate the separate rates of heat deposition for protons and electrons. When comparing with other regions of the heliosphere, the fast solar wind has the lowest density and the least frequent Coulomb collisions. This makes the fast wind an optimal testing ground for studies of collisionless kinetic processes associated with the dissipation of plasma turbulence. Data from the Helios and Ulysses plasma instruments were collected to determine mean radial trends in the temperatures and the electron heat conduction flux between 0.29 and 5.4 AU. The derived heating rates apply specifically for these mean plasma properties and not for the full range of measured values around the mean. We found that the protons receive about 60% of the total plasma heating in the inner heliosphere, and that this fraction increases to approximately 80% by the orbit of Jupiter. A major factor affecting the uncertainty in this fraction is the uncertainty in the measured radial gradient of the electron heat conduction flux. The empirically derived partitioning of heat between protons and electrons is in rough agreement with theoretical predictions from a model of linear Vlasov wave damping. For a modeled power spectrum consisting only of Alfvenic fluctuations, the best agreement was found for a distribution of wavenumber vectors that evolves toward isotropy as distance increases.

preprint2009arXiv

Testing Models of Accretion-driven Coronal Heating and Stellar Wind Acceleration for T Tauri Stars

Classical T Tauri stars are pre-main-sequence objects that undergo simultaneous accretion, wind outflow, and coronal X-ray emission. The impact of plasma on the stellar surface from magnetospheric accretion streams is likely to be a dominant source of energy and momentum in the upper atmospheres of these stars. This paper presents a set of models for the dynamics and heating of three distinct regions on T Tauri stars that are affected by accretion: (1) the shocked plasmas directly beneath the magnetospheric accretion streams, (2) stellar winds that are accelerated along open magnetic flux tubes, and (3) closed magnetic loops that resemble the Sun's coronal active regions. For the loops, a self-consistent model of coronal heating was derived from numerical simulations of solar field-line tangling and turbulent dissipation. Individual models are constructed for the properties of 14 well-observed stars in the Taurus-Auriga star-forming region. Predictions for the wind mass loss rates are, on average, slightly lower than the observations, which suggests that disk winds or X-winds may also contribute to the measured outflows. For some of the stars, however, the modeled stellar winds do appear to contribute significantly to the measured mass fluxes. Predictions for X-ray luminosities from the shocks and loops are in general agreement with existing observations. The stars with the highest accretion rates tend to have X-ray luminosities dominated by the high-temperature (5-10 MK) loops. The X-ray luminosities for the stars having lower accretion rates are dominated by the cooler accretion shocks.

preprint2008arXiv

Ion Temperatures in the Low Solar Corona: Polar Coronal Holes at Solar Minimum

In the present work we use a deep-exposure spectrum taken by the SUMER spectrometer in a polar coronal hole in 1996 to measure the ion temperatures of a large number of ions at many different heights above the limb between 0.03 and 0.17 solar radii. We find that the measured ion temperatures are almost always larger than the electron temperatures and exhibit a non-monotonic dependence on the charge-to-mass ratio. We use these measurements to provide empirical constraints to a theoretical model of ion heating and acceleration based on gradually replenished ion-cyclotron waves. We compare the wave power required to heat the ions to the observed levels to a prediction based on a model of anisotropic magnetohydrodynamic turbulence. We find that the empirical heating model and the turbulent cascade model agree with one another, and explain the measured ion temperatures, for charge-to-mass ratios smaller than about 0.25. However, ions with charge-to-mass ratios exceeding 0.25 disagree with the model; the wave power they require to be heated to the measured ion temperatures shows an increase with charge-to-mass ratio (i.e., with increasing frequency) that cannot be explained by a traditional cascade model. We discuss possible additional processes that might be responsible for the inferred surplus of wave power.