The Coronal Global Evolutionary Model (CGEM) provides data-driven simulations of the magnetic field in the solar corona to better understand the build-up of magnetic energy that leads to eruptive events. The CGEM project has developed six capabilities. CGEM modules (1) prepare time series of full-disk vector magnetic field observations to (2) derive the changing electric field in the solar photosphere over active-region scales. This local electric field is (3) incorporated into a surface flux transport model that reconstructs a global electric field that evolves magnetic flux in a consistent way. These electric fields drive a (4) 3D spherical magneto-frictional (SMF) model, either at high-resolution over a restricted range of solid angle or at lower resolution over a global domain, to determine the magnetic field and current density in the low corona. An SMF-generated initial field above an active region and the evolving electric field at the photosphere are used to drive (5) detailed magneto-hydrodynamic (MHD) simulations of active regions in the low corona. SMF or MHD solutions are then used to compute emissivity proxies that can be compared with coronal observations. Finally, a lower-resolution SMF magnetic field is used to initialize (6) a global MHD model that is driven by an SMF electric-field time series to simulate the outer corona and heliosphere, ultimately connecting Sun to Earth. As a demonstration, this report features results of CGEM applied to observations of the evolution of NOAA Active Region 11158 in February 2011.

We describe the PDFI_SS software library, which is designed to find the electric field at the Sun's photosphere from a sequence of vector magnetogram and Doppler velocity measurements, and estimates of horizontal velocities obtained from local correlation tracking using the recently upgraded FLCT code. The library, a collection of Fortran subroutines, uses the "PDFI" technique described by Kazachenko et al. (2014), but modified for use in spherical, Plate-Carr\'ee geometry on a staggered grid. The domain over which solutions are found is a subset of the global spherical surface, defined by user-specified limits of colatitude and longitude. Our staggered-grid approach, based on that of Yee (1966), is more conservative and self-consistent compared to the centered, Cartesian grid used by Kazachenko et al. (2014). The library can be used to compute an end-to-end solution for electric fields from data taken by the HMI instrument aboard NASA's SDO Mission. This capability has been incorporated into the HMI pipeline processing system operating at SDO's JSOC. The library is written in a general and modular way so that the calculations can be customized to modify or delete electric field contributions, or used with other data sets. Other applications include "nudging" numerical models of the solar atmosphere to facilitate assimilative simulations. The library includes an ability to compute "global" (whole-Sun) electric field solutions. The library also includes an ability to compute Potential Magnetic Field solutions in spherical coordinates. This distribution includes a number of test programs which allow the user to test the software.

Solar coronal dimmings have been observed extensively in the past two decades and are believed to have close association with coronal mass ejections (CMEs). Recent study found that coronal dimming is the only signature that could differentiate powerful ares that have CMEs from those that do not. Therefore, dimming might be one of the best candidates to observe the stellar CMEs on distant Sun-like stars. In this study, we investigate the possibility of using coronal dimming as a proxy to diagnose stellar CMEs. By simulating a realistic solar CME event and corresponding coronal dimming using a global magnetohydrodynamics model (AWSoM: Alfven-wave Solar Model), we first demonstrate the capability of the model to reproduce solar observations. We then extend the model for simulating stellar CMEs by modifying the input magnetic flux density as well as the initial magnetic energy of the CME flux rope. Our result suggests that with improved instrument sensitivity, it is possible to detect the coronal dimming signals induced by the stellar CMEs.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

The most energetic solar flares are typically associated with the ejection of a cloud of coronal material into the heliosphere in the form of a coronal mass ejection (CME). However, there exist large flares which are not accompanied by a CME. The existence of these non-eruptive flares raises the question of whether such flares suffer from a lack of access to nearby open fields in the vicinity above the flare (reconnection) site. In this study, we use a sample of 56 flares from Sunspot Cycles 23 and 24 to test whether active regions that produce eruptive X-class flares are preferentially located near coronal magnetic field domains that are open to the heliosphere, as inferred from a potential field source surface model. The study shows that X-class flares having access to open fields are eruptive at a higher rate than those for which access is lacking. The significance of this result should be moderated due to the small number of non-eruptive X-class flares in the sample, based on the associated Bayes factor.

The heliospheric magnetic field is of pivotal importance in solar and space physics. The field is rooted in the Sun's photosphere, where it has been observed for many years. Global maps of the solar magnetic field based on full disk magnetograms are commonly used as boundary conditions for coronal and solar wind models. Two primary observational constraints on the models are (1) the open field regions in the model should approximately correspond to coronal holes observed in emission, and (2) the magnitude of the open magnetic flux in the model should match that inferred from in situ spacecraft measurements. In this study, we calculate both MHD and PFSS solutions using fourteen different magnetic maps produced from five different types of observatory magnetograms, for the time period surrounding July, 2010. We have found that for all of the model/map combinations, models that have coronal hole areas close to observations underestimate the interplanetary magnetic flux, or, conversely, for models to match the interplanetary flux, the modeled open field regions are larger than coronal holes observed in EUV emission. In an alternative approach, we estimate the open magnetic flux entirely from solar observations by combining automatically detected coronal holes for Carrington rotation 2098 with observatory synoptic magnetic maps. This approach also underestimates the interplanetary magnetic flux. Our results imply that either typical observatory maps underestimate the Sun's magnetic flux, or a significant portion of the open magnetic flux is not rooted in regions that are obviously dark in EUV and X-ray emission.

With the global view and high-cadence observations from SDO/AIA and STEREO, many spatially separated solar eruptive events appear to be coupled. However, the mechanisms for "sympathetic" events are still largely unknown. In this study, we investigate the impact of an erupting flux rope on surrounding solar structures through large-scale magnetic coupling. We build a realistic environment of the solar corona on 2011 February 15 using a global magnetohydrodynamics (MHD) model and initiate coronal mass ejections (CMEs) in active region (AR) 11158 by inserting Gibson-Low analytical flux ropes. We show that a CME's impact on the surrounding structures depends not only on the magnetic strength of these structures and their distance to the source region, but also on the interaction between the CME with the large-scale magnetic field. Within the CME expansion domain where the flux rope field directly interacts with the solar structures, expansion-induced reconnection often modifies the overlying field, thereby increasing the decay index. This effect may provide a primary coupling mechanism underlying the sympathetic eruptions. The magnitude of the impact is found to depend on the orientation of the erupting flux rope, with the largest impacts occurring when the flux rope is favorably oriented for reconnecting with the surrounding regions. Outside the CME expansion domain, the influence of the CME is mainly through field line compression or post-eruption relaxation. Based on our numerical experiments, we discuss a way to quantify the eruption impact, which could be useful for forecasting purposes.

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.

The most violent space weather events (eruptive solar flares and coronal mass ejections) are driven by the release of free magnetic energy stored in the solar corona. Energy can build up on timescales of hours to days, and then may be suddenly released in the form of a magnetic eruption, which then propagates through interplanetary space, possibly impacting the Earth's space environment. Can we use the observed evolution of the magnetic and velocity fields in the solar photosphere to model the evolution of the overlying solar coronal field, including the storage and release of magnetic energy in such eruptions? The objective of CGEM, the Coronal Global Evolutionary Model, funded by the NASA/NSF Space Weather Modeling program, is to develop and evaluate such a model for the evolution of the coronal magnetic field. The evolving coronal magnetic field can then be used as a starting point for magnetohydrodynamic (MHD) models of the corona, which can then be used to drive models of heliospheric evolution and predictions of magnetic field and plasma density conditions at 1AU.

Seeking to establish whether active region upflow material contributes to the slow solar wind, we examine in detail the plasma upflows from Active Region (AR)10978, which crossed the Sun's disc in the interval 8 to 16 December, 2007 during Carrington rotation (CR)2064. In previous work, using data from the Hinode/EUV Imaging Spectrometer, upflow velocity evolution was extensively studied as the region crossed the disc while a linear force-free magnetic extrapolation was used to confirm aspects of the velocity evolution and to establish the presence of quasi-separatrix layers at the upflow source areas. The plasma properties, temperature, density and first ionisation potential bias (FIP-bias) were measured with the spectrometer during the disc passage of the active region. Global potential field source surface (PFSS) models showed that AR 10978 was completely covered by the closed field of a helmet streamer that is part of the streamer belt. Thus it is not clear how any of the upflowing AR-associated plasma could reach the source surface at 2.5 R(Sun)and contribute to the slow solar wind. However a detailed examination of solar-wind in-situ data obtained by the Advanced Composition Explorer (ACE) spacecraft at the L1 point shows that the increase in O^7+/O^6+, C^6+/C^5+ and Fe/O - a FIP-bias proxy - are present before the heliospheric current sheet crossing. These increases, along with an accompanying reduction in proton velocity and an increase in density are characteristic of both AR and slow-wind plasma. Finally we describe a two-step reconnection process by which some of the upflowing plasma from the AR could reach the heliosphere.

The shapes of solar coronal loops are sensitive to the presence of electrical currents that are the carriers of the nonpotential energy available for impulsive activity. We use this information in a new method for modeling the coronal magnetic field of AR 11158 as a nonlinear force-free field (NLFFF). The observations used are coronal images around time of major flare activity on 2011/02/15, together with the surface line-of-sight magnetic field measurements. The data are from the Helioseismic and Magnetic Imager and Atmospheric Imaging Assembly (HMI and AIA, respectively) onboard the Solar Dynamics Observatory (SDO). The model fields are constrained to approximate the coronal loop configurations as closely as possible, while also subject to the force-free constraints. The method does not use transverse photospheric magnetic field components as input, and is thereby distinct from methods for modeling NLFFFs based on photospheric vector magnetograms. We validate the method using observations of AR 11158 at a time well before major flaring, and subsequently review the field evolution just prior to and following an X2.2 flare and associated eruption. The models indicate that the energy released during the instability is about $1\times10^{32}$ erg, consistent with what is needed to power such a large eruptive flare. Immediately prior to the eruption the model field contains a compact sigmoid bundle of twisted flux that is not present in the post-eruption models, which is consistent with the observations. The core of that model structure is twisted by $\approx0.9$ full turns about its axis.

The high-cadence, comprehensive view of the solar corona by SDO/AIA shows many events that are widely separated in space while occurring close together in time. In some cases, sets of coronal events are evidently causally related, while in many other instances indirect evidence can be found. We present case studies to highlight a variety of coupling processes involved in coronal events. We find that physical linkages between events do occur, but concur with earlier studies that these couplings appear to be crucial to understanding the initiation of major eruptive or explosive phenomena relatively infrequently. We note that the post-eruption reconfiguration time scale of the large-scale corona, estimated from the EUV afterglow, is on average longer than the mean time between CMEs, so that many CMEs originate from a corona that is still adjusting from a previous event. We argue that the coronal field is intrinsically global: current systems build up over days to months, the relaxation after eruptions continues over many hours, and evolving connections easily span much of a hemisphere. This needs to be reflected in our modeling of the connections from the solar surface into the heliosphere to properly model the solar wind, its perturbations, and the generation and propagation of solar energetic particles. However, the large-scale field cannot be constructed reliably by currently available observational resources. We assess the potential of high-quality observations from beyond Earth's perspective and advanced global modeling to understand the couplings between coronal events in the context of CMEs and solar energetic particle events.

We test a nonlinear force-free field (NLFFF) optimization code in spherical geometry using an analytical solution from Low and Lou. Several tests are run, ranging from idealized cases where exact vector field data are provided on all boundaries, to cases where noisy vector data are provided on only the lower boundary (approximating the solar problem). Analytical tests also show that the NLFFF code in the spherical geometry performs better than that in the Cartesian one when the field of view of the bottom boundary is large, say, $20^\circ \times 20^\circ$. Additionally, We apply the NLFFF model to an active region observed by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) both before and after an M8.7 flare. For each observation time, we initialize the models using potential field source surface (PFSS) extrapolations based on either a synoptic chart or a flux-dispersal model, and compare the resulting NLFFF models. The results show that NLFFF extrapolations using the flux-dispersal model as the boundary condition have slightly lower, therefore better, force-free and divergence-free metrics, and contain larger free magnetic energy. By comparing the extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the Atmospheric Imaging Assembly (AIA) on board SDO, we find that the NLFFF performs better than the PFSS not only for the core field of the flare productive region, but also for large EUV loops higher than 50 Mm.

We present a method for performing data-driven simulations of solar active region formation and evolution. The approach is based on magnetofriction, which evolves the induction equation assuming the plasma velocity is proportional to the Lorentz force. The simulations of active region coronal field are driven by temporal sequences of photospheric magnetograms from the Helioseismic Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO). Under certain conditions, the data-driven simulations produce flux ropes that are ejected from the modeled active region due to loss of equilibrium. Following the ejection of flux ropes, we find an enhancement of the photospheric horizontal field near the polarity inversion line. We also present a method for the synthesis of mock coronal images based on a proxy emissivity calculated from the current density distribution in the model. This method yields mock coronal images that are somewhat reminiscent of images of active regions taken by instruments such as SDO's Atmospheric Imaging Assembly (AIA) at extreme ultraviolet wavelengths.

The variable magnetic field of the solar photosphere exhibits periodic reversals as a result of dynamo activity occurring within the solar interior. We decompose the surface field as observed by both the Wilcox Solar Observatory and the Michelson Doppler Imager into its harmonic constituents, and present the time evolution of the mode coefficients for the past three sunspot cycles. The interplay between the various modes is then interpreted from the perspective of general dynamo theory, where the coupling between the primary and secondary families of modes is found to correlate with large-scale polarity reversals for many examples of cyclic dynamos. Mean-field dynamos based on the solar parameter regime are then used to explore how such couplings may result in the various long-term trends in the surface magnetic field observed to occur in the solar case.

The three-dimensional (3D) coordinates of stereoscopically triangulated loops provide strong constraints for magnetic field models of active regions in the solar corona. Here we use STEREO/A and B data from some 500 stereoscopically triangulated loops observed in four active regions (2007 Apr 30, May 9, May 19, Dec 11), together with SOHO/MDI line-of-sight magnetograms. We measure the average misalignment angle between the stereoscopic loops and theoretical magnetic field models, finding a mismatch of $\mu=19^\circ-46^\circ$ for a potential field model, which is reduced to $\mu=14^\circ-19^\circ$ for a non-potential field model parameterized by twist parameters. The residual error is commensurable with stereoscopic measurement errors ($\mu_{SE} \approx 8^\circ-12^\circ$). We developed a potential field code that deconvolves a line-of-sight magnetogram into three magnetic field components $(B_x, B_y, B_z)$, as well as a non-potential field forward-fitting code that determines the full length of twisted loops ($L \approx 50-300$ Mm), the number of twist turns (median $N_{twist}=0.06$), the nonlinear force-free $\alpha$-parameter (median $\alpha \approx 4 \times 10^{-11}$ cm$^{-1}$), and the current density (median $j_z \approx 1500$ Mx cm$^{-2}$ s$^{-1}$). All twisted loops are found to be far below the critical value for kink instability, and Joule dissipation of their currents is found be be far below the coronal heating requirement. The algorithm developed here, based on an analytical solution of nonlinear force-free fields that is accurate to second order (in the force-free parameter $\alpha$), represents the first code that enables fast forward-fitting to photospheric magnetograms and stereoscopically triangulated loops in the solar corona.

Ground level events (GLEs) occupy the high-energy end of gradual solar energetic particle (SEP) events. They are associated with coronal mass ejections (CMEs) and solar flares, but we still do not clearly understand the special conditions that produce these rare events. During Solar Cycle 23, a total of 16 GLEs were registered, using ground-based neutron monitor data. We first ask if these GLEs are clearly distinguishable from other SEP events observed from space. Setting aside possible difficulties in identifying all GLEs consistently, we then try to find observables which may unmistakably isolate these GLEs by studying the basic properties of the associated eruptions and the active regions (ARs) that produced them. It is found that neither the magnitudes of the CMEs and flares nor the complexities of the ARs give sufficient conditions for GLEs. It is possible to find CMEs, flares or ARs that are not associated with GLEs but that have more extreme properties than those associated with GLEs. We also try to evaluate the importance of magnetic field connection of the AR with Earth on the detection of GLEs and their onset times. Using the potential field source surface (PFSS) model, a half of the GLEs are found to be well-connected. However, the GLE onset time with respect to the onset of the associated flare and CME does not strongly depend on how well-connected the AR is. The GLE onset behavior may be largely determined by when and where the CME-driven shock develops. We could not relate the shocks responsible for the onsets of past GLEs with features in solar images, but the combined data from the Solar TErrestrial RElations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) have the potential to change this for GLEs that may occur in the rising phase of Solar Cycle 24.

Presently, many models of the coronal magnetic field rely on photospheric vector magnetograms but these data have been shown to be problematic as the sole boundary information for nonlinear force-free field (NLFFF) extrapolations. Magnetic fields in the corona manifest themselves in high-energy images (X-rays and EUV) in the shapes of coronal loops, providing an additional constraint that at present is not used due to the mathematical complications of incorporating such input into numerical models. Projection effects and the limited number of usable loops further complicate the use of coronal information. We develop and test an algorithm to use images showing coronal loops in the modeling of the solar coronal magnetic field. We first fit projected field lines with field lines of constant-\als force-free fields to approximate the three-dimensional distribution of currents in the corona along a sparse set of trajectories. We then apply a Grad-Rubin-like iterative technique to obtain a volume-filling nonlinear force-free model of the magnetic field, modifying method presented in \citet{Wheatland2007}. We thoroughly test the technique on known analytical and solar-like model magnetic fields previously used for comparing different extrapolation techniques \citep{Schrijver2006, Schrijver2008} and compare the results with those obtained by presently available methods that rely only on the photospheric data. We conclude that we have developed a functioning method of modeling the coronal magnetic field by combining the line-of-sight component of photospheric magnetic field with information from coronal images. Vector magnetograms over the full or partial photospheric boundary of the numerical domain could optionally be used.

We examine the role of small-scale granulation in helping to drive supergranulation and even larger scales of convection. The granulation is modeled as localized cooling events introduced at the upper boundary of a 3-D simulation of compressible convection in a rotating spherical shell segment. With a sufficient number of stochastic cooling events compared to uniform cooling, we find that supergranular scales are realized, along with a differential rotation that becomes increasingly solar-like.

The coming decade will see the routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness. To capitalize on the new (or soon to be available) facilities such as SDO, ATST and FASR, and the challenges they present in the visualization and synthesis of multi-wavelength datasets, we propose that realistic, sophisticated, 3D active region and flare modeling is timely and critical, and will be a forefront of coronal studies over the coming decade. To make such modeling a reality, a broad, concerted effort is needed to capture the wealth of information resulting from the data, develop a synergistic modeling effort, and generate the necessary visualization, interpretation and model-data comparison tools to accurately extract the key physics.

We constrain the velocity spectral distribution of global-scale solar convective cells at depth using techniques of local helioseismology. We calibrate the sensitivity of helioseismic waves to large-scale convective cells in the interior by analyzing simulations of waves propagating through a velocity snapshot of global solar convection via methods of time-distance helioseismology. Applying identical analysis techniques to observations of the Sun, we are able to bound from above the magnitudes of solar convective cells as a function of spatial convective scale. We find that convection at a depth of $r/R_\odot = 0.95$ with spatial extent $\ell <20$, where $\ell$ is the spherical harmonic degree, comprise weak flow systems, on the order of 15 m/s or less. Convective features deeper than $r/R_\odot = 0.95$ are more difficult to image due to the rapidly decreasing sensitivity of helioseismic waves.

The magnetic field of the solar corona has a large-scale dipole character, which maps into the bipolar field in the solar wind. Using standard representations of the coronal field, we show that high-energy ions can be trapped stably in these large-scale closed fields. The drift shells that describe the conservation of the third adiabatic invariant may have complicated geometries. Particles trapped in these zones would resemble the Van Allen Belts and could have detectable consequences. We discuss potential sources of trapped particles.

Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models have been found to have markedly different field line configurations and to provide widely varying estimates of the magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, force-free vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic field data covering larger areas are needed so that more electric currents associated with the full active regions of interest are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems.

The solar magnetic field is key to understanding the physical processes in the solar atmosphere. Nonlinear force-free codes have been shown to be useful in extrapolating the coronal field upward from underlying vector boundary data. However, we can only measure the magnetic field vector routinely with high accuracy in the photosphere, and unfortunately these data do not fulfill the force-free condition. We must therefore apply some transformations to these data before nonlinear force-free extrapolation codes can be self-consistently applied. To this end, we have developed a minimization procedure that yields a more chromosphere-like field, using the measured photospheric field vectors as input. The procedure includes force-free consistency integrals, spatial smoothing, and -- newly included in the version presented here -- an improved match to the field direction as inferred from fibrils as can be observed in, e.g., chromospheric H$\alpha$ images. We test the procedure using a model active-region field that included buoyancy forces at the photospheric level. The proposed preprocessing method allows us to approximate the chromospheric vector field to within a few degrees and the free energy in the coronal field to within one percent.

The solar magnetic field is key to understanding the physical processes in the solar atmosphere. Nonlinear force-free codes have been shown to be useful in extrapolating the coronal field from underlying vector boundary data [see Schrijver et al. 2006 for an overview]. However, we can only measure the magnetic field vector routinely with high accuracy in the photosphere with, e.g., Hinode/SOT, and unfortunately these data do not fulfill the force-free consistency condition as defined by Aly (1989). We must therefore apply some transformations to these data before nonlinear force-free extrapolation codes can be legitimately applied. To this end, we have developed a minimization procedure that uses the measured photospheric field vectors as input to approximate a more chromospheric like field The method was dubbed preprocessing. See Wiegelmann et al. 2006 for details]. The procedure includes force-free consistency integrals and spatial smoothing. The method has been intensively tested with model active regions [see Metcalf et al. 2008] and been applied to several ground based vector magnetogram data before. Here we apply the preprocessing program to photospheric magnetic field measurements with the Hinode/SOT instrument.

Solar flares and coronal mass ejections are associated with rapid changes in field connectivity and powered by the partial dissipation of electrical currents in the solar atmosphere. A critical unanswered question is whether the currents involved are induced by the motion of pre-existing atmospheric magnetic flux subject to surface plasma flows, or whether these currents are associated with the emergence of flux from within the solar convective zone. We address this problem by applying state-of-the-art nonlinear force-free field (NLFFF) modeling to the highest resolution and quality vector-magnetographic data observed by the recently launched Hinode satellite on NOAA Active Region 10930 around the time of a powerful X3.4 flare. We compute 14 NLFFF models with 4 different codes and a variety of boundary conditions. We find that the model fields differ markedly in geometry, energy content, and force-freeness. We discuss the relative merits of these models in a general critique of present abilities to model the coronal magnetic field based on surface vector field measurements. For our application in particular, we find a fair agreement of the best-fit model field with the observed coronal configuration, and argue (1) that strong electrical currents emerge together with magnetic flux preceding the flare, (2) that these currents are carried in an ensemble of thin strands, (3) that the global pattern of these currents and of field lines are compatible with a large-scale twisted flux rope topology, and (4) that the ~10^32 erg change in energy associated with the coronal electrical currents suffices to power the flare and its associated coronal mass ejection.

The global scales of solar convection are studied through three-dimensional simulations of compressible convection carried out in spherical shells of rotating fluid which extend from the base of the convection zone to within 15 Mm of the photosphere. Such modelling at the highest spatial resolution to date allows study of distinctly turbulent convection, revealing that coherent downflow structures associated with giant cells continue to play a significant role in maintaining the strong differential rotation that is achieved. These giant cells at lower latitudes exhibit prograde propagation relative to the mean zonal flow, or differential rotation, that they establish, and retrograde propagation of more isotropic structures with vortical character at mid and high latitudes. The interstices of the downflow networks often possess strong and compact cyclonic flows. The evolving giant-cell downflow systems can be partly masked by the intense smaller scales of convection driven closer to the surface, yet they are likely to be detectable with the helioseismic probing that is now becoming available. Indeed, the meandering streams and varying cellular subsurface flows revealed by helioseismology must be sampling contributions from the giant cells, yet it is difficult to separate out these signals from those attributed to the faster horizontal flows of supergranulation. To aid in such detection, we use our simulations to describe how the properties of giant cells may be expected to vary with depth, how their patterns evolve in time, and analyze the statistical features of correlations within these complex flow fields.

The "Frequency-Agile Solar Radiotelescope" (FASR) litteraly opens up a new dimension in addition to the 3D Euclidian geometry: the frequency dimension. The 3D geometry is degenerated to 2D in all images from astronomical telescopes, but the additional frequency dimension allows us to retrieve the missing third dimension by means of physical modeling. We call this type of 3D reconstruction "Frequency Tomography". In this study we simulate a realistic 3D model of an active region, composed of 500 coronal loops with the 3D geometry [x(s),y(s),z(s)] constrained by magnetic field extrapolations and the physical parameters of the density n_e(s) and temperature T_e(s) given by hydrostatic solutions. We simulate a series of 20 radio images in a frequency range of f=0.1-10 GHz, anticipating the capabilities of FASR, and investigate what physical information can be retrieved from such a dataset. We discuss also forward-modeling of the chromospheric and Quiet Sun density and temperature structure, another primary goal of future FASR science.

The differential rotation of the sun, as deduced from helioseismology, exhibits a prominent radial shear layer near the top of the convection zone wherein negative radial gradients of angular velocity are evident in the low- and mid-latitude regions spanning the outer 5% of the solar radius. Supergranulation and related scales of turbulent convection are likely to play a significant role in the maintenance of such radial gradients, and may influence dynamics on a global scale in ways that are not yet understood. To investigate such dynamics, we have constructed a series of three-dimensional numerical simulations of turbulent compressible convection within spherical shells, dealing with shallow domains to make such modeling computationally tractable. These simulations are the first models of solar convection in a spherical geometry that can explicitly resolve both the largest dynamical scales of the system (of order the solar radius) as well as smaller-scale convective overturning motions comparable in size to solar supergranulation (20--40 Mm). We find that convection within these simulations spans a large range of horizontal scales, and that the radial angular velocity gradient in these models is typically negative, especially in low- and mid-latitude regions. Analyses of the angular momentum transport indicates that such gradients are maintained by Reynolds stresses associated with the convection, transporting angular momentum inward to balance the outward transport achieved by viscous diffusion and large-scale flows in the meridional plane. We suggest that similar mechanisms associated with smaller-scale convection in the sun may contribute to the maintenance of the observed radial shear layer located immediately below the solar photosphere.