The National Science Foundation (NSF) Daniel K. Inouye Solar Telescope (DKIST) is about to start operations at the summit of Haleakala (Hawaii). DKIST will join the early science phases of the NASA and ESA Parker Solar Probe and Solar Orbiter encounter missions. By combining in-situ measurements of the near-sun plasma environment and detail remote observations of multiple layers of the Sun, the three observatories form an unprecedented multi-messenger constellation to study the magnetic connectivity inside the solar system. This white paper outlines the synergistic science that this multi-messenger suite enables.

Gradient- and curl-type or E- and B-type polarizations have been routinely analyzed to study the physics contributing to the cosmic microwave background polarization and galactic foregrounds. They characterize the parity-even and parity-odd properties of the underlying physical mechanisms, for example hydromagnetic turbulence in the case of dust polarization. Here we study spectral correlation functions characterizing the parity-even and parity-odd parts of linear polarization for homogeneous and inhomogeneous turbulence to show that only the inhomogeneous helical case can give rise to a parity-odd polarization signal. We also study nonhelical turbulence and suggest that a strong nonvanishing (here negative) skewness of the E polarization is responsible for an enhanced ratio of the EE to the BB (quadratic) correlation in both helical and nonhelical cases. This could explain the enhanced EE/BB ratio observed recently for dust polarization. We close with a preliminary assessment of using linear polarization of the Sun to characterize its helical turbulence without being subjected to the pi ambiguity that magnetic inversion techniques have to address.

Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.

Accurate prediction of solar activity calls for precise calibration of solar cycle models. Consequently we aim to find optimal parameters for models which describe the physical processes on the solar surface, which in turn act as proxies for what occurs in the interior and provide source terms for coronal models. We use a genetic algorithm to optimize surface flux transport models using National Solar Observatory (NSO) magnetogram data for Solar Cycle 23. This is applied to both a 1D model that inserts new magnetic flux in the form of idealized bipolar magnetic regions, and also to a 2D model that assimilates specific shapes of real active regions. The genetic algorithm searches for parameter sets (meridional flow speed and profile, supergranular diffusivity, initial magnetic field, and radial decay time) that produce the best fit between observed and simulated butterfly diagrams, weighted by a latitude-dependent error structure which reflects uncertainty in observations. Due to the easily adaptable nature of the 2D model, the optimization process is repeated for Cycles 21, 22, and 24 in order to analyse cycle-to-cycle variation of the optimal solution. We find that the ranges and optimal solutions for the various regimes are in reasonable agreement with results from the literature, both theoretical and observational. The optimal meridional flow profiles for each regime are almost entirely within observational bounds determined by magnetic feature tracking, with the 2D model being able to accommodate the mean observed profile more successfully. Differences between models appear to be important in deciding values for the diffusive and decay terms. In like fashion, differences in the behaviours of different solar cycles lead to contrasts in parameters defining the meridional flow and initial field strength.

We develop a two-scale formalism to determine global magnetic helicity spectra in systems where the local magnetic helicity has opposite signs on both sides of the equator, giving rise to cancelation with conventional methods. We verify this approach using first a synthetic one-dimensional magnetic field and then two-dimensional slices from a three-dimensional alpha effect-type dynamo-generated magnetic field with forced turbulence of opposite helicity above and below the midplane of the domain. We then apply this formalism to global solar synoptic vector magnetograms. To improve the statistics, data from three consecutive Carrington rotations (2161--2163) are combined into a single map. We find that the spectral magnetic helicity representative of the northern hemisphere is negative at all wavenumbers and peaks at ~ 0.06 Mm^{-1} (scales around 100 Mm). There is no evidence of bihelical magnetic fields that are found in three-dimensional turbulence simulations of helicity-driven alpha effect-type dynamos.

Recently, the National Solar Observatory (NSO) Solar-atmosphere Pipeline Working Group has undertaken the production of synoptic maps from Helioseismic and Magnetic Imager (HMI) magnetograms. A set of maps has been processed spanning the data available for 2010-2015 using twice daily images (taken at UT midnight and noon) and running them through the same algorithms used to produce SOLIS/VSM 6302l mean-magnetic and spatial-variance maps. The contents of this document provide an overview of what these maps look like, and the processing steps used to generate them from the original HMI input data.

Analysis of synoptic data from the Vector Stokes Magnetograph (VSM) of the Synoptic Optical Long-term Investigations of the Sun (SOLIS) and the NASA/NSO Spectromagnetograph (SPM) at the NSO/Kitt Peak Vacuum Telescope facility shows that the reversals of solar polar magnetic fields exhibit elements of a stochastic process, which may include the development of specific patterns of emerging magnetic flux, and the asymmetry in activity between northern and southern hemispheres. The presence of such irregularities makes the modeling and prediction of polar field reversals extremely hard if possible. In a classical model of solar activity cycle, the unipolar magnetic regions (UMRs) of predominantly following polarity fields are transported polewards due to meridional flows and diffusion. The UMRs gradually cancel out the polar magnetic field of the previous cycle, and re-build the polar field of opposite polarity setting the stage for the next cycle. We show, however, that this deterministic picture can be easily altered by the developing of a strong center of activity, or by the emergence of an extremely large active region, or by a "strategically placed" coronal hole. We demonstrate that the activity occurring during the current cycle 24 may be the result of this randomness in the evolution of the solar surface magnetic field.

Coronal mass ejections (CMEs) with angular width $> 30^{\circ}$ have been observed to occur at a higher rate during solar cycle 24 compared to cycle 23, per sunspot number. This result is supported by data from three independent databases constructed using Large Angle and Spectrometric Coronagraph Experiment (LASCO) coronagraph images, two employing automated detection techniques and one compiled manually by human observers. According to the two databases that cover a larger field of view, the enhanced CME rate actually began shortly after the cycle 23 polar field reversal, in 2004, when the polar fields returned with a 40\% reduction in strength and interplanetary radial magnetic field became $\approx 30\%$ weaker. This result is consistent with the link between anomalous CME expansion and heliospheric total pressure decrease recently reported by Gopalswamy et al.

We study the correlation between abrupt permanent changes of magnetic field during X-class flares observed by the GONG and HMI instruments, and the hard X-ray (HXR) emission observed by RHESSI, to relate the photospheric field changes to the coronal restructuring and investigate the origin of the field changes. We find that spatially the early RHESSI emission corresponds well to locations of the strong field changes. The field changes occur predominantly in the regions of strong magnetic field near the polarity inversion line (PIL). The later RHESSI emission does not correspond to significant field changes as the flare footpoints are moving away from the PIL. Most of the field changes start before or around the start time of the detectable HXR signal, and they end at about the same time or later than the detectable HXR flare emission. Some of the field changes propagate with speed close to that of the HXR footpoint at a later phase of the flare. The propagation of the field changes often takes place after the strongest peak in the HXR signal when the footpoints start moving away from the PIL, i.e. the field changes follow the same trajectory as the HXR footpoint, but at an earlier time. Thus, the field changes and HXR emission are spatio-temporally related but not co-spatial nor simultaneous. We also find that in the strongest X-class flares the amplitudes of the field changes peak a few minutes earlier than the peak of the HXR signal. We briefly discuss this observed time delay in terms of the formation of current sheets during eruptions.

We compare the average photospheric current helicity $H_c$, photospheric twist parameter $\alpha$ (a well-known proxy for the full relative magnetic helicity), and subsurface kinetic helicity $H_k$ for 194 active regions observed between 2006-2013. We use 2440 Hinode photospheric vector magnetograms, and the corresponding subsurface fluid velocity data derived from GONG (2006-2012) and HMI (2010-2013) Dopplergrams. We find a significant hemispheric bias in all three parameters. The subsurface kinetic helicity is preferentially positive in the southern hemisphere and negative in the northern hemisphere. The photospheric current helicity and the $\alpha$ parameter have the same bias for strong fields ($|B|>1000$ G) and no significant bias for weak fields (100 G $<|B|<500$ G). We find no significant region-by-region correlation between the subsurface kinetic helicity and either the strong-field current helicity or $\alpha$. Subsurface fluid motions of a given handedness correspond to photospheric helicities of both signs in approximately equal numbers. However, common variations appear in annual averages of these quantities over all regions. Furthermore, in a subset of 77 regions we find significant correlations between the temporal profiles of the subsurface and photospheric helicities. In these cases, the sign of the linear correlation coefficient matches the sign relationship between the helicities, indicating that the photospheric magnetic field twist is sensitive to the twisting motions below the surface.

It is shown that expressions for the global Lorentz force associated with a flaring active region derived by Fisher {\it et al.} (2012) can be used to estimate the Lorentz force changes for strong fields in large structures over photospheric subdomains within active regions. Gary's~(2001) model for the stratified solar atmosphere is used to demonstrate that in large-scale structures with typical horizontal magnetic length scale $\gg 300$~km and with strong magnetic fields ($\ge 1$~kG at the $\tau =1$ opacity layer at 5000~\AA), the Lorentz force acting on the photosphere may be approximated by a surface integral based on photospheric boundary data alone. These conditions cover many of the sunspot fields and major neutral lines that have been studied using Fisher {\it et al.}'s (2012) expressions over the past few years. The method gives a reasonable estimate of flare-related Lorentz force changes based on photospheric magnetogram observations provided that the Lorentz force changes associated with the flare have a lasting effect on the observed fields, and are not immediately erased by post-flare equilibration processes.

Magnetic flux synoptic charts are critical for reliable modeling of the corona and heliosphere. Until now, however, these charts were provided without any estimate of uncertainties. The uncertainties are due to instrumental noise in the measurements and to the spatial variance of the magnetic flux distribution that contributes to each bin in the synoptic chart. We describe here a simple method to compute synoptic magnetic flux maps and their corresponding magnetic flux spatial variance charts that can be used to estimate the uncertainty in the results of coronal models. We have tested this approach by computing a potential-field-source-surface model of the coronal field for a Monte Carlo simulation of Carrington synoptic magnetic flux maps generated from the variance map. We show that these uncertainties affect both the locations of source-surface neutral lines and the distributions of coronal holes in the models.

We study the relationship between decaying active region magnetic fields, coronal holes and the global coronal magnetic structure using Global Oscillations Network Group (GONG) synoptic magnetograms, Solar Terrestrial RElations Observatory (STEREO) extreme ultra-violet (EUV) synoptic maps and coronal potential-field source-surface (PFSS) models. We analyze 14 decaying regions and associated coronal holes occurring between early 2007 and late 2010, four from cycle 23 and 10 from cycle 24. We investigate the relationship between asymmetries in active regions' positive and negative magnetic intensities, asymmetric magnetic decay rates, flux imbalances, global field structure and coronal hole formation. Whereas new emerging active regions caused changes in the large-scale coronal field, the coronal fields of the 14 decaying active regions only opened under the condition that the global coronal structure remained almost unchanged. This was because the dominant slowly-varying, low-order multipoles prevented opposing-polarity fields from opening and the remnant active-region flux preserved the regions' low-order multipole moments long after the regions had decayed. Thus the polarity of each coronal hole necessarily matched the polar field on the side of the streamer belt where the corresponding active region decayed. For magnetically isolated active regions initially located within the streamer belt, the more intense polarity generally survived to form the hole. For non-isolated regions, flux imbalance and topological asymmetry prompted the opposite to occur in some cases.

We study the evolution of the observed photospheric magnetic field and the modeled global coronal magnetic field during the past 3 1/2 solar activity cycles observed since the mid-1970s. We use synoptic magnetograms and extrapolated potential-field models based on longitudinal full-disk photospheric magnetograms from the NSO's three magnetographs at Kitt Peak, the Synoptic Optical Long-term Investigations of the Sun (SOLIS) vector spectro-magnetograph (VSM), the spectro-magnetograph and the 512-channel magnetograph instruments, and from the U. Stanford's Wilcox Solar Observatory. The associated multipole field components are used to study the dominant length scales and symmetries of the coronal field. Polar field changes are found to be well correlated with active fields over most of the period studied, except between 2003-6 when the active fields did not produce significant polar field changes. Of the axisymmetric multipoles, only the dipole and octupole follow the poles whereas the higher orders follow the activity cycle. All non-axisymmetric multipole strengths are well correlated with the activity cycle. The tilt of the solar dipole is therefore almost entirely due to active-region fields. The axial dipole and octupole are the largest contributors to the global field except while the polar fields are reversing. This influence of the polar fields extends to modulating eruption rates. According to the Computer Aided CME Tracking (CACTus), Solar Eruptive Event Detection System (SEEDS), and Nobeyama radioheliograph prominence eruption catalogs, the rate of solar eruptions is found to be systematically higher for active years between 2003-2012 than for those between 1997-2002. This behavior appears to be connected with the weakness of the late-cycle 23 polar fields as suggested by Luhmann. We see evidence that the process of cycle 24 field reversal is well advanced at both poles.

We perform a statistical study of permanent changes in longitudinal fields associated with solar flares by tracking magnetic features. The YAFTA feature tracking algorithm is applied to GONG++ one-minute magnetograms for 77 X-class and M-class flares to analyze the evolution and interaction of the magnetic features and to estimate the amount of canceled magnetic flux. We find that significantly more magnetic flux decreases than increases occurred during the flares, consistent with a model of collapsing loop structure for flares. Correlations between both total (unsigned) and net (signed) flux changes and the GOES peak X-ray flux are dominated by X-class flares at limb locations. The flux changes were accompanied in most cases by significant cancellation, most of which occurred during the flares. We find that the field strength and complexity near the polarity inversion line are approximately equally important in the flux cancellation processes that accompany the flares. We do not find a correlation between the flux cancellation events and the stepwise changes in the magnetic flux in the region.

We analyze the spatial and temporal variations of the abrupt photospheric magnetic changes associated with six major flares using 12-minute, 0.5-arcsec vector magnetograms from NASA's Helioseismic and Magnetic Imager instrument (HMI) on the Solar Dynamics Observatory satellite. The six major flares occurred near the main magnetic neutral lines of four active regions, NOAA 11158, 11166, 11283 and 11429. During all six flares the neutral line field vectors became stronger and more horizontal, in each case almost entirely due to strengthening of the horizontal field components parallel to the neutral line. In all six cases the neutral line pre-flare fields were more vertical than the reference potential fields, and collapsed abruptly and permanently closer to potential field tilt angles during every flare, implying that the relaxation of magnetic stress associated with non-potential tilt angles plays a major role during major flares. The shear angle with respect to the reference potential field did not show such a pattern, demonstrating that flare processes do not generally relieve magnetic stresses associated with photospheric magnetic shear. The horizontal fields became significantly and permanently more aligned with the neutral line during the four largest flares, suggesting that the collapsing field is on average more aligned with the neutral line than the pre-flare neutral line field. The vertical Lorentz force had a large, abrupt, permanent downward change during each of the flares, consistent with loop collapse. The horizontal Lorentz force changes acted mostly parallel to the neutral line in opposite directions on each side, a signature of the fields contracting during the flare, pulling the two sides of the neutral line towards each other. The greater effect of the flares on field tilt than on shear may be explained by photospheric line-tying.

We have used Transition Region and Coronal Explorer (TRACE) 1600 \AA images and Global Oscillation Network Group (GONG) magnetograms to compare ultraviolet (UV) emissions from the chromosphere to longitudinal magnetic field changes in the photosphere during four X-class solar flares. An abrupt, significant, and persistent change in the magnetic field occurred across more than ten pixels in the GONG magnetograms for each flare. These magnetic changes lagged the GOES flare start times in all cases, showing that they were consequences and not causes of the flares. Ultraviolet emissions were spatially coincident with the field changes. The UV emissions tended to lag the GOES start times for the flares, and led the changes in the magnetic field in all pixels except one. The UV emissions led the photospheric field changes by 4 minutes on average with the longest lead being 9 minutes, however, the UV emissions continued for tens of minutes, and more than an hour in some cases, after the field changes were complete. The observations are consistent with the picture in which an Alfv\'{e}n wave from the field reconnection site in the corona propagates field changes outward in all directions near the onset of the impulsive phase, including downwards through the chromosphere and into the photosphere, causing the photospheric field changes, whereas the chromosphere emits in the UV in the form of flare kernels, ribbons and sequential chromospheric brightenings during all phases of the flare.

The evolution of the photospheric magnetic field during the declining phase and minimum of Cycle 23 and the recent rise of Cycle 24 are compared with the behavior during previous cycles. We used longitudinal full-disk magnetograms from the NSO's three magnetographs at Kitt Peak, the Synoptic Optical Long-term Investigations of the Sun (SOLIS) Vector Spectro-Magnetograph (VSM), the Spectromagnetograph and the 512-Channel Magnetograph instruments, and longitudinal full-disk magnetograms from the Mt. Wilson 150-foot tower. We analyzed 37 years of observations from these two observatories that have been observing daily, weather permitting, since 1974, offering an opportunity to study the evolving relationship between the active region and polar fields in some detail over several solar cycles. It is found that the annual averages of a proxy for the active region poloidal magnetic field strength, the magnetic field strength of the high-latitude poleward streams, and the time derivative of the polar field strength are all well correlated in each hemisphere. These results are based on statistically significant cyclical patterns in the active region fields and are consistent with the Babcock-Leighton phenomenological model for the solar activity cycle. There was more hemispheric asymmetry in the activity level, as measured by total and maximum active region flux, during late Cycle 23 (after around 2004), when the southern hemisphere was more active, and Cycle 24 up to the present, when the northern hemisphere has been more active, than at any other time since 1974. The active region net proxy poloidal fields effectively disappeared in both hemispheres around 2004, and the polar fields did not become significantly stronger after this time. We see evidence that the process of Cycle 24 field reversal has begun at both poles.

The active region NOAA 11158 produced the first X-class flare of Solar Cycle 24, an X2.2 flare at 01:44 UT on 2011 February 15. Here we analyze SDO/HMI magnetograms covering a 12-hour interval centered at the time of this flare. We describe the spatial distributions of the photospheric magnetic changes associated with this flare, including the abrupt changes in the field vector, vertical electric current and Lorentz force vector. We also trace these parameters' temporal evolution. The abrupt magnetic changes were concentrated near the neutral line and in two neighboring sunspots. Near the neutral line, the field vectors became stronger and more horizontal during the flare and the shear increased. This was due to an increase in strength of the horizontal field components near the neutral line, most significant in the horizontal component parallel to the neutral line but the perpendicular component also increased in strength. The vertical component did not show a significant, permanent overall change at the neutral line. The increase in total flux at the neutral line was accompanied by a compensating flux decrease in the surrounding volume. In the two sunspots near the neutral line the azimuthal flux abruptly decreased during the flare but this change was permanent in only one of the spots. There was a large, abrupt, downward vertical Lorentz force change during the flare, consistent with results of past analyses and recent theoretical work. The horizontal Lorentz force acted in opposite directions along each side of neutral line, with the two sunspots at each end subject to abrupt torsional forces. The shearing forces were consistent with field contraction and decrease of shear near the neutral line, whereas the field itself became more sheared as a result of the flux collapsing towards the neutral line from the surrounding volume.

Between 2008/3/24 and 2008/4/2, the three active regions NOAA active regions 10987, 10988 and 10989 were observed daily by the Synoptic Optical Long-term Investigations of the Sun (SOLIS) Vector Spectro-Magnetograph (VSM) while they traversed the solar disk. We use these measurements and the nonlinear force-free magnetic field code XTRAPOL to reconstruct the coronal magnetic field for each active region and compare model field lines with images from the Solar Terrestrial RElations Observatory (STEREO) and Hinode X-ray Telescope (XRT) telescopes. Synoptic maps made from continuous, round-the-clock Global Oscillations Network Group (GONG) magnetograms provide information on the global photospheric field and potential-field source-surface models based on these maps describe the global coronal field during the Whole Heliospheric Interval (WHI) and its neighboring rotations. Features of the modeled global field, such as the coronal holes and streamer belt locations, are discussed in comparison with extreme ultra-violet and coronagraph observations from STEREO. The global field is found to be far from a minimum, dipolar state. From the nonlinear models we compute physical quantities for the active regions such as the photospheric magnetic and electric current fluxes, the free magnetic energy and the relative helicity for each region each day where observations permit. The interconnectivity of the three regions is addressed in the context of the potential-field source-surface model. Using local and global quantities derived from the models, we briefly discuss the different observed activity levels of the regions.

Besides their own intrinsic interest, correct interpretation of solar surface magnetic field observations is crucial to our ability to describe the global magnetic structure of the solar atmosphere. Photospheric magnetograms are often used as lower boundary conditions in models of the corona, but not data from the nearly force-free chromosphere. National Solar Observatory's (NSO) Synoptic Optical Long-term Investigations of the Sun VSM (Vector Spectromagnetograph) produces full-disk line-of-sight magnetic flux images deriving from both photospheric and chromospheric layers on a daily basis. In this paper, we investigate key properties of the magnetic field in these two layers using more than five years of VSM data. We find from near-equatorial measurements that the east-west inclination angle of most photospheric fields is less than about 12{\deg}, while chromospheric fields expand in all directions to a significant degree. Using a simple stereoscopic inversion, we find evidence that photospheric polar fields are also nearly radial but that during 2008 the chromospheric field in the south pole was expanding superradially. We obtain a spatially resolved polar photospheric flux distribution up to 80{\deg} latitude whose strength increases poleward approximately as cosine(colatitude) to the power 9-10. This distribution would give a polar field strength of 5-6 G. We briefly discuss implications for future synoptic map construction and modeling.

We characterize the changes in the longitudinal photospheric magnetic field during 38 X-class and 39 M-class flares within $65^{\circ}$ of disk-center using 1-minute GONG magnetograms. In all 77 cases we identify at least one site in the flaring active region where clear, permanent, stepwise field changes occurred. The median duration of the field changes was about 15 minutes and was approximately equal for X-class and for M-class flares. The absolute values of the field changes ranged from the detection limit of $\sim\!\!10$~G to as high as $\sim\!\!450$~G in two exceptional cases. The median value was 69~G. Field changes were significantly stronger for X-class than for M-class flares and for limb flares than for disk-center flares. Longitudinal field changes less than 100~G tended to decrease longitudinal field strengths, both close to disk-center and close to the limb, while field changes greater than 100~G showed no such pattern. Likewise, longitudinal flux strengths tended to decrease during flares. Flux changes, particularly net flux changes near disk-center, correlated better than local field changes with GOES peak X-ray flux. The strongest longitudinal field and flux changes occurred in flares observed close to the limb. We estimate the change of Lorentz force associated with each flare and find that this is large enough in some cases to power seismic waves. We find that longitudinal field decreases would likely outnumber increases at all parts of the solar disk within $65^{\circ}$ of disk-center, as in our observations, if photospheric field tilts increase during flares as predicted by Hudson et al.

Observations from the Hinode/XRT telescope and STEREO/SECCHI/EUVI are utilized to study polar coronal jets and plumes. The study focuses on the temporal evolution of both structures and their relationship. The data sample, spanning April 7-8 2007, shows that over 90% of the 28 observed jet events are associated with polar plumes. EUV images (STEREO/SECCHI) show plume haze rising from the location of approximately 70% of the polar X-ray (Hinode/XRT) and EUV jets, with the plume haze appearing minutes to hours after the jet was observed. The remaining jets occurred in areas where plume material previously existed causing a brightness enhancement of the latter after the jet event. Short-lived, jet-like events and small transient bright points are seen (one at a time) at different locations within the base of pre-existing long-lived plumes. X-ray images also show instances (at least two events) of collimated-thin jets rapidly evolving into significantly wider plume-like structures that are followed by the delayed appearance of plume haze in the EUV. These observations provide evidence that X-ray jets are precursors of polar plumes, and in some cases cause brightenings of plumes. Possible mechanisms to explain the observed jet and plume relationship are discussed.

Coronal loop emission profiles are often of remarkably constant width along their entire lengths, contradicting expectations based on model coronal magnetic field strengths decreasing with height. Meanwhile Paul Bellan has produced a theoretical model in which an initially empty, twisted force-free loop, on being filled with plasma via upflow at each foot point, in the absence of significant gravitational effects, forms a narrow, filamentary loop of constant cross-section. In this paper, we focus on equilibrium states that include stratification by uniform gravity while retaining the effects of magnetic field twist. Comparing these with related force-free equilibria, it is found that injection of low-$\beta $ plasma under coronal conditions is not likely to change the shape of a loop significantly. These linear equilibria apply to the interiors and boundaries of loops only, with external influences modeled by boundary total pressures. The effects of total pressure balance with surroundings and of gravitational stratification are to inhibit the pinching of a loop to a constant cross-section. Only if the plasma $\beta$ were high enough for the plasma to reconfigure the external field and the hydrostatic scale height much greater than the loop size could the final state have nearly constant cross section. We do not expect this to occur in the corona.

We present numerical magnetohydrostatic solutions describing the gravitationally stratified, bulk equilibrium of cool, dense prominence plasma embedded in a near-potential coronal field. These solutions are calculated using the FINESSE magnetohydrodynamics equilibrium solver and describe the morphologies of magnetic field distributions in and around prominences and the cool prominence plasma that these fields support. The equilibrium condition for this class of problem is usually different in distinct subdomains, separated by free boundaries, across which solutions are matched by suitable continuity or jump conditions describing force balance. We employ our precise finite element elliptic solver to calculate solutions not accessible by previous analytical techniques with temperature or entropy prescribed as free functions of the magnetic flux function, including a range of values of the polytropic index, temperature variations mainly across magnetic field lines and photospheric field profiles sheared close to the polarity inversion line. Out of the many examples computed here, perhaps the most noteworthy is one which reproduces precisely the three-part structure often encountered in observations: a cool dense prominence within a cavity/flux rope embedded in a hot corona. The stability properties of these new equilibria, which may be relevant to solar eruptions, can be determined in the form of a full resistive MHD spectrum using a companion hyperbolic stability solver.

We consider the potential magnetic field associated with a helical electric line current flow, idealizing the near-potential coronal field within which a highly localized twisted current structure is embedded. It is found that this field has a significant axial component off the helical magnetic axis where there is no current flow, such that the flux winds around the axis. The helical line current field, in including the effects of flux rope writhe, is therefore more topologically complex than straight line and ring current fields sometimes used in solar flux rope models. The axial flux in magnetic fields around confined current structures may be affected by the writhe of these current structures such that the field twists preferentially with the same handedness as the writhe. This property of fields around confined current structures with writhe may be relevant to classes of coronal magnetic flux rope, including structures observed to have sigmoidal forms in soft X-rays and prominence magnetic fields. For example, ``bald patches'' and the associated heating by Parker current sheet dissipation seem likely. Thus some measurements of flux rope magnetic helicities may derive from external, near-potential fields. The predicted hemispheric preference for positive and negative magnetic helicities is consistent with observational results for prominences and sigmoids and past theoretical results for flux rope internal fields.

The scale heights of stratification and the widths of steady solar coronal loops exhibit properties unexplained by standard theory: observed scale heights are often much greater than static theory predicts, while the nearly-constant widths of loop emission signatures defy theoretical expectations for large flux tubes in stratified media. In this work we relate the cross-sectional profile of a coronal flux tube to its density scale height in steady-state plasma flow regimes. Steady flows may shorten or lengthen the scale height according to how the tube cross-sectional area varies with arclength. In a near-potential corona the flux tubes are expected to be sufficiently expansive in many active regions for scale heights to be increased by steady flows. On the other hand, cases where scale lengths are actually increased to observed sizes form a small part of the solution space, close to regimes where density profiles reverse. Therefore, although steady flows are the only steady process known to be capable of extending scale heights significantly, they are not expected to be not responsible for the majority of extended active region scale heights.

We analyzed a coronal loop observed with the Normal Incidence Spectrometer (NIS), which is part of the Coronal Diagnostic Spectrometer (CDS) on board the Solar and Heliospheric Observatory (SOHO). The measured Doppler shifts and proper motions along the selected loop strongly indicate unidirectional flows. Analysing the Emission Measure Curves of the observed spectral lines, we estimated that the temperature along the loop was about 380000 K. We adapted a solution of the ideal MHD steady equations to our set of measurements. The derived energy balance along the loop, as well as the advantages/disadvantages of this MHD model for understanding the characteristics of solar coronal loops are discussed.

An analytical MHD model of a normal-polarity prominence with compressible flow is presented. The exact solution is constructed via a systematic nonlinear separation of variables method used to calculate several classes of MHD equilibria in Cartesian geometry and uniform gravity. Although the model is 2D, a third magnetic/velocity vector field component is included and the highly sheared fields observed in prominences are reproduced. A description is given of the balance of gas pressure gradient with gravity and the Lorentz or inertial forces acting along and across the prominence. It is found that the flow does not significantly influence the heating profile. The analyzed model has dimensions, plasma density, temperature and velocity profiles which agree with those in the observations literature.

An analytical MHD model of coronal loops with compressible flows and including heating is compared to observational data. The model is constructed via a systematic nonlinear separation of the variables method used to calculate several classes of exact MHD equilibria in Cartesian geometry and uniform gravity. By choosing a particularly versatile solution class with a large parameter space we are able to calculate models whose loop length, shape, plasma density, temperature and velocity profiles are fitted to loops observed with TRACE, SoHO/CDS and SoHO/SUMER. Synthetic emission profiles are also calculated and fitted to the observed emission patterns. An analytical discussion is given of the two-dimenional balance of the Lorentz force and the gas pressure gradient, gravity and inertial forces acting along and across the loop. These models are the first to include a fully consistent description of the magnetic field, 2D geometry, plasma density and temperature, flow velocity and thermodynamics of loops. The consistently calculated heating profiles which are largely dominated by radiative losses are influenced by the flow and are asymmetrical being concentrated at the upflow footpoint.

We present a systematic method for constructing two-dimensional magnetohydrodynamic equilibria with compressible flow in Cartesian geometry. This systematic method has already been developed in spherical geometry and applied in modelling solar and stellar winds and outflows (Vlahakis & Tsinganos,1998) but is derived here in Cartesian geometry in the context of the solar atmosphere for the first time. Using the method we find several new classes of solutions, some of which generalise known solutions, including the Kippenhahn & Schl\"uter (1957) and Hood & Anzer (1990) solar prominence models and the Tsinganos, Surlantzis & Priest (1993) coronal loop model with flow, and some of which are completely new. Having developed the method in full and summarised the several classes of solutions, we explore in a some detail one of the classes to illustrate the general construction method. From one of the new classes of solutions we calculate two loop-like solutions, one of which is the first exact two-dimensional magnetohydrodynamic equilibrium with trans-Alfv\'enic flow.