Oscillatory Modes of a Prominence – PCTR – Corona Slab Model

Oscillations of magnetic structures in the solar corona have often been interpreted in terms of magnetohydrodynamic waves. We study the adiabatic magnetoacoustic modes of a prominence plasma slab with a uniform longitudinal magnetic field, surrounded by a prominence – corona transition region (PCTR) and a coronal medium. Considering linear small-amplitude oscillations, we deduce the dispersion relation for the magnetoacoustic slow and fast modes by assuming evanescentlike perturbations in the coronal medium. In the system without PCTR, a classification of the oscillatory modes according to the polarisation of their eigenfunctions is made to distinguish modes with fastlike or slowlike properties. Internal and external slow modes are governed by the prominence and coronal properties, respectively, and fast modes are mostly dominated by prominence conditions for the observed wavelengths. In addition, the inclusion of an isothermal PCTR does not substantially influence the mode frequencies, but new solutions (PCTR slow modes) are present.     

Coronal Mass Ejections,Magnetic Fields,and the Green Corona in Cycle 23

We study a time – latitudinal distribution of CMEs observed by the SOHO spacecraft, their projected speeds and associated magnetic fields, as well as the north – south (N – S) asymmetry of solar surface magnetic fields, and the coronal green line intensities. We have found that (a) there exists an intricate relation between the average projected velocity of CMEs and the mean value of large-scale magnetic fields; (b) there exists a pronounced N – S asymmetry in both the distribution and the number of CMEs; (c) this asymmetry is in favor of the northern hemisphere at the beginning of the cycle, and of the southern hemisphere from 2001 onward, being, in fact, (d) closely related with the N – S asymmetry in the distribution of large-scale magnetic fields and the coronal green line intensities.     

Damped Oscillations of Coronal Loops

A mechanism of damped oscillations of a coronal loop is investigated. The loop is treated as a thin toroidal flux rope with two stationary photospheric footpoints, carrying both toroidal and poloidal currents. The forces and the flux-rope dynamics are described within the framework of ideal magnetohydrodynamics (MHD). The main features of the theory are the following: i) Oscillatory motions are determined by the Lorentz force that acts on curved current-carrying plasma structures and ii) damping is caused by drag that provides the momentum coupling between the flux rope and the ambient coronal plasma. The oscillation is restricted to the vertical plane of the flux rope. The initial equilibrium flux rope is set into oscillation by a pulse of upflow of the ambient plasma. The theory is applied to two events of oscillating loops observed by the Transition Region and Coronal Explorer (TRACE). It is shown that the Lorentz force and drag with a reasonable value of the coupling coefficient (c _d ) and without anomalous dissipation are able to accurately account for the observed damped oscillations. The analysis shows that the variations in the observed intensity can be explained by the minor radial expansion and contraction. For the two events, the values of the drag coefficient consistent with the observed damping times are in the range c _d ≈2 – 5, with specific values being dependent on parameters such as the loop density, ambient magnetic field, and the loop geometry. This range is consistent with a previous MHD simulation study and with values used to reproduce the observed trajectories of coronal mass ejections (CMEs).     

Nonlinear excitation of small-scale Alfvén waves by fast waves and plasma heating in the solar atmosphere

We study a nonlinear mechanism for the excitation of kinetic Alfvén waves (KAWs) by fast magneto-acoustic waves (FWs) in the solar atmosphere. Our focus is on the excitation of KAWs that have very small wavelengths in the direction perpendicular to the background magnetic field. Because of their small perpendicular length scales, these waves are very efficient in the energy exchange with plasmas and other waves. We show that the nonlinear coupling of the energy of the finite-amplitude FWs to the small-scale KAWs can be much faster than other dissipation mechanisms for fast wave, such as electron viscous damping, Landau damping, and modulational instability. The nonlinear damping of the FWs due to decay FW = KAW + KAW places a limit on the amplitude of the magnetic field in the fast waves in the solar corona and solar-wind at the level B/B ₀∼10⁻². In turn, the nonlinearly excited small-scale KAWs undergo strong dissipation due to resistive or Landau damping and can provide coronal and solar-wind heating. The transient coronal heating observed by Yohkoh and SOHO may be produced by the kinetic Alfvén waves that are excited by parametric decay of fast waves propagating from the reconnection sites.     

MHD waves in coronal loops with a shell

We consider a model of a coronal loop in the form of a cord surrounded by a coaxial shell. Two slow magnetosonic waves longitudinally propagate within a thin flux tube on the m=0 cylindrical mode with velocities close to the tube velocities in the cord and the shell. One wave propagates inside the cord, while the other propagates inside the shell. A peculiar feature of the second wave is that the plasma in the cord and the shell oscillates with opposite phases. There are two fast magnetosonic waves on each of the cylindrical modes with m>0. If the plasma density in the shell is lower than that in the surrounding corona, then one of the waves is radiated into the corona, which causes the loop oscillations to be damped, while the other wave is trapped by the cord, but can also be radiated out under certain conditions. If the plasma density in the shell is higher than that in the cord, then one of the waves is trapped by the shell, while the other wave can also be trapped by the shell under certain conditions. In the wave trapped by the shell and the wave radiated by the tube, the plasma in the cord and the shell oscillates with opposite phases.     

Space Weather Application Using Projected Velocity Asymmetry of Halo CMEs

Halo coronal mass ejections (HCMEs) originating from regions close to the center of the Sun are likely to be responsible for severe geomagnetic storms. It is important to predict geoeffectiveness of HCMEs by using observations when they are still near the Sun. Unfortunately, coronagraphic observations do not provide true speeds of CMEs because of projection effects. In the present paper, we present a new technique to allow estimates of the space speed and approximate source location using projected speeds measured at different position angles for a given HCME (velocity asymmetry). We apply this technique to HCMEs observed during 2001 – 2002 and find that the improved speeds are better correlated with the travel times of HCMEs to Earth and with the magnitudes of ensuing geomagnetic storms.     

Sunspot Group Decay

We examine daily records of sunspot group areas (measured in millionths of a solar hemisphere or μHem) for the last 130 years to determine the rate of decay of sunspot group areas. We exclude observations of groups when they are more than 60° in longitude from the central meridian and only include data when at least three days of observations are available following the date of maximum area for a group’s disk passage. This leaves data for over 18 000 measurements of sunspot group decay. We find that the decay rate increases linearly from 28 μHem day⁻¹ to about 140 μHem day⁻¹ for groups with areas increasing from 35 μHem to 1000 μHem. The decay rate tends to level off for groups with areas larger than 1000 μHem. This behavior is very similar to the increase in the number of sunspots per group as the area of the group increases. Calculating the decay rate per individual sunspot gives a decay rate of about 3.65 μHem day⁻¹ with little dependence upon the area of the group. This suggests that sunspots decay by a Fickian diffusion process with a diffusion coefficient of about 10 km² s⁻¹. Although the 18 000 decay rate measurements are lognormally distributed, this can be attributed to the lognormal distribution of sunspot group areas and the linear relationship between area and decay rate for the vast majority of groups. We find weak evidence for variations in decay rates from one solar cycle to another and for different phases of each sunspot cycle. However, the strongest evidence for variations is with latitude and the variations with cycle and phase of each cycle can be attributed to this variation. High latitude spots tend to decay faster than low latitude spots.     

Analysis of Two Coronal Loops with Combined TRACE and SOHO/CDS Data

We use an innovative research technique to analyze combined images from the Coronal Diagnostic Spectrometer (CDS) on the Solar and Heliospheric Observatory (SOHO) and the Transition Region and Coronal Explorer (TRACE). We produce a high spatial and temporal resolution simulated CDS raster or “composite” map from TRACE data and use this composite map to jointly analyze data from both instruments. We show some of the advantages of using the “composite” map method for coronal loop studies. We investigate two postflare loop structures. We find cool material (250 000 K) concentrated at the tips or apex of the loops. This material is found to be above its scale height and therefore not in hydrostatic equilibrium. The exposure times of the composite map and TRACE images are used to give an estimate of another loop’s cooling time. The contribution to the emission in the TRACE images for the spectral lines present in its narrow passband is estimated by using the CDS spectral data and CHIANTI to derive synthetic spectra. We obtain cospatial and cotemporal data collected by both instruments in SOHO Joint Observations Program (JOP) 146 and show how the combination of these data can be utilized to obtain more accurate measurements of coronal plasmas than if analyzed individually. Electronic Supplementary Material  The online version of this article () contains supplementary material, which is available to authorized users.     

Dissipation of Standing Slow Magnetoacoustic Waves in Hot Coronal Loops

The damping of standing slow waves in hot (T>6 MK) coronal loops of semicircular shape is revisited in both the linear and nonlinear regimes. Dissipation by thermal conduction, compressive viscosity, radiative cooling, and heating are examined for nonstratified and stratified loops. We find that for typical conditions of hot SUMER loops, thermal conduction increases the period of damped oscillations over the sound-crossing time, whereas the decay times are mostly shaped by compressive viscosity. Damping from optically thin radiation is negligible. We also find that thermal conduction alone results in slower damping of the density and velocity waves compared to the observations. Only when compressive viscosity is added do these waves damp out at the same rate as the observed rapidly decaying modes of hot SUMER loop oscillations, in contrast to most current work, which has pointed to thermal conduction as the dominant mechanism. We compare the linear predictions with numerical hydrodynamic calculations. Under the effects of gravity, nonlinear viscous dissipation leads to a reduction of the decay time compared to the homogeneous case. In contrast, the linear results predict that the damping rates are barely affected by gravity.     

The Effect of Abnormal Granulation on Acoustic Wave Travel Times and Mode Frequencies

Observations indicate that in plage areas (i.e. in active regions outside sunspots) acoustic waves travel faster than in the quiet Sun, leading to shortened travel times and higher p-mode frequencies. Coupled with the 11-year variation of solar activity, this may also explain the solar cycle variation of oscillation frequencies. While it is clear that the ultimate cause of any difference between the quiet Sun and plage is the presence of magnetic fields of order 100 G in the latter, the mechanism by which the magnetic field exerts its influence has not yet been conclusively identified. One possible such mechanism is suggested by the observation that granular motions in plage areas tend to be slightly “abnormal”, dampened compared to the quiet Sun. In this paper we consider the effect that abnormal granulation observed in active regions should have on the propagation of acoustic waves. Any such effect is found to be limited to a shallow surface layer where sound waves propagate nearly vertically. The magnetically suppressed turbulence implies higher sound speeds, leading to shorter travel times. This time shift Δ τ is independent of the travel distance, while it shows a characteristic dependence on the assumed plage field strength. As a consequence of the variation of the acoustic cutoff with height, Δ τ is expected to be significantly higher for higher frequency waves within the observed regime of 3 – 5 mHz. The lower group velocity near the upper reflection point further leads to an increased envelope time shift, as compared to the phase shift. p-mode frequencies in plage areas are increased by a corresponding amount, Δ ν/ν=ν Δ τ. These characteristics of the time and frequency shifts are in accordance with observations. The calculated overall amplitudes of the time and frequency shifts are comparable to, but still significantly less than (by a factor of 2 to 5), those suggested by measurements.     

Intensity and Magnetic Field Distribution of Sunspots

We study the relationship between the brightness (I) and magnetic field (B) distributions of sunspots using 272 samples observed at the San Fernando Observatory and the National Solar Observatory, Kitt Peak, whose characteristics varied widely. We find that the I – B relationship has a quadratic form for the spots with magnetic field less than about 2000 G. The slope of the linear part of the I – B curve varies by about a factor of three for different types of spots. In general the slope increases as the spot approaches disk center. The I – B slope does not have a clear dependency on the spot size but the lower limit appears to increase as a function of the ratio of umbra and penumbra area. The I – B slope changes as a function of age of the sunspots. We discuss various sunspot models using these results.     

A Study of Surges: II. On the Relationship between Chromospheric Surges and Coronal Mass Ejections

Liu et al. (Astrophys. J. 628, 1056, 2005a) described one surge – coronal mass ejection (CME) event showing a close relationship between solar chromospheric surge ejection and CME that had not been noted before. In this work, large Hα surges (>72 Mm, or 100 arcsec) are studied. Eight of these were associated with CMEs. According to their distinct morphological features, Hα surges can be classified into three types: jetlike, diffuse, and closed loop. It was found that all of the jetlike surges were associated with jetlike CMEs (with angular widths ≤30 degrees); the diffuse surges were all associated with wide-angle CMEs (e.g., halo); the closed-loop surges were not associated with CMEs. The exclusive relation between Hα surges and CMEs indicates difference in magnetic field configurations. The jetlike surges and related narrow CMEs propagate along coronal fields that are originally open. The unusual transverse mass motions in the diffuse surges are suggested to be due to magnetic reconnections in the corona that produce wide-angle CMEs. For the closed-loop surges, their paths are just outlining stable closed loops close to the solar surface. Thus no CMEs are associated with them.     

MHD waves at a spherical interface modelling EIT waves

Eruptive events such as flares and coronal mass ejections (CMEs) are known to generate global waves propagating over distances comparable to the solar radius in different layers of the solar atmosphere. Here we investigate the propagation of coronal EIT waves, modelled as fast magnetoacoustic modes propagating at a spherical interface in the presence of a purely radial magnetic field. Based on a simplified equilibrium we derive the dispersion relation of the waves. The generation and propagation of EIT waves at the spherical interface is studied numerically for different values of spherical degree and preliminary conclusions are reached regarding the properties of EIT waves. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Steps for Building a Calibrated Flux-Transport Dynamo for the Sun

Given the complexity involved in a flux-transport-type dynamo driven by both Babcock – Leighton and tachocline α effects, we present here a step-by-step procedure for building a flux-transport dynamo model calibrated to the Sun as a guide for anyone who wishes to build this kind of model. We show that a plausible sequence of steps to reach a converged solution in such a dynamo consists of (i) numerical integration of a classical α – ω dynamo driven by a tachocline α effect, (ii) continued integration with inclusion of meridional circulation to convert the model into a flux-transport dynamo driven by only a tachocline α effect, (iii) final integration with inclusion of a Babcock – Leighton surface α effect, resulting in a flux-transport dynamo that can be calibrated to obtain a close fit of model output with solar observations.     

Dynamical processes in flux tabes and their role in chromospheric heating

We model the dynamical interaction between magnetic flux tubes and granules in the solar photosphere which leads to the excitation of transverse (kink) and longitudinal (sausage) tube waves. The investigation is motivated by the interpretation of network oscillations in terms of flux tube waves. The calculations show that for magnetic field strengths typical of the network, the energy flux in transverse waves is higher than in longitudinal waves by an order of magnitude. But for weaker fields, such as those that might be found in internetwork regions, the energy fluxes in the two modes are comparable. Using observations of footpoint motions, the energy flux in transverse waves is calculated and the implications for chromospheric heating are pointed out.     

Solar Cycle Variation of the Interplanetary Forward Shock Drivers Observed at 1 AU

Forecasting space weather more accurately from solar observations requires an understanding of the variations in physical properties of interplanetary (IP) shocks as solar activity changes. We examined the characteristics (occurrence rate, physical parameters, and types of shock driver) of IP shocks. During the period of 1995 – 2001, a total of 249 forward IP shocks were observed. In calculating the shock parameters, we used the solar wind data from Wind at the solar minimum period (1995 – 1997) and from ACE since 1998 including the solar maximum period (1999 – 2001). Most of IP shocks (68%) are concentrated in the solar maximum period. The values of physical quantities of IP shocks, such as the shock speed, the sonic Mach number, and the ratio of plasma density compression, are larger at solar maximum than at solar minimum. However, the ratio of IMF compression is larger at solar minimum. The IP shock drivers are classified into four groups: magnetic clouds (MCs), ejecta, high speed streams (HSSs), and unidentified drivers. The MC is the most dominant and strong shock driver and 150 out of total 249 IP shocks are driven by MCs. The MC is a principal and very effective shock driver not only at solar maximum but also at solar minimum, in contrast to results from previous studies, where the HSS is considered as the dominant IP shock driver.     

RHESSI Microflares: II. Implications for Loop Structure and Evolution

We present simple analytic models which predict the peak X-ray emission measure and temperature attained in flares in which the chromospheric evaporation process takes place either in a single ‘monolithic’ loop or in a loop consisting of several filaments that are created successively as the energy release process proceeds in time. As possible mechanisms driving chromospheric evaporation we consider both classical heat conduction from the loop top and non-thermal electron beams. The model predictions are tested for a set of 18 well studied RHESSI microflares. The results suggest beam driven evaporation in filamented loops as being capable of accounting for the observed emission measures and temperatures though there are issues with the very high beam densities needed. On the other hand, estimates of the emission measures achieved by conductive evaporation which are derived by using the Rosner – Tucker – Vaiana (RTV) scaling law are much larger than the observed ones. Possible reasons for this discrepancy are discussed.     

Positively Drifting Structures During the 18 March 2003 Solar Flare

We analyze the 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. The drifting structures take place before the soft X-ray maximum, almost at the maximum of hard X-ray flux at 25 – 50 keV. For the first time, the positive drift in this kind of radio structures is detected in such a high frequency range. Their global drifting rate is roughly estimated as 3.6 GHz s⁻¹. They appear in four groups, lasting in total for less than 6 s, and have a broad bandwidth of more than 2 GHz but a smaller ratio of the bandwidth of the drifting structures to mean frequency than that of the lower frequency range. The lifetime of each individual burst in this event can be derived by using the high temporal resolution of the spectrometer at PMO and has an average value of 36.3 ms. Since the negative drifting structures observed in the 0.6 – 4.5 GHz frequency range were interpreted to be a radio signature of a plasmoid ejected upward (moving out of the Sun), the present observation may imply that it is possible for a plasmoid to move downward during a solar flare. However, for a confirmation of this suggestion direct radio imaging observation would be needed.     

High-frequency modes in solar-like stars

p-mode oscillations in solar-like stars are excited by the outer convection zone in these stars and reflected close to the surface. The p modes are trapped inside an acoustic cavity, but the modes only stay trapped up to a given frequency [known as the acoustic cut-off frequency  (ν_ac)  ] as modes with larger frequencies are generally not reflected at the surface. This means that modes with frequency larger than the acoustic cut-off frequency must be travelling waves. The high-frequency modes may provide information about the physics in the outer layers of the stars and the excitation source and are therefore highly interesting as it is the estimation of these two phenomena that cause some of the largest uncertainties when calculating stellar oscillations.
High-frequency modes have been detected in the Sun, in β Hydri and in α Cen A and α Cen B by smoothing the so-called echelle diagram and the large frequency separation as a function of frequency has been estimated. The large frequency separation has been compared with a simple model of the acoustic cavity which suggests that the reflectivity of the photosphere is larger at high frequency than predicted by standard models of the solar atmosphere and that the depth of the excitation source is larger than what has been estimated by other models and might depend on the order n and degree l of the modes.     

The Butterfly Diagram in the Eighteenth Century

Digitized images of the drawings by J.C. Staudacher were used to determine sunspot positions for the period 1749 – 1796. From the entire set of drawings, 6285 sunspot positions were obtained for a total of 999 days. Various methods have been applied to find the orientation of the solar disk, which is not given for the vast majority of the drawings by Staudacher. Heliographic latitudes and longitudes in the Carrington rotation frame were determined. The resulting butterfly diagram shows a highly-populated Equator during the first two cycles (cycles 0 and 1 in the usual counting since 1749). An intermediate period is cycle 2, whereas cycles 3 and 4 show a typical butterfly shape. A tentative explanation may be the transient dominance of a quadrupolar magnetic field during the first two cycles.