Solar energetic particles (SEPs) are an important product of solar activity. They are connected to solar active regions and flares, coronal mass ejections (CMEs), EUV waves, shocks, Type II and III radio emissions, and X-ray bursts. These phenomena are major probes of the partition of energy in solar eruptions, as well as for the organization, dynamics, and relaxation of coronal and interplanetary magnetic fields. Many of these phenomena cause terrestrial space weather, posing multiple hazards for humans and their technology from space to the ground. Since particular flares, shocks, CMEs, and EUV waves produce SEP events but others do not, since propagation effects from the low corona to 1 AU appear important for some events but not others, and since Type II and III radio emissions and X-ray bursts are sometimes produced by energetic particles leaving these acceleration sites, it is necessary to study the whole system with a multi-frequency and multi-instrument perspective that combines both in-situ and remote observations with detailed modelling of phenomena. This article demonstrates this comprehensive approach, and shows its necessity, by analysing a trio of unusual and striking solar eruptions, radio and X-ray bursts, and SEP events that occurred on 4 November 2015. These events show both strong similarities and differences from standard events and each other, despite having very similar interplanetary conditions and only two are sites and CME genesis regions. They are therefore major targets for further in-depth observational studies, and for testing both existing and new theories and models. Based on the very limited modelling available we identify the aspects that are and are not understood, and we discuss ideas that may lead to improved understanding of the SEP, radio, and space-weather events.

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories.

We present spectropolarimetric imaging observations of the solar corona at low frequencies (80 - 240 MHz) using the Murchison Widefield Array (MWA). These images are the first of their kind, and we introduce an algorithm to mitigate an instrumental artefact by which the total intensity signal contaminates the polarimetric images due to calibration errors. We then survey the range of circular polarization (Stokes V) features detected in over 100 observing runs near solar maximum during quiescent periods. First, we detect around 700 compact polarized sources across our dataset with polarization fractions ranging from less than 0.5% to nearly 100%. These sources exhibit a positive correlation between polarization fraction and total intensity, and we interpret them as a continuum of plasma emission noise storm (Type I burst) continua sources associated with active regions. Second, we report a characteristic "bullseye" structure observed for many low-latitude coronal holes in which a central polarized component is surrounded by a ring of the opposite sense. The central component does not match the sign expected from thermal bremsstrahlung emission, and we speculate that propagation effects or an alternative emission mechanism may be responsible. Third, we show that the large-scale polarimetric structure at our lowest frequencies is reasonably well-correlated with the line-of-sight (LOS) magnetic field component inferred from a global potential field source surface (PFSS) model. The boundaries between opposite circular polarization signs are generally aligned with polarity inversion lines in the model at a height roughly corresponding to that of the radio limb. This is not true at our highest frequencies, however, where the LOS magnetic field direction and polarization sign are often not straightforwardly correlated.

Using a refined setup process, we simulated the propagation of six observed Coronal Mass Ejections (CMEs) with the 2012 Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme (BATS-R-US) code from the Sun to the Earth or STEREO A and compared the outputs with observations. A linear relation between the average CME speed below 6 solar radii and the flux rope current is demonstrated and used to tune the simulations. The simulations correctly predict if and when an observable CME shock reaches one astronomical unit (AU). The arrival time predictions of the CME shocks at 1 AU have an accuracy of 0.9 $\pm$ 1.9 hours. The simulated initial CME speeds and average accelerations are close to the model and data of Gopalswamy et al., 2000. The approach shows promise for predicting the sense of the predominant shock-associated change in the magnetic field component $B_z$. However, the magnetic fields and plasma conditions in the solar wind and CME are not predicted well quantitatively.

Solar radio emission, especially at metre-wavelengths, is well known to vary over small spectral ($\lesssim$100\,kHz) and temporal ($<1$\,s) spans. It is comparatively recently, with the advent of a new generation of instruments, that it has become possible to capture data with sufficient resolution (temporal, spectral and angular) that one can begin to characterize the solar morphology simultaneously along the axes of time and frequency. This ability is naturally accompanied by an enormous increase in data volumes and computational burden, a problem which will only become more acute with the next generation of instruments such as the Square Kilometre Array (SKA). The usual approach, which requires manual guidance of the calibration process, is impractical. Here we present the "Automated Imaging Routine for Compact Arrays for the Radio Sun (AIRCARS)", an end-to-end imaging pipeline optimized for solar imaging with arrays with a compact core. We have used AIRCARS so far on data from the Murchison Widefield Array (MWA) Phase-I. The dynamic range of the images is routinely from a few hundred to a few thousand. In the few cases, where we have pushed AIRCARS to its limits, the dynamic range can go as high as $\sim$75000. The images made represent a substantial improvement in the state-of-the-art in terms of imaging fidelity and dynamic range. This has the potential to transform the multi-petabyte MWA solar archive from raw visibilities into science-ready images. AIRCARS can also be tuned to upcoming telescopes like the SKA, making it a very useful tool for the heliophysics community.

We describe the motivation and design details of the "Phase II" upgrade of the Murchison Widefield Array (MWA) radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the MWA in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing MWA core. These new tiles enhance the surface brightness sensitivity of the MWA and will improve the ability of the MWA to estimate the slope of the Epoch of Reionisation power spectrum by a factor of ~3.5. The remaining 56 tiles are deployed on long baselines, doubling the maximum baseline of the array and improving the array u,v coverage. The improved imaging capabilities will provide an order of magnitude improvement in the noise floor of MWA continuum images. The upgrade retains all of the features that have underpinned the MWA's success (large field-of-view, snapshot image quality, pointing agility) and boosts the scientific potential with enhanced imaging capabilities and by enabling new calibration strategies.

We present coronal density profiles derived from low-frequency (80-240 MHz) imaging of three type III solar radio bursts observed at the limb by the Murchison Widefield Array (MWA). Each event is associated with a white light streamer at larger heights and is plausibly associated with thin extreme ultraviolet rays at lower heights. Assuming harmonic plasma emission, we find average electron densities of 1.8 x10^8 cm^-3 down to 0.20 x10^8 cm^-3 at heights of 1.3 to 1.9 solar radii. These values represent roughly 2.4-5.4x enhancements over canonical background levels and are comparable to the highest streamer densities obtained from data at other wavelengths. Assuming fundamental emission instead would increase the densities by a factor of 4. High densities inferred from type III source heights can be explained by assuming that the exciting electron beams travel along overdense fibers or by radio propagation effects that may cause a source to appear at a larger height than the true emission site. We review the arguments for both scenarios in light of recent results. We compare the extent of the quiescent corona to model predictions to estimate the impact of propagation effects, which we conclude can only partially explain the apparent density enhancements. Finally, we use the time- and frequency-varying source positions to estimate electron beam speeds of between 0.24 and 0.60 c.

Low frequency imaging radio arrays such as MWA, LWA and LOFAR have been recently commissioned, and significantly more advanced and flexible arrays are planned for the near term. These powerful instruments offer new opportunities for direct solar imaging at high time and frequency resolution. They can also probe large volumes of the heliosphere simultaneously, by virtue of very large fields of view. They allow highly detailed, spatially resolved study of solar and heliospheric radio bursts, which are complemented by heliospheric propagation studies using both background astronomical radio emissions as well as the bursts themselves. In this paper, the state of the art in such wide field solar and heliospheric radio studies is summarized, including recent results from the Murchison Widefield Array (MWA). The prospects for major advances in observational capabilities in the near future are reviewed, with particular emphasis on the RAPID system developed at Haystack Observatory

We present low-frequency (80-240 MHz) radio imaging of type III solar radio bursts observed by the Murchison Widefield Array (MWA) on 2015/09/21. The source region for each burst splits from one dominant component at higher frequencies into two increasingly-separated components at lower frequencies. For channels below ~132 MHz, the two components repetitively diverge at high speeds (0.1-0.4 c) along directions tangent to the limb, with each episode lasting just ~2 s. We argue that both effects result from the strong magnetic field connectivity gradient that the burst-driving electron beams move into. Persistence mapping of extreme ultraviolet (EUV) jets observed by the Solar Dynamics Observatory reveals quasi-separatrix layers (QSLs) associated with coronal null points, including separatrix dome, spine, and curtain structures. Electrons are accelerated at the flare site toward an open QSL, where the beams follow diverging field lines to produce the source splitting, with larger separations at larger heights (lower frequencies). The splitting motion within individual frequency bands is interpreted as a projected time-of-flight effect, whereby electrons traveling along the outer field lines take slightly longer to excite emission at adjacent positions. Given this interpretation, we estimate an average beam speed of 0.2 c. We also qualitatively describe the quiescent corona, noting in particular that a disk-center coronal hole transitions from being dark at higher frequencies to bright at lower frequencies, turning over around 120 MHz. These observations are compared to synthetic images based on the Magnetohydrodynamic Algorithm outside a Sphere (MAS) model, which we use to flux-calibrate the burst data.

We present techniques developed to calibrate and correct Murchison Widefield Array (MWA) low frequency (72-300 MHz) radio observations for polarimetry. The extremely wide field-of-view, excellent instantaneous (u, v)-coverage and sensitivity to degree-scale structure that the MWA provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the MWA, we discuss future directions for polarimetric science at low frequencies to answer outstanding questions relating to polarised source counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the nature of the interstellar and intergalactic media, and the solar environment.

An intense type IV radio burst was observed by the STEREO Behind (STB) spacecraft located about 144 degres behind Earth. The burst was associated with a large solar eruption that occurred on the backside of the Sun (N05E151) close to the disk center in the STB view. The eruption was also observed by the STEREO Ahead (STA) spacecraft (located at 149 degrees ahead of Earth) as an eruption close to the west limb (N05W60) in that view. The type IV burst was complete in STB observations in that the envelope reached the lowest frequency and then receded to higher frequencies. The burst was partial viewed from STA, revealing only the edge coming down to the lowest frequency. The type IV burst was not observed at all near Earth because the source was 61 degrees behind the east limb. The eruption was associated with a low-frequency type II burst observed in all three views, although it was not very intense. Solar energetic particles were also observed at both STEREOs and at SOHO, suggesting that the shock was much extended, consistent with the very high speed of the CME (about 2048 km/s). These observations suggest that the type IV emission is directed along a narrow cone above the flare site. We confirm this result statistically using the type IV bursts of solar cycle 23.

Electron density irregularities in the ionosphere are known to be magnetically anisotropic, preferentially elongated along the lines of force. While many studies of their morphology have been undertaken by topside sounding and whistler measurements, it is only recently that detailed regional-scale reconstructions have become possible, enabled by the advent of widefield radio telescopes. Here we present a new approach for visualising and studying field-aligned irregularities (FAIs), which involves transforming interferometric measurements of TEC gradients onto a magnetic shell tangent plane. This removes the perspective distortion associated with the oblique viewing angle of the irregularities from the ground, facilitating the decomposition of dynamics along and across magnetic field lines. We apply this transformation to the dataset of Loi et al. [2015a], obtained on 15 October 2013 by the Murchison Widefield Array (MWA) radio telescope and displaying prominent FAIs. We study these FAIs in the new reference frame, quantifying field-aligned and field-transverse behaviour, examining time and altitude dependencies, and extending the analysis to FAIs on sub-array scales. We show that the inclination of the plane can be derived solely from the data, and verify that the best-fit value is consistent with the known magnetic inclination. The ability of the model to concentrate the fluctuations along a single spatial direction may find practical application to future calibration strategies for widefield interferometry, by providing a compact representation of FAI-induced distortions.

Geomagnetically-aligned density structures with a range of sizes exist in the near-Earth plasma environment, including 10-100 km-wide VLF/HF wave-ducting structures. Their small diameters and modest density enhancements make them difficult to observe, and there is limited evidence for any of the several formation mechanisms proposed to date. We present a case study of an event on 26 August 2014 where a travelling ionospheric disturbance (TID) shortly precedes the formation of a complex collection of field-aligned ducts, using data obtained by the Murchison Widefield Array (MWA) radio telescope. Their spatiotemporal proximity leads us to suggest a causal interpretation. Geomagnetic conditions were quiet at the time, and no obvious triggers were noted. Growth of the structures proceeds rapidly, within 0.5 hr of the passage of the TID, attaining their peak prominence 1-2 hr later and persisting for several more hours until observations ended at local dawn. Analyses of the next two days show field-aligned structures to be preferentially detectable under quiet rather than active geomagnetic conditions. We used a raster scanning strategy facilitated by the speed of electronic beamforming to expand the quasi-instantaneous field of view of the MWA by a factor of three. These observations represent the broadest angular coverage of the ionosphere by a radio telescope to date.

We present observations of high-amplitude rapid (2 s) variability toward two bright, compact extragalactic radio sources out of several hundred of the brightest radio sources in one of the 30x30 deg MWA Epoch of Reionization fields using the Murchison Widefield Array (MWA) at 155 MHz. After rejecting intrinsic, instrumental, and ionospheric origins we consider the most likely explanation for this variability to be interplanetary scintillation (IPS), likely the result of a large coronal mass ejection propagating from the Sun. This is confirmed by roughly contemporaneous observations with the Ooty Radio Telescope. We see evidence for structure on spatial scales ranging from <1000 km to >1e6 km. The serendipitous night-time nature of these detections illustrates the new regime that the MWA has opened for IPS studies with sensitive night-time, wide-field, low-frequency observations. This regime complements traditional dedicated strategies for observing IPS and can be utilized in real-time to facilitate dedicated follow-up observations. At the same time, it allows large-scale surveys for compact (arcsec) structures in low-frequency radio sources despite the 2 arcmin resolution of the array.

The fields of solar radiophysics and solar system radio physics, or radio heliophysics, will benefit immensely from an instrument with the capabilities projected for SKA. Potential applications include interplanetary scintillation (IPS), radio-burst tracking, and solar spectral radio imaging with a superior sensitivity. These will provide breakthrough new insights and results in topics of fundamental importance, such as the physics of impulsive energy releases, magnetohydrodynamic oscillations and turbulence, the dynamics of post-eruptive processes, energetic particle acceleration, the structure of the solar wind and the development and evolution of solar wind transients at distances up to and beyond the orbit of the Earth. The combination of the high spectral, time and spatial resolution and the unprecedented sensitivity of the SKA will radically advance our understanding of basic physical processes operating in solar and heliospheric plasmas and provide a solid foundation for the forecasting of space weather events.

Low-frequency, wide field-of-view (FoV) radio telescopes such as the Murchison Widefield Array (MWA) enable the ionosphere to be sampled at high spatial completeness. We present the results of the first power spectrum analysis of ionospheric fluctuations in MWA data, where we examined the position offsets of radio sources appearing in two datasets. The refractive shifts in the positions of celestial sources are proportional to spatial gradients in the electron column density transverse to the line of sight. These can be used to probe plasma structures and waves in the ionosphere. The regional (10-100 km) scales probed by the MWA, determined by the size of its FoV and the spatial density of radio sources (typically thousands in a single FoV), complement the global (100-1000 km) scales of GPS studies and local (0.01-1 km) scales of radar scattering measurements. Our data exhibit a range of complex structures and waves. Some fluctuations have the characteristics of travelling ionospheric disturbances (TIDs), while others take the form of narrow, slowly-drifting bands aligned along the Earth's magnetic field.

Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct wide-angle observation of an extensive array of field-aligned ducts bridging the upper ionosphere and inner plasmasphere, using a novel ground-based imaging technique. We establish their heights and motions by feature-tracking and parallax analysis. The structures are strikingly organized, appearing as regularly-spaced, alternating tubes of overdensities and underdensities strongly aligned with the Earth's magnetic field. These findings represent the first direct visual evidence for the existence of such structures.

The amplitude of density turbulence in the extended solar corona, especially near the dissipation scale, impinges on several problems of current interest. Radio sources observed through the turbulent solar wind are broadened due to refraction by and scattering off density inhomogeneities, and observations of scatter broadening are often employed to constrain the turbulence amplitude. The extent of such scatter broadening is usually computed using the structure function, which gives a measure of the spatial correlation measured by an interferometer. Most such treatments have employed analytical approximations to the structure function that are valid in the asymptotic limits $s \gg l_{i}$ or $s \ll l_{i}$, where $s$ is the interferometer spacing and $l_{i}$ is the inner scale of the density turbulence spectrum. We instead use a general structure function (GSF) that straddles these regimes, and quantify the errors introduced by the use of these approximations. We have included the effects of anisotropic scattering for distant cosmic sources viewed through the solar wind at small elongations. We show that the regimes where the GSF predictions are more accurate than those of the asymptotic expressions are not only of practical relevance, but are where inner scale effects influence estimates of scatter broadening. Taken together, we argue that the GSF should henceforth be used for scatter broadening calculations and estimates of turbulence amplitudes in the solar corona and solar wind.

The angular broadening of compact radio sources observed through a medium having turbulent density irregularities is usually estimated using the phase structure function. We employ an exact formulation for the phase structure function that helps in obtaining an accurate estimate of angular broadening when the baseline lengths are comparable to the inner scale of the turbulent spectrum.

The Sun has remained a difficult source to image for radio telescopes, especially at the low radio frequencies. Its morphologically complex emission features span a large range of angular scales, emission mechanisms involved and brightness temperatures. In addition, time and frequency synthesis, the key tool used by most radio interferometers to build up information about the source being imaged is not effective for solar imaging, because many of the features of interest are short lived and change dramatically over small fractional bandwidths. Building on the advances in radio frequency technology, digital signal processing and computing, the kind of instruments needed to simultaneously capture the evolution of solar emission in time, frequency, morphology and polarization over a large spectral span with the requisite imaging fidelity, and time and frequency resolution have only recently begun to appear. Of this class of instruments, the Murchison Widefield Array (MWA) is best suited for solar observations. The MWA has now entered a routine observing phase and here we present some early examples from MWA observations.

Significant new opportunities for astrophysics and cosmology have been identified at low radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency. The telescope will enable new advances along four key science themes, including searching for redshifted 21 cm emission from the epoch of reionisation in the early Universe; Galactic and extragalactic all-sky southern hemisphere surveys; time-domain astrophysics; and solar, heliospheric, and ionospheric science and space weather. The Murchison Widefield Array is located in Western Australia at the site of the planned Square Kilometre Array (SKA) low-band telescope and is the only low-frequency SKA precursor facility. In this paper, we review the performance properties of the Murchison Widefield Array and describe its primary scientific objectives.

The Murchison Widefield Array is a low frequency (80 - 300 MHz) SKA Precursor, comprising 128 aperture array elements (known as tiles) distributed over an area of 3 km diameter. The MWA is located at the extraordinarily radio quiet Murchison Radioastronomy Observatory in the mid-west of Western Australia, the selected home for the Phase 1 and Phase 2 SKA low frequency arrays. The MWA science goals include: 1) detection of fluctuations in the brightness temperature of the diffuse redshifted 21 cm line of neutral hydrogen from the epoch of reionisation; 2) studies of Galactic and extragalactic processes based on deep, confusion-limited surveys of the full sky visible to the array; 3) time domain astrophysics through exploration of the variable radio sky; and 4) solar imaging and characterisation of the heliosphere and ionosphere via propagation effects on background radio source emission. This paper concentrates on the capabilities of the MWA for solar science and summarises some of the solar science results to date, in advance of the initial operation of the final instrument in 2013.

We seek to reconcile observations of small source sizes in the solar corona at 327 MHz with predictions of scattering models that incorporate refractive index effects, inner scale effects and a spherically diverging wavefront. We use an empirical prescription for the turbulence amplitude $C_{N}^{2}(R)$ based on VLBI observations by Spangler and coworkers of compact radio sources against the solar wind for heliocentric distances $R \approx$ 10--50 $R_{\odot}$. We use the Coles & Harmon model for the inner scale $l_{i}(R)$, that is presumed to arise from cyclotron damping. In view of the prevalent uncertainty in the power law index that characterizes solar wind turbulence at various heliocentric distances, we retain this index as a free parameter. We find that the inclusion of spherical divergence effects suppresses the predicted source size substantially. We also find that inner scale effects significantly reduce the predicted source size. An important general finding for solar sources is that the calculations substantially underpredict the observed source size. Three possible, non-exclusive, interpretations of this general result are proposed. First and simplest, future observations with better angular resolution will detect much smaller sources. Consistent with this, previous observations of small sources in the corona at metric wavelengths are limited by the instrument resolution. Second, the spatially-varying level of turbulence $C_{N}^{2}(R)$ is much larger in the inner corona than predicted by straightforward extrapolation Sunwards of the empirical prescription, which was based on observations between 10--50 $R_{\odot}$. Either the functional form or the constant of proportionality could be different. Third, perhaps the inner scale is smaller than the model, leading to increased scattering.

We present the first spectroscopic images of solar radio transients from the prototype for the Murchison Widefield Array (MWA), observed on 2010 March 27. Our observations span the instantaneous frequency band 170.9-201.6 MHz. Though our observing period is characterized as a period of `low' to `medium' activity, one broadband emission feature and numerous short-lived, narrowband, non-thermal emission features are evident. Our data represent a significant advance in low radio frequency solar imaging, enabling us to follow the spatial, spectral, and temporal evolution of events simultaneously and in unprecedented detail. The rich variety of features seen here reaffirms the coronal diagnostic capability of low radio frequency emission and provides an early glimpse of the nature of radio observations that will become available as the next generation of low frequency radio interferometers come on-line over the next few years.

Giant pulses and giant micropulses from pulsars are distinguished from normal pulsed emission by their large fluxes, rarity, approximately power-law distribution of fluxes and, typically, occurrence in restricted phase windows. Here existing observations of flux distributions are manipulated into a common format and interpreted in terms of theories for wave growth in inhomogeneous media, with the aim of constraining the emission mechanism and source physics for giant pulses and micropulses. Giant micropulses near 2 GHz (PSRs B8033-45 and B1706-44) and 0.4 GHz (PSR B0950+08) have indices $\alpha = 6.5 \pm 0.7$ for the probability distribution $P(E)$ of the electric field $E$, with $P(E) \propto E^{-\alpha}$. Giant pulses (PSRs B0531+24, B1937+214, and B1821-24) have $\alpha$ ranging from $4.6 \pm 0.2$ to $9 \pm 2$, possibly increasing with frequency. These are similar enough to regard giant micropulses and pulses as a single phenomenon with a common physical explanation. The power-law functional form and values of $\alpha$ observed are consistent with predictions for nonlinear wave collapse, but inconsistent with known self-organized critical systems, nonlinear decay processes, and elementary burst theory. While relativistic beaming may be important, its statistics are yet to be predicted theoretically and collapse is currently the favored interpretation. Other possibilities remain, including stochastic growth theory (consistent with normal pulse emission) and, less plausibly, refractive lensing. Unresolved issues remain for all four interpretations and suggestions for further work are given. The differences between normal and giant pulse emission suggest they have distinct source regions and emission processes.

Individual pulses from pulsars have intensity-phase profiles that differ widely from pulse to pulse, from the average profile, and from phase to phase within a pulse. Widely accepted explanations do not exist for this variability or for the mechanism producing the radiation. The variability corresponds to the field statistics, particularly the distribution of wave field amplitudes, which are predicted by theories for wave growth in inhomogeneous media. This paper shows that the field statistics of the Vela pulsar (PSR B0833-45) are well-defined and vary as a function of pulse phase, evolving from Gaussian intensity statistics off-pulse to approximately power-law and then lognormal distributions near the pulse peak to approximately power-law and eventually Gaussian statistics off-pulse again. Detailed single-component fits confirm that the variability corresponds to lognormal statistics near the peak of the pulse profile and Gaussian intensity statistics off-pulse. The lognormal field statistics observed are consistent with the prediction of stochastic growth theory (SGT) for a purely linear system close to marginal stability. The simplest interpretations are that the pulsar's variability is a direct manifestation of an SGT state and the emission mechanism is linear (either direct or indirect), with no evidence for nonlinear mechanisms like modulational instability and wave collapse which produce power-law field statistics. Stringent constraints are placed on nonlinear mechanisms: they must produce lognormal statistics when suitably ensemble-averaged. Field statistics are thus a powerful, potentially widely applicable tool for understanding variability and constraining mechanisms and source characteristics of coherent astrophysical and space emissions.

The variability of the Vela pulsar (PSR B0833-45) corresponds to well-defined field statistics that vary with pulsar phase, ranging from Gaussian intensity statistics off-pulse to approximately power-law statistics in a transition region and then lognormal statistics on-pulse, excluding giant micropulses. These data are analyzed here in terms of two superposed wave populations, using a new calculation for the amplitude statistics of two vectorially-combined transverse fields. Detailed analyses show that the approximately power-law and lognormal distributions observed are fitted well at essentially all on-pulse phases by Gaussian-lognormal and double-lognormal combinations, respectively. These good fits, plus the smooth but significant variations in fit parameters across the source, provide strong evidence that the approximately power-law statistics observed in the transition region are not intrinsic. Instead, the data are consistent with normal pulsar emission having lognormal statistics at all phases. This is consistent with generation in an inhomogeneous source obeying stochastic growth theory (SGT) and with the emission mechanism being purely linear (either direct or indirect). A nonlinear mechanism is viable only if it produces lognormal statistics when suitably ensemble-averaged. Variations in the SGT fit parameters with phase imply that the radiation is relatively more variable near the pulse edges than near the center, as found in earlier work. In contrast, Vela's giant micropulses come from a very restricted phase range and have power-law statistics with indices ($6.7 \pm 0.6$) not inconsistent with nonlinear wave collapse. These results imply that normal pulses have a different source and generation mechanism than giant micropulses, as suggested previously on other grounds.

Individual pulses from pulsars have intensity-phase profiles that differ widely from pulse to pulse, from the average profile, and from phase to phase within a pulse. Widely accepted explanations for pulsar radio emission and its time variability do not exist. Here, by analysing data near the peak of the Vela pulsar's average profile, we show that Vela's variability corresponds to lognormal field statistics, consistent with the prediction of stochastic growth theory (SGT) for a purely linear system close to marginal stability. Vela's variability is therefore a direct manifestation of an SGT state and the field statistics constrain the emission mechanism to be linear (either direct or indirect), ruling out nonlinear mechanisms like wave collapse. Field statistics are thus a powerful, potentially widely applicable tool for understanding variability and constraining mechanisms and source characteristics of coherent astrophysical and space emissions.