Topological changes of the photospheric magnetic field inside active regions: A prelude to flares?

The detection of magnetic field variations as a signature of flaring activity is one of the main goals in solar physics. Past efforts gave apparently no unambiguous observations of systematic changes. In the present study, we discuss recent results from observations that scaling laws of turbulent current helicity inside a given flaring active region change in response to large flares in that active region. Such changes can be related to the evolution of current structures by a simple geometrical argument, which has been tested using high Reynolds number direct numerical simulations of the MHD equations. Interpretation of the observed data within this picture indicates that the change in scaling behavior of the current helicity seems to be associated with a topological reorganization of the footpoint of the magnetic field loops, namely with the dissipation of small scales structures in turbulent media.     

Transient Magnetic and Doppler Features Related to the White-Light Flares in NOAA 10486

Rapidly moving transient features have been detected in magnetic and Doppler images of super-active region NOAA 10486 during the X17/4B flare of 28 October 2003 and the X10/2B flare of 29 October 2003. Both these flares were extremely energetic white-light events. The transient features appeared during impulsive phases of the flares and moved with speeds ranging from 30 to 50 km?s^?1. These features were located near the previously reported compact acoustic (Donea and Lindsey, Astrophys. J. 630, 1168, 2005) and seismic sources (Zharkova and Zharkov, Astrophys. J. 664, 573, 2007). We examine the origin of these features and their relationship with various aspects of the flares, viz., hard X-ray emission sources and flare kernels observed at different layers: i) photosphere (white-light continuum), ii) chromosphere (Hα 6563 Å), iii) temperature minimum region (UV 1600 Å), and iv) transition region (UV 284 Å).     

Calculating and Testing Nonlinear Force-Free Fields

Improvements to an existing method for calculating nonlinear force-free magnetic fields (Wheatland, Solar Phys. 238, 29, 2006) are described. In particular a solution of the 3-D Poisson equation using 2-D Fourier transforms is presented. The improved nonlinear force-free method is demonstrated in application to linear force-free test cases with localized nonzero values of the normal component of the field in the boundary. These fields provide suitable test cases for nonlinear force-free calculations because the boundary conditions involve localized nonzero values of the normal components of the field and of the current density, and because (being linear force-free fields) they have more direct numerical solutions. Despite their simplicity, fields of this kind have not been recognized as test cases for nonlinear methods before. The examples illustrate the treatment of the boundary conditions on current in the nonlinear force-free method, and in particular the limitations imposed by field lines that connect outside of the boundary region.     

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.     

Determination of the True Shape of Coronal Loops

Using line of sight velocity measurements from the SUMER and CDS instruments aboard SOHO, in conjunction with a simple geometrical model, we reconstructed the true, 3D shape and the velocity of plasma flow along coronal loops. The projection of the loop on the sky and the position of the footpoints define a family of curves. Assuming that the loop is located on a plane, the line of sight velocity can be used to select the most plausible solution. For two loops, observed in the Ne viii? 770 Å and O v? 630 Å spectral lines, we find asymmetric, subsonic uni-directional flows, with velocity maxima of ≈?80 km?s^?1 near the footpoints. The loops are highly inclined with respect to the vertical, by 55^○ and 70^○, respectively; thus the true height of the loop tops from the photospheric level is ≈?20^′′, comparable to the isothermal scale height.     

Two-Dimensional Spectroscopy of Photospheric Shear Flows in a Small <Emphasis Type="Italic">δ</Emphasis> Spot

In recent high-resolution observations of complex active regions, long-lasting and well-defined regions of strong flows were identified in major flares and associated with bright kernels of visible, near-infrared, and X-ray radiation. These flows, which occurred in the proximity of the magnetic neutral line, significantly contributed to the generation of magnetic shear. Signatures of these shear flows are strongly curved penumbral filaments, which are almost tangential to sunspot umbrae rather than exhibiting the typical radial filamentary structure. Solar active region NOAA 10756 was a moderately complex β δ sunspot group, which provided an opportunity to extend previous studies of such shear flows to quieter settings. We conclude that shear flows are a common phenomenon in complex active regions and δ spots. However, they are not necessarily a prerequisite condition for flaring. Indeed, in the present observations, the photospheric shear flows along the magnetic neutral line are not related to any change of the local magnetic shear. We present high-resolution observations of NOAA 10756 obtained with the 65-cm vacuum reflector at Big Bear Solar Observatory (BBSO). Time series of speckle-reconstructed white-light images and two-dimensional spectroscopic data were combined to study the temporal evolution of the three-dimensional vector flow field in the β δ sunspot group. An hour-long data set of consistent high quality was obtained, which had a cadence of better than 30 seconds and subarcsecond spatial resolution.     

Low-Energy Cutoffs in Electron Spectra of Solar Flares: Statistical Survey

The Reuven Ramaty High Energy Spectroscopic Imager (RHESSI) X-ray data base (February 2002 – May 2006) has been searched to find solar flares with weak thermal components and flat photon spectra. Using a regularized inversion technique, we determine the mean electron flux distribution from count spectra for a selection of events with flat photon spectra in the 15 – 20 keV energy range. Such spectral behavior is expected for photon spectra either affected by photospheric albedo or produced by electron spectra with an absence of electrons in a given energy range (e.g., a low-energy cutoff in the mean electron spectra of nonthemal particles). We have found 18 cases that exhibit a statistically significant local minimum (a dip) in the range of 13 – 19 keV. The positions and spectral indices of events with low-energy cutoff indicate that such features are likely to be the result of photospheric albedo. It is shown that if the isotropic albedo correction is applied, all low-energy cutoffs in the mean electron spectrum are removed, and hence the low-energy cutoffs in the mean electron spectrum of solar flares above ∼ 12 keV cannot be viewed as real features. If low-energy cutoffs exist in the mean electron spectra, their energies should be less than ∼ 12 keV.     

RHESSI Observations of the Neupert Effect in Three Solar Flares

Previous observations show that in many solar flares there is a causal correlation between the hard X-ray flux and the derivative of the soft X-ray flux. This so-called Neupert effect is indicative of a strong link between the primary energy release to accelerate particles and plasma heating. It suggests a flare model in which the hard X-rays are electron – ion bremsstrahlung produced by energetic electrons as they lose their energy in the lower corona and chromosphere and the soft X-rays are thermal bremsstrahlung from the “chromospheric evaporation” plasma heated by those same electrons. Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observes in a broad energy band and its high spectral resolution and coverage of the low-energy range allow us to separate the thermal continuum from the nonthermal component, which gives us an opportunity to investigate the Neupert effect. In this paper, we use the parameters derived from RHESSI observations to trace the primary energy release and the plasma response: The hard X-ray flux or spectral hardness is compared with the derivative of plasma thermal energy in three impulsive flares on 10 November 2002 and on 3 and 25 August 2005. High correlations show that the Neupert effect does hold for the two hard X-ray peaks of the 10 November 2002 flare, for the first peaks of the 3 August 2005 flare, and for the beginning period of the 25 August 2005 flare.     

Frequency Distributions of Microwave Pulses for the 18 March 2003 Solar Flare

We analyze the pulses in high-frequency drift radio structures observed by the spectrometer at Purple Mountain Observatory (PMO) over the frequency range of 4.5 – 7.5 GHz during the 18 March 2003 solar flare. A number of individual pulses are determined from the drifting radio structures after the detected gradual component subtraction. The frequency distributions of microwave pulse occurrence as functions of peak flux, duration, bandwidth, and time interval between two adjacent pulses exhibit a power-law behavior, i.e. . From regression fitting in log-log space, we obtain the power-law indexes, α _P=7.38±0.40 for the peak flux, α _D=5.39±0.86 for the duration, and α _B=6.35±0.56 for the bandwidth. We find that the frequency distribution for the time interval displays a broken power law. The break occurs at about 500 ms, and their indexes are α _W1=1.56±0.08 and α _W2=3.19±0.12, respectively. Our results are consistent with the previous findings of hard X-ray pulses, type III bursts, and decimetric millisecond spikes.     

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.     

A Comparison of Feature Classification Methods for Modeling Solar Irradiance Variation

Physical understanding of total and spectral solar irradiance variation depends upon establishing a connection between the temporal variability of spatially resolved solar structures and spacecraft observations of irradiance. One difficulty in comparing models derived from different data sets is that the many ways for identifying solar features such as faculae, sunspots, quiet Sun, and various types of “network” are not necessarily consistent. To learn more about classification differences and how they affect irradiance models, feature “masks” are compared as derived from five current methods: multidimensional histogram analysis of NASA/National Solar Observatory/Kitt Peak spectromagnetograph data, statistical pattern recognition applied to SOHO/Michelson Doppler Imager photograms and magnetograms, threshold masks allowing for influence of spatial surroundings applied to NSO magnetograms, and “one-trigger” and “three-trigger” algorithms applied to California State University at Northridge Cartesian Full Disk Telescope intensity observations. In general all of the methods point to the same areas of the Sun for labeling sunspots and active-region faculae, and available time series of area measurements from the methods correlate well with each other and with solar irradiance. However, some methods include larger label sets, and there are important differences in detail, with measurements of sunspot area differing by as much as a factor of two. The methods differ substantially regarding inclusion of fine spatial scale in the feature definitions. The implications of these differences for modeling solar irradiance variation are discussed. K.L. Harvey and S.R. Walton are deseased, to whom this paper is dedicated.     

原子钟噪声的蒙特卡络模拟方法

讨论了原子钟噪声的蒙特卡络模拟方法．具体过程是：先产生模拟所需的正态分布随机数，从实际钟的阿仑方差估算出噪声电平，然后用１组递点函数产生高质量的模拟钟时差的随机序列．在对１台商品铯钟进行具体模拟后，给出了真实数据与模拟值的阿仑方差之问的相符程度，最后对实际问题进行了分析．     

Monte Carlo computations of the escape of atomic nitrogen from Mars

Monte Carlo simulations were carried out to compute the escape flux of atomic nitrogen for the low and high solar activity martian thermospheres. The total escape of atomic nitrogen at low and high solar activities was found to be 3.03×¹⁰⁵ and , respectively. The escape flux of atomic nitrogen at low and high solar activities from photodissociation of N₂ was found to be 2.75×¹⁰⁵ and , respectively. The remainder of the contribution is from dissociative recombination, which is only important at high solar activity were it comprises about 25% of the total escape. The relative contributions to the total N escape flux from thermal motion of the background atmosphere, winds and co-rotation, and photoionization and subsequent solar wind pickup are also considered here. We find that the total predicted escape fluxes are observed to increase by 20 and 25% at low and high solar activities owing to thermal motion of the background atmosphere. At low and high solar activities, we find that the co-rotation and wind velocities combined translate to a maximum transferable energy of ∼0.0103 and 0.0181 eV, respectively, and that the total escape flux contribution from winds and co-rotation is negligible. Photoionization was found to be a minor process only impacting those source atoms produced with energies close to the escape energy, between 1.5 and 2 eV. The contributions to the total escape fluxes at low and high solar activities from photoionization and subsequent solar wind pickup are found to be about 8 and 13%, respectively.     

Monte Carlo simulations of the water vapor plumes on Enceladus

Monte Carlo simulations are used to model the July 14, 2005 UVIS stellar occultation observations of the water vapor plumes on Enceladus. These simulations indicate that the observations can be best fit if the water molecules ejected along the Tiger Stripes in the South Polar region of Enceladus have a vertical surface velocity of 300-500 m/s at the surface. The high surface velocity suggests that the plumes on Enceladus originate from some depth beneath the surface. The total escape rate of water molecules is 4-6×¹⁰²⁷ s⁻¹, or 120-180 kg/s, consistent with previous works, and more than 100 times the estimated mass escape rate for ice particles. The average deposition rate in the South Polar region is on the order of 10¹¹ cm⁻² s⁻¹, yielding a resurfacing rate as high as 3×¹⁰⁻⁴ cm/yr. The globally averaged deposition rate of water molecules is about one order of magnitude lower.     

Adiabatic Deceleration of Solar Energetic Particles as Deduced from Monte Carlo Simulations of Interplanetary Transport

Monte Carlo simulations of interplanetary transport are employed to study adiabatic energy losses of solar protons during propagation in the interplanetary medium. We consider four models. The first model is based on the diffusion-convection equation. Three other models employ the focused transport approach. In the focused transport models, we simulate elastic scattering in the local solar wind frame and magnetic focusing. We adopt three methods to treat scattering. In two models, we simulate a pitch-angle diffusion as successive isotropic or anisotropic small-angle scatterings. The third model treats large-angle scatterings as numerous small-chance isotropizations. The deduced intensity–time profiles are compared with each other, with Monte Carlo solutions to the diffusion-convection equation, and with results of the finite-difference scheme by Ruffolo (1995). A numerical agreement of our Monte Carlo simulations with results of the finite-difference scheme is good. For the period shortly after the maximum intensity time, including deceleration can increase the decay rate of the near-Earth intensity essentially more than would be expected based on advection from higher momenta. We, however, find that the excess in the exponential-decay rate is time dependent. Being averaged over a reasonably long period, the decay rate of the near-Earth intensity turns out to be close to that expected based on diffusion, convection, and advection from higher momenta. We highlight a variance of the near-Earth energy which is not small in comparison with the energy lost. It leads to blurring of any fine details in the accelerated particle spectra. We study the impact of realistic spatial dependencies of the mean free path on adiabatic deceleration and on the near-Earth intensity magnitude. We find that this impact is essential whenever adiabatic deceleration itself is important. It is also found that the initial angular distribution of particles near the Sun can markedly affect MeV-proton energy losses and intensities observed at 1 AU. Computations invoked during the study are described in detail.     

RHESSI耀斑中非热电子低能截止的研究

太阳耀斑中硬X射线(HXR)光子谱的低能变平过去一般认为是由于耀斑中非热电子的低能截止造成的,但现在也有作者认为耀斑光子与下层大气的逆康普顿散射(albedo效应)或者其他作用也能够使得HXR光子谱出现低能变平的情形.采用Gan etal.(2001,2002)中提出的求非热电子低能截止的方法,统计分析了Ramaty High EnergySolar Spectroscopy Imager(RHESSI)卫星在2002--2005年间观测的100个耀斑,发现经albedo校正,有18个耀斑的HXR光子谱可以利用单幂律谱来拟合,在80个可以用双幂律谱来拟合HXR光子谱的耀斑中,有21个耀斑可以直接用单幂律电子谱加一个低能截止来解释.低能截止范围为20-50keV,平均值约为30keV.同时也分析了耀斑光子谱特征的其他可能解释.     

银河系球状星团轨道参数分布的Monte Carlo模拟

选取前文中所列出的29个累积光谱型为F型的球状星团中的3个作为样本，深入研究了初始观测资料的不确定性和选用不同的银河系引力势模型，对样本星团轨道参数的影响。首先采用Monte Carlo方法产生3个样本球状星团的模拟初始观测数据，而后，以这些模拟数据为初始条件，在3种不同的银河系引力势模型下进行轨道计算，得到此3个样本的模拟轨道参数。模拟计算的结果表明：根据模拟初始数据生成的样本轨道参数分布形态大致可分为高斯分布、准高斯分布和非高斯分布等3类；初始观测数据的不确定性对样本轨道参数分布的影响，与样本星团的选择和轨道参数的类型有关；选用不同的银河系引力势模型，对3个样本星团的各个轨道参数的分布和形态结构也会产生不同程度的影响。该工作的结果，可供深入研究球状星团的整体运动和动力学性质等问题参考。