A phenomenological theory describes radial evolution of plasma turbulence in the solar wind from 1 to 50 astronomical units. The theory includes a simple closure for local anisotropic magnetohydrodynamic turbulence, spatial transport, and driving by large-scale shear and pickup ions. Results compare well to plasma and magnetic field data from the Voyager 2 spacecraft, providing a basis for a concise, tractable description of turbulent energy transport in a variety of astrophysical plasmas. [S0031-9007(99) Low-frequency fluctuations in the solar wind plasma represent perhaps the most extensively studied type of magnetohydrodynamic (MHD) turbulence, having been observed by spacecraft instruments for more than thirty years [1][2][3]. The observed turbulence displays properties expected of both hydrodynamic and MHD theories, including distinctive spectra and correlations [3,4]. Solar wind turbulence is a crucial element in coupling the lower corona plasma and the earth's magnetosphere, and in the transport of energetic charged particles throughout the solar-terrestrial environment. It is also a prototype for understanding stellar and galactic winds and astrophysical plasma flows in general. There has been notable progress in understanding the cascade process [5][6][7][8][9][10][11][12] that accompanies solar wind turbulence. So far, however, no single quantitative model has explained how turbulent energy flows from the largest interacting structures to the smallest dissipative scales where it is deposited as heat. In this Letter we present such a theory, based upon the dynamics of large-scale "eddies," which, controlled by a single similarity scale, drives a cascade that supplies thermal energy to the fluid plasma. The theoretical results compare well with measurements by the Voyager 2 spacecraft at heliocentric radial distances r from 1 astronomical units (AU) to beyond 30 AU. This motivates the development of similar phenomenological turbulence theories for nonlinear MHD flows in a variety of astrophysical plasmas.From the Helios and Mariner missions reaching inside 0.3 AU, to the Voyager and Pioneer explorations beyond 50 AU, spacecraft instruments have returned magnetic field data and plasma data (proton temperature, velocity, and density) that reveal the organized large-scale structure of the heliospheric plasma, along with transient mesoscale features such as coronal mass ejections and an ubiquitous but nonuniform admixture of fluctuations. Substantial fluctuation energy resides in an inferred range of spatial scales between the ion inertial scale (ഠ10 6 cm at 1 AU) and the observed correlation scale l (ഠ6 3 10 11 cm at 1 AU). The radial dependence of fluctuations in the low latitude solar wind is illustrated using Voyager 2 data in Figs. 1-3, which portray magnetic field variance (energy density in the turbulent magnetic field), correlation length, and proton temperature, from 1 AU to beyond 30 AU. To simultaneously explain these three data sets is a significant challenge. The main objective of this Letter is...