We model the contribution of the nearest young supernova remannt Vela to the local cosmic ray flux taking into account both the influence of the Local Superbubble and the effect of anisotropic diffusion. The dominant contribution of this source in the energy region around the cosmic ray knee can naturally explain the observed fluxes of individual groups of nuclei and their total flux. Adding the CR flux from a 2-3 Myr old local CR source suggested earlier, the CR spectra in the whole energy range between 200 GeV and the transition to extragalactic CRs are described well by the combined fluxes from these two local Galactic sources.Keywords: High energy cosmic rays, Vela supernova remnant, Local Superbubble, Galactic magnetic field.The measured energy spectrum of cosmic rays (CRs) extends smoothly over more than 11 decades as a nearly featureless power law, I(E) ∝ E −β . One of its most prominent features is the knee, a break in the all-particle energy spectrum at the energy E k 4 PeV, which was discovered by Kulikov and Khristiansen in the data of the MSU experiment already in 1958 [1]. The second knee corresponds to a change in the spectral slope of the all-particle energy spectrum at 5 × 10 17 eV where the slope hardens by ∆β 0.2. There is a general consensus that the knee in the total CR spectrum at E k 4 PeV coincides with a suppression of the primary proton and/or helium flux, and that the composition becomes increasingly heavier in the energy range between the knee and 10 17 eV [2][3][4][5].Explanations for the origin of the knee fall in two main categories, connecting it either with a change in the propagation or the injection of CRs. In the first case, the knee energy may either corresponds to the rigidity at which the CR Larmor radius R L is of the order of the coherence length l c of the turbulent magnetic field in the Galactic disk [6,7]. Alternatively, the knee corresponds to a transition between the dominance of pitch angle scattering to Hall diffusion or drift along the regular field [8][9][10]. In both cases, the energy dependence of the confinement time changes which in turn induces a steepening of the CR spectrum [6][7][8][9][10][11]. In the second class of models, the knee is connected to properties in the injection spectrum of the Galactic CR sources. For instance, the knee might correspond to the maximal rigidity to which CRs can be accelerated by the population of Galactic CR sources dominating the CR flux below PeV [12][13][14]. Alternatively, the knee may be caused by a break in the source CR energy spectrum at this rigidity [15,16]. A variant of this model is the suggestion that the spectrum below the knee is dominated by a single, nearby source and that the knee correspond to the maximal energy of this specific source [17,18]. All these models lead to a sequence of knees at ZE k , a behaviour first suggested by Peters [19].In the isotropic diffusion approximation one defines a scalar diffusion coefficient which depends on energy as D(E) = D 0 (E/E 0 ) δ . Measurements of the Boron and Carb...