The RHESSI Microflare Height Distribution

We present the first in-depth statistical survey of flare source heights observed by RHESSI. Flares were found using a flare-finding algorithm designed to search the 6 – 10 keV count-rate when RHESSI’s full sensitivity was available in order to find the smallest events (Christe et al. in Astrophys. J. 677, 1385, 2008). Between March 2002 and March 2007, a total of 25 006 events were found. Source locations were determined in the 4 – 10 keV, 10 – 15 keV, and 15 – 30 keV energy ranges for each event. In order to extract the height distribution from the observed projected source positions, a forward-fit model was developed with an assumed source height distribution where height is measured from the photosphere. We find that the best flare height distribution is given by g(h)∝exp (−h/λ) where λ=6.1±0.3 Mm is the scale height. A power-law height distribution with a negative power-law index, γ=3.1±0.1 is also consistent with the data. Interpreted as thermal loop-top sources, these heights are compared to loops generated by a potential-field model (PFSS). The measured flare heights distribution are found to be much steeper than the potential-field loop height distribution, which may be a signature of the flare energization process.     

Imaging Observations of Quasi-Periodic Pulsatory Nonthermal Emission in Two-Ribbon Solar Flares

Using RHESSI and some auxiliary observations we examine possible connections between the spatial and temporal structure of nonthermal hard X-ray (HXR) emission sources from the two-ribbon flares of 29 May 2003 and 19 January 2005. In each of these events quasi-periodic pulsations (QPP) with time period of 1 – 3 minutes are evident in both hard X rays and microwaves. The sources of nonthermal HXR emission are situated mainly at the footpoints of the flare arcade loops observed by TRACE and the SOHO/EIT instrument in the EUV range. At least one of the sources moves systematically during and after the QPP phase in each flare. The sources move predominantly parallel to the magnetic inversion line during the 29 May flare and along flare ribbons during the QPP phase of both flares. By contrast, the sources start to show movement perpendicular to the flare ribbons with velocity comparable to that along the ribbons’ movement after the QPP phase. The sources of each pulse are localized in distinct parts of the ribbon during the QPP phase. The measured velocity of the sources and the estimated energy release rate do not correlate well with the flux of the HXR emission calculated from these sources. The sources of microwaves and thermal HXRs are situated near the apex of the flare loop arcade and are not stationary either. Almost all of the QPP as well as some pulses of nonthermal HXR emission during the post-QPP phase reveal soft – hard – soft spectral behavior, indicating separate acts of electron acceleration and injection. In our opinion at least two different flare scenarios based on the Nakariakov et al. (2006, Astron. Astrophys. 452, 343) model and on the idea of current-carrying loop coalescence are suitable for interpreting the observations. However, it is currently not possible to choose between them owing to observational limitations.     

Time-Dependent Stochastic Modeling of Solar Active Region Energy

A time-dependent model for the energy of a flaring solar active region is presented based on an existing stochastic jump-transition model (Wheatland and Glukhov in Astrophys. J. 494, 858, 1998; Wheatland in Astrophys. J. 679, 1621, 2008 and Solar Phys. 255, 211, 2009). The magnetic free energy of an active region is assumed to vary in time due to a prescribed (deterministic) rate of energy input and prescribed (random) jumps downwards in energy due to flares. The existing model reproduces observed flare statistics, in particular flare frequency – size and waiting-time distributions, but modeling presented to date has considered only the time-independent choices of constant energy input and constant flare-transition rates with a power-law distribution in energy. These choices may be appropriate for a solar active region producing a constant mean rate of flares. However, many solar active regions exhibit time variation in their flare productivity, as exemplified by NOAA active region (AR) 11029, observed during October – November 2009 (Wheatland in Astrophys. J. 710, 1324, 2010). Time variation is incorporated into the jump-transition model for two cases: (1) a step change in the rates of flare transitions, and (2) a step change in the rate of energy supply to the system. Analytic arguments are presented describing the qualitative behavior of the system in the two cases. In each case the system adjusts by shifting to a new stationary state over a relaxation time which is estimated analytically. The model exhibits flare-like event statistics. In each case the frequency – energy distribution is a power law for flare energies less than a time-dependent rollover set by the largest energy the system is likely to attain at a given time. The rollover is not observed if the mean free energy of the system is sufficiently large. For Case 1, the model exhibits a double exponential waiting-time distribution, corresponding to flaring at a constant mean rate during two intervals (before and after the step change), if the average energy of the system is large. For Case 2 the waiting-time distribution is a simple exponential, again provided the average energy of the system is large. Monte Carlo simulations of Case 1 are presented which confirm the estimate for the relaxation time and the expected forms of the frequency – energy and waiting-time distributions. The simulation results provide a qualitative model for observed flare statistics in AR 11029.     

Seismic Emission from A M9.5-Class Solar Flare

Following the discovery of a few significant seismic sources at 6.0 mHz from the large solar flares of October 28 and 29, 2003, we have extended SOHO/MDI helioseismic observations to moderate M-class flares. We report the detection of seismic waves emitted from the β γ δ active region NOAA 9608 on September 9, 2001. A quite impulsive solar flare of type M9.5 occurred from 20:40 to 20:48 UT. We used helioseismic holography to image seismic emission from this flare into the solar interior and computed time series of egression power maps in 2.0 mHz bands centered at 3.0 and 6.0 mHz. The 6.0 mHz images show an acoustic source associated with the flare some 30 Mm across in the East – West direction and 15 Mm in the North – South direction nestled in the southern penumbra of the main sunspot of AR 9608. This coincides closely with three white-light flare kernels that appear in the sunspot penumbra. The close spatial correspondence between white-light and acoustic emission adds considerable weight to the hypothesis that the acoustic emission is driven by heating of the lower photosphere. This is further supported by a rough hydromechanical model of an acoustic transient driven by sudden heating of the low photosphere. Where direct heating of the low photosphere by protons or high-energy electrons is unrealistic, the strong association between the acoustic source and co-spatial continuum emission can be regarded as evidence supporting the back-warming hypothesis, in which the low photosphere is heated by radiation from the overlying chromosphere. This is to say that a seismic source coincident with strong, sudden radiative emission in the visible continuum spectrum indicates a photosphere sufficiently heated so as to contribute significantly to the continuum emission observed.     

Comparisons of Supergranule Characteristics During the Solar Minima of Cycles 22/23 and 23/24

Supergranulation is a component of solar convection that manifests itself on the photosphere as a cellular network of around 35 Mm across, with a turnover lifetime of 1 – 2 days. It is strongly linked to the structure of the magnetic field. The horizontal, divergent flows within supergranule cells carry local field lines to the cell boundaries, while the rotational properties of supergranule upflows may contribute to the restoration of the poloidal field as part of the dynamo mechanism, which controls the solar cycle. The solar minimum at the transition from cycle 23 to 24 was notable for its low level of activity and its extended length. It is of interest to study whether the convective phenomena that influence the solar magnetic field during this time differed in character from periods of previous minima. This study investigates three characteristics (velocity components, sizes and lifetimes) of solar supergranulation. Comparisons of these characteristics are made between the minima of cycles 22/23 and 23/24 using MDI Doppler data from 1996 and 2008, respectively. It is found that whereas the lifetimes are equal during both epochs (around 18 h), the sizes are larger in 1996 (35.9 ± 0.3 Mm) than in 2008 (35.0 ± 0.3 Mm), while the dominant horizontal velocity flows are weaker (139 ± 1 m s⁻¹ in 1996; 141 ± 1 m s⁻¹ in 2008). Although numerical differences are seen, they are not conclusive proof of the most recent minimum being inherently unusual.     

Energy-Dependent Timing of Thermal Emission in Solar Flares

We report solar flare plasma to be multi-thermal in nature based on the theoretical model and study of the energy-dependent timing of thermal emission in ten M-class flares. We employ high-resolution X-ray spectra observed by the Si detector of the “Solar X-ray Spectrometer” (SOXS). The SOXS onboard the Indian GSAT-2 spacecraft was launched by the GSLV-D2 rocket on 8 May 2003. Firstly we model the spectral evolution of the X-ray line and continuum emission flux F(ε) from the flare by integrating a series of isothermal plasma flux. We find that the multi-temperature integrated flux F(ε) is a power-law function of ε with a spectral index (γ)≈−4.65. Next, based on spectral-temporal evolution of the flares we find that the emission in the energy range E=4 – 15 keV is dominated by temperatures of T=12 – 50 MK, while the multi-thermal power-law DEM index (δ) varies in the range of −4.4 and −5.7. The temporal evolution of the X-ray flux F(ε,t) assuming a multi-temperature plasma governed by thermal conduction cooling reveals that the temperature-dependent cooling time varies between 296 and 4640 s and the electron density (n _e) varies in the range of n _e=(1.77 – 29.3)×10¹⁰ cm⁻³. Employing temporal evolution technique in the current study as an alternative method for separating thermal from nonthermal components in the energy spectra, we measure the break-energy point, ranging between 14 and 21±1.0 keV.     

Evaluation of a Selected Case of the Minimum Dissipative Rate Method for Non-Force-Free Solar Magnetic Field Extrapolations

The minimum dissipative rate (MDR) method for deriving a coronal non-force-free magnetic field solution is partially evaluated. These magnetic field solutions employ a combination of three linear (constant-α) force-free-field solutions with one being a potential field (i.e., α=0). The particular case of the solutions where the other two α’s are of equal magnitude but of opposite sign is examined. This is motivated by studying the SOLIS (Synoptic Optical Long-term Investigation of the Sun (SOLIS), a National Solar Observatory facility) vector magnetograms of AR 10987, which show a global α value consistent with an α=0 value as evaluated by (∇×B) _z /B _z over the region. Typical of the current state of the observing technology, there is no definitive twist for input into the general MDR method. This suggests that the special α case, of two α’s with equal magnitudes and opposite signs, is appropriate given the data. Only for an extensively twisted active region does a dominant, nonzero α normally emerge from a distribution of local values. For a special set of conditions, is it found that (i) the resulting magnetic field is a vertically inflated magnetic field resulting from the electric currents being parallel to the photosphere, similar to the results of Gary and Alexander (Solar Phys. 186:123, 1999), and (ii) for α≈(α _max/2), the Lorentz force per unit volume normalized by the square of the magnetic field is on the order of 1.4×10⁻¹⁰ cm⁻¹. The Lorentz force (F _L) is a factor of ten higher than that of the magnetic force d(B ²/8π)/dz, a component of F _L. The calculated photospheric electric current densities are an order of magnitude smaller than the maximum observed in all active regions. Hence both the Lorentz force density and the generated electric current density seem to be physically consistent with possible solar dynamics. The results imply that the field could be inflated with an overpressure along the neutral line. However, the implementation of this or any other extrapolation method using the electric current density as a lower boundary condition must be done cautiously, with the current magnetography.     

Granular-Scale Elementary Flux Emergence Episodes in a Solar Active Region

We analyse data from Hinode spacecraft taken over two 54-minute periods during the emergence of AR 11024. We focus on small-scale portions within the observed solar active region and discover the appearance of very distinctive small-scale and short-lived dark features in Ca ii H chromospheric filtergrams and Stokes I images. The features appear in regions with close-to-zero longitudinal magnetic field, and are observed to increase in length before they eventually disappear. Energy release in the low chromospheric line is detected while the dark features are fading. Three complete series of these events are detected with remarkably similar properties, i.e. lifetime of ≈ 12 min, maximum length and area of 2 – 4 Mm and 1.6 – 4 Mm², respectively, and all with associated brightenings. In time series of magnetograms a diverging bipolar configuration is observed accompanying the appearance of the dark features and the brightenings. The observed phenomena are explained as evidencing elementary flux emergence in the solar atmosphere, i.e. small-scale arch filament systems rising up from the photosphere to the lower chromosphere with a length scale of a few solar granules. Brightenings are explained as being the signatures of chromospheric heating triggered by reconnection of the rising loops (once they have reached chromospheric heights) with pre-existing magnetic fields, as well as being due to reconnection/cancellation events in U-loop segments of emerging serpentine fields. The characteristic length scale, area and lifetime of these elementary flux emergence events agree well with those of the serpentine field observed in emerging active regions. We study the temporal evolution and dynamics of the events and compare them with the emergence of magnetic loops detected in quiet Sun regions and serpentine flux emergence signatures in active regions. The physical processes of the emergence of granular-scale magnetic loops seem to be the same in the quiet Sun and active regions. The difference is the reduced chromospheric emission in the quiet Sun attributed to the fact that loops are emerging in a region of lower ambient magnetic field density, making interactions and reconnection less likely to occur. Incorporating the novel features of granular-scale flux emergence presented in this study, we advance the scenario for serpentine flux emergence.     

Two-component photosphere models of a 2N/M2-class solar flare

Physical state of the photosphere during a 2N/M2 solar flare on July 18, 2000, was studied. We used Echelle Zeeman spectrograms obtained by V. G. Lozitsky in orthogonal circular polarizations with a solar spectrograph. Semiempirical photospheric models were constructed for three moments in time in the initial and main phases of the flare using the SIR code applied to Stokes I and V profiles of seven iron and chromium lines. The photospheric model of the flare contains two components: a magnetic-field component and nonmagnetic environment. The height distributions of the temperature, magnetic field, and line-of-sight velocity were derived. The temperature in the nonmagnetic component had a nonmonotonous run with height. The models include layers in the middle and upper photosphere in which temperature is enhanced relative to an unperturbed photosphere model. As the flare developed, the temperature in the lower layers was increasing by 500–800 K. The magnetic field increased by 0.05 T and 0.08–0.1 T in the lower and upper photosphere during the flare, respectively, with the vertical temperature gradient decreasing from 0.0012 to 0.0008 T/km. The model for the onset phase of the flare indicates that there were upflows and downflows of substance in the lower and upper photosphere, respectively. The flow velocities decreased appreciably in the main phase of the flare. The model parameters of the nonmagnetic environment were only slightly different from those of the unperturbed photosphere.     

Microwave Neutral Line Associated Source and a Current Sheet

Neutral Line associated Sources (NLSs) are quasi-stationary microwave sources projected onto vicinities of the neutral line of the photospheric magnetic field. NLSs are often precursors of powerful flares, but their nature is unclear. We endeavor to reveal the structure of an NLS and to analyze a physical connection between such a source with a site of energy release in the corona above NOAA 10488 (October/November 2003). Evolution of this AR includes emergence and collision of two bipolar magnetic structures, rise of the main magnetic separator, and the appearance of an NLS underneath. The NLS appears at a contact site of colliding sunspots, whose relative motion goes on, resulting in a large shear along a tangent. Then the nascent NLS becomes the main source of microwave fluctuations in the AR. The NLS emission at 17 GHz is dominated by either footpoints or the top of a loop-like structure, an NLS loop, which connects two colliding sunspots. During a considerable amount of time, the emission dominates over that footpoint of the NLS loop, where the magnetic field is stronger. At that time, the NLS resembles a usual sunspot–associated radio source, whose brightness center is displaced towards the periphery of a sunspot. Microwave emission of an X2.7 flare is mainly concentrated in an ascending flare loop, initially coinciding with the NLS loop. The top of this loop is located at the base of a non-uniform bar-like structure visible in soft X-rays and at 34 GHz at the flare onset. We reveal i) upward lengthening of this bar before the flare onset, ii) the motion of the top of an apparently ascending flare loop along the axis of this bar, and iii) a non-thermal microwave source, whose descent along the bar was associated with the launching of a coronal ejection. We connect the bar with a probable position of a nearly vertical diffusion region, a site of maximal energy release inside an extended pre-flare current sheet. The top of the NLS loop is located at the bottom of this region. A combination of the NLS loop and diffusion region constitutes the skeleton of a quasi-stationary microwave NLS.     

An Extreme Solar Event of 20 January 2005: Properties of the Flare and the Origin of Energetic Particles

The famous extreme solar and particle event of 20 January 2005 is analyzed from two perspectives. Firstly, using multi-spectral data, we study temporal, spectral, and spatial features of the main phase of the flare, when the strongest emissions from microwaves up to 200 MeV gamma-rays were observed. Secondly, we relate our results to a long-standing controversy on the origin of solar energetic particles (SEP) arriving at Earth, i.e., acceleration in flares, or shocks ahead of coronal mass ejections (CMEs). Our analysis shows that all electromagnetic emissions from microwaves up to 2.22 MeV line gamma-rays during the main flare phase originated within a compact structure located just above sunspot umbrae. In particular, a huge (≈ 10⁵ sfu) radio burst with a high frequency maximum at 30 GHz was observed, indicating the presence of a large number of energetic electrons in very strong magnetic fields. Thus, protons and electrons responsible for various flare emissions during its main phase were accelerated within the magnetic field of the active region. The leading, impulsive parts of the ground-level enhancement (GLE), and highest-energy gamma-rays identified with π ⁰-decay emission, are similar and closely correspond in time. The origin of the π ⁰-decay gamma-rays is argued to be the same as that of lower-energy emissions, although this is not proven. On the other hand, we estimate the sky-plane speed of the CME to be 2 000 – 2 600 km s⁻¹, i.e., high, but of the same order as preceding non-GLE-related CMEs from the same active region. Hence, the flare itself rather than the CME appears to determine the extreme nature of this event. We therefore conclude that the acceleration, at least, to sub-relativistic energies, of electrons and protons, responsible for both the major flare emissions and the leading spike of SEP/GLE by 07 UT, are likely to have occurred nearly simultaneously within the flare region. However, our analysis does not rule out a probable contribution from particles accelerated in the CME-driven shock for the leading GLE spike, which seemed to dominate at later stages of the SEP event. S.N. Kuznetsov deceased 17 May 2007.     

Chromospheric Height and Density Measurements in a Solar Flare Observed with RHESSI – II. Data Analysis

We present an analysis of hard X-ray imaging observations from one of the first solar flares observed with the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) spacecraft, launched on 5 February 2002. The data were obtained from the 22 February 2002, 11:06 UT flare, which occurred close to the northwest limb. Thanks to the high energy resolution of the germanium-cooled hard X-ray detectors on RHESSI we can measure the flare source positions with a high accuracy as a function of energy. Using a forward-fitting algorithm for image reconstruction, we find a systematic decrease in the altitudes of the source centroids z(ε) as a function of increasing hard X-ray energy ε, as expected in the thick-target bremsstrahlung model of Brown. The altitude of hard X-ray emission as a function of photon energy ε can be characterized by a power-law function in the ε=15–50 keV energy range, viz., z(ε)≈2.3(ε/20 keV)^−1.3 Mm. Based on a purely collisional 1-D thick-target model, this height dependence can be inverted into a chromospheric density model n(z), as derived in Paper I, which follows the power-law function n _e(z)=1.25×10¹³(z/1 Mm)^−2.5 cm⁻³. This density is comparable with models based on optical/UV spectrometry in the chromospheric height range of h≲1000 km, suggesting that the collisional thick-target model is a reasonable first approximation to hard X-ray footpoint sources. At h≈1000–2500 km, the hard X-ray based density model, however, is more consistent with the `spicular extended-chromosphere model' inferred from radio sub-mm observations, than with standard models based on hydrostatic equilibrium. At coronal heights, h≈2.5–12.4 Mm, the average flare loop density inferred from RHESSI is comparable with values from hydrodynamic simulations of flare chromospheric evaporation, soft X-ray, and radio-based measurements, but below the upper limits set by filling-factor insensitive iron line pairs.     

Global Forces in Eruptive Solar Flares: The Lorentz Force Acting on the Solar Atmosphere and the Solar Interior

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 (Astron. Soc. Pac. CS-383, 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 (toward 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.     

On 17 – 22 January 2005 Events in Space Weather

This contribution is a follow-up to the recent paper of Kuznetsov et al. (Contrib. Astron. Obs. Skalnaté Pleso 36, 85, 2006) on the ground level enhancement (GLE) on 20 January 2005. We focused on a study of Forbush decrease (FD) of 17 – 18 and 21 – 22 January 2005, respectively. The data from the neutron monitor at Lomnicky Štít (1 min counts) and from the Geomagnetic Observatory in Hurbanovo, both in Slovakia, were used as the basis for our investigation. The data on magnetic field and solar wind from GOES 10 and 12, SOHO-CELIAS, ACE and WIND satellites were used for better understanding of the global evolution of the event. The magnetic field is transformed to the RTN (Radial – Tangential – Normal) system where only the disturbed part of the field is compared, i.e., daily variations and a constant part are subtracted. The field reduction method is described. Our results are temporal vector diagrams of variation of all parameters at all positions from where we used the data. The amplitudes of |B| exceed 100 nT and variations during the arrival of the wavefront of CME take place simultaneously at the ground-based station and at GOES satellites. The character of the variations is as if there would be regions with the dominant electric charge of opposite signs, or electric currents with different orientations in the CME. On the basis of the values v _p and n _p and using certain assumptions we determined the mass of CME on 17 January and 21 January, respectively, of 10¹² kg. A decrease of the cosmic ray level runs suddenly (during 10 minutes), starting, however, about two hours after a sudden change of the magnetic field.     

Do EUV Nanoflares Account for Coronal Heating?

Recent observations with EUV imaging instruments such as SOHO/EIT and TRACE have shown evidence for flare-like processes at the bottom end of the energy scale, in the range of E _th≈10²⁴–10²⁷ erg. Here we compare these EUV nanoflares with soft X-ray microflares and hard X-ray flares across the entire energy range. From the observations we establish empirical scaling laws for the flare loop length, L(T)∼T, the electron density, n _e(T)∼T ², from which we derive scaling laws for the loop pressure, p(T)∼T ³, and the thermal energy, E _th∼T ⁶. Extrapolating these scaling laws into the picoflare regime we find that the pressure conditions in the chromosphere constrain a height level for flare loop footpoints, which scales with h _eq(T)∼T ^−0.5. Based on this chromospheric pressure limit we predict a lower cutoff of flare loop sizes at L _∖min≲5 Mm and flare energies E _∖min≲10²⁴ erg. We show evidence for such a rollover in the flare energy size distribution from recent TRACE EUV data. Based on this energy cutoff imposed by the chromospheric boundary condition we find that the energy content of the heated plasma observed in EUV, SXR, and HXR flares is insufficient (by 2–3 orders of magnitude) to account for coronal heating.     

Evaluating Mean Magnetic Field in Flare Loops

We analyze multiple-wavelength observations of a two-ribbon flare exhibiting apparent expansion motion of the flare ribbons in the lower atmosphere and rising motion of X-ray emission at the top of newly-formed flare loops. We evaluate magnetic reconnection rate in terms of V _r B _r by measuring the ribbon-expansion velocity (V _r) and the chromospheric magnetic field (B _r) swept by the ribbons. We also measure the velocity (V _t) of the apparent rising motion of the loop-top X-ray source, and estimate the mean magnetic field (B _t) at the top of newly-formed flare loops using the relation 〈V _t B _t〉≈〈V _r B _r〉, namely, conservation of reconnection flux along flare loops. For this flare, B _t is found to be 120 and 60 G, respectively, during two emission peaks five minutes apart in the impulsive phase. An estimate of the magnetic field in flare loops is also achieved by analyzing the microwave and hard X-ray spectral observations, yielding B=250 and 120 G at the two emission peaks, respectively. The measured B from the microwave spectrum is an appropriately-weighted value of magnetic field from the loop top to the loop leg. Therefore, the two methods to evaluate coronal magnetic field in flaring loops produce fully-consistent results in this event.     

Oscillatory Response of the 3D Solar Atmosphere to the Leakage of Photospheric Motion

The direct propagation of acoustic waves, driven harmonically at the solar photosphere, into the three-dimensional solar atmosphere is examined numerically in the framework of ideal magnetohydrodynamics. It is of particular interest to study the leakage of 5-minute global solar acoustic oscillations into the upper, gravitationally stratified and magnetised atmosphere, where the modelled solar atmosphere possesses realistic temperature and density stratification. This work aims to complement and bring further into the 3D domain our previous efforts (by Erdélyi et al., 2007, Astron. Astrophys. 467, 1299) on the leakage of photospheric motions and running magnetic-field-aligned waves excited by these global oscillations. The constructed model atmosphere, most suitable perhaps for quiet Sun regions, is a VAL IIIC derivative in which a uniform magnetic field is embedded. The response of the atmosphere to a range of periodic velocity drivers is numerically investigated in the hydrodynamic and magnetohydrodynamic approximations. Among others the following results are discussed in detail: i) High-frequency waves are shown to propagate from the lower atmosphere across the transition region, experiencing relatively low reflection, and transmitting most of their energy into the corona; ii) the thin transition region becomes a wave guide for horizontally propagating surface waves for a wide range of driver periods, and particularly at those periods that support chromospheric standing waves; iii) the magnetic field acts as a waveguide for both high- and low-frequency waves originating from the photosphere and propagating through the transition region into the solar corona. Electronic Supplementary Material  The online version of this article () contains supplementary material, which is available to authorized users.     

Solar Intranetwork Magnetic Elements: Bipolar Flux Appearance

This study aims to quantify characteristic features of the bipolar flux appearance of solar intranetwork (IN) magnetic elements. To attack this problem, we use the Narrowband Filter Imager (NFI) magnetograms from the Solar Optical Telescope (SOT) on board Hinode; these data are from quiet and enhanced network areas. Cluster emergence of mixed polarities and IN ephemeral regions (ERs) are the most conspicuous forms of bipolar flux appearance within the network. Each of the clusters is characterized by a few well-developed ERs that are partially or fully coaligned in magnetic axis orientation. On average, the sampled IN ERs have a total maximum unsigned flux of several 10¹⁷ Mx, a separation of 3 – 4 arcsec, and a lifetime of 10 – 15 minutes. The smallest IN ERs have a maximum unsigned flux of several 10¹⁶ Mx, separations of less than 1 arcsec, and lifetimes as short as 5 minutes. Most IN ERs exhibit a rotation of their magnetic axis of more than 10 degrees during flux emergence. Peculiar flux appearance, e.g., bipole shrinkage followed by growth or the reverse, is not unusual. A few examples show repeated shrinkage–growth or growth–shrinkage, like magnetic floats in the dynamic photosphere. The observed bipolar behavior seems to carry rich information on magnetoconvection in the subphotospheric layer.     

Magnetic Field Structures in a Facular Region Observed by THEMIS and <Emphasis Type="Italic">Hinode</Emphasis>

We study the space – time distributions of intensity fluctuations in 2 – 3 hour sequences of multi-spectral, high-resolution, high-cadence, broad-band filtergram images of the Sun made by the SOT – FG system aboard the Hinode spacecraft. In the frequency range 5.5<f<8.0 mHz both G-band and Ca ii H-line oscillations are suppressed in the presence of magnetic fields, but the suppression disappears for f>10 mHz. By looking at G-band frequencies above 10 mHz we find that the oscillatory power, both at these frequencies and at lower frequencies, lies in a mesh pattern with cell scale 2 – 3 Mm, clearly larger than normal granulation, and with correlation times on the order of hours. The mesh pattern lies in the dark lanes between stable cells found in time-integrated G-band intensity images. It also underlies part of the bright pattern in time-integrated H-line emission. This discovery may reflect dynamical constraints on the sizes of rising granular convection cells together with the turbulence created in strong intercellular downflows.