The plasma of the lower solar atmosphere consists of mostly neutral particles, whereas the upper solar atmosphere is mostly ionised particles and electrons. A shock that propagates upwards in the solar atmosphere therefore undergoes a transition where the dominant fluid is either neutral or ionised. An upwards propagating shock also passes a point where the sound and Alfv\'en speed are equal. At this point the energy of the acoustic shock can separated into fast and slow components. How the energy is distributed between the two modes depends on the angle of magnetic field. The separation of neutral and ionised species in a gravitationally stratified atmosphere is investigated. The role of two-fluid effects on the structure of the shocks post-mode-conversion and the frictional heating is quantified for different levels of collisional coupling. Two-fluid numerical simulations are performed using the (P\underline{I}P) code of a wave steepening into a shock in an isothermal, partially-ionised atmosphere. The collisional coefficient is varied to investigate the regimes where the plasma and neutral species are weakly, strongly and finitely coupled. The propagation speeds of the compressional waves hosted by neutral and ionised species vary, therefore velocity drift between the two species is produced as the plasma attempts to propagate faster than the neutrals. This is most extreme for a fast-mode shock. We find that the collisional coefficient drastically changes the features present in the system, specifically the mode conversion height, type of shocks present, and the finite shock widths created by the two-fluid effects. In the finitely-coupled regime fast-mode shock widths can exceed the pressure scale height leading to a new potential observable of two-fluid effects in the lower solar atmosphere.

Recent observations show cool, oscillating prominence threads fading when observed in cool spectral lines and appearing in warm spectral lines. A proposed mechanism to explain this evolution is that the threads were heated by turbulence driven by the Kelvin--Helmholtz instability that developed as a result of wave-driven shear flows on the surface of the thread. As the Kelvin--Helmholtz instability is an instability that works to mix the fluids, in the solar corona it can be expected to work by mixing the cool prominence material with that of the hot corona to form a warm boundary layer. In this paper we develop a simple phenomenological model of nonlinear Kelvin--Helmholtz mixing, using it to determine the characteristic density and temperature of the mixing layer, which for the case under study with constant pressure across the two fluids are $\rho_{\rm mixed}=\sqrt{\rho_1\rho_2}$ and $T_{\rm mixed}=\sqrt{T_1T_2}$. One result from the model is that it provides an accurate, as determined by comparison with simulation results, determination of the kinetic energy in the mean velocity field. A consequence of this is that the magnitude of turbulence, and with it the energy that can be dissipated on fast time-scales, as driven by this instability can be determined. For the prominence-corona system, the mean temperature rise possible from turbulent heating is estimated to be less than 1\% of the characteristic temperature (which is found to be $10^5$\,K). These results highlight that mixing, and not heating, are likely to be the cause of the observed transition between cool to warm material in Okamoto et. al (2015). One consequence of this result is that the mixing creates a region with higher radiative loss rates on average than either of the original fluids, meaning that this instability could contribute a net loss of thermal energy from the corona, i.e. coronal cooling.

The nonlinear magnetic Kelvin-Helmholtz instability (KHi), and the turbulence it creates, appears in many astrophysical systems. This includes those systems where the local plasma conditions are such that the plasma is not fully ionised, for example in the lower solar atmosphere and molecular clouds. In a partially ionised system, the fluids couple via collisions which occur at characteristic frequencies, therefore neutral and plasma species become decoupled for sufficiently high-frequency dynamics. Here we present high-resolution 2D two-fluid simulations of the nonlinear KHi for a system that traverses the dynamic scales between decoupled fluids and coupled dynamics. We discover some interesting phenomena, including the presence of a density coupling that is independent of the velocity coupling. Using these simulations we analyse the heating rate, and two regimes appear. The first is a regime where the neutral flow is decoupled from the magnetic field that is characterised with a constant heating rate, then at larger scales the strong coupling approximation holds and the heating rate. At large scales with the KHi layer width to the $-2$ power. There is an energy cascade in the simulation, but the nature of the frictional heating means the heating rate is determined by the largest scale of the turbulent motions, a fact that has consequences for understanding turbulent dissipation in multi-fluid systems.

Slow-mode shocks are important in understanding fast magnetic reconnection, jet formation and heating in the solar atmosphere, and other astrophysical systems. The atmospheric conditions in the solar chromosphere allow both ionised and neutral particles to exist and interact. Under such conditions, fine substructures exist within slow-mode shocks due to the decoupling and recoupling of the plasma and neutral species. We study numerically the fine substructure within slow-mode shocks in a partially ionised plasma, in particular, analysing the formation of an intermediate transition within the slow-mode shock. High-resolution 1D numerical simulations are performed using the (P\underline{I}P) code using a two-fluid approach. We discover that long-lived intermediate (Alfven) shocks can form within the slow-mode shock, where there is a shock transition from above to below the Alfven speed and a reversal of the magnetic field across the shock front. The collisional coupling provides frictional heating to the neutral fluid, resulting in a Sedov-Taylor-like expansion with overshoots in the neutral velocity and neutral density. The increase in density results in a decrease of the Alfven speed and with this the plasma inflow is accelerated to above the Alfven speed within the finite width of the shock leading to the intermediate transition. This process occurs for a wide range of physical parameters and an intermediate shock is present for all investigated values of plasma-$\beta$, neutral fraction, and magnetic angle. As time advances the magnitude of the magnetic field reversal decreases since the neutral pressure cannot balance the Lorentz force. The intermediate shock is long-lived enough to be considered a physical structure, independent of the initial conditions.

Magnetic reconnection may be the fundamental process allowing energy stored in magnetic fields to be released abruptly, solar flares and coronal mass ejection (CME) being archetypal natural plasma examples. Magnetic reconnection is much too slow a process to be efficient on the large scales, but accelerates once small enough scales are formed in the system. For this reason, the fractal reconnection scenario was introduced (Shibata and Tanuma 2001) to explain explosive events in the solar atmosphere: it was based on the recursive triggering and collapse via tearing instability of a current sheet originally thinned during the rise of a filament in the solar corona. Here we compare the different fractal reconnection scenarios that have been proposed, and derive generalized scaling relations for the recursive triggering of fast, `ideal' - i.e. Lundquist number independent - tearing in collapsing current sheet configurations with arbitrary current profile shapes. An important result is that the Sweet-Parker scaling with Lundquist number, if interpreted as the aspect ratio of the singular layer in an ideally unstable sheet, is universal and does not depend on the details of the current profile in the sheet. Such a scaling however must not be interpreted in terms of stationary reconnection, rather it defines a step in the accelerating sequence of events of the ideal tearing mediated fractal cascade. We calculate scalings for the expected number of plasmoids for such generic profiles and realistic Lundquist numbers.

The Kelvin-Helmholtz instability has been proposed as a mechanism to extract energy from magnetohydrodynamic (MHD) kink waves in flux tubes, and to drive dissipation of this wave energy through turbulence. It is therefore a potentially important process in heating the solar corona. However, it is unclear how the instability is influenced by the oscillatory shear flow associated with an MHD wave. We investigate the linear stability of a discontinuous oscillatory shear flow in the presence of a horizontal magnetic field within a Cartesian framework that captures the essential features of MHD oscillations in flux tubes. We derive a Mathieu equation for the Lagrangian displacement of the interface and analyse its properties, identifying two different instabilities: a Kelvin-Helmholtz instability and a parametric instability involving resonance between the oscillatory shear flow and two surface Alfv\'{e}n waves. The latter occurs when the system is Kelvin-Helmholtz stable, thus favouring modes that vary along the flux tube, and as a consequence provides an important and additional mechanism to extract energy. When applied to flows with the characteristic properties of kink waves in the solar corona, both instabilities can grow, with the parametric instability capable of generating smaller scale disturbances along the magnetic field than possible via the Kelvin-Helmholtz instability. The characteristic time-scale for these instabilities is $\sim 100$ s, for wavelengths of $200$ km. The parametric instability is more likely to occur for smaller density contrasts and larger velocity shears, making its development more likely on coronal loops than on prominence threads.

Prominences are incredibly dynamic across the whole range of their observable spatial scales, with observations revealing gravity-driven fluid instabilities, waves, and turbulence. With all these complex motions, it would be expected that instabilities driven by shear in the internal fluid motions would develop. However, evidence of these have been lacking. Here we present the discovery in a prominence, using observations from the Interface Region Imaging Spectrograph (IRIS), of a shear flow instability, the Kelvin-Helmholtz sinusoidal-mode of a fluid channel, driven by flows in the prominence body. This finding presents a new mechanism through which we can create turbulent motions from the flows observed in quiescent prominences. The observation of this instability in a prominence highlights their great value as a laboratory for understanding the complex interplay between magnetic fields and fluid flows that play a crucial role in a vast range of astrophysical systems.

We analyze solar quiescent prominence bubble characteristics and instability dynamics using Hinode/Solar Optical Telescope (SOT) data. We measure bubble expansion rate, prominence downflows, and the profile of the boundary layer brightness and thickness as a function of time. The largest bubble analyzed rises into the prominence with a speed of about 1.3 km/s until it is destabilized by a localized shear flow on the boundary. Boundary layer thickness grows gradually as prominence downflows deposit plasma onto the bubble with characteristic speeds of 20 - 35 km/s. Lateral downflows initiate from the thickened boundary layer with characteristic speeds of 25 - 50 km/s, "draining" the layer of plasma. Strong shear flow across one bubble boundary leads to an apparent coupled Kelvin-Helmholtz Rayleigh-Taylor (KH-RT) instability. We measure shear flow speeds above the bubble of 10 km/s and infer interior bubble flow speeds on the order of 100 km/s. Comparing the measured growth rate of the instability to analytic expressions, we infer a magnetic flux density across the bubble boundary of ~10^{-3} T (10 gauss) at an angle of ~70 degrees to the prominence plane. The results are consistent with the hypothesis that prominence bubbles are caused by magnetic flux that emerges below a prominence, setting up the conditions for RT, or combined KH-RT, instability flows that transport flux, helicity, and hot plasma upward into the overlying coronal magnetic flux rope.

We investigate the onset of 2D time-dependent magnetic reconnection that is triggered using an external velocity driver located away from, and perpendicular to, an equilibrium Harris current sheet. Previous studies have typically utilised an internal trigger to initiate reconnection, e.g. initial conditions centred on the current sheet. Numerical simulations solving the compressible, resistive magnetohydrodynamics equations were performed to investigate the reconnection onset within different atmospheric layers of the Sun, namely the corona, chromosphere and photosphere. A reconnecting state is reached for all atmospheric heights considered, with the dominant physics being highly dependent on atmospheric conditions. The coronal case achieves a sharp rise in electric field for a range of velocity drivers. For the chromosphere, we find a larger velocity amplitude is required to trigger reconnection. For the photospheric environment, the electric field is highly dependent on the inflow speed; a sharp increase in electric field is obtained only as the velocity entering the reconnection region approaches the Alfven speed. Additionally, the role of ambipolar diffusion is investigated for the chromospheric case and we find that the ambipolar diffusion alters the structure of the current density in the inflow region. The rate at which flux enters the reconnection region is controlled by the inflow velocity. This determines all aspects of the reconnection start-up process, i.e. the early onset of reconnection is dominated by the advection term in Ohms law in all atmospheric layers. A lower plasma-$\beta$ enhances reconnection and creates a large change in the electric field. A high plasma-$\beta$ hinders the reconnection, yielding a sharp rise in the electric field only when the velocity flowing into the reconnection region approaches the local Alfven speed.

This work examines the effect of the embedded magnetic field strength on the non-linear development of the magnetic Rayleigh-Taylor Instability (RTI) (with a field-aligned interface) in an ideal gas close to the incompressible limit in three dimensions. Numerical experiments are conducted in a domain sufficiently large so as to allow the predicted critical modes to develop in a physically realistic manner. The ratio between gravity, which drives the instability in this case (as well as in several of the corresponding observations), and magnetic field strength is taken up to a ratio which accurately reflects that of observed astrophysical plasma, in order to allow comparison between the results of the simulations and the observational data which served as inspiration for this work. This study finds reduced non-linear growth of the rising bubbles of the RTI for stronger magnetic fields, and that this is directly due to the change in magnetic field strength, rather than the indirect effect of altering characteristic length scales with respect to domain size. By examining the growth of the falling spikes, the growth rate appears to be enhanced for the strongest magnetic field strengths, suggesting that rather than affecting the development of the system as a whole, increased magnetic field strengths in fact introduce an asymmetry to the system. Further investigation of this effect also revealed that the greater this asymmetry, the less efficiently the gravitational energy is released. By better understanding the under-studied regime of such a major phenomenon in astrophysics, deeper explanations for observations may be sought, and this work illustrates that the strength of magnetic fields in astrophysical plasmas influences observed RTI in subtle and complex ways.

In astrophysical systems with partially ionized plasma the motion of ions is governed by the magnetic field while the neutral particles can only feel the magnetic field's Lorentz force indirectly through collisions with ions. The drift in the velocity between ionized and neutral species plays a key role in modifying important physical processes like magnetic reconnection, damping of magnetohydrodynamic waves, transport of angular momentum in plasma through the magnetic field, and heating. This paper investigates the differences between Doppler velocities of calcium ions and neutral hydrogen in a solar prominence to look for velocity differences between the neutral and ionized species. We simultaneously observed spectra of a prominence over an active region in H I 397 nm, H I 434 nm, Ca II 397 nm, and Ca II 854 nm using a high dispersion spectrograph of the Domeless Solar Telescope at Hida observatory, and compared the Doppler velocities, derived from the shift of the peak of the spectral lines presumably emitted from optically-thin plasma. There are instances when the difference in velocities between neutral atoms and ions is significant, e.g. 1433 events (~ 3 % of sets of compared profiles) with a difference in velocity between neutral hydrogen atoms and calcium ions greater than 3sigma of the measurement error. However, we also found significant differences between the Doppler velocities of two spectral lines emitted from the same species, and the probability density functions of velocity difference between the same species is not significantly different from those between neutral atoms and ions. We interpreted the difference of Doppler velocities as a result of motions of different components in the prominence along the line of sight, rather than the decoupling of neutral atoms from plasma.

The role of slow-mode MHD shocks in magnetic reconnection is one of great importance for energy conversion and transport, but in many astrophysical plasmas the plasma is not fully ionised. In this paper, we investigate, using numerical simulations, the role of collisional coupling between a proton-electron charge-neutral fluid and a neutral hydrogen fluid for the 1D Riemann problem initiated in a constant pressure and density background state by a discontinuity in the magnetic field. This system, in the MHD limit, is characterised by two waves: a fast-mode rarefaction wave that drives a flow towards a slow-mode MHD shock. The system evolves through four stage: initiation, weak coupling, intermediate coupling and a quasi steady state. The initial stages are characterised by an over-pressured neutral region that expands with characteristics of a blast wave. In the later stages, the system tends towards a self-similar solution where the main drift velocity is concentrated in the thin region of the shock front. Due to the nature of the system, the neutral fluid is overpressured by the shock when compared to a purely hydrodynamic shock which results in the neutral fluid expanding to form the shock precursor. The thickness of the shockfront once it has formed proportional to the ionisation fraction to the power -1.2, which is a smaller exponent than would be naively expected from simple scaling arguments. One interesting result is that the shock front is a continuous transition of the physical variables for sub-sonic velocity upstream of the shock front (a c-shock) to a sharp jump in the physical variables followed by a relaxation to the downstream values for supersonic upstream velocity (a j-shock). The frictional heating that results from the velocity drift across the shock front can amount to approximately two per cent of the reference magnetic energy.

The magnetic Rayleigh-Taylor instability has been shown to play a key role in many astrophysical systems. The equation for the growth rate of this instability in the incompressible limit, and the most-unstable mode that can be derived from it, are often used to estimate the strength of the magnetic field that is associated with the observed dynamics. However, there are some issues with the interpretations given. Here we show that the class of most unstable modes $k_u$ for a given $\theta$, the class of modes often used to estimate the strength of the magnetic field from observations, for the system leads to the instability growing as $\sigma^2=1/2 A g k_u$, a growth rate which is independent of the strength of the magnetic field and which highlights that small scales are preferred by the system, but not does not give the fastest growing mode for that given $k$. We also highlight that outside of the interchange ($\mathbf{k}\cdot\mathbf{B}=0$) and undular ($\mathbf{k}$ parallel to $\mathbf{B}$) modes, all the other modes have a perturbation pair of the same wavenumber and growth rate that when excited in the linear regime can result in an interference pattern that gives field aligned filamentary structure often seen in 3D simulations. The analysis was extended to a sheared magnetic field, where it was found that it was possible to extend the results for a non-sheared field to this case. We suggest that without magnetic shear it is too simplistic to be used to infer magnetic field strengths in astrophysical systems.

Quiescent prominences host a diverse range of flows, including Rayleigh-Taylor instability driven upflows and impulsive downflows, and so it is no surprise that turbulent motions also exist. As prominences are believed to have a mean horizontal guide field, investigating any turbulence they host could shed light on the nature of MHD turbulence in a wide range of astrophysical systems. In this paper we have investigated the nature of the turbulent prominence motions using structure function analysis on the velocity increments estimated from H$\alpha$ Dopplergrams constructed with observational data from Hinode SOT. The pdf of the velocity increments shows that as we look at increasingly small spatial separations the distribution displays greater departure from a reference Gaussian distribution, hinting at intermittency in the velocity field. Analysis of the even order structure functions for both the horizontal and vertical separations showed the existence of two distinct regions displaying different exponents of the power law with the break in the power law at approximately 2000km. We hypothesise this to be a result of internal turbulence excited in the prominence by the dynamic flows of the system found at this spatial scale. We found that the scaling exponents of the p-th order structure functions for these two regions generally followed the p/2 (smaller scales) and p/4 (larger scales) laws that are the same as those predicted for weak MHD turbulence and Kraichnan-Iroshnikov turbulence respectively. However, the existence of the p/4 scaling at larger scales than the p/2 scaling is inconsistent with the increasing nonlinearity expected in MHD turbulence. Estimating the heating from the turbulent energy dissipation showed that this turbulent heating would be very inefficient, but that the mass diffusion through turbulence driven reconnection was of the order of $10^{10}$cm$^2$/s.

The recent observations of Singh et al. (2012) have shown multiple plasma ejections and the intermittent nature of magnetic reconnection in the solar chromosphere, highlighting the need for fast reconnection to occur in highly collisional plasma. However, the physical process through which fast magnetic reconnection occurs in partially ionized plasma, like the solar chromosphere, is still poorly understood. It has been shown that for sufficiently high magnetic Reynolds numbers, Sweet-Parker current sheets can become unstable leading to tearing mode instability and plasmoid formation, but when dealing with a partially ionized plasma the strength of coupling between the ions and neutrals plays a fundamental role in determining the dynamics of the system. We propose that as the reconnecting current sheet thins and the tearing instability develops, plasmoid formation passes through strongly, intermediately, and weakly coupled (or decoupled) regimes, with the time scale for the tearing mode instability depending on the frictional coupling between ions and neutrals. We present calculations for the relevant time scales for fractal tearing in all three regimes. We show that as a result of the tearing mode instability and the subsequent non-linear instability due to the plasmoid-dominated reconnection, the Sweet-Parker current sheet tends to have a fractal-like structure, and when the chromospheric magnetic field is sufficiently strong the tearing instability can reach down to kinetic scales, which are hypothesized to be necessary for fast reconnection.

Kepler data from G-, K- and M-type stars are used to study conditions that lead to superflares with energies above $10^{34} {\rm erg}$. From the 117,661 stars included, 380 show superflares with a total of 1690 such events. We study whether parameters, like effective temperature or the rotation rate, have any effect on the superflare occurrence rate or energy. With increasing effective temperature we observe a decrease in the superflare rate, which is analogous to the previous findings of a decrease in dynamo activity with increasing effective temperature. For slowly rotating stars, we find a quadratic increase of the mean occurrence rate with the rotation rate up to a critical point, after which the rate decreases linearly. Motivated by standard dynamo theory, we study the behavior of the relative starspot coverage, approximated as the relative brightness variation. For faster rotating stars, an increased fraction of stars shows higher spot coverage, which leads to higher superflare rates. A turbulent dynamo is used to study the dependence of the Ohmic dissipation as a proxy of the flare energy on the differential rotation or shear rate. The resulting statistics of the dissipation energy as a function of dynamo number is similar to the observed flare statistics as a function of the inverse Rossby number and shows similarly strong fluctuations. This supports the idea that superflares might well be possible for solar-type G stars.

The source of the non-thermal energy required for the heating of the upper solar atmosphere to temperatures in excess of a million degrees and the acceleration of the solar wind to hundreds of kilometres per second is still unclear. One such mechanism for providing the required energy flux is incompressible torsional Alfv\'en and kink magnetohydrodynamic (MHD) waves, which are magnetically dominated waves supported by the Sun's pervasive and complex magnetic field. In particular, propagating MHD kink waves have recently been observed to be ubiquitous throughout the solar atmosphere, but, until now, critical details of the transport of the kink wave energy throughout the Sun's atmosphere were unclear. Here, the ubiquity of the waves is exploited for statistical studies in the highly dynamic solar chromosphere. This large-scale investigation allows for the determination of the chromospheric kink wave velocity power spectra, a missing link necessary for determining the energy transport between the photosphere and corona. Crucially, the power spectra contains evidence for horizontal photospheric motions being the main mechanism for kink wave generation in the quiescent Sun. In addition, a comparison to measured coronal power spectra is provided, revealing frequency-dependent transmission profiles suggesting there is enhanced damping of kink waves in the lower corona.

This work examines infalling matter following an enormous Coronal Mass Ejection (CME) on 2011 June 7. The material formed discrete concentrations, or blobs, in the corona and fell back to the surface, appearing as dark clouds against the bright corona. In this work we examined the density and dynamic evolution of these blobs in order to formally assess the intriguing morphology displayed throughout their descent. The blobs were studied in five wavelengths (94, 131, 171, 193 and 211 \AA) using the Solar Dynamics Observatory Atmospheric Imaging Assembly (SDO/AIA), comparing background emission to attenuated emission as a function of wavelength to calculate column densities across the descent of four separate blobs. We found the material to have a column density of hydrogen of approximately 2 $\times$ 10$^{19}$ cm$^{-2}$, which is comparable with typical pre-eruption filament column densities. Repeated splitting of the returning material is seen in a manner consistent with the Rayleigh-Taylor instability. Furthermore, the observed distribution of density and its evolution are also a signature of this instability. By approximating the three-dimensional geometry (with data from STEREO-A), volumetric densities were found to be approximately 2 $\times$ 10$^{-14}$ g cm$^{-3}$, and this, along with observed dominant length-scales of the instability, was used to infer a magnetic field of the order 1 G associated with the descending blobs.

The launch of the Hinode satellite has allowed for seeing-free observations at high-resolution and high-cadence making it well suited to study the dynamics of quiescent prominences. In recent years it has become clear that quiescent prominences support small-amplitude transverse oscillations, however, sample sizes are usually too small for general conclusions to be drawn. We remedy this by providing a statistical study of transverse oscillations in vertical prominence threads. Over a three-hour period of observations it was possible to measure the properties of 3436 waves, finding periods from 50 to 6000 s with typical velocity amplitudes ranging between 0.2 to 23 km s$^{-1}$. The large number of observed waves allows the determination of the frequency dependence of the wave properties and derivation of the velocity power spectrum for the transverse waves. For frequencies less than 7 mHz, the frequency-dependence of the velocity power is consistent with the velocity power spectra generated from observations of the horizontal motions of magnetic elements in the photosphere, suggesting that the prominence transverse waves are driven by photospheric motions. However, at higher frequencies the two distributions significantly diverge, with relatively more power found at higher frequencies in the prominence oscillations. These results highlight that waves over a large frequency range are ubiquitous in prominences, and that a significant amount of the wave energy is found at higher frequency.

The dense prominence material is believed to be supported against gravity through the magnetic tension of dipped coronal magnetic field. For quiescent prominences, which exhibit many gravity-driven flows, hydrodynamic forces are likely to play an important role in the determination of both the large and small scale magnetic field distributions. In this study, we present the first steps toward creating three-dimensional magneto-hydrostatic prominence model where the prominence is formed in the dips of a coronal flux tube. Here 2.5D equilibria are created by adding mass to an initially force-free magnetic field, then performing a secondary magnetohydrodynamic relaxation. Two inverse polarity magnetic field configurations are studied in detail, a simple o-point configuration with a ratio of the horizontal field (B_x) to the axial field (B_y) of 1:2 and a more complex model that also has an x-point with a ratio of 1:11. The models show that support against gravity is either by total pressure or tension, with only tension support resembling observed quiescent prominences. The o-point of the coronal flux tube was pulled down by the prominence material, leading to compression of the magnetic field at the base of the prominence. Therefore tension support comes from the small curvature of the compressed magnetic field at the bottom and the larger curvature of the stretched magnetic field at the top of the prominence. It was found that this method does not guarantee convergence to a prominence-like equilibrium in the case where an x-point exists below the prominence flux tube. The results imply that a plasma beta of ~ 0.1 is necessary to support prominences through magnetic tension.

Recent observations of solar type stars with the Kepler satellite by Maehara et al. have revealed the existence of superflares (with energy of 10^33 - 10^35 erg) on Sun-like stars, which are similar to our Sun in their surface temperature (5600 K - 6000 K) and slow rotation (rotational period > 10 days). From the statistical analysis of these superflares, it was found that superflares with energy 10^34 erg occur once in 800 years and superflares with 10^35 erg occur once in 5000 years on Sun-like stars. In this paper, we examine whether superflares with energy of 10^33 - 10^35 erg could occur on the present Sun through the use of simple order-of-magnitude estimates based on current ideas relating to the mechanisms of the solar dynamo.

The launch of Hinode satellite led to the discovery of rising plumes, dark in chromospheric lines, in quiescent prominences that propagate from large (~10 Mm) bubbles that form at the base of the prominences. These plumes present a very interesting opportunity to study Magnetohydrodynamic (MHD) phenomena in quiescent prominences, but obstacles still remain. One of the biggest issues is that of the magnetic field strength, which is not easily measurable in prominences. In this paper we present a method that may be used to determine a prominence's plasma \beta when rising plumes are observed. Using the classic fluid dynamic solution for flow around a circular cylinder with an MHD correction, the compression of the prominence material can be estimated. This has been successfully confirmed through simulations; application to a prominence gave an estimate of the plasma beta as \beta=0.47 to 1.13 with an error of 0.080 for the range \gamma=1.4 to 1.7. Using this method it may be possible to estimate the plasma beta of observed prominences, therefore helping our understanding of a prominence's dynamics in terms of MHD phenomena.

The launch of the Hinode satellite has allowed unprecedented high resolution, stable images of solar quiescent prominences to be taken over extended periods of time. These new images led to the discovery of dark upflows that propagated from the base of prominences developing highly turbulent profiles. As yet, how these flows are driven is not fully understood. To study the physics behind this phenomena we use 3-D magnetohydrodynamic (MHD) simulations to investigate the nonlinear stability of the Kippenhahn-Shl\"{u}ter prominence model to the magnetic Rayleigh-Taylor instability. The model simulates the rise of a buoyant tube inside a quiescent prominence, where the upper boundary between the tube and prominence model is perturbed to excite the interchange of magnetic field lines. We found upflows of constant velocity (maximum found 6\,km\,s$^{-1}$) and a maximum plume width $\approx 1500$\,km which propagate through a height of approximately 6\,Mm in the no guide field case. The case with the strong guide field (initially $B_y=2B_x$) results in a large plume that rises though the prominence model at $\sim 5$\,km\,s$^{-1}$ with width $\sim 900$\,km (resulting in width of 2400\,km when viewed along the axis of the prominence) reaching a height of $\sim 3.1$\,Mm. In both cases nonlinear processes were important for determining plume dynamics.

We report findings from 0.2" resolution observations of the 2007 October 03 quiescent prominence observed with the Solar Optical Telescope on the Hinode satellite. The observations show clear ejections from the top of the quiescent prominence of plasma blobs. The ejections, originating from the top of rising prominence threads, are impulsively accelerated to Alfvenic velocities and then undergo ballistic motion. The ejections have a characteristic size between ~ 1000 - 2000 km. These characteristics are similar to downwardly propagating knots (typical size ~ 700 km) that have been observed in prominence threads, we suggest that the plasma blob ejections could be the upward moving counterpart to the downwardly propagating knots. We discuss the tearing instability as a possible mechanism to explain the ejections.

We present the results from 1.5D diffusion simulations of the Kippenhahn-Schlueter prominence model magnetic field evolution under the influence of the ambipolar terms of Cowling resistivity. We show that initially the evolution is determined by the ratio of the horizontal and vertical magnetic fields, which gives current sheet thinning (thickening) when this ratio is large (small) and a marginal case where a new characteristic current sheet length scale is formed. After a timespan greater than the Cowling resistivity time, the current sheet thickens as a power law of $t$ independent of the ratio of the field strengths. These results imply that when Cowling resistivity is included in the model, the tearing instability time scale is reduced by more than one order of magnitude when the ratio of the horizontal field to the vertical field is 20\% or less. These results imply that, over the course of its lifetime, the structure of the prominence can be significantly altered by Cowling resistivity, and in some cases will allow the tearing instability to occur.

We studied spicular jets over a plage area and derived their dynamic characteristics using Hinode Solar Optical Telescope (SOT) high-resolution images. The target plage region was near the west limb of the solar disk. This location permitted us to study the dynamics of spicular jets without the overlapping effect of spicular structures along the line of sight. In this work, to increase the ease with which we can identify spicules on the disk, we applied the image processing method `MadMax' developed by Koutchmy et al. (1989). It enhances fine, slender structures (like jets), over a diffuse background. We identified 169 spicules over the target plage. This sample permits us to derive statistically reliable results regarding spicular dynamics. The properties of plage spicules can be summarized as follows: (1) In a plage area, we clearly identified spicular jet features. (2) They were shorter in length than the quiet region limb spicules, and followed ballistic motion under constant deceleration. (3) The majority (80%) of the plage spicules showed the cycle of rise and retreat, while 10% of them faded out without a complete retreat phase. (4) The deceleration of the spicule was proportional to the velocity of ejection (i.e. the initial velocity).