Two of the central problems in our understanding of the solar chromosphere are how the upper chromosphere is heated and what drives spicules. Estimates of the required chromospheric heating, based on radiative and conductive losses suggest a rate of ${\sim} 0.1 \mathrm{\:erg\:cm^{-3}\:s^{-1}}$ in the lower chromosphere dropping to ${\sim} 10^{-3} \mathrm{\:erg\:cm^{-3}\:s^{-1}}$ in the upper chromosphere (\citet{Avrett1981}). The chromosphere is also permeated by spicules, higher density plasma from the lower atmosphere propelled upwards at speeds of ${\sim} 10-20 \mathrm{\:km\:s^{-1}}$, for so called Type-I spicules (\citet{Pereira2012,Zhang2012}), reaching heights of ${\sim} 3000-5000 \mathrm{\:km}$ above the photosphere. A clearer understanding of chromospheric dynamics, its heating and the formation of spicules, is thus of central importance to solar atmospheric science. For over thirty years it has been proposed that photospheric driving of MHD waves may be responsible for both heating and spicule formation. This paper presents results from a high-resolution MHD treatment of photospheric driven Alfv\'en and kink waves propagating upwards into an expanding flux tube embedded in a model chromospheric atmosphere. We show that the ponderomotive coupling from Alfv\'en and kink waves into slow modes generates shocks which both heat the upper chromosphere and drive spicules. These simulations show that wave driving of the solar chromosphere can give a local heating rate which matches observations and drive spicules consistent with Type-I observations all within a single coherent model.

Physical processes which may lead to solar chromospheric heating are analyzed using high-resolution 1.5D non-ideal MHD modelling. We demonstrate that it is possible to heat the chromospheric plasma by direct resistive dissipation of high-frequency Alfv\'en waves through Pedersen resistivity. However this is unlikely to be sufficient to balance radiative and conductive losses unless unrealistic field strengths or photospheric velocities are used. The precise heating profile is determined by the input driving spectrum since in 1.5D there is no possibility of Alfv\'en wave turbulence. The inclusion of the Hall term does not affect the heating rates. If plasma compressibility is taken into account, shocks are produced through the ponderomotive coupling of Alfv\'en waves to slow modes and shock heating dominates the resistive dissipation. In 1.5D shock coalescence amplifies the effects of shocks and for compressible simulations with realistic driver spectra the heating rate exceeds that required to match radiative and conductive losses. Thus while the heating rates for these 1.5D simulations are an overestimate they do show that ponderomotive coupling of Alfv\'en waves to sound waves is more important in chromospheric heating than Pedersen dissipation through ion-neutral collisions.

Aims: We perform an advanced study of the fast magnetoacoustic sausage oscillations of coronal loops in the context of MHD coronal seismology to establish the dependence of the sausage mode period and cut-off wavenumber on the plasma-beta of the loop-filling plasma. A parametric study of the ratios for different harmonics of the mode is also carried out. Methods: Full magnetohydrodynamic numerical simulations were performed using Lare2d, simulating hot, dense loops in a magnetic slab environment. The symmetric Epstein profile and a simple step-function profile were both used to model the density structure of the simulated loops. Analytical expressions for the cut-off wavenumber and the harmonic ratio between the second longitudinal harmonic and the fundamental were also examined. Results: It was established that the period of the global sausage mode is only very weakly dependent on the value of the plasma-beta inside a coronal loop, which justifies the application of this model to hot flaring loops. The cut-off wavenumber kc for the global mode was found to be dependent on both internal and external values of the plasma-beta, again only weakly. By far the most important factor in this case was the value of the density contrast ratio between the loop and the surroundings. Finally, the deviation of the harmonic ratio P1/2P2 from the ideal non-dispersive case was shown to be considerable at low k, again strongly dependent on plasma density. Quantifying the behaviour of the cut-off wavenumber and the harmonic ratio has significant applications to the field of coronal seismology.

This paper extends the models of Craig & McClymont (1991) and McLaughlin & Hood (2004) to include finite $\beta$ and nonlinear effects. We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 2001) and consider an initial condition in $ {\bf{v}}\times{\bf{B}} \cdot {\hat{\bf{z}}}$ (a natural variable of the system). We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a 'cusp-like' point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of the final current density left in the system. The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime.