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.

In this work we study how the input data cadence affects the photospheric energy and helicity injection estimates in eruptive NOAA active region 11158. We sample the novel 2.25-minute vector magnetogram and Dopplergram data from the \emph{Helioseismic and Magnetic Imager} (HMI) instrument onboard the \emph{Solar Dynamics Observatory} (SDO) spacecraft to create input datasets of variable cadences ranging from 2.25 minutes to 24 hours. We employ state-of-the-art data processing, velocity and electric field inversion methods for deriving estimates of the energy and helicity injections from these datasets. We find that the electric field inversion methods that reproduce the observed magnetic field evolution through the use of Faraday's law are more stable against variable cadence: the PDFI (PTD-Doppler-FLCT-Ideal) electric field inversion method produces consistent injection estimates for cadences from 2.25 minutes up to 2 hours, implying that the photospheric processes acting on time scales below 2 hours contribute little to the injections, or that they are below the sensitivity of the input data and the PDFI method. On other hand, the electric field estimate derived from the output of DAVE4VM (Differential Affine Velocity Estimator for Vector Magnetograms), which does not fulfil Faraday's law exactly, produces significant variations in the energy and helicity injection estimates in the 2.25-minute to 2-hour cadence range. We present also a third, novel DAVE4VM-based electric field estimate, which corrects the poor inductivity of the raw DAVE4VM estimate. This method is less sensitive to the changes of cadence, but still faces significant issues for the lowest of considered cadences ($\geq$2 hours). We find several potential problems in both PDFI- and DAVE4VM-based injection estimates and conclude that the quality of both should be surveyed further in controlled environments.

Stealth coronal mass ejections (CMEs) are events in which there are almost no observable signatures of the CME eruption in the low corona but often a well-resolved slow flux rope CME observed in the coronagraph data. We present results from a three-dimensional numerical magnetohydrodynamics (MHD) simulation of the 1--2 June 2008 slow streamer blowout CME that Robbrecht et al. (2009) called "the CME from nowhere." We model the global coronal structure using a 1.4 MK isothermal solar wind and a low-order potential field source surface representation of the Carrington Rotation 2070 magnetogram synoptic map. The bipolar streamer belt arcade is energized by simple shearing flows applied in the vicinity of the helmet streamer's polarity inversion line. The flows are large scale and impart a shear typical of that expected from the differential rotation. The slow expansion of the energized helmet streamer arcade results in the formation of a radial current sheet. The subsequent onset of expansion-induced flare reconnection initiates the stealth CME while gradually releasing the stored magnetic energy. We present favorable comparisons between our simulation results and the multiviewpoint SOHO-LASCO (Large Angle and Spectrometric Coronagraph) and STEREO-SECCHI (Sun Earth Connection Coronal and Heliospheric Investigation) coronagraph observations of the preeruption streamer structure and the initiation and evolution of the stealth streamer blowout CME.

The minimum-energy configuration for the magnetic field above the solar photosphere is curl-free (hence, by Ampere's law, also current-free), so can be represented as the gradient of a scalar potential. Since magnetic fields are divergence free, this scalar potential obeys Laplace's equation, given an appropriate boundary condition (BC). With measurements of the full magnetic vector at the photosphere, it is possible to employ either Neumann or Dirichlet BCs there. Historically, the Neumann BC was used with available line-of-sight magnetic field measurements, which approximate the radial field needed for the Neumann BC. Since each BC fully determines the 3D vector magnetic field, either choice will, in general, be inconsistent with some aspect of the observed field on the boundary, due to the presence of both currents and noise in the observed field. We present a method to combine solutions from both Dirichlet and Neumann BCs to determine a hybrid, "least-squares" potential field, which minimizes the integrated square of the residual between the potential and actual fields. This has advantages in both not overfitting the radial field used for the Neumann BC, and maximizing consistency with the observations. We demonstrate our methods with SDO/HMI vector magnetic field observations of AR 11158, and find that residual discrepancies between the observed and potential fields are significant, and are consistent with nonzero horizontal photospheric currents. We also analyze potential fields for two other active regions observed with two different vector magnetographs, and find that hybrid potential fields have significantly less energy than the Neumann fields in every case --- by more than 10^(32) erg in some cases. This has major implications for estimates of free magnetic energy in coronal field models, e.g., non-linear force-free field extrapolations.

How much electromagnetic energy crosses the photosphere in evolving solar active regions? With the advent of high-cadence vector magnetic field observations, addressing this fundamental question has become tractable. In this paper, we apply the "PTD-Doppler-FLCT-Ideal" (PDFI) electric field inversion technique of Kazachenko et al. (2014) to a 6-day HMI/SDO vector magnetogram and Doppler velocity sequence, to find the electric field and Poynting flux evolution in active region NOAA 11158, which produced an X2.2 flare early on 2011 February 15. We find photospheric electric fields ranging up to $2$ V/cm. The Poynting fluxes range from $[-0.6$ to $2.3]\times10^{10}$ ergs$\cdot$cm$^{-2}$s$^{-1}$, mostly positive, with the largest contribution to the energy budget in the range of $[10^9$-$10^{10}]$ ergs$\cdot$cm$^{-2}$s$^{-1}$. Integrating the instantaneous energy flux over space and time, we find that the total magnetic energy accumulated above the photosphere from the initial emergence to the moment before the X2.2 flare to be $E=10.6\times10^{32}$ ergs, which is partitioned as $2.0$ and $8.6\times10^{32}$ ergs, respectively, between free and potential energies. Those estimates are consistent with estimates from preflare non-linear force-free field (NLFFF) extrapolations and the Minimum Current Corona estimates (MCC), in spite of our very different approach. This study of photospheric electric fields demonstrates the potential of the PDFI approach for estimating Poynting fluxes and opens the door to more quantitative studies of the solar photosphere and more realistic data-driven simulations of coronal magnetic field evolution.

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.

We model the evolution of buoyant magnetic flux tubes in the Sun's convection zone. A flux tube is assumed to lie initially near the top of the stably stratified radiative core below the convection zone, but a segment of it is perturbed into the convection zone by gradual heating and convective overshoot motions. The ends ("footpoints") of the segment remain anchored at the base of the convection zone, and if the segment is sufficiently long, it may be buoyantly unstable, rising through the convection zone in a short time. The length of the flux tube determines the ratio of buoyancy to magnetic tension: short loops of flux are arrested before reaching the top of the convection zone, while longer loops emerge to erupt through the photosphere. Using Spruit's convection zone model, we compute the minimum footpoint separation $L_c$ required for erupting flux tubes. We explore the dependence of $L_c$ on the initial thermal state of the perturbed flux tube segment and on its initial magnetic field strength. Following an investigation of thermal diffusion time scales and the dynamic rise times of unstable flux tube segments, we conclude that the most likely origin for magnetic flux which erupts to the surface is from short length scale perturbations ($L < L_c$) which are initially stable, but which are subsequently destabilized either by diffusion of heat into the tube or by stretching of the anchor points until $L$ just exceeds $L_c$. In either case, the separation of the anchor points of the emergent tube should lie between the critical distance for a tube in mechanical equilibrium and one in thermal equilibrium. Finally, after comparing the dispersion of dynamic rise times with the much shorter observed active region formation time scales, we conclude that active regions form from the emergence of a single flux tube segment.

Solar active region (AR) 12192 of October 2014 hosts the largest sunspot group in 24 years. It is the most prolific flaring site of Cycle 24, but surprisingly produced no coronal mass ejection (CME) from the core region during its disk passage. Here, we study the magnetic conditions that prevented eruption and the consequences that ensued. We find AR 12192 to be "big but mild"; its core region exhibits weaker non-potentiality, stronger overlying field, and smaller flare-related field changes compared to two other major flare-CME-productive ARs (11429 and 11158). These differences are present in the intensive-type indices (e.g., means) but generally not the extensive ones (e.g., totals). AR 12192's large amount of magnetic free energy does not translate into CME productivity. The unexpected behavior suggests that AR eruptiveness is limited by some relative measure of magnetic non-potentiality over the restriction of background field, and that confined flares may leave weaker photospheric and coronal imprints compared to their eruptive counterparts.

Solar Active Region NOAA 11158 has hosted a number of strong flares, including one X2.2 event. The complexity of current density and current helicity are studied through cancellation analysis of their sign-singular measure, which features power-law scaling. Spectral analysis is also performed, revealing the presence of two separate scaling ranges with different spectral index. The time evolution of parameters is discussed. Sudden changes of the cancellation exponents at the time of large flares, and the presence of correlation with EUV and X-ray flux, suggest that eruption of large flares can be linked to the small scale properties of the current structures.

Photospheric electric fields, estimated from sequences of vector magnetic field and Doppler measurements, can be used to estimate the flux of magnetic energy (the Poynting flux) into the corona and as time-dependent boundary conditions for dynamic models of the coronal magnetic field. We have modified and extended an existing method to estimate photospheric electric fields that combines a poloidal-toroidal (PTD) decomposition of the evolving magnetic field vector with Doppler and horizontal plasma velocities. Our current, more comprehensive method, which we dub the "{\bf P}TD-{\bf D}oppler-{\bf F}LCT {\bf I}deal" (PDFI) technique, can now incorporate Doppler velocities from non-normal viewing angles. It uses the \texttt{FISHPACK} software package to solve several two-dimensional Poisson equations, a faster and more robust approach than our previous implementations. Here, we describe systematic, quantitative tests of the accuracy and robustness of the PDFI technique using synthetic data from anelastic MHD (\texttt{ANMHD}) simulations, which have been used in similar tests in the past. We find that the PDFI method has less than $1%$ error in the total Poynting flux and a $10%$ error in the helicity flux rate at a normal viewing angle $(\theta=0$) and less than $25%$ and $10%$ errors respectively at large viewing angles ($\theta<60^\circ$). We compare our results with other inversion methods at zero viewing angle, and find that our method's estimates of the fluxes of magnetic energy and helicity are comparable to or more accurate than other methods. We also discuss the limitations of the PDFI method and its uncertainties.

The zero point of measured photospheric Doppler shifts is uncertain for at least two reasons: instrumental variations (from, e.g., thermal drifts), and the convective blueshift, a known correlation between intensity and upflows. Calibrated Doppler velocities would be useful for (i) improving estimates of the Poynting flux of magnetic energy across the photosphere, and (ii) constraining processes underlying flux cancellation, the mutual apparent loss of magnetic flux in closely spaced, opposite-polarity magnetogram features. We present a method to absolutely calibrate line-of-sight (LOS) velocities in solar active regions (ARs) near disk center using three successive vector magnetograms and one Dopplergram coincident with the central magnetogram. It exploits the fact that Doppler shifts measured along polarity inversion lines (PILs) of the LOS magnetic field determine one component of the velocity perpendicular to the magnetic field, and optimizes consistency between changes in LOS flux near PILs and the transport of transverse magnetic flux by LOS velocities, assuming ideal electric fields govern the magnetic evolution. We apply our method to vector magnetograms of AR 11158, observed by the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory, and find clear evidence of offsets in the Doppler zero point, in the range of 50 -- 550 m s$^{-1}$. In addition, we note that a simpler calibration can be determined from an LOS magnetogram and Dopplergram pair from the median Doppler velocity among all near-disk-center PIL pixels. We briefly discuss shortcomings in our initial implementation, and suggest ways to address these. In addition, as a step in our data reduction, we discuss use of temporal continuity in the transverse magnetic field direction to correct apparently spurious fluctuations in resolution of the 180$^\circ$ ambiguity.

We compute the change in the Lorentz force integrated over the outer solar atmosphere implied by observed changes in vector magnetograms that occur during large, eruptive solar flares. This force perturbation should be balanced by an equal and opposite force perturbation acting on the solar photosphere and solar interior. The resulting expression for the estimated force change in the solar interior generalizes the earlier expression presented by Hudson, Fisher and Welsch (CS-383, ASP, 221, 2008), providing horizontal as well as vertical force components, and provides a more accurate result for the vertical component of the perturbed force. We show that magnetic eruptions should result in the magnetic field at the photosphere becoming more horizontal, and hence should result in a downward (towards the solar interior) force change acting on the photosphere and solar interior, as recently argued from an analysis of magnetogram data by Wang and Liu (Astrophys. J. Lett. 716, L195, 2010). We suggest the existence of an observational relationship between the force change computed from changes in the vector magnetograms, the outward momentum carried by the ejecta from the flare, and the properties of the helioseismic disturbance driven by the downward force change. We use the impulse driven by the Lorentz-force change in the outer solar atmosphere to derive an upper limit to the mass of erupting plasma that can escape from the Sun. Finally, we compare the expected Lorentz-force change at the photosphere with simple estimates from flare-driven gasdynamic disturbances and from an estimate of the perturbed pressure from radiative backwarming of the photosphere in flaring conditions.

We present a series of numerical simulations of the quiet Sun plasma threaded by magnetic fields that extend from the upper convection zone into the low corona. We discuss an efficient, simplified approximation to the physics of optically thick radiative transport through the surface layers, and investigate the effects of convective turbulence on the magnetic structure of the Sun's atmosphere in an initially unipolar (open field) region. We find that the net Poynting flux below the surface is on average directed toward the interior, while in the photosphere and chromosphere the net flow of electromagnetic energy is outward into the solar corona. Overturning convective motions between these layers driven by rapid radiative cooling appears to be the source of energy for the oppositely directed fluxes of electromagnetic energy.

The availability of vector magnetogram sequences with sufficient accuracy and cadence to estimate the time derivative of the magnetic field allows us to use Faraday's law to find an approximate solution for the electric field in the photosphere, using a Poloidal-Toroidal Decomposition (PTD) of the magnetic field and its partial time derivative. Without additional information, however, the electric field found from this technique is under-determined -- Faraday's law provides no information about the electric field that can be derived the gradient of a scalar potential. Here, we show how additional information in the form of line-of-sight Doppler flow measurements, and motions transverse to the line-of-sight determined with ad-hoc methods such as local correlation tracking, can be combined with the PTD solutions to provide much more accurate solutions for the solar electric field, and therefore the Poynting flux of electromagnetic energy in the solar photosphere. Reliable, accurate maps of the Poynting flux are essential for quantitative studies of the buildup of magnetic energy before flares and coronal mass ejections.

Solar flares are currently understood as the explosive release of energy stored in the form of stressed magnetic fields. In many cases, the released energy seems to take the form of large numbers of electrons accelerated to high energies or alternatively plasma heated to very high temperatures. The transport of this energy into the remaining portion of the atmosphere results in violent mass motion and strong emission across the electromagnetic spectrum. One important phenomenon observed during flares is the appearance in coronal magnetic loops of large amounts of upflowing, soft X-ray emitting plasma. It is believed that this is due to chromospheric evaporation, the process of heating cool chromospheric material beyond its ability to radiate. The pressure increase in the evaporated plasma leads to a number of interesting phenomena in the flare chromosphere. The sudden pressure increase initiates a downward moving "chromospheric condensation", an overdense region which gradually decelerates as it accretes material and propagates into the gravitationally stratified chromosphere. Solutions to an equation of motion for this condensation shows that its motion decays after about one minute of propagation into the chromosphere. When the front of this downflowing region is supersonic relative to the atmosphere ahead of it, a radiating shock will form. If the downflow is rapid enough, the shock strength should be sufficient to excite UV radiation normally associated with the transition region, and furthermore, the radiating shock will be brighter than the transition region. These results lead to a number of observationally testable relationships between the optical and ultraviolet spectra from the condensation and radiating shock.

Determining the electric field (E-field) distribution on the Sun's photosphere is essential for quantitative studies of how energy flows from the Sun's photosphere, through the corona, and into the heliosphere. This E-field also provides valuable input for data-driven models of the solar atmosphere and the Sun-Earth system. We show how Faraday's Law can be used with observed vector magnetogram time series to estimate the photospheric E-field, an ill-posed inversion problem. Our method uses a "poloidal-toroidal decomposition" (PTD) of the time derivative of the vector magnetic field. The PTD solutions are not unique; the gradient of a scalar potential can be added to the PTD E-field without affecting consistency with Faraday's Law. We present an iterative technique to determine a potential function consistent with ideal MHD evolution; but this E-field is also not a unique solution to Faraday's Law. Finally, we explore a variational approach that minimizes an energy functional to determine a unique E-field, similar to Longcope's "Minimum Energy Fit". The PTD technique, the iterative technique, and the variational technique are used to estimate E-fields from a pair of synthetic vector magnetograms taken from an MHD simulation; and these E-fields are compared with the simulation's known electric fields. These three techniques are then applied to a pair of vector magnetograms of solar active region NOAA AR8210, to demonstrate the methods with real data.

We introduce two new methods that are designed to improve the realism and utility of large, active region-scale 3D MHD models of the solar atmosphere. We apply these methods to RADMHD, a code capable of modeling the Sun's upper convection zone, photosphere, chromosphere, transition region, and corona within a single computational volume. We first present a way to approximate the physics of optically-thick radiative transfer without having to take the computationally expensive step of solving the radiative transfer equation in detail. We then briefly describe a rudimentary assimilative technique that allows a time series of vector magnetograms to be directly incorporated into the MHD system.

Where does solar flare energy come from? More specifically, assuming that the ultimate source of flare energy is mechanical energy in the convection zone, how is this translated into energy dissipated or stored in the corona? This question appears to have been given relatively little thought, as attention has been focussed predominantly on mechanisms for the rapid dissipation of coronal magnetic energy by way of MHD instabilities and plasma micro instabilities. We consider three types of flare theory: the steady state "photospheric dynamo" model in which flare power represents coronal dissipation of currents generated simultaneously by sub-photospheric flows; the "magnetic energy storage" model where sub-photospheric flows again induce coronal currents but which in this case are built up over a longer period before being released suddenly; and "emerging flux" models, in which new magnetic flux rising to the photosphere already contains free energy, and does not require subsequent stressing by photospheric motions. We conclude that photospheric dynamos can power only very minor flares; that coronal energy storage can in principle meet the requirements of a major flare, although perhaps not the very largest flares, but that difficulties in coupling efficiently to the energy source may limit this mechanism to moderate sized flares; and that emerging magnetic flux tubes, generated in the solar interior, can carry sufficient free energy to power even the largest flares ever observed.

We estimated photospheric velocities by separately applying the Fourier Local Correlation Tracking (FLCT) and Differential Affine Velocity Estimator (DAVE) methods to 2708 co-registered pairs of SOHO/MDI magnetograms, with nominal 96-minute cadence and ~2" pixels, from 46 active regions (ARs) from 1996-1998 over the time interval t45 when each AR was within 45^o of disk center. For each magnetogram pair, we computed the average estimated radial magnetic field, B; and each tracking method produced an independently estimated flow field, u. We then quantitatively characterized these magnetic and flow fields by computing several extensive and intensive properties of each; extensive properties scale with AR size, while intensive properties do not depend directly on AR size. Intensive flow properties included moments of speeds, horizontal divergences, and radial curls; extensive flow properties included sums of these properties over each AR, and a crude proxy for the ideal Poynting flux, the total |u| B^2. Several magnetic quantities were also computed, including: total unsigned flux; a measure of the amount of unsigned flux near strong-field polarity inversion lines, R; and the total B^2. Next, using correlation and discriminant analysis, we investigated the associations between these properties and flares from the GOES flare catalog, when averaged over both t45 and shorter time windows, of 6 and 24 hours. We found R and total |u| B^2 to be most strongly associated with flares; no intensive flow properties were strongly associated with flares.

We describe the computational techniques employed in the recently updated Fourier local correlation tracking (FLCT) method. The FLCT code is then evaluated using a series of simple, 2D, known flow patterns that test its accuracy and characterize its errors.

Estimates of velocities from time series of photospheric and/or chromospheric vector magnetograms can be used to determine fluxes of magnetic energy (the Poynting flux) and helicity across the magnetogram layer, and to provide time-dependent boundary conditions for data-driven simulations of the solar atmosphere above this layer. Velocity components perpendicular to the magnetic field are necessary both to compute these transport rates and to derive model boundary conditions. Here, we discuss some possible approaches to estimating perpendicular flows from magnetograms. Since Doppler shifts contain contributions from flows parallel to the magnetic field, perpendicular velocities are not generally recoverable from Doppler shifts alone. The induction equation's vertical component relates evolution in $B_z$ to the perpendicular flow field, but has a finite null space, meaning some ``null'' flows, e.g., motions along contours of normal field, do not affect $B_z$. Consequently, additional information is required to accurately specify the perpendicular flow field. Tracking methods, which analyze $\partial_t B_z$ in a neighborhood, have a long heritage, but other approaches have recently been developed. In a recent paper, several such techniques were tested using synthetic magnetograms from MHD simulations. Here, we use the same test data to characterize: 1) the ability of the induction equation's normal component, by itself, to estimate flows; and 2) a tracking method's ability to recover flow components that are perpendicular to $\mathbf{B}$ and parallel to contours of $B_z$. This work has been supported by NASA Heliophysics Theory grant NNG05G144G.

The objective of this paper is to investigate whether a convective dynamo can account quantitatively for the observed lower limit of X-ray surface flux in solar-type main sequence stars. Our approach is to use 3D numerical simulations of a turbulent dynamo driven by convection to characterize the dynamic behavior, magnetic field strengths, and filling factors in a non-rotating stratified medium, and to predict these magnetic properties at the surface of cool stars. We use simple applications of stellar structure theory for the convective envelopes of main-sequence stars to scale our simulations to the outer layers of stars in the F0--M0 spectral range, which allows us to estimate the unsigned magnetic flux on the surface of non-rotating reference stars. With these estimates we use the recent results of \citet{Pevtsov03} to predict the level of X-ray emission from such a turbulent dynamo, and find that our results compare well with observed lower limits of surface X-ray flux. If we scale our predicted X-ray fluxes to \ion{Mg}{2} fluxes we also find good agreement with the observed lower limit of chromospheric emission in K dwarfs. This suggests that dynamo action from a convecting, non-rotating plasma is a viable alternative to acoustic heating models as an explanation for the basal emission level seen in chromospheric, transition region, and coronal diagnostics from late-type stars.

We present three-dimensional MHD simulations of buoyant magnetic flux tubes that rise through a stratified model convection zone in the presence of solar rotation. The equations of MHD are solved in the anelastic approximation, and the results are used to determine the effects of solar rotation on the dynamic evolution an Omega-loop. We find that the Coriolis force significantly suppresses the degree of fragmentation at the apex of the loop during its ascent toward the photosphere. If the initial axial field strength of the tube is reduced, then, in the absence of forces due to convective motions, the degree of apex fragmentation is also reduced. We show that the Coriolis force slows the rise of the tube, and induces a retrograde flow in both the magnetized and unmagnetized plasma of an emerging active region. Observationally, we predict that this flow will appear to originate at the leading polarity, and will terminate at the trailing polarity.

We present three-dimensional numerical simulations of the rise and fragmentation of twisted, initially horizontal magnetic flux tubes which evolve into emerging Omega-loops. The flux tubes rise buoyantly through an adiabatically stratified plasma that represents the solar convection zone. The MHD equations are solved in the anelastic approximation, and the results are compared with studies of flux tube fragmentation in two dimensions. We find that if the initial amount of field line twist is below a critical value, the degree of fragmentation at the apex of a rising Omega-loop depends on its three-dimensional geometry: the greater the apex curvature of a given Omega-loop, the lesser the degree of fragmentation of the loop as it approaches the photosphere. Thus, the amount of initial twist necessary for the loop to retain its cohesion can be reduced substantially from the two-dimensional limit. The simulations also suggest that as a fragmented flux tube emerges through a relatively quiet portion of the solar disk, extended crescent-shaped magnetic features of opposite polarity should form and steadily recede from one another. These features eventually coalesce after the fragmented portion of the Omega-loop emerges through the photosphere.