Modeling the (upper) solar atmosphere including the magnetic field

The atmosphere of the Sun is highly structured and dynamic in nature. From the photosphere and chromosphere into the transition region and the corona plasma-β changes from above to below one, i.e., while in the lower atmosphere the energy density of the plasma dominates, in the upper atmosphere the magnetic field plays the governing role – one might speak of a “magnetic transition”. Therefore the dynamics of the overshooting convection in the photosphere, the granulation, is shuffling the magnetic field around in the photosphere. This leads not only to a (re-)structuring of the magnetic field in the upper atmosphere, but induces also the dynamic reaction of the coronal plasma, e.g., due to reconnection events. Therefore the (complex) structure and the interaction of various magnetic patches is crucial to understand the structure, dynamics and heating of coronal plasma as well as its acceleration into the solar wind.

The present article will emphasize the need for three-dimensional modeling accounting for the complexity of the solar atmosphere to understand these processes. Some advances on 3D modeling of the upper solar atmosphere in magnetically closed as well as open regions will be presented together with diagnostic tools to compare these models to observations. This highlights the recent success of these models which in many respects closely match the observations.     

Asymmetric precipitation in a coronal loop as explanation of a singular observed spectrum

Almost 10 years of solar submillimeter observations have shown new aspects of solar activity, such as the presence of rapid solar spikes associated with the launch of coronal mass ejections and an increasing submillimeter spectral component in flares. We analyse the singular microwave–submillimeter spectrum of an M class solar flare on 20 December, 2002. Flux density observations measured by Sun patrol telescopes and the Solar Submillimeter Telescope are used to build the radio spectrum, which is fitted using Ramaty’s code. At submillimeter frequencies the spectrum shows a component different from the microwave classical burst. The fitting is achieved proposing two homogeneous sources of emission. This theoretical fitting is in agreement with differential precipitation through a magnetically asymmetric loop or set of loops. From a coronal magnetic field model we infer an asymmetric magnetic structure at the flare location. The model proposed to quantify the differential precipitation rates due to the asymmetry results in a total precipitation ratio Q₂/Q₁≈10⁴–10⁵, where Q₁(Q₂) represents the total precipitation in the loop foot with the high (low) magnetic field intensity. This ratio agrees with the electron total number ratio of the two sources proposed to fit the radio spectrum.     

Quasi-biennial oscillations in the cross-correlation of properties of macrospicules

Jets, whatever small (e.g. spicules) or large (e.g. macrospicules) their size, may play a key role in momentum and energy transport from photosphere to chromosphere and at least to the low corona. Here, we investigate the properties of abundant, large-scale dynamic jets observable in the solar atmosphere: the macrospicules (MS). These jets are observationally more distinct phenomena than their little, and perhaps more ubiquitous, cousins, the spicules. Investigation of long-term variation of the properties of macrospicules may help to a better understanding of their underlying physics of generation and role in coronal heating. Taking advantage of the high temporal and spatial resolution of the Solar Dynamics Observatory, a new dataset, with several hundreds of macrospicules, was constructed encompassing a period of observations over six years. Here, we analyse the measured properties and relations between these properties of macrospicules as function of time during the observed time interval. We found that cross-correlations of several of these macrospicule properties display a strong oscillatory pattern. Next, wavelet analysis is used to provide more detailed information about the temporal behaviour of the various properties of MS. For coronal hole macrospicules, a significant peak is found at around 2-year period. This peak also exists partially or is shifted to longer period, in the case of quiet Sun macrospicules. These observed findings may be rooted in the underlying mechanism generating the solar magnetic field, i.e. the global solar dynamo.     

Three-dimensional global simulation of multiple ICMEs’ interaction and propagation from the Sun to the heliosphere following the 25–28 October 2003 solar events

This study performs simulations of interplanetary coronal mass ejection (ICME) propagation in a realistic three-dimensional (3D) solar wind structure from the Sun to the Earth by using the newly developed hybrid code, HAFv.2+3DMHD. This model combines two simulation codes, Hakamada–Akasofu–Fry code version 2 (HAFv.2) and a fully 3D, time-dependent MHD simulation code. The solar wind structure is simulated out to 0.08 AU (18 Rs) from source surface maps using the HAFv.2 code. The outputs at 0.08 AU are then used to provide inputs for the lower boundary, at that location, of the 3D MHD code to calculate solar wind and its evolution to 1 AU and beyond. A dynamic disturbance, mimicking a particular flare’s energy output, is delivered to this non-uniform structure to model the evolution and interplanetary propagation of ICMEs (including their shocks). We then show the interaction between two ICMEs and the dynamic process during the overtaking of one shock by the other. The results show that both CMEs and heliosphere current sheet/plasma sheet were deformed by interacting with each other.     

Source investigation of impulsive He-rich particle events

We have investigated the source characteristic and coronal magnetic field structure of six impulsive solar energetic particle (SEP) events selected from Wang et al. [Wang, Y.-M., Pick, M., Mason, G.M. Coronal holes, jets, and the origin of ³He-rich particle events. ApJ 639, 495, 2006] and Pick et al. [Pick, M., Mason, G.M., Wang, Y.-M., Tan, C., Wang, L. Solar source regions for ³He-rich solar energetic particle events identified using imaging radio, optical, and energetic particle observations. ApJ 648, 1247, 2006]. Some results are obtained: first, 2 events are associated with wide (≈100°) CMEs (hereafter wide CME events), another 4 events are associated with narrow (?40°) CMEs (hereafter narrow CME events); second, the coronal magnetic field configuration of narrow CME events appear more simple than that of the wide CME events; third, the photospheric magnetic field evolutions of all these events show new emergence of fluxes, while one case also shows magnetic flux cancellation; fourth, the EUV jets usually occurred very close to the footpoint of the magnetic field loop, while meter type III bursts occurred near or at the top of the loop and higher than EUV jets. Furthermore, the heights of type III bursts are estimated from the result of the coronal magnetic field extrapolations.     

From solar nanoflares to stellar giant flares: Scaling laws and non-implications for coronal heating

In this study we explore physical scaling laws applied to solar nanoflares, microflares, and large flares, as well as to stellar giant flares. Solar flare phenomena exhibit a fractal volume scaling, V(L)  L^1.9, with L being the flare loop length scale, which explains the observed correlation between the total emission measure EM_p and flare peak temperature T_p in both solar and stellar flares. However, the detected stellar flares have higher emission measures EM_p than solar flares at the same flare peak temperature T_p, which can be explained by a higher electron density that is caused by shorter heating scale height ratios s_H/L ≈ 0.04–0.1. Using these scaling laws we calculate the total radiated flare energies E_X and thermal flare energies E_T and find that the total counts C are a good proxy for both parameters. Comparing the energies of solar and stellar flares we find that even the smallest observed stellar flares exceed the largest solar flares, and thus their observed frequency distributions are hypothetically affected by an upper cutoff caused by the maximum active region size limit. The powerlaw slopes fitted near the upper cutoff can then not reliably be extrapolated to the microflare regime to evaluate their contribution to coronal heating.     

Magnetic coupling of the Sun–Earth system – The view from STEREO

The STEREO mission, launched on October 25 2006, will provide the first stereoscopic view of the Sun and its magnetic environment. A suite of identical instruments on two continuously separating spacecraft will monitor the onset of solar eruptive phenomena, and track them as the shocks and ejected material propagate through the interplanetary medium (IPM). The combination of remote sensing and in situ instrumentation will provide new insights into the onset of eruptions, the extent of their effects on the global structure of the low corona, and their subsequent evolution through the IPM. These observations will provide new constraints on the processes involved and allow us to distinguish between competing models of eruptive solar phenomena.     

Probing the first stars and black holes in the early Universe with the Dark Ages Radio Explorer (DARE)

A concept for a new space-based cosmology mission called the Dark Ages Radio Explorer (DARE) is presented in this paper. DARE’s science objectives include: (1) When did the first stars form? (2) When did the first accreting black holes form? (3) When did Reionization begin? (4) What surprises does the end of the Dark Ages hold (e.g., Dark Matter decay)? DARE will use the highly-redshifted hyperfine 21-cm transition from neutral hydrogen to track the formation of the first luminous objects by their impact on the intergalactic medium during the end of the Dark Ages and during Cosmic Dawn (redshifts z = 11–35). It will measure the sky-averaged spin temperature of neutral hydrogen at the unexplored epoch 80–420 million years after the Big Bang, providing the first evidence of the earliest stars and galaxies to illuminate the cosmos and testing our models of galaxy formation. DARE’s approach is to measure the expected spectral features in the sky-averaged, redshifted 21-cm signal over a radio bandpass of 40–120 MHz. DARE orbits the Moon for a mission lifetime of 3 years and takes data above the lunar farside, the only location in the inner solar system proven to be free of human-generated radio frequency interference and any significant ionosphere. The science instrument is composed of a low frequency radiometer, including electrically-short, tapered, bi-conical dipole antennas, a receiver, and a digital spectrometer. The smooth frequency response of the antennas and the differential spectral calibration approach using a Markov Chain Monte Carlo technique will be applied to detect the weak cosmic 21-cm signal in the presence of the intense solar system and Galactic foreground emissions.     

Variations of the solar acoustic high-degree mode frequencies over solar cycle 23

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present variations of the solar acoustic mode frequencies caused by the solar activity cycle. High-degree (100 < ? < 900) solar acoustic modes were analyzed using global helioseismology analysis techniques over most of solar cycle 23. We followed the methodology described in details in [Korzennik, S.G., Rabello-Soares, M.C., Schou, J. On the determination of Michelson Doppler Imager high-degree mode frequencies. ApJ 602, 481–515, 2004] to infer unbiased estimates of high-degree mode parameters ([see also Rabello-Soares, M.C., Korzennik, S.G., Schou, J. High-degree mode frequencies: changes with solar cycle. ESA SP-624, 2006]). We have removed most of the known instrumental and observational effects that affect specifically high-degree modes. We show that the high-degree changes are in good agreement with the medium-degree results, except for years when the instrument was highly defocused. We analyzed and discuss the effect of defocusing on high-degree estimation. Our results for high-degree modes confirm that the frequency shift scaled by the relative mode inertia is a function of frequency and it is independent of degree.     

Flaring loop parameters estimated from solar decimeter type U-like and type J-like fine structures

This work presents the analysis of five fine structures in the solar radio emission, observed between June 2000 and October 2001 by the Brazilian Solar Spectroscope (BSS), in the decimeter frequency band of 950–2500 MHz. Based on their morphological characteristics identified in the dynamic spectra, the fine structures had been classified as type U-like or type J-like bursts. Such emissions are variants of the type III bursts. They support the hypothesis of generation by plasma emission mechanism, from interaction of electron beams accelerated during solar flares, propagating along closed magnetic structures, within the trapped plasma of the solar corona. The spectral and temporal characteristics of the five fine structures had been obtained from the dynamic spectra and the parameters of the agent and the emitting source have been determined, assuming both fundamental and harmonic emissions. The analysis revealed the flux density of the structures is less than 20–80 s.f.u. For assumption of harmonic emission, the interval of values for the source parameters estimated are: the loop size is (0.3–5.1) × 10¹⁰ cm; the electron beam velocity is in the range of 0.16–0.53 c; the temperature of coronal loop top is of the order of (0.25–1.55) × 10⁷ K; and the low limit for the magnetic field is of 7–26 G. These results are in agreement with previous determinations reported in the literature.     

The dynamic chromosphere: Pushing the boundaries of observations and models

The interface between the bright solar surface and the million-degree corona continues to hold the key to many unsolved problems in solar physics. Advances in instrumentation now allow us to observe the dynamic structures of the solar chromosphere down to less than $0_{.}^{″}1$ with cadences of just a few seconds and in multiple polarisation states. Such observational progress has been matched by the ever-increasing sophistication of numerical models, which have become necessary to interpret the complex observations. With an emphasis on the quiet Sun, I will review recent progress in the observation and modelling of the chromosphere. Models have come a long way from 1D static atmospheres, but their predictions still fail to reproduce several key observed features. Nevertheless, they have given us invaluable insight into the physical processes that energise the atmosphere. With more physics being added to models, the gap between predictions and observations is narrowing. With the next generation of solar observatories just around the corner, the big question is: will they close the gap?     

第21,22周峰年太阳射电爆发中毫秒级精细结构的分析

本文对第21周的2840MHz及22周峰年的2545MHz，2645MHz和2840MHz波段上所观测到的285个快速事件进行分析，发现两个峰年的快速事件具有多种多样的快速活动（快速精细结构）．它们与耀斑的对应关系也做了统计分析．     

A MHD-turbulence model for solar corona

The disposition of energy in the solar corona has always been a problem of great interest. It remains an open question how the low temperature photosphere supports the occurence of solar extreme phenomena. In this work, a turbulent heating mechanism for the solar corona through the framework of reduced magnetohydrodynamics (RMHD) is proposed. Two-dimensional incompressible long time simulations of the average energy disposition have been carried out with the aim to reveal the characteristics of the long time statistical behavior of a two-dimensional cross-section of a coronal loop and the importance of the photospheric time scales in the understanding of the underlying mechanisms. It was found that for a slow, shear type photospheric driving the magnetic field in the loop self-organizes at large scales via an inverse MHD cascade. The system undergoes three distinct evolutionary phases. The initial forcing conditions are quickly “forgotten” giving way to an inverse cascade accompanied with and ending up to electric current dissipation. Scaling laws are being proposed in order to quantify the nonlinearity of the system response which seems to become more impulsive for decreasing resistivity. It is also shown that few, if any, qualitative changes in the above results occur by increasing spatial resolution.     

Using LASCO-C1 spectroscopy for coronal diagnostics

The LASCO-C1 telescope was designed to perform spectral analysis of coronal structures by means of a tunable Fabry–Pérot interferometer acquiring images at different wavelengths. Results from spectral scans of the Fe xiv 5303 Å green coronal emission line are presented. Physical quantities like the ion temperature (line widths), and the flow velocity along the line of sight (Doppler shifts) are obtained over the entire corona.     

Statistical ionospheric E layer properties measured with the Cyprus Digisonde and comparisons with IRI predictions

This paper reports the diurnal, seasonal, and long term variability of the E layer critical frequency (foE) and peak height (hmE) derived from Digisonde measurements from 2009 to 2016 at the low-middle latitude European station of Nicosia, Cyprus (geographical coordinates: 35°N, 33°E, geomagnetic lat. 29.38°N, I = 51.7°). Manually scaled monthly median values of foE and hmE are compared with IRI-2012 predictions with a view to assess the predictability of IRI. Results show that in general, IRI slightly overestimates foE values both at low and high solar activity. At low solar activity, overestimations are mostly limited to 0.25?MHz (equivalent electron density, 0.775?×?10³?el/m^?3) but can go as high as 0.5?MHz (equivalent electron density, 3.1?×?10³?el/m^?3, during noon) around equinox. In some months, underestimations, though sporadic in nature, up to 0.25?MHz are noted (mostly during sunrise and sunset). At high solar activity, a similar pattern of over-/underestimation is evident. During the entire period of study, over-/under estimations are mostly limited to 0.25?MHz. In very few cases, these exceed 0.25?MHz but are limited to 0.5?MHz. Analysis of hmE reveals that: (1) hmE remains almost constant during ±2 to ±4?h around local noon, (2) hmE values are higher in winter than in spring, summer and autumn, (3) there are two maxima near sunrise and sunset with a noontime minimum in between. During the entire period of study, significant differences between observed hmE and the IRI predictions have been noted. IRI fails to predict hmE and outputs a constant value of 110?km, which is higher than most of the observed values. Over- and under estimations range from 3 to 13?km and from 0 to 3?km respectively.     

A statistical study of CME-associated flare during the solar cycle 24

We investigate on the relationship between flares and coronal mass ejections (CMEs) in which a flare started before and after the CME events which differ in their physical properties, indicating potentially different initiation mechanisms. The physical properties of two types flare-correlated CME remain an interesting and important question in space weather. We study the relationship between flares and CMEs using a different approach requiring both temporal and spatial constraints during the period from December 1, 2008 to April 30, 2017 in which the CMEs data were acquired by SOHO/LASCO (Solar and Heliospheric Observatory/Large Angle Spectrometric Coronagraph) over the solar cycle 24. The soft X-ray flare flux data, such as flare class, location, onset time and integrated flux, are collected from Geostationary Environmental satellite (GOES) and XRT Flare catalogs. We selected 307 CMEs-flares pairs applying simultaneously temporal and spatial constraints in all events for the distinguish between two associated CME-flare types. We study the correlated properties of coincident flares and CMEs during this period, specifically separating the sample into two types: flares that precede a CME and flares that follow a CME. We found an opposite correlation relationship between the acceleration and velocity of CMEs in the After- and Before-CMEs events. We found a log-log relation between the width and mass of CMEs in the two associated types. The CMEs and flares properties show that there were significant differences in all physical parameters such as (mass, angular width, kinetic energy, speed and acceleration) between two flare-associated CME types.     

Observational study of magnetic chirality of solar active regions

We present the evolution of magnetic field and relationship with the magnetic (current) helicity in solar active regions from a series of photospheric vector magnetograms obtained at Huairou Solar Observing Station near Beijing, and also longitudinal magnetograms by MDI of SOHO, white light and 171 Å images by TRACE and soft X-ray images by Yohkoh.The conclusions in the analysis of the formation process of complex and delta magnetic configuration in some super active regions are the following: (1) The magnetic shear and gradient provide the non-potentiality of the magnetic field of active regions reflecting the existence of electric current. (2) Some of large-scale delta active regions could be due to the emergence of highly sheared non-potential magnetic flux bundles from the subatmosphere with amount of magnetic helicity, in addition to the emergence of twisted magnetic ropes. (3) We also present some results on the study of the magnetic (current) helicity in solar active regions.     

Synoptic magnetic field in cycle 23 and in the beginning of the cycle 24

The SOHO/MDI data provide the uniform time series of the synoptic magnetic maps which cover the period of the cycle 23 and the beginning of the cycle 24. It is very interesting period because of the long and deep solar minimum between the cycles 23 and 24. Synoptic structure of the solar magnetic field shows variability during solar cycles. It is known that the magnetic activity contributes to the solar irradiance. The axisymmetrical distribution of the magnetic flux (Fig. 3c) is closely associated with the ‘butterfly’ diagram in the EUV emission (Benevolenskaya et al., 2001). And, also, the magnetic field (B_∥) shows the non-uniform distributions of the solar activity with longitude, so-called ‘active zones’, and ‘coronal holes’ in the mid-latitude. Polar coronal holes are forming after the solar maxima and they persist during the solar minima. SOHO/EIT data in the emission of Fe XII (195 Å) could be a proxy for the coronal holes tracking. The active longitudinal zones or active longitude exist due to the reappearance of the activity and it is clearly seen in the synoptic structure of the solar cycle. On the descending branch of the solar cycle 23 active zones are less pronounced comparing with previous cycles 20, 21 and 22. Moreover, the weak polar magnetic field precedes the long and deep solar minimum. In this paper we have discussed the development of solar cycles 23 and 24 in details.     

Quiet sun XUV and EUV spectroscopy

The two XUV–EUV spectrometers on SOHO have collected a large amount of data in the 6000–10⁶ K solar plasma temperature range. These data have allowed us to greatly enhance our knowledge of the processes acting in the solar atmosphere, from the chromosphere to the corona. Some results on the quiet Sun structure (network, quiet Sun versus coronal hole), on the dynamics (velocities, waves, transient events), and the main characteristics of the quiet Sun atmosphere are presented and discussed.