BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. S. Wheatland

M. S. Wheatland appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-phphysics.space-phApplicationsastro-ph.HEphysics.class-phphysics.ed-phphysics.plasm-phphysics.pop-ph

Catalog footprint

What is connected

17works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

17 published item(s)

preprint2019arXiv

Mapping Magnetic Field Lines for an Accelerating Solar Wind

Mapping of magnetic field lines is important for studies of the solar wind and the sources and propagation of energetic particles between the Sun and observers. A recently developed mapping approach is generalised to use a more advanced solar wind model that includes the effects of solar wind acceleration, non-radial intrinsic magnetic fields and flows at the source surface/inner boundary, and conservation of angular momentum. The field lines are mapped by stepping along local magnetic field $\bm{B}$ and via a Runge-Kutta algorithm, leading to essentially identical maps. The new model's maps for Carrington rotation CR 1895 near solar minimum (19 April to 15 May 1995) and a solar rotation between CR 2145 and CR 2146 near solar maximum (14 January to 9 February 2014) are compared with the published maps for a constant solar wind model. The two maps are very similar on a large scale near both solar minimum and solar maximum, meaning that the field line orientations, winding angles, and connectivity generally agree very well. However, close inspection shows that the field lines have notable small-scale structural differences. An interpretation is that inclusion of the acceleration and intrinsic azimuthal velocity has significant effects on the local structure of the magnetic field lines. Interestingly, the field lines are more azimuthal for the accelerating solar wind model for both intervals. In addition, predictions for the pitch angle distributions (PADs) for suprathermal electrons agree at the $90$ -- $95\%$ level with observations for both solar wind models for both intervals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Is Cyclotron Maser Emission in Solar Flares Driven by a Horseshoe Distribution?

Since the early 1980s, decimetric spike bursts have been attributed to electron cyclotron maser emission (ECME) by the electrons that produce hard X-ray bursts as they precipitate into the chromosphere in the impulsive phase of a solar flare. Spike bursts are regarded as analogous to the auroral kilometric radiation (AKR), which is associated with the precipitation of auroral electrons in a geomagnetic substorm. Originally, a loss-cone-driven version of ECME, developed for AKR, was applied to spike bursts, but it is now widely accepted that a different, horseshoe-driven, version of EMCE applies to AKR. We explore the implications of the assumption that horseshoe-driven ECME also applies to spike bursts. We develop a 1D model for the acceleration of the electrons by a parallel electric field, and show that under plausible assumptions it leads to a horseshoe distribution of electrons in a solar flare. A second requirement for horseshoe-driven ECME is an extremely low plasma density, referred to as a density cavity. We argue that a coronal density cavity should develop in association with a hard X-ray burst, and that such a density cavity can overcome a long-standing problem with the escape of ECME through the second-harmonic absorption layer. Both the horseshoe distribution and the associated coronal density cavity are highly localized, and could not be resolved in the statistically large number of local precipitation regions needed to explain a hard X-ray burst. The model highlights the "number problem" in the supply of the electrons needed to explain a hard X-ray burst.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Estimating electric current densities in solar active regions

Electric currents in solar active regions are thought to provide the energy released via magnetic reconnection in solar flares. Vertical electric current densities $J_z$ at the photosphere may be estimated from vector magnetogram data, subject to substantial uncertainties. The values provide boundary conditions for nonlinear force- free modelling of active region magnetic fields. A method is presented for estimating values of $J_z$ taking into account uncertainties in vector magnetogram field values, and minimizing $J_z^2$ across the active region. The method is demonstrated using the boundary values of the field for a force-free twisted bipole, with the addition of noise at randomly chosen locations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modeling sunspot and starspot decay by turbulent erosion

Disintegration of sunspots (and starspots) by fluxtube erosion, originally proposed by Simon and Leighton, is considered. A moving boundary problem is formulated for a nonlinear diffusion equation that describes the sunspot magnetic field profile. Explicit expressions for the sunspot decay rate and lifetime by turbulent erosion are derived analytically and verified numerically. A parabolic decay law for the sunspot area is obtained. For moderate sunspot magnetic field strengths, the predicted decay rate agrees with the results obtained by Petrovay and Moreno-Insertis. The new analytical and numerical solutions significantly improve the quantitative description of sunspot and starspot decay by turbulent erosion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Influence of Spatial Resolution on Nonlinear Force-Free Modeling

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bulk Energization of Electrons in Solar Flares by Alfvén Waves

Bulk energization of electrons to $10\,-\,20\,$keV in solar flares is attributed to dissipation of Alfvén waves that transport energy and potential downward to an acceleration region near the chromosphere. The acceleration involves the parallel electric field that develops in the limit of inertial Alfvén waves (IAWs). A two-potential model for IAWs is used to relate the parallel potential to the cross-field potential transported by the waves. We identify a maximum parallel potential in terms of a maximum current density that corresponds to the threshold for the onset of anomalous resistivity. This maximum is of order $10\,$kV when the threshold is that for the Buneman instability. We argue that this restricts the cross-field potential in an Alfvén wave to about $10\,$kV. Effective dissipation requires a large number of up- and down-current paths associated with multiple Alfvén waves. The electron acceleration occurs in localized, transient, anomalously-conducting regions (LTACRs) and is associated with the parallel electric field determined by Ohm's law with an anomalous resistivity. We introduce an idealized model in which the LTACRs are (upward-)current sheets, a few skin depths in thickness, separated by much-larger regions of weaker return current. We show that this model can account semi-quantitatively for bulk energization.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear force-free modeling of the corona in spherical coordinates

We present a code for solving the nonlinear force-free equations in spherical polar geometry, with the motivation of modeling the magnetic field in the corona. The code is an implementation of the Grad-Rubin method. Our method is applicable to a spherical domain of arbitrary angular size. The implementation is based on a global spectral representation for the magnetic field which makes no explicit assumptions about the form of the magnetic field at the transverse boundaries of the domain. We apply the code to a bipolar test case with analytic boundary conditions, and we demonstrate the convergence of the Grad-Rubin method, and the self-consistency of the resulting numerical solution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transfer of Energy, Potential, and Current by Alfvén Waves in Solar Flares

Alfvén waves play three related roles in the impulsive phase of a solar flare: they transport energy from a generator region to an acceleration region; they map the cross-field potential (associated with the driven energy release) from the generator region onto the acceleration region; and within the acceleration region they damp by setting up a parallel electric field that accelerates electrons and transfers the wave energy to them. The Alfvén waves may also be regarded as setting up new closed current loops, with field-aligned currents that close across field lines at boundaries. A model is developed for large-amplitude Alfvén waves that shows how Alfvén waves play these roles in solar flares. A picket-fence structure for the current flow is incorporated into the model to account for the "number problem" and the energy of the accelerated electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Magnetostatic Grad-Rubin Code for Coronal Magnetic Field Extrapolations

The coronal magnetic field cannot be directly observed, but in principle it can be reconstructed from the comparatively well observed photospheric magnetic field. A popular approach uses a nonlinear force-free model. Non-magnetic forces at the photosphere are significant meaning the photospheric data are inconsistent with the force-free model, and this causes problems with the modeling (De Rosa et al., Astrophys. J. 696, 1780, 2009). In this paper we present a numerical implementation of the Grad-Rubin method for reconstructing the coronal magnetic field using a magnetostatic model. This model includes a pressure force and a non-zero magnetic Lorentz force. We demonstrate our implementation on a simple analytic test case and obtain the speed and numerical error scaling as a function of the grid size.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling a falling slinky

A slinky is an example of a tension spring: in an unstretched state a slinky is collapsed, with turns touching, and a finite tension is required to separate the turns from this state. If a slinky is suspended from its top and stretched under gravity and then released, the bottom of the slinky does not begin to fall until the top section of the slinky, which collapses turn by turn from the top, collides with the bottom. The total collapse time t_c (typically ~0.3 s for real slinkies) corresponds to the time required for a wave front to propagate down the slinky to communicate the release of the top end. We present a modification to an existing model for a falling tension spring (Calkin 1993) and apply it to data from filmed drops of two real slinkies. The modification of the model is the inclusion of a finite time for collapse of the turns of the slinky behind the collapse front propagating down the slinky during the fall. The new finite-collapse time model achieves a good qualitative fit to the observed positions of the top of the real slinkies during the measured drops. The spring constant k for each slinky is taken to be a free parameter in the model. The best-fit model values for k for each slinky are approximately consistent with values obtained from measured periods of oscillation of the slinkies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The free energy of NOAA active region AR 11029

The NOAA active region AR 11029 was a small but highly active sunspot region which produced 73 GOES soft X-ray flares. The flares appear to show a departure from the well known power-law frequency-size distribution. Specifically, too few GOES C-class and no M-class flares were observed by comparison with a power-law distribution (Wheatland in Astrophys. J. 710, 1324, 2010). This was conjectured to be due to the region having insufficient magnetic energy to power large events. We construct nonlinear force-free extrapolations of the coronal magnetic field of active region AR 11029 using data taken on 24 October by the SOLIS Vector-SpectroMagnetograph (SOLIS/VSM), and data taken on 27 October by the Hinode Solar Optical Telescope SpectroPolarimeter (Hinode/SP). Force-free modeling with photospheric magnetogram data encounters problems because the magnetogram data are inconsistent with a force-free model, and we employ a recently developed `self-consistency' procedure which addresses this and accommodates uncertainties in the boundary data (Wheatland and Regnier in Astrophys. J. 700, L88, 2009). We calculate the total energy and free energy of the self-consistent solution and find that the free energy was 4x10^29 erg on 24 October, and 7x10^31 erg on 27 October. An order of magnitude scaling between RHESSI non-thermal energy and GOES peak X-ray flux is established from a sample of flares from the literature and is used to estimate flare energies from observed GOES peak X-ray flux. Based on the scaling, we conclude that the estimated free energy of AR 11029 on 27 October when the flaring rate peaked is sufficient to power M-class or X-class flares, and hence the modeling does not appear to support the hypothesis that the absence of large flares is due to the region having limited energy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Achieving Self-Consistent Nonlinear Force-free Modeling of Solar Active Regions

A nonlinear force-free solution is constructed for the coronal magnetic field in NOAA solar active region AR 10953 based on a photospheric vector magnetogram derived from Hinode satellite observations on 30 April 2007, taking into account uncertainties in the boundary data and using improved methods for merging multiple-instrument data. The solution demonstrates the "self-consistency" procedure of Wheatland & Regnier (2009), for the first time including uncertainties. The self-consistency procedure addresses the problem that photospheric vector magnetogram data are inconsistent with the force-free model, and in particular that the boundary conditions on vertical electric current density are over-specified and permit the construction of two different nonlinear force-free solutions. The procedure modifies the boundary conditions on current density during a sequence of cycles until the two nonlinear force-free solutions agree. It hence constructs an accurate single solution to the force-free model, with boundary values close, but not matched exactly, to the vector magnetogram data. The inclusion of uncertainties preserves the boundary conditions more closely at points with smaller uncertainties. The self-consistent solution obtained for active region AR 10953 is significantly non-potential, with magnetic energy E/E_0 = 1.08, where E_0 is the energy of the reference potential (current-free) magnetic field. The self-consistent solution is shown to be robust against changes in the details of the construction of the two force-free models at each cycle. This suggests that reliable nonlinear force-free modeling of active regions is possible if uncertainties in vector magnetogram boundary data are included.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evidence for departure from a power-law flare size distribution for a small solar active region

Active region 11029 was a small, highly flare-productive solar active region observed at a time of extremely low solar activity. The region produced only small flares: the largest of the $>70$ Geostationary Observational Environmental Satellite (GOES) events for the region has a peak 1--$8Å$ flux of $2.2\times 10^{-6} {\rm W} {\rm m}^{-2}$ (GOES C2.2). The background-subtracted GOES peak-flux distribution suggests departure from power-law behavior above $10^{-6} {\rm W} {\rm m}^{-2}$, and a Bayesian model comparison strongly favors a power-law plus rollover model for the distribution over a simple power-law model. The departure from the power law is attributed to this small active region having a finite amount of energy. The rate of flaring in the region varies with time, becoming very high for two days coinciding with the onset of an increase in complexity of the photospheric magnetic field. The observed waiting-time distribution for events is consistent with a piecewise-constant Poisson model. These results present challenges for models of flare statistics and of energy balance in solar active regions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Time-dependent Stochastic Modeling of Solar Active Region Energy

A time-dependent model for the energy of a flaring solar active region is presented based on a stochastic jump-transition model (Wheatland and Glukhov 1998; Wheatland 2008; Wheatland 2009). The magnetic free energy of the model active region varies in time due to a prescribed (deterministic) rate of energy input and prescribed (random) flare jumps downwards in energy. The model has been shown to reproduce observed flare statistics, for specific time-independent choices for the energy input and flare transition rates. However, many solar active regions exhibit time variation in flare productivity, as exemplified by NOAA active region AR 11029 (Wheatland 2010). In this case a time-dependent model is needed. Time variation is incorporated for two cases: 1. a step change in the rates of flare jumps; and 2. a step change in the rate of energy supply to the system. Analytic arguments are presented describing the qualitative behavior of the system in the two cases. In each case the system adjusts by shifting to a new stationary state over a relaxation time which is estimated analytically. The new model retains flare-like event statistics. In each case the frequency-energy distribution is a power law for flare energies less than a time-dependent rollover set by the largest energy the system is likely to attain at a given time. For Case 1, the model exhibits a double exponential waiting-time distribution, corresponding to flaring at a constant mean rate during two intervals (before and after the step change), if the average energy of the system is large. For Case 2 the waiting-time distribution is a simple exponential, again provided the average energy of the system is large. Monte Carlo simulations of Case~1 are presented which confirm the analytic estimates. The simulation results provide a qualitative model for observed flare statistics in active region AR 11029.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Brightness and Waiting-time Distributions of a Type III Radio Storm observed by STEREO/WAVES

Type III solar radio storms, observed at frequencies below approximately 16 MHz by space borne radio experiments, correspond to the quasi-continuous, bursty emission of electron beams onto open field lines above active regions. The mechanisms by which a storm can persist in some cases for more than a solar rotation whilst exhibiting considerable radio activity are poorly understood. To address this issue, the statistical properties of a type III storm observed by the STEREO/WAVES radio experiment are presented, examining both the brightness distribution and (for the first time) the waiting-time distribution. Single power law behavior is observed in the number distribution as a function of brightness; the power law index is approximately 2.1 and is largely independent of frequency. The waiting-time distribution is found to be consistent with a piecewise-constant Poisson process. This indicates that during the storm individual type III bursts occur independently and suggests that the storm dynamics are consistent with avalanche type behavior in the underlying active region.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A Statistical Solar Flare Forecast Method

A Bayesian approach to solar flare prediction has been developed, which uses only the event statistics of flares already observed. The method is simple, objective, and makes few ad hoc assumptions. It is argued that this approach should be used to provide a baseline prediction for certain space weather purposes, upon which other methods, incorporating additional information, can improve. A practical implementation of the method for whole-Sun prediction of Geostationary Observational Environment Satellite (GOES) events is described in detail, and is demonstrated for 4 November 2003, the day of the largest recorded GOES flare. A test of the method is described based on the historical record of GOES events (1975-2003), and a detailed comparison is made with US National Oceanic and Atmospheric Administration (NOAA) predictions for 1987-2003. Although the NOAA forecasts incorporate a variety of other information, the present method out-performs the NOAA method in predicting mean numbers of event days, for both M-X and X events. Skill scores and other measures show that the present method is slightly less accurate at predicting M-X events than the NOAA method, but substantially more accurate at predicting X events, which are important contributors to space weather.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The coronal mass ejection waiting-time distribution

The distribution of times $Δt$ between coronal mass ejections (CMEs) in the Large Angle and Spectrometric Coronagraph (LASCO) CME catalog for the years 1996-2001 is examined. The distribution exhibits a power-law tail $\propto (Δt)^γ$ with an index $γ\approx -2.36\pm 0.11$ for large waiting times ($Δt>10 {\rm hours}$). The power-law index of the waiting-time distribution varies with the solar cycle: for the years 1996-1998 (a period of low activity), the power-law index is $γ\approx-1.86\pm 0.14$, and for the years 1999-2001 (a period of higher activity), the index is $γ\approx-2.98\pm 0.20$. The observed CME waiting-time distribution, and its variation with the cycle, may be understood in terms of CMEs occurring as a time-dependent Poisson process. The CME waiting-time distribution is compared with that for greater than C1 class solar flares in the Geostationary Operational Environmental Satellite (GOES) catalog for the same years. The flare and CME waiting-time distributions exhibit power-law tails with very similar indices and time variation.

0 citations0 reviewsNo code linkedOpen access