NASA's WB-57 High Altitude Research Program provides a deployable, mobile, stratospheric platform for scientific research. Airborne platforms are of particular value for making coronal observations during total solar eclipses because of their ability both to follow the Moon's shadow and to get above most of the atmospheric airmass that can interfere with astronomical observations. We used the 2017 Aug 21 eclipse as a pathfinding mission for high-altitude airborne solar astronomy, using the existing high-speed visible-light and near-/mid-wave infrared imaging suite mounted in the WB-57 nose cone. In this paper, we describe the aircraft, the instrument, and the 2017 mission; operations and data acquisition; and preliminary analysis of data quality from the existing instrument suite. We describe benefits and technical limitations of this platform for solar and other astronomical observations. We present a preliminary analysis of the visible-light data quality and discuss the limiting factors that must be overcome with future instrumentation. We conclude with a discussion of lessons learned from this pathfinding mission and prospects for future research at upcoming eclipses, as well as an evaluation of the capabilities of the WB-57 platform for future solar astronomy and general astronomical observation.

The Spectral Imaging of the Coronal Environment (SPICE) instrument is a high-resolution imaging spectrometer operating at extreme ultraviolet (EUV) wavelengths. In this paper, we present the concept, design, and pre-launch performance of this facility instrument on the ESA/NASA Solar Orbiter mission. The goal of this paper is to give prospective users a better understanding of the possible types of observations, the data acquisition, and the sources that contribute to the instrument's signal. The paper discusses the science objectives, with a focus on the SPICE-specific aspects, before presenting the instrument's design, including optical, mechanical, thermal, and electronics aspects. This is followed by a characterisation and calibration of the instrument's performance. The paper concludes with descriptions of the operations concept and data processing. The performance measurements of the various instrument parameters meet the requirements derived from the mission's science objectives. The SPICE instrument is ready to perform measurements that will provide vital contributions to the scientific success of the Solar Orbiter mission.

Jets (transient/collimated plasma ejections) occur frequently throughout the solar corona and contribute mass/energy to the corona and solar wind. By combining numerical simulations and high-resolution observations, we have made substantial progress recently on determining the energy buildup and release processes in these jets. Here we describe a study of 27 equatorial coronal-hole jets using Solar Dynamics Observatory/AIA and HMI observations on 2013 June 27-28 and 2014 January 8-10. Out of 27 jets, 18 (67%) are associated with mini-filament ejections; the other 9 (33%) do not show mini-filament eruptions but do exhibit mini-flare arcades and other eruptive signatures. This indicates that every jet in our sample involved a filament-channel eruption. From the complete set of events, 6 jets (22%) are apparently associated with tiny flux-cancellation events at the polarity inversion line, and 2 jets (7%) are associated with sympathetic eruptions of filaments from neighboring bright points. Potential-field extrapolations of the source-region photospheric magnetic fields reveal that all jets originated in the fan-spine topology of an embedded bipole associated with an extreme ultraviolet coronal bright point. Hence, all our jets are in agreement with the breakout model of solar eruptions. We present selected examples and discuss the implications for the jet energy build-up and initiation mechanisms.

Small, impulsive jets commonly occur throughout the solar corona, but are especially visible in coronal holes. Evidence is mounting that jets are part of a continuum of eruptions that extends to much larger coronal mass ejections and eruptive flares. Because coronal-hole jets originate in relatively simple magnetic structures, they offer an ideal testbed for theories of energy buildup and release in the full range of solar eruptions. We analyzed an equatorial coronal-hole jet observed by SDO/AIA on 09 January 2014, in which the magnetic-field structure was consistent with the embedded-bipole topology that we identified and modeled previously as an origin of coronal jets. In addition, this event contained a mini-filament, which led to important insights into the energy storage and release mechanisms. SDO/HMI magnetograms revealed footpoint motions in the primary minority-polarity region at the eruption site, but show negligible flux emergence or cancellation for at least 16 hours before the eruption. Therefore, the free energy powering this jet probably came from magnetic shear concentrated at the polarity inversion line within the embedded bipole. We find that the observed activity sequence and its interpretation closely match the predictions of the breakout jet model, strongly supporting the hypothesis that the breakout model can explain solar eruptions on a wide range of scales.

I present a family of algorithms to reduce noise in astrophysical im- ages and image sequences, preserving more information from the original data than is retained by conventional techniques. The family uses locally adaptive filters ("noise gates") in the Fourier domain, to separate coherent image structure from background noise based on the statistics of local neighborhoods in the image. Processing of solar data limited by simple shot noise or by additive noise reveals image structure not easily visible in the originals, preserves photometry of observable features, and reduces shot noise by a factor of 10 or more with little to no apparent loss of resolution, revealing faint features that were either not directly discernible or not sufficiently strongly detected for quantitative analysis. The method works best on image sequences containing related subjects, for example movies of solar evolution, but is also applicable to single images provided that there are enough pixels. The adaptive filter uses the statistical properties of noise and of local neighborhoods in the data, to discriminate between coherent features and incoherent noise without reference to the specific shape or evolution of the those features. The technique can potentially be modified in a straightforward way to exploit additional a priori knowledge about the functional form of the noise.

We examine the long-term time evolution (1965-2015) of the relationships between solar wind proton temperature (Tp) and speed (Vp) and between the proton density (np) and speed using OMNI solar wind observations taken near Earth. We find a long-term decrease in the proton temperature-speed (Tp-Vp) slope that lasted from 1972 to 2010, but has been trending upward since 2010. Since the solar wind proton density-speed (np-Vp) relationship is not linear like the Tp-Vp relationship, we perform power law fits for np-Vp. The exponent (steepness in the np-Vp relationship) is correlated with the solar cycle. This exponent has a stronger correlation with current sheet tilt angle than with sunspot number because the sunspot number maxima vary considerably from cycle to cycle and the tilt angle maxima do not. To understand this finding, we examined the average np for different speed ranges, and found that for the slow wind np is highly correlated with the sunspot number with a lag of ~4 years. The fast wind np variation was less, but in phase with the cycle. This phase difference may contribute to the np-Vp exponent correlation with the solar cycle. These long-term trends are important since empirical formulas based on fits to Tp and Vp data are commonly used to identify ICMEs, but these formulas do not include any time dependence. Changes in the solar wind density over a solar cycle will create corresponding changes in the near Earth space environment and the overall extent of the heliosphere.

Above the top of the solar corona, the young slow solar wind transitions from low-beta, magnetically structured flow dominated by radial structures, to high-beta, less structured flow dominated by hydrodynamics. This transition, long inferred via theory, is readily apparent in the sky region close to 10 degrees from the Sun, in processed, background-subtracted solar wind images. We present image sequences collected by the STEREO/HI1 instrument in 2008 Dec, covering apparent distances from approximately 4 to 24 degrees from the center of the Sun and spanning this transition in large-scale morphology of the wind. We describe the observation and novel techniques to extract evolving image structure from the images, and we use those data and techniques to present and quantify the clear textural shift in the apparent structure of the corona and solar wind in this altitude range. We demonstrate that the change in apparent texture is due both to anomalous fading of the radial striae that characterize the corona, and to anomalous relative brightening of locally dense puffs of solar wind that we term "flocculae;" and show that these phenomena are inconsistent with smooth radial flow, but consistent with onset of hydrodynamic or MHD instabilities leading to a turbulent cascade in the young solar wind.

Chromospheric rapid blueshifted excursions (RBEs) are suggested to be the disk counterparts of type II spicules at the limb and believed to contribute to the coronal heating process. Previous identification of RBEs was mainly based on feature detection using Dopplergrams. In this paper, we study RBEs on 2011 October 21 in a very quiet region at the disk center, which were observed with the high-cadence imaging spectroscopy of the Ca II 8542 A line from the Interferometric Bidimensional Spectrometer (IBIS). By using an automatic spectral analysis algorithm, a total of 98 RBEs are identified during a 11 minute period. Most of these RBEs have either a round or elongated shape, with an average area of 1.2 arcsec^2. The detailed temporal evolution of spectra from IBIS makes possible a quantitative determination of the velocity (~16 km/s) and acceleration (~400 m/s^2) of Ca II 8542 RBEs, and reveal an additional deceleration (~-160 m/s^2) phase that usually follows the initial acceleration. In addition, we also investigate the association of RBEs with the concomitant photospheric magnetic field evolution, using coordinated high-resolution and high-sensitivity magnetograms made by Hinode. Clear examples are found where RBEs appear to be associated with the preceding magnetic flux emergence and/or the subsequent flux cancellation. However, a further analysis with the aid of the Southwest Automatic Magnetic Identification Suite does not yield a significant statistical association between these RBEs and magnetic field evolution. We discuss the implications of our results in the context of understanding the driving mechanism of RBEs.

The tenuous supersonic solar wind that streams from the top of the corona passes through a natural boundary -- the Alfv\'en surface -- that marks the causal disconnection of individual packets of plasma and magnetic flux from the Sun itself. The Alfv\'en surface is the locus where the radial motion of the accelerating solar wind passes the radial Alfv\'en speed, and therefore any displacement of material cannot carry information back down into the corona. It is thus the natural outer boundary of the solar corona, and the inner boundary of interplanetary space. Using a new and unique motion analysis to separate inbound and outbound motions in synoptic visible-light image sequences from the COR2 coronagraph on board the STEREO-A spacecraft, we have identified inbound wave motion in the outer corona beyond 6 solar radii for the first time, and used it to determine that the Alfv\'en surface is at least 12.5 solar radii from the Sun over the polar coronal holes and 17 solar radii in the streamer belt, well beyond the distance planned for NASA's upcoming Solar Probe Plus mission. To our knowledge this is the first measurement of inbound waves in the outer solar corona, and the first direct measurement of lower bounds for the Alfv\'en surface.

We explore the nature of the small-scale solar dynamo by tracking magnetic features. We investigate two previously-explored categories of the small-scale solar dynamo: shallow and deep. Recent modeling work on the shallow dynamo has produced a number of scenarios for how a strong network concentration can influence the formation and polarity of nearby small-scale magnetic features. These scenarios have measurable signatures, which we test for here using magnetograms from the Narrowband Filter Imager (NFI) on Hinode. We find no statistical tendency for newly-formed magnetic features to cluster around or away from network concentrations, nor do we find any statistical relationship between their polarities. We conclude that there is no shallow or "surface" dynamo on the spatial scales observable by Hinode/NFI. In light of these results, we offer a scenario in which the sub-surface field in a deep solar dynamo is stretched and distorted via turbulence, allowing the field to emerge at random locations on the photosphere.

Observing the Sun at high time and spatial scales is a step towards understanding the finest and fundamental scales of heating events in the solar corona. The Hi-C instrument has provided the highest spatial and temporal resolution images of the solar corona in the EUV wavelength range to date. Hi-C observed an active region on 11 July 2012, which exhibits several interesting features in the EUV line at 193\AA: one of them is the existence of short, small brightenings ``sparkling" at the edge of the active region; we call these EUV Bright Dots (EBDs). Individual EBDs have a characteristic duration of 25s with a characteristic length of 680 km. These brightenings are not fully resolved by the SDO/AIA instrument at the same wavelength, however, they can be identified with respect to the Hi-C location of the EBDs. In addition, EBDs are seen in other chromospheric/coronal channels of SDO/AIA suggesting a temperature between 0.5 and 1.5 MK. Based on their frequency in the Hi-C time series, we define four different categories of EBDs: single peak, double peak, long duration, and bursty EBDs. Based on a potential field extrapolation from an SDO/HMI magnetogram, the EBDs appear at the footpoints of large-scale trans-equatorial coronal loops. The Hi-C observations provide the first evidence of small-scale EUV heating events at the base of these coronal loops, which have a free magnetic energy of the order of 10$^{26}$ erg.

The removal of magnetic flux from the quiet-sun photosphere is important for maintaining the statistical steady-state of the magnetic field there, for determining the magnetic flux budget of the Sun, and for estimating the rate of energy injected into the upper solar atmosphere. Magnetic feature death is a measurable proxy for the removal of detectable flux. We used the SWAMIS feature tracking code to understand how nearly 20000 detected magnetic features die in an hour-long sequence of Hinode/SOT/NFI magnetograms of a region of quiet Sun. Of the feature deaths that remove visible magnetic flux from the photosphere, the vast majority do so by a process that merely disperses the previously-detected flux so that it is too small and too weak to be detected. The behavior of the ensemble average of these dispersals is not consistent with a model of simple planar diffusion, suggesting that the dispersal is constrained by the evolving photospheric velocity field. We introduce the concept of the partial lifetime of magnetic features, and show that the partial lifetime due to Cancellation of magnetic flux, 22 h, is 3 times slower than previous measurements of the flux turnover time. This indicates that prior feature-based estimates of the flux replacement time may be too short, in contrast with the tendency for this quantity to decrease as resolution and instrumentation have improved. This suggests that dispersal of flux to smaller scales is more important for the replacement of magnetic fields in the quiet Sun than observed bipolar cancellation. We conclude that processes on spatial scales smaller than those visible to Hinode dominate the processes of flux emergence and cancellation, and therefore also the quantity of magnetic flux that threads the photosphere.

Plasma flows within prominences/filaments have been observed for many years and hold valuable clues concerning the mass and energy balance within these structures. Previous observations of these flows primarily come from H-alpha and cool EUV lines (e.g., 304A) where estimates of the size of the prominence threads has been limited by the resolution of the available instrumentation. Evidence of `counter-steaming' flows has previously been inferred from these cool plasma observations but now, for the first time, these flows have been directly imaged along fundamental filament threads within the million degree corona (at 193A). In this work we present observations of an active region filament observed with Hi-C that exhibits anti-parallel flows along adjacent filament threads. Complementary data from SDO/AIA and HMI are presented. The ultra-high spatial and temporal resolution of Hi-C allow the anti-parallel flow velocities to be measured (70-80 km/s) and gives an indication of the resolvable thickness of the individual strands (0.8'' +/- 0.1''). The temperature distribution of the plasma flows was estimated to be log T(K) = 5.45 +/- 0.10 using EM loci analysis. We find that SDO/AIA cannot clearly observe these anti-parallel flows nor measure their velocity or thread width due to its larger pixel size. We suggest that anti-parallel/counter-streaming flows are likely commonplace within all filaments and are currently not observed in EUV due to current instrument spatial resolution.

The High-resolution Coronal Imager (Hi-C) has provided Fe XII 193A images of the upper transition region moss at an unprecedented spatial (~0.3-0.4 arcsec) and temporal (5.5s) resolution. The Hi-C observations show in some moss regions variability on timescales down to ~15s, significantly shorter than the minute scale variability typically found in previous observations of moss, therefore challenging the conclusion of moss being heated in a mostly steady manner. These rapid variability moss regions are located at the footpoints of bright hot coronal loops observed by SDO/AIA in the 94A channel, and by Hinode/XRT. The configuration of these loops is highly dynamic, and suggestive of slipping reconnection. We interpret these events as signatures of heating events associated with reconnection occurring in the overlying hot coronal loops, i.e., coronal nanoflares. We estimate the order of magnitude of the energy in these events to be of at least a few $10^{23}rg, also supporting the nanoflare scenario. These Hi-C observations suggest that future observations at comparable high spatial and temporal resolution, with more extensive temperature coverage are required to determine the exact characteristics of the heating mechanism(s).

Using simultaneous observations of the MDI and EIT instruments on board the SoHO spacecraft, we determined in flight the plate scale of the EIT. We found a value of 2.629+-0.001 arc seconds per pixel, in fair agreement with the 2.627+-0.001 arc seconds per pixel value deduced from recent laboratory measurements of the focal length, and much higher by 7 sigma than the 2.622 arc seconds per pixel value of the pre-flight calibrations. The plate scale is found to be constant across the field of view, confirming the negligible distortion level predicted by the theoretical models of the EIT. Furthermore, the 2 sigma difference between our results and the latest laboratory measurements, although statistically small, may confirm a recent work suggesting that the solar photospheric radius may be 0.5 Mm lower than the classically adopted value of 695.99 Mm.

The solar corona and heliosphere are visible via sunlight that is Thomson-scattered off of free electrons, yielding a radiance against the celestial sphere. In this second part of a three-article series, we discuss linear polarization of this scattered light parallel and perpendicular to the plane of scatter in the context of heliopheric imaging far from the Sun. The difference between these two radiances, (pB), varies quite differently with scattering angle, compared to the sum that would be detected in unpolarized light (B). The difference between these two quantities has long been used in a coronagraphic context for background subtraction and to extract some three-dimensional information about the corona; we explore how these effects differ in the wider-field heliospheric imaging case where small-angle approximations do not apply. We develop an appropriately-simplified theory of polarized Thomson scattering in the heliosphere, discuss signal-to-noise considerations, invert the scattering equations analytically to solve the three dimensional object location problem for small objects, discuss exploiting polarization for background subtraction, and generate simple forward models of several classes of heliospheric feature. We conclude that pB measurements of heliospheric material are much more localized to the Thomson surface than are B measurements, that the ratio pB/B can be used to track solar wind features in three dimensions for scientific and space weather applications better in the heliosphere than corona; and that, by providing an independent measurement of background signal, pB measurements may be used to reduce the effect of background radiances including the stably polarized zodiacal light.

Disconnection of open magnetic flux by reconnection is required to balance the injection of open flux by CMEs and other eruptive events. Making use of recent advances in heliospheric background subtraction, we have imaged many abrupt disconnection events. These events produce dense plasma clouds whose distinctie shape can now be traced from the corona across the inner solar system via heliospheric imaging. The morphology of each initial event is characteristic of magnetic reconnection across a current sheet, and the newly-disconnected flux takes the form of a "U"-shaped loop that moves outward, accreting coronal and solar wind material. We analyzed one such event on 2008 December 18 as it formed and accelerated at 20 m/s^2 to 320 km/s, expanding self-similarly until it exited our field of view 1.2 AU from the Sun. From acceleration and photometric mass estimates we derive the coronal magnetic field strength to be 8uT, 6 Rs above the photosphere, and the entrained flux to be 1.6x10^11 Wb (1.6x10^19 Mx). We model the feature's propagation by balancing inferred magnetic tension force against accretion drag. This model is consistent with the feature's behavior and accepted solar wind parameters. By counting events over a 36 day window, we estimate a global event rate of 1/day and a global solar minimum unsigned flux disconnection rate of 6x10^13 Wb/y (6x10^21 Mx/y) by this mechanism. That rate corresponds to ~0.2 nT/y change in the radial heliospheric field at 1 AU, indicating that the mechanism is important to the heliospheric flux balance.

Heliospheric imagers offer the promise of remote sensing of large-scale structures present in the solar wind. The STEREO/HI-2 imagers, in particular, offer high resolution, very low noise observations of the inner heliosphere but have not yet been exploited to their full potential. This is in part because the signal of interest, Thomson scattered sunlight from free electrons, is ~1000 times fainter than the background visual field in the images, making background subtraction challenging. We have developed a procedure for separating the Thomson-scattered signal from the other background/foreground sources in the HI-2 data. Using only the Level 1 data from STEREO/HI-2, we are able to generate calibrated imaging data of the solar wind with sensitivity of a few times 1e-17 Bsun, compared to the background signal of a few times 1e-13 Bsun. These images reveal detailed spatial structure in CMEs and the solar wind at projected solar distances in excess of 1 AU, at the instrumental motion-blur resolution limit of 1-3 degree. CME features visible in the newly reprocessed data from December 2008 include leading-edge pileup, interior voids, filamentary structure, and rear cusps. "Quiet" solar wind features include V shaped structure centered on the heliospheric current sheet, plasmoids, and "puffs" that correspond to the density fluctuations observed in-situ. We compare many of these structures with in-situ features detected near 1 AU. The reprocessed data demonstrate that it is possible to perform detailed structural analyses of heliospheric features with visible light imagery, at distances from the Sun of at least 1 AU.

Magnetic null points can be located numerically in a potential field extrapolation or their average density can be estimated from the Fourier spectrum of a magnetogram. We use both methods to compute the null point density from a quiet Sun magnetogram made with Hinode's NFI and from magnetograms from SOHO's MDI in both its high-resolution and low-resolution modes. All estimates of the super-chromospheric column density (z>1.5 Mm) agree with one another and with the previous measurements: 0.003 null points per square Mm of solar surface.

We demonstrate that magnetic reconnection is not necessary to initiate fast CMEs. The Aly-Sturrock conjecture states that the magnetic energy of a given force free boundary field is maximized when the field is open. This is problematic for CME initiation because it leaves little or no magnetic energy to drive the eruption, unless reconnection is present to allow some of the field to escape without opening. Thus, it has been thought that reconnection must be present to initiate CMEs. This theory has not been subject to rigorous numerical testing because conventional MHD numerical models contain numerical diffusion, which introduces uncontrolled numerical reconnection. We use a quasi-Lagrangian simulation technique to run the first controlled experiments of CME initiation in the complete lack of reconnection. We find that a flux rope confined by an arcade, when twisted beyond a critical amount, can escape to an open state, allowing some of the surrounding arcade to shrink and releasing magnetic energy from the global field. This mechanism includes a true ideal MHD instability. We conclude that reconnection is not a necessary trigger for fast CME eruptions.

Using the 2004 Venus transit of the Sun to constrain a semi-empirical point-spread function for the TRACE EUV solar telescope, we have measured the effect of stray light in that telescope. We find that 43% of 171A EUV light that enters TRACE is scattered, either through diffraction off the entrance filter grid or through other nonspecular effects. We carry this result forward, via known-PSF deconvolution of TRACE images, to identify its effect on analysis of TRACE data. Known-PSF deconvolution by this derived PSF greatly reduces the effect of visible haze in the TRACE 171A images, enhances bright features, and reveals that the smooth background component of the corona is considerably less bright (and hence much more rarefied) than commonly supposed. Deconvolution reveals that some prior conlclusions about the Sun appear to have been based on stray light in the images. In particular, the diffuse background "quiet corona" becomes consistent with hydrostatic support of the coronal plasma; feature contrast is greatly increased, possibly affecting derived parameters such as the form of the coronal heating function; and essentially all existing differential emission measure studies of small features appear to be affected by contamination from nearby features. We speculate on further implications of stray light for interpretation of EUV images from TRACE and similar instruments, and advocate deconvolution as a standard tool for image analysis with future instruments such as SDO/AIA.

We examine the nature of large-scale, coronal, propagating wave fronts (``EIT waves'') and find they are incongruous with solutions using fast-mode MHD plane-wave theory. Specifically, we consider the following properties: non-dispersive single pulse manifestions, observed velocities below the local Alfven speed, and different pulses which travel at any number of constant velocities, rather than at the ``predicted'' fast-mode speed. We discuss the possibility of a soliton-like explanation for these phenomena, and show how it is consistent with the above-mentioned aspects.

Feature tracking and recognition are increasingly common tools for data analysis, but are typically implemented on an ad-hoc basis by individual research groups, limiting the usefulness of derived results when selection effects and algorithmic differences are not controlled. Specific results that are affected include the solar magnetic turnover time, the distributions of sizes, strengths, and lifetimes of magnetic features, and the physics of both small scale flux emergence and the small-scale dynamo. In this paper, we present the results of a detailed comparison between four tracking codes applied to a single set of data from SOHO/MDI, describe the interplay between desired tracking behavior and parameterization of tracking algorithms, and make recommendations for feature selection and tracking practice in future work.

Stereoscopic spectral imaging is an observing technique that affords rapid acquisition of limited spectral information over an entire image plane simultaneously. Light from a telescope is dispersed into multiple spectral orders, which are imaged separately, and two or more of the dispersed images are combined using an analogy between the (x,y,\lambda) spectral data space and conventional (x,y,z) three-space. Because no photons are deliberately destroyed during image acquisition, the technique is much more photon-efficient in some observing regimes than existing techniques such as scanned-filtergraph or scanned-slit spectral imaging. Hybrid differential stereoscopy, which uses a combination of conventional cross-correlation stereoscopy and linear approximation theory to extract the central wavelength of a spectral line, has been used to produce solar Stokes-V (line-of-sight) magnetograms in the 617.34 nm Fe I line, and more sophisticated inversion techniques are currently being used to derive Doppler and line separation data from EUV images of the solar corona collected in the neighboring lines of He-II and Si-XI at 30.4 nm. In this paper we develop an analytic a priori treatment of noise in the line shift signal derived from hybrid differential stereoscopy. We use the analysis to estimate the noise level and measurement precision in a high resolution solar magnetograph based on stereoscopic spectral imaging, compare those estimates to a test observation made in 2003, and discuss implications for future instruments.

Fine-scale structure in the corona appears not to be well resolved by current imaging instruments. Assuming this to be true offers a simple geometric explanation for several current puzzles in coronal physics, including: the apparent uniform cross-section of bright threadlike structures in the corona; the low EUV contrast (long apparent scale height) between the top and bottom of active region loops; and the inconsistency between loop densities derived by spectral and photometric means. Treating coronal loops as a mixture of diffuse background and very dense, unresolved filamentary structures address these problems with a combination of high plasma density within the structures, which greatly increases the emissivity of the structures, and geometric effects that attenuate the apparent brightness of the feature at low altitudes. It also suggests a possible explanation for both the surprisingly high contrast of EUV coronal loops against the coronal background, and the uniform ``typical'' height of the bright portion of the corona (about 0.3 solar radii) in full-disk EUV images. Some ramifications of this picture are discussed, including an estimate (10-100 km) of the fundamental scale of strong heating events in the corona.

We have developed a new, quasi-Lagrangian approach for numerical modeling of magnetohydrodynamics in low to moderate $\beta$ plasmas such as the solar corona. We introduce the concept of a ``fluxon'', a discretized field line. Fluxon models represent the magnetic field as a skeleton of such discrete field lines, and interpolate field values from the geometry of the skeleton where needed, reversing the usual direction of the field line transform. The fluxon skeleton forms the grid for a collection of 1-D Eulerian models of plasma along individual flux tubes. Fluxon models have no numerical resistivity, because they preserve topology explicitly. Our prototype code, \emph{FLUX}, is currently able to find 3-D nonlinear force-free field solutions with a specified field topology, and work is ongoing to validate and extend the code to full magnetohydrodynamics. FLUX has significant scaling advantages over conventional models: for ``magnetic carpet'' models, with photospheric line-tied boundary conditions, FLUX simulations scale in complexity like a conventional 2-D grid although the full 3-D field is represented. The code is free software and is available online. In this current paper we introduce fluxons and our prototype code, and describe the course of future work with the code.