The power spectrum of the evolving solar wind shows evidence of a spectral break between an inertial range of turbulent fluctuations at higher frequencies and a "$1/f$" like region at lower frequencies. In the ecliptic plane at $\sim 1$ AU, this break occurs approximately at timescales of a few hours, and is observed in the power spectra of components of velocity and magnetic field. The "$1/f$" energy range is of more direct coronal origin than the inertial range, and carries signatures of the complex magnetic field structure of the solar corona, and of footpoint stirring in the solar photosphere. To quantify the scaling properties we use generic statistical methods such as generalised structure functions and PDFs, focusing on solar cycle dependence and on anisotropy with respect to the background magnetic field. We present structure function analysis of magnetic and velocity field fluctuations, using a novel technique to decompose the fluctuations into directions parallel and perpendicular to the mean local background magnetic field. Whilst the magnetic field is close to "$1/f$", we show that the velocity field is "$1/f^{\alpha}$" with $\alpha\neq1$. For the velocity, the value of $\alpha$ varies between parallel and perpendicular fluctuations and with the solar cycle. There is also variation in $\alpha$ with solar wind speed. We have examined the PDFs in the fast, quiet solar wind and intriguingly, whilst parallel and perpendicular are distinct, both the ${\boldmath$B$}$ field and velocity show the same PDF of their perpendicular flucutations, which is close to gamma or inverse Gumbel. These results point to distinct physical processes in the corona, and to their mapping out into the solar wind. The scaling exponents obtained constrain the models for these processes.

Solar wind fluctuations, such as magnetic field or velocity, show power law power spectra suggestive both of an inertial range of intermittent turbulence (with $\sim -5/3$ exponent) and at lower frequencies, of fluctuations of coronal origin (with $\sim -1$ exponent). The ULYSSES spacecraft spent many months in the quiet fast solar wind above the Sun's polar coronal holes in a highly ordered magnetic field. We use statistical analysis methods such as the generalized structure function (GSF) and extended self-similarity (ESS) to quantify the scaling of the moments of the probability density function of fluctuations in the magnetic field. The GSFs give power law scaling in the ``$f^{-1}$'' range of the form $<| y(t+\tau)-y(t)|^{m}>\sim\tau^{\zeta(m)}$, but ESS is required to reveal scaling in the inertial range, which is of the form $<| y(t+\tau)-y(t)|^{m}>\sim [g(\tau)]^{\zeta(m)}$. We find that $g(\tau)$ is independent of spacecraft position and $g(\tau)\sim\tau^{-log_{10}(\tilde{\lambda}\tau)}$. The ``$f^{-1}$'' scaling fluctuates with radial spacecraft position. This confirms that, whereas the ``$f^{-1}$'' fluctuations are directly influenced by the corona, the inertial range fluctuations are consistent with locally evolving turbulence, but with an ``envelope'' $g(\tau)$, which captures the formation of the quiet fast solar wind.

We investigate the spatial correlation properties of the solar wind using simultaneous observations by the ACE and WIND spacecraft. We use mutual information as a nonlinear measure of correlation and compare this to linear correlation. We find that the correlation lengthscales of fluctuations in density and magnetic field magnitude vary strongly with the solar cycle, whereas correlation lengths of fluctuations in B field components do not. We find the correlation length of |B| ~ 120 Re at solar minimum and ~ 270 Re at maximum and the correlation length of density ~ 75 Re at minimum and ~ 170 Re at minimum. The components of the B field have correlation lengths ~ correlation length |B| at minimum.

Particle acceleration by means of non-linear plasma wave interactions is of great topical interest. Accordingly, in this paper we focus on the electron surfing process. Self-consistent kinetic simulations, using both relativistic Vlasov and PIC (Particle In Cell) approaches, show here that electrons can be accelerated to highly relativistic energies (up to 100 m_e c^2) if the phase speed of the electrostatic wave is mildly relativistic (0.6c to 0.9c for the magnetic field strengths considered). The acceleration is strong because of relativistic stabilisation of the nonlinearly saturated electrostatic wave, seen in both relativistic Vlasov and PIC simulations. An inverse power law momentum distribution can arise for the most strongly accelerated electrons. These results are of relevance to observed rapid changes in the radio synchrotron emission intensities from microquasars, gamma ray bursts and other astrophysical objects that require rapid acceleration mechanisms for electrons.

The determination of delay times of solar wind conditions at the sunward libration point to effects on Earth is investigated using mutual information. This measures the amount of information shared between two timeseries. We consider the mutual information content of solar wind observations, from WIND, and the geomagnetic indices. The success of five commonly used schemes for estimating interplanetary propagation times is examined. Propagation assuming a fixed plane normal at 45 degrees to the GSE x-axis (i.e. the Parker Spiral estimate) is found to give optimal mutual information. The mutual information depends on the point in space chosen as the target for the propagation estimate, and we find that it is maximized by choosing a point in the nightside rather than dayside magnetosphere. In addition, we employ recurrence plot analysis to visualize contributions to the mutual information, this suggests that it appears on timescales of hours rather than minutes.

The identification of pre-acceleration mechanisms for cosmic ray ions in supernova remnant shocks is an important problem in astrophysics. Recent particle-in-cell (PIC) shock simulations have shown that inclusion of the full electron kinetics yields non-time-stationary solutions, in contrast to previous hybrid (kinetic ions, fluid electrons) simulations. Here, by running a PIC code at high phase space resolution, ion acceleration mechanisms associated with the time dependence of a supercritical collisionless perpendicular shock are examined. In particular the components of $\int \mathbf{F} \cdot \mathbf{v} dt$ are analysed along trajectories for ions that reach both high and low energies. Selection mechanisms for the ions that reach high energies are also examined. In contrast to quasi-stationary shock solutions, the suprathermal protons are selected from the background population on the basis of the time at which they arrive at the shock, and thus are generated in bursts.

We investigate the distribution of fluctuations in solar irradiance when integrated over the full disk, obtained using extreme ultraviolet/soft X-ray observations from the SOHO CELIAS/SEM instrument. This time series sums over both the contributions of single distinguishable flares, and of many other processes. By detrending we select events with timescales of less than a few hours such as waves, slow flows, and CMEs. The statistics generated by this range of phenomena can be characterised by power-law-tailed distributions. We show that (i) during the high-activity period 2000 Jan-June the tail exponent a(T)=1.5 +/- 0.1; (ii) during the low-activity period 1996 Jan-June a(T)=3.0 +/- 0.2; and (iii) in general a(T) decreases with increasing activity.

Measurement of the energy emitted from accreting astrophysical systems provides an observational constraint on the plasma processes that may be operating within the disk. Here we examine the continual time variation over the past six years of the total X-ray flux from the microquasar GRS 1915+105. The application of differencing and rescaling techniques to RXTE/ASM data shows that the small amplitude fluctuations scale up to 12-17 days. A 17-day timescale in the X-ray fluctuations corresponds to half the measured binary orbital period of this system (33.5 +/- 1.5 days). While this may be coincidental, it is possible that these two timescales may be linked by, for example, a turbulent cascade in the accretion disk driven by a tidally-induced two-armed spiral shock corotating with the binary system. Temporal scaling is found only in the ever-present small fluctuations, and not in the intermittent larger-amplitude fluctuations. This is consistent with the basic model for this source which consists of a steady, cold outer disk and an unstable inner disk.

Whilst direct observations of internal transport in accretion disks are not yet possible, measurement of the energy emitted from accreting astrophysical systems can provide useful information on the physical mechanisms at work. Here we examine the unbroken multi-year time variation of the total X-ray flux from three sources: Cygnus X-1, the microquasar GRS1915+105, and for comparison the nonaccreting Crab nebula. To complement previous analyses, we demonstrate that the application of advanced statistical methods to these observational time-series reveals important contrasts in the nature and scaling properties of the transport processes operating within these sources. We find the Crab signal resembles Gaussian noise; the Cygnus X-1 signal is a leptokurtic random walk whose self-similar properties persist on timescales up to three years; and the GRS1915+105 signal is similar to that from Cygnus X-1, but with self-similarity extending possibly to only a few days. This evidence of self-similarity provides a robust quantitative characterisation of anomalous transport occuring within the systems.

Motivated by recent unsuccessful efforts to detect the predicted flux of TeV gamma-rays from supernova remnants, we present a critical examination of the theory on which these predictions are based. Three crucial problems are identified: injection, maximum achievable particle energy and spectral index. In each case significant new advances in understanding have been achieved, which cast doubt on prevailing paradigms such as Bohm diffusion and single-fluid MHD. This indicates that more realistic analytical models, backed by more sophisticated numerical techniques should be employed to obtain reliable predictions. Preliminary work on incorporating the effects of anomalous transport suggest that the resulting spectrum should be significantly softer than that predicted by conventional theory.

As noted by Chang, the hypothesis of Self-Organised Criticality provides a theoretical framework in which the low dimensionality seen in magnetospheric indices can be combined with the scaling seen in their power spectra and the recently-observed plasma bursty bulk flows. As such, it has considerable appeal, describing the aspects of the magnetospheric fuelling:storage:release cycle which are generic to slowly-driven, interaction-dominated, thresholded systems rather than unique to the magnetosphere. In consequence, several recent numerical "sandpile" algorithms have been used with a view to comparison with magnetospheric observables. However, demonstration of SOC in the magnetosphere will require further work in the definition of a set of observable properties which are the unique "fingerprint" of SOC. This is because, for example, a scale-free power spectrum admits several possible explanations other than SOC. A more subtle problem is important for both simulations and data analysis when dealing with multiscale and hence broadband phenomena such as SOC. This is that finite length systems such as the magnetosphere or magnetotail will by definition give information over a small range of orders of magnitude, and so scaling will tend to be narrowband. Here we develop a simple framework in which previous descriptions of magnetospheric dynamics can be described and contrasted. We then review existing observations which are indicative of SOC, and ask if they are sufficient to demonstrate it unambiguously, and if not, what new observations need to be made?

Observations of synchrotron radiation across a wide range of wavelengths provide clear evidence that electrons are accelerated to relativistic energies in supernova remnants (SNRs). However, a viable mechanism for the pre-acceleration of such electrons to mildly relativistic energies has not yet been established. In this paper an electromagnetic particle-in-cell (PIC) code is used to simulate acceleration of electrons from background energies to tens of keV at perpendicular collisionless shocks associated with SNRs. Free energy for electron energization is provided by ions reflected from the shock front, with speeds greater than the upstream electron thermal speed. The PIC simulation results contain several new features, including: the acceleration, rather than heating, of electrons via the Buneman instability; the acceleration of electrons to speeds exceeding those of the shock-reflected ions producing the instability; and strong acceleration of electrons perpendicular to the magnetic field. Electron energization takes place through a variety of resonant and non-resonant processes, of which the strongest involves stochastic wave-particle interactions. In SNRs the diffusive shock process could then supply the final step required for the production of fully relativistic electrons. The mechanisms identified in this paper thus provide a possible solution to the electron pre-acceleration problem.

Self-organised criticality (SOC) has been suggested as a potentially powerful unifying paradigm for interpreting the structure of, and signals from, accretion systems. After reviewing the most promising sites where SOC might be observable, we consider the theoretical arguments for supposing that SOC can occur in accretion discs. Perhaps the most rigorous evidence is provided by numerical modelling of energy dissipation due to magnetohydrodynamic turbulence in accretion discs by G Geertsema & A Achterberg (A&A {\bf 255}, 427 (1992)); we investigate how "sandpile"-type dynamics arise in this model. It is concluded that the potential sites for SOC in accretion systems are numerous and observationally accessible, and that theoretical support for the possible occurrence of SOC can be derived from first principles.

The theory of first order Fermi acceleration at shocks assumes that particles diffuse due to scattering off slow-moving magnetic irregularities. However, cosmic rays are closely tied to magnetic field lines, and the transport process, particularly across the direction of the field, is likely to be more complicated. To describe cross field transport we employ recent extensions of the Rechester-Rosenbluth theory in localised stochastic regions of magnetic field. During acceleration at a shock and when the motion along the field is diffusive, there is a transition at a critical energy from \lq\lq sub-diffusive\rq\rq\ motion, where the mean square displacement of a particle increases with time as $t^{1/2}$, to compound diffusion, a combined process involving diffusion along a magnetic field which is itself wandering. Requiring this critical energy to be less than the system cut-off in a SNR of radius $R$ places an upper limit on the coherence length, $\lambda$, of the magnetic field for which diffusive shock acceleration can occur at quasi-perpendicular shocks when field line wandering dominates the effective transport of particles: $\lambda<R(U/v)^{1/2}$ where $U$ is the shock speed and $v$ the particle speed. The acceleration rate for particles in the sub-diffusive regime is derived and its implication for the origin of high energy cosmic rays is discussed.