Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: this is the second in a series of papers in which we study high - resolution images of bl lacertae made at 15 ghz with the vlba , under the mojave program ( monitoring of jets in active galactic nuclei with vlba experiments , lister et al . , 2009 ) . in ( * ? ? ? * hereafter paper i ) we investigated a quasi - stationary bright radio feature ( component ) in the jet located 0.26 mas from the core , ( 0.34 pc , projected ) and identified it as a recollimation shock ( rcs ) . numerous components appear to emanate from this shock , or pass through it . they propagate superluminally downstream , and their tracks cluster around an axis that connects the core and the rcs . this behavior is highly similar to the results of numerical modeling @xcite , in which mhd waves or shocks are emitted by an rcs . in the simulations , the jet has a magnetic field that dominates the dynamics , and is in the form of a helix with a high pitch angle , @xmath10 . in bl lac the motions of the components are similar to those in the numerical models , and in addition the electric vector position angle ( evpa ) is longitudinal ; i.e. , parallel to the jet axis . for a jet dominated by helical field , this indicates that the toroidal component is substantial ( @xmath11 ) , a necessary condition for the comparison of the observations with the numerical simulations . hence , in paper i , we assumed that the superluminal components in bl lac are compressions in the beam established by slow-  and/or fast-  mode magnetosonic waves or shocks traveling downstream on a helical field . it has been common to assume that the evpa is perpendicular to the projection of the magnetic field vector b that is in the synchrotron emission region . this is correct in the frame of an optically - thin emission region , but may well be incorrect in the frame of the observer if the beam is moving relativistically @xcite . lyutikov , et al . show that if the jet is cylindrical and not resolved transversely , and if the b field has a helical form , then the evpa will be either longitudinal or perpendicular to the jet , depending on the pitch angle . this is partly seen in the polarization survey results of @xcite , where the bl lac objects tend to have longitudinal evpa in the inner jet , whereas the quasars have a broad distribution of evpa , relative to the jet direction . this suggests that in bl lacs the field may be helical , with pitch angles large enough to produce longitudinal evpa , although strong transverse shocks in a largely tangled field are also a possibility ( e.g. * ? ? ? the wide distribution of evpa values in quasars suggests that oblique shocks , rather than helical structures , might dominate the field order . however , a distribution of helical pitch angles could also explain the evpas in quasars , if symmetry is broken between the near and far sides of the jet . it has been suggested @xcite that this difference in the magnetic field is fundamental to the generic differences between quasars and bl lacs and , by inference , between fanaroff & riley class ii and i sources , respectively @xcite . bl lacs often show a bend in the jet , and the literature contains examples showing that in some cases the evpa stays longitudinal around the bend ; e.g. , 1803 + 784 ; @xcite , 1749 + 701 ; @xcite , and bl lac itself ; @xcite . in these examples the fractional polarization @xmath12 rises smoothly along the jet to values as high as @xmath13 . the field must be well - ordered for the polarization to be that high . in this paper we assume that the field is in a rather tight helix ( in the beam frame ) and that the moving patterns ( the transvese disturbances ) are alfvn waves propagating along the longitudinal component of the field .    in a plasma dominated by the magnetic field , alfvn waves are transverse displacements of the field ( and , perforce , of the plasma ) , analogous to waves on a whip . the tension is provided by the magnetic field @xmath14 , and the wave velocity is proportional to the square root of the tension divided by the ( relativistic ) mass density . alfvn waves have been employed in various astronomical contexts , including the acceleration of cosmic rays @xcite , the solar wind @xcite , the jupiter - io system @xcite , turbulence in the ism @xcite , the bow shock of mars @xcite , and the solar atmosphere@xcite . in our case they are transverse waves on a relativistically - moving beam of plasma threaded with a helical magnetic field . the appropriate formulas for the phase speeds of the mhd waves are given in the appendix of paper i.    changes in the ridge lines of bl lacs are also seen frequently . @xcite showed that in 1.4 years the bl lac object 0735 + 178 changed from having a `` staircase '' structure to being straight , and that there were prominent transverse motions . @xcite also studied 1803 + 784 and described various models that might explain the structure . @xcite studied the ridge line in 0836 + 710 at several frequencies and over a range of epochs . they showed that the ridge line corresponds to the maximum pressure in the jet . they discussed the concept of transverse velocity , and concluded that their measured transverse motions are likely to be caused by a `` moving wave pattern '' ; this was elaborated in @xcite . in our work here on bl lac we also see transverse motions , but their patterns move longitudinally and we identify them as alfvn waves . we calculate the resulting transverse velocity of the wave motion and show that it is non - relativistic . it has been more customary to discuss the fast radio components in a relativistic jet in hydrodynamic ( hd ) terms . we note here only a few examples of this . the shock - in - jet model @xcite was used by hughes , aller , & aller ( 1989a , 1989b , 1991 ) to develop models of several sources , including bl lac @xcite and 3c279 @xcite . @xcite recognized two threads of emission in 3c273 that they explained with kelvin - helmholtz instabilities , and this was developed more by @xcite . @xcite discussed the patterns and motions in 3c120 in terms of helical instability modes . in all these studies the magnetic field is needed of course for the synchrotron radiation , but it also is explicitly used to explain observed polarization changes as due to compression of the transverse components of magnetic field by the hd shock . but in these works the magnetic field has no dynamical role in the jet . on the contrary , in this paper , as in paper i , we assume that the dynamics in the jet are dominated by the magnetic field . the plan for this paper is as follows . in section [ sec : observations ] we briefly describe the observations . the definition of the ridge line of a jet is considered in section  [ sec : ridgelines ] , and the transverse waves and their velocities , including the behavioral change in 2010 , are presented and discussed in section  [ sec : waves ] . excitation of the waves by changes in the p.a . of the rcs is considered in section  [ sec : excitation ] . in section  [ xverse_mhd ] we identify the waves as alfvn waves , discuss their properties , and present simple models of the system .    for bl lac @xmath15 , and the linear scale is @xmath16 . an apparent speed of @xmath17 corresponds to @xmath18 . for this study of bl lac we use 114 epochs of high - resolution observations made with the vlba at 15 ghz , between 1995.27 and 2012.98 . most of the observations ( 75/114 ) were made under the mojave program @xcite , a few were taken from our earlier 2-cm program on the vlba @xcite , and the rest were taken from the vlba archive . the data were all reduced by the mojave team , using standard calibration programs @xcite . following the reduction to fringe visibilities we calculated three main products at nearly every epoch :    ( 1 ) : :    an image , consisting of a large number of `` clean delta functions ''    produced by the algorithm used for deconvolution , convolved with a    `` median restoring beam '' , defined in section  [ sec : ridgelines ] . ( 2 ) : :    a model , consisting of a set of gaussian `` components '' found by    model - fitting in the visibility plane ; each component has a centroid ,    an ellipticity , a size ( fwhm ) , and a flux density . the gaussians are    circular when possible . the total set of components sums to the image ,    but in this paper we only use components that have been reliably    measured at four or more epochs , have flux density    @xmath19 , and can be tracked unambiguously from epoch    to epoch . a typical epoch shows 4 - 6 of these `` robust '' components .    the rcs is a permanent component and , together with the core , usually    produces more than half of the total flux density from the jet . the    centroids of the robust components for each epoch are plotted on the    images in figure  [ fig1 ] . +    the centroid locations are measured relative to the core , which we    take to be the bright spot at the north end of the source ; it usually    is regarded as the optically - thick ( @xmath20 ) region of    the jet . in principle , the core can move on the sky . we considered    this in paper i , and concluded that any motions are less than    @xmath21 in a few years , and they were ignored . our    positional accuracy is conservatively estimated as    @xmath22 , and in this paper we again ignore any    possible core motions . ( 3 ) : :    the ridge line , shown in figure 1 and discussed in    section  [ sec : ridgelines ] . the image , the components , and the ridge line are not independent , but each is advantageous when discussing different aspects of the source . in most cases the ridge line runs down the smallest gradient from the peak of the image , and the centroids of the components lie on the ridge line . however , when the jet has a sharp bend the algorithm can fail , as in figure  [ fig1]c . this is discussed in section  [ sec : ridgelines ] . the components move in a roughly radial direction , and plots of @xmath23 as well as the sky ( ra  dec ) tracks are shown in paper i and in @xcite . the tracks cluster around an axis at @xmath24 and appear to emanate from a strong quasi - stationary component , c7 , that we identified as a recollimation shock ( rcs ) in paper i. the moving components have superluminal speeds ; the fastest has @xmath25 in units of the speed of light . we are dealing with moving patterns on the jet of bl lac , and in order to quantify them we first need to define the ridge line of a jet . at least four definitions have been used previously . @xcite used the line that connects the components at a single epoch , in studying 1803 + 784 . @xcite investigated three methods of finding the ridge line : at each radius making a transverse gaussian fit and connecting the maxima of the fits , using the geometrical center , and using the line of maximum emission . they found no significant differences among these procedures , for the case they studied , 0836 + 710 . they showed that the intensity ridge line is a robust structure , and that it corresponds to the pressure maximum in the jet .    to quantify a ridge line we start with the image as in figure  [ fig1 ] , which is the convolution of the `` clean delta functions '' with a smoothing beam . since we are comparing ridge lines from different epochs , we have used a constant `` median beam '' for smoothing , and not the individual ( `` native '' ) smoothing beams . the latter vary a little according to the observing circumstances for each epoch , and their use would effectively introduce `` instrumental errors '' into the ridge lines . the median beam is a gaussian with @xmath26 and @xmath27 each of the three parameters is the median of the corresponding parameters for all the epochs . the algorithm for the ridge line starts at the core , and at successive steps ( 0.1 mas ) down the image finds the midpoint , where the integral of the intensity across the jet , along a circular arc centered on the core , is equal on the two sides of the arc . the successive midpoints are then smoothed with a third - order spline . ridge lines are shown on the three images in figure  [ fig1 ] . in figure [ fig1]a the bends in the jet are gradual and the algorithm works very well , as indeed would any of the methods mentioned above . in figure  [ fig1]b there are two sharp bends and our algorithm makes a smooth line that misses the corners of the bends . in this case connecting the components would be better , if the modelling procedure actually put components at the corners . in figure [ fig1]c the jet appears to bifurcate , and our algorithm picks the west track . in this case a visual inspection of the image is required to see what is going on . in fact there is another problem with figure  [ fig1]c . the image has a step to the east ( looking downstream ) about 1 mas from the core , where a short ew section connects two longer ns sections . since the restoring beam is nearly ns the details of this step can not be reconstructed . the calculated ridge line in figure  [ fig1]c does not reproduce the step , but makes a smooth track . figure  [ fig : ridge_all ] shows nearly all the ridge lines that we consider in this paper ; a few are not shown because they occur very close in time to another one . in all cases the rcs is located close to the semi - circle , drawn 0.25 mas from the core . successive panels are adjacent in time , although there is a 1-yr gap in the data between panels ( d ) and ( e ) . the only other substantial data gap is seen in panel ( a ) , from 1998.18 to 1999.04 . in figure [ fig : ridge_all ] the epochs are set nearly equally among the panels , with the separations picked to emphasize the various waves that are discussed below . it is important to establish the reliability of the ridge lines because our analysis rests on them , and some of the structures that we interpret as waves are smaller than the synthesized vlba beam . we first note that as with all vlbi our sampling of the @xmath28 plane is sparse , and different samplings can produce different ridge lines . to see how strong this effect is , we emulated an observation with missing antennas by analyzing a data set with and without one and two antennas , and we did this analysis both with the native restoring beams and the median restoring beam described above . the results for 2005 - 09 - 16 are shown in figure  [ fig : nativevsmedian ] ; they are similar to the results we obtained for two other epochs . in figure [ fig : nativevsmedian]a we show two ridge lines , the solid one is calculated with the full data set and the dashed line is obtained when data from the sc and hn antennas are omitted . the latter calculation does not use many of the baselines , including the longest ones . the chief effect is a shift of the pattern downstream , by roughly 0.1 mas . this shift is not a statistical effect , but is mainly due to the different smoothing beams that were used for the two cases . we found that the differences in the ridge lines increased with increasing difference in the p.a.s of the smoothing beams . in figure [ fig : nativevsmedian]a the difference in p.a of the smoothing beams is 17 .    in figure [ fig : nativevsmedian]b we used the median beam . in this case the curves are close with differences of typically @xmath29as out to 4 mas , where the surface brightness becomes low . beyond 4 mas the differences rise to @xmath30as . another way to investigate the reliability of the ridge lines is to examine pairs of ridge lines measured independently but close together in time . the full data set contains 10 pairs where the separation is no more than 10 days , and these are all shown in figure  [ fig : pairs ] . they are calculated with the median restoring beam . note that the bottom three panels have a different vertical scale than the others . in general the comparison is very good within 4 mas of the core . panel ( i ) contains one ridge line that stops at 3.6 mas because the brightness at the ridge becomes too low ; this limit also can be seen in a few places in the other figures . panel ( i ) contains the only pair that has a continuous offset , @xmath31 . these data were taken during an exceptional flux outburst at 15 ghz in bl lac , seen in the mojave data ( unpublished ) , and roughly coincident with outbursts seen at shorter wavelengths @xcite . an extra coreshift leading to a position offset is expected with such an event @xcite . in any event , this pair appears to be different from the others , and we do not include it in the statistics . figure  [ fig : histogram ] shows the histogram of separations between the paired ridge lines , after excluding those in panel ( i ) of figure  [ fig : pairs ] . in forming the ridge lines a 3-pixel smoothing was used , and for the histogram we have used every third point . the median separation is @xmath32 . thus the repeatability of the ridge lines is accurate to about @xmath32 . the reliability also depends on the effect discussed in connection with figure  [ fig1 ] , that the ridge - finding algorithm can smooth around a corner , and can be in error by perhaps @xmath33 . however , the error is roughly constant over short time spans , as in figure  [ fig : pairs ] panel ( e ) where the sharp bend at @xmath34 is smoothed the same in the two curves . this smoothing will have little effect on calculations of wave velocity , which is our main quantitative use of the ridge lines . we ignore the smoothing in this paper .    from this investigation we conclude that caution must be taken in interpreting the ridge lines , especially when comparing ridge lines obtained at different epochs , or with different frequencies . the details of the restoring beam can have a noticeable effect on the ridge line , and to avoid misinterpretation the restoring beam should be the same for all the ridge lines that are being intercompared . when considering these ridge lines it is important to keep the geometry in mind : the jet has a small angle to the line - of - sight ( los ) , and the foreshortening is about a factor of 10 ( paper i ) . also , the projected images in figure  [ fig1 ] can hide three - dimensional motions . to work with skew and non - planar disturbances , we use the coordinate systems shown in figure  [ fig : coords ] . east , north , and the los form the left - hand system ( x , y , z ) and the jet lies at angle @xmath35 from the los in the @xmath36 formed by the los and the mean jet axis . this plane is perpendicular to the sky plane and is at angle p.a . from the @xmath37 axis . the rotated system @xmath38 is used to describe transverse motions : @xmath39 is in the sagittal plane , @xmath40 is perpendicular to it , and @xmath41 is along the jet . by `` transverse motion '' we mean that a point on the beam has a motion in the @xmath42 plane : @xmath43 . the component @xmath44 lies in the sagittal plane and its projection on the sky is along the projection of the jet . this component therefore is not visible , although a bright feature moving in the @xmath39 direction might be seen as moving slowly along the jet . however , the @xmath45 component remains perpendicular to the los as @xmath35 or p.a . changes , and its full magnitude is always seen . thus a measured transverse motion is a lower limit . if the beam is relativistic then time compression of the forward motion must be added ; see section  [ sec : transverse ] .    ] some of the panels in figure  [ fig : ridge_all ] show disturbances that appear to move down the jet , and at other epochs the jet is fairly quiet . we now consider several of the disturbances in detail , starting with the structures seen in figure  [ fig : ridge_all ] , panel ( b ) . figure  [ fig : ridge9900 ] is an expanded view of figure  [ fig : ridge_all ] , panel ( b ) . it includes ridge lines for 14 consecutive epochs over a period of about 1.6 yr . beyond 1 mas the early epochs ( solid lines ) show the jet bending to the se . later epochs show the bend farther downstream , and at 2000.31 and later the jet bends to the sw before bending se . we anticipate a result from section  [ sec : velocity ] and draw vector a at @xmath46 across the tracks . the intersections of vector a with the tracks are shown in the inset in figure  [ fig : ridge9900 ] . the velocity implied by the line in the inset is close to @xmath17 or @xmath47 . the pattern on the ridge line is moving superluminally downstream at nearly constant velocity . we consider three possible explanations for this . * 1 ) * we see the projection of a conical pattern due to a ballistic flow from a swinging nozzle , like water from a hose . the argument against this is that line b in figure  [ fig : ridge9900 ] is parallel to vector a and approximately tangent to the western crest ; this feature of the ridge lines is not radial from the core as it would be if it were a ballistic flow . in figure [ fig : ridge_all ] all the panels except ( a ) , ( b ) , and ( e ) show clearly that the excursions of the ridge lines are constrained to lie in a cylinder , not a cone . * 2 ) * the moving pattern is due to a helical kink instability that is advected downstream with the flow . in the kink the field would be stretched out and become largely parallel to the observed bends in the jet that , in this case , seem to be transverse waves @xcite . this would produce an evpa normal to the wave crest in figure  [ fig : ridge9900 ] rather than longitudinal . but in bl lac the evpa tends to be longitudinal , even along the bends . in figure [ fig : pol ] we show the polarization image for 2005 - 09 - 23 , taken from the mojave website@xmath48 . similar polarization images for bl lac , at several wavelengths , are shown in osullivan & gabuzda ( 2009 , figure 19 ) for epoch 2006 - 07 - 02 . both of these epochs are part of the large wave d shown later in figure  [ fig : ridge0506 ] . in these polarization images the evpa is nearly parallel to the jet out to about 5 mas and @xmath12 is high on the ridge , indicating that the magnetic field remains in a relatively tightly - coiled helix around the bend and is not nearly parallel to the axis , as it should be for an advected kink instability . wave d is the largest wave in the bl lac data , and seems to have the cleanest longitudinal polarization . at other epochs the evpa tends to be longitudinal , but can be off by up to 40 . we have only one epoch of polarization data for wave a , but that one does show an evpa that is tightly longitudinal in the bend . thus we believe that the evpa results preclude the identification of the structures seen in figure  [ fig : ridge9900 ] as due to a kink instability .    * 3 ) * the moving patterns are transverse mhd waves ; i.e. , alfvn waves . for this to be possible the plasma must be dynamically dominated by a helical magnetic field . this condition for the jet of a bl lac has been suggested many times ; see e.g. , @xcite , @xcite . note that we implicitly assumed the helical , strong - field case in discussing the kink instability , in the preceding paragraph , and we also assumed it in paper i. thus , we assume that the moving pattern under vector a in figure  [ fig : ridge9900 ] is an alfvn wave , with velocity @xmath49 .    in figure [ fig : ridge9900 ] a second wave is seen between @xmath50 and @xmath51 mas , where the ridge lines for epochs 2000.31 and later bend to the sw . the two waves in figure  [ fig : ridge9900 ] can be thought of as one wave with a crest to the west . this wave is generated by a swing of the nozzle to the west followed by a swing back to the east about 2 years later , as discussed below in section  [ sec : excitation ] . the 1999 - 2000 wave is displayed in a different form in figure  [ fig : waves9901 ] , which shows the ridge lines from 1999.37 to 2001.97 . vertical spacing is proportional to epoch , and the axes have been rotated by 13 ; arrows at top show north and east . tick marks on the right are 0.1 mas apart . the dots show the points described later in section  [ sec : velocity ] , where the slope changes , and the solid line a is a linear fit through the points , with speed @xmath52 . this wave is prominent until 2000.99 . in 2001.22 the structure has changed . there are alternate possibilities to explain this new structure , b. it may be a new wave , with the crests connected with line b ( drawn with the same slope as line a ) . in this case the wave must have been excited somehow far from the rcs . the fit of line b to the wave crests is poor and would be improved if acceleration were included , but there is not enough data for that . alternatively , structure b may simply be a relic of the trailing side of wave a , perhaps relativistically boosted by the changing geometry ( the bend ) seen in figure  [ fig : ridge_all ] panel ( c ) . a third wave c is shown by the dashed line that again is drawn with the same slope .    . the three lines are parallel and all have slope 0.92 mas @xmath53 . [ fig : waves9901],scaledwidth=49.0% ]    panel ( c ) of figure  [ fig : ridge_all ] shows the ridge lines projected on the sky for 2001  2002 . wave b from figure  [ fig : waves9901 ] is seen as the bump to the east at @xmath51 mas , which moves downstream at succeeding epochs . the projected axis of the jet is curved at these epochs , and the possible acceleration noted above for wave b may simply be a relativistic effect inherent in the changing geometry . wave a in figure  [ fig : waves9901 ] is barely visible in figure  [ fig : ridge_all ] panel ( a ) as a gentle bump in 1999.04 , so it is first apparent in early 1999 at a distance @xmath54 mas from the core . this is reminiscent of the behavior of the components discussed in paper i ; figure 3 of that paper shows that most of the components first become visible near @xmath50 mas . wave c also appears to start near @xmath54 mas .    in figure [ fig : ridge9900 ] the short arrow c shows an eastward swing of the inner jet between 2000.01 and 2000.31 . this is seen in figure  [ fig : waves9901 ] in the ridge line for 2000.31 , which shows a new inner p.a . the effect of these p.a . swings on the beam is discussed in section  [ sec : excitation ] . the different panels in figure  [ fig : ridge_all ] show that the jet can be bent , and even when relatively straight , can lie at different p.a.s . hence there is no unique rotation angle for the ridge lines in a plot such as that in figure  [ fig : waves9901 ] . the rotation angle used in figure  [ fig : waves9901 ] was found by the velocity algorithm described in section  [ sec : velocity ] for wave a.    further examples of waves are shown in @xmath55 , omitting the extraneous ridge lines to avoid confusion . the wave motions are indicated by the arrows , which are propagation vectors derived in section  [ sec : velocity ] . table  [ tbl : waves ] lists the details for these waves . @xmath56 is the measured proper motion , @xmath0 is the apparent speed in units of @xmath1 , @xmath57 the wave speed in the coordinate frame of the galaxy , assuming @xmath58 , and p.a . is the projected direction of the propagation vector . the amplitude is an estimate of the projected distance ( in mas ) across the wave , perpendicular to the propagation vector . wave d is the largest such feature seen in the data . unfortunately , there was an 11-month data gap prior to 2005.71 , and the wave can not be seen at earlier times . the amplitudes of the larger waves appear to be comparable with the wavelength , as suggested for example by the inclination angle @xmath59 shown in figure  [ fig : ridge0506 ] : @xmath60 . but this is an illusion caused by the foreshortening , which is approximately a factor of 10 ( paper i ) , so the deprojected value of @xmath59 is about 5 . note that this is a lower limit , since the transverse motion can have a component in the @xmath61 plane in figure  [ fig : coords ] . figure  [ fig : movie ] contains one frame of a movie of bl lac showing the jet motions and ridge line fits at 15  ghz . the full movie is available in the electronic version of this paper . ccccccccc epoch & n & @xmath56 & & @xmath62&@xmath63 & p.a . & & amplitude + & & ( mas y@xmath64 ) & & & & ( deg ) & & ( mas ) + a 1999.37 - 2000.99&14 & 0.92 & @xmath65 & 3.9 & 0.979 & @xmath66 & @xmath67 & 0.5 + d 2005.71 - 2006.86 & 5 & 1.25 & @xmath68 & 5.6 & 0.987 & @xmath69 & @xmath70 & 0.9 + e 2008.33 - 2008.88 & 8 & 3.01 & @xmath71 & 13.5 & 0.998 & @xmath72 & @xmath73 & 0.3 + f 2009.33 - 2009.96 & 6 & 1.11 & @xmath74 & 5.0 & 0.985 & @xmath75 & @xmath76 & 0.2 +    notes . columns are as follows : ( 1 ) wave label , ( 2 ) inclusive range of epochs , ( 3 ) number of epochs , ( 4 ) apparent speed , ( 5 ) error , ( 6 ) apparent speed in units of @xmath1 , ( 7 ) speed in galaxy frame , assuming @xmath58 , ( 8) p.a . of the wave , ( 9 ) error , ( 10 ) estimated amplitude . mas , with a major axis position angle at @xmath77 , as indicated in the lower left corner of the frame . the fitted ridge line is shown as a dashed line in the image , and again as a solid line to the left of the image . these have also been linearly interpolated between the individual vlba epochs . the points of changing slope ( see section  [ sec : velocity ] ) at individual vlba epochs are shown as the small symbols . at left the ridge lines are shown with different colors for the various waves . the yellow @xmath78 line on the right is advancing at the speed of light ( @xmath79 ) and is included for reference . the entire movie is available in the web version of this paper . [ fig : movie],scaledwidth=47.0% ]      in figure  [ fig : ridge_all ] panels ( g ) and ( h ) we see that by 2010 the earlier transverse wave activity in the jet has subsided , and that after 2010.5 the jet is well - aligned at @xmath80 with a weak wiggle . but the wiggle is not stationary . figure  [ fig : ridge1013 ] shows the ridge lines plotted on axes rotated by 95 , and spaced proportionately to epoch . most of the ridge lines have a quasi - sinusoidal form . almost all the epochs show a negative peak in the inner jet , with a minimum near @xmath81 . this is a quasi - standing feature , of variable amplitude . at most epochs there is a positive peak near @xmath82 . this also is a quasi - standing feature , but less distinct than the inner one .    what is causing the quasi - standing features ? the patterns can hardly be true standing waves because that requires a reflection region . a rotating helix would project as a traveling wave , as on a barber pole , so a simple barber - pole model is excluded . possible motions of the core are only about @xmath21 ( paper i ) , so any registration errors due to core motion are much smaller than the observed changes , which are up to @xmath33 . there is little indication of wave motion in figure  [ fig : ridge1013 ] , at least not at the speeds seen in figure  [ fig : ridge_all ] . it appears then , that during the period 2010 - 2013 , the jet was essentially straight but with a set of weak quasi - stationary patterns , with variable amplitude . we estimated the velocity of wave a in figure  [ fig : ridge9900 ] in two independent ways . in the first we assume that there is a constant propagation vector , and we shift and superpose the ridge lines on a grid of ( @xmath56 , p.a . ) where @xmath56 is the speed of the wave and p.a . is its propagation direction . if the ridge lines form a simple wave , then the solution is found when the lines lie on top of each other . this is shown in figure  [ fig : superpose ] , where a reasonable fit can be selected by eye . the result is @xmath83 at @xmath84 . this solution is somewhat subjective and the quoted errors do not have the usual statistical significance .    as an alternative procedure to visually aligning the ridge lines , we developed a method of identifying a characteristic point on the wave , just downstream of the crest , where the wave amplitude has begun to decrease . define the slope of the ridge line as @xmath85 in pixels , where in figure  [ fig : waves9901 ] , @xmath86 and @xmath37 are rotated ra and dec , and take the first downstream location where the slope exceeds @xmath87 . this point is marked with the dot @xmath88 on the ridge line for 2000.57 in figure  [ fig : ridge9900 ] . the @xmath86 and @xmath37 positions vs time for these locations are then fit independently using the same methods as described in @xcite to extract a vector proper motion for this characteristic point on the wave . the two methods agree well and the analytic solution is @xmath89 at @xmath90 , and the apparent speed is @xmath91 . the propagation vector is shown in figure  [ fig : ridge9900 ] and the speed and direction of the wave are listed in table  [ tbl : waves ] . the table also includes @xmath63 the speed of the wave in the galaxy frame , assuming @xmath58 . this calculation assumes that the ridge lines lie in a plane ; i.e. , are not twisted . this is not neccessarily the case . rather , since the inner jet , near the accretion disk , may wobble in 3 dimensions , ( mckinney et al , 2013 ) it seems likely that the rcs may execute 3-dimensional motion and that the downstream jet will also . see section  [ sec : excitation ] . note that the p.a . of the first and last propagation vectors in table  [ tbl : waves ] ( -1670 , -1671 ) is the same ( to within the uncertainties ) as the p.a . of the axis ( -1666 ) defined in paper i as the line connecting the core with the mean position of the recollimation shock . in the context to be developed later , the jet acts as a whip being shaken rapidly at the rcs , and tension in the whip continually pulls it towards the mean pa .    in table  [ tbl : waves ] the speeds for the first , second and fourth waves are all similar at @xmath92 , but wave e ( 2008 ) is 3 times faster . wave e has @xmath93 , which is comparable to the speed for the fastest component in bl lac , @xmath94 , although the components speeds vary widely , from @xmath95 to @xmath96 @xcite . wave e is also distinguished by its polarization ; the evpa is transverse not longitudinal like the others . we defer further discussion of wave e to another paper . the ridge waves are relativistic transverse waves with apparent speeds @xmath0 from 3.9 to 13.5 times the speed of light , and we assume that they have a small amplitude . from the usual formula for apparent speed , @xmath97 and taking values of @xmath98 from table  [ tbl : waves ] and using @xmath58 , we find @xmath99 for the speed of the waves in the frame of the host galaxy . we now discuss the jet motion in terms of the coordinate system @xmath100 shown in figure  [ fig : coords ] . consider a transverse motion that is in the ( @xmath101 ) plane . let the beam contain a co - moving beacon that is at the origin and emits a pulse at time @xmath102 , where @xmath103 is in the coordinate frame of the galaxy . when @xmath104 yr the signal from the origin will have traveled @xmath105 down the z axis , towards the observer . also at @xmath104 the beacon has moved from the origin to the point @xmath106 where @xmath107 is the transverse speed , and @xmath108 is the longitudinal speed of the beam , both in the frame of the galaxy . at this point the beacon emits a second signal that also travels at the speed of light . in the z - direction , this signal trails the first one by ( @xmath109 ) years . the apparent transverse speed of the beacon in the direction perpendicular to the jet , in the galaxy frame , is then @xmath110 and is to be differentiated from the apparent speed @xmath111 commonly used in studies of superluminal motion , which is the apparent speed _ along _ the jet . note the close relation between equations  [ vapp ] and  [ vtrans ] . equation  [ vtrans ] can be inverted to find @xmath112 , a lower limit to the transverse speed .    for wave a in figure  [ fig : ridge9900 ] we obtain an estimate for the transverse speed at @xmath113  mas by taking the transverse motion as 0.5 mas and the time interval as @xmath114 yr , giving @xmath115 and @xmath116 and , from equation  [ vtrans ] with @xmath58 and @xmath117 ( paper i ) , @xmath118 . this is a model - dependent rough value , but it shows that the transverse speed is non - relativistic . this is necessary for consistency , since the derivation of the relativistic form of the mhd wave speeds shown in paper i assumes that the velocity perturbation is small . we suggested in paper i that component 7 is a recollimation shock , and that the fast components emanate from it . if this is correct , then the rcs should be a nozzle and its orientation should dictate the direction of the jet . in this section we investigate this possibility . we first note that it is not possible to make a detailed mapping between the p.a . of the rcs and the later wave shape , for two reasons . first , the algorithm for the ridge line smooths over 3 pixels ( 0.3 mas ) , and thus smooths over any sharp features in the advected pattern . the second reason is more speculative . our conjecture is that the wave is launched by plasma flowing through the nozzle and moving close to ballistically until its direction is changed by a swing in the p.a . of the nozzle . but magnetic tension in the jet continually pulls it towards the axis , and this means that it will bend , and that small - scale features will be stretched out and made smooth . we start by comparing the p.a . of the rcs with the p.a . of the downstream ridge line at @xmath50 mas . figure  [ fig : pa ] shows the p.a . of the rcs measured at 15 ghz and at 43 ghz . the latter is calculated from data kindly provided by the boston university vlbi group . we used the result found in paper i , that the 15 ghz core is a blend of the first two 43 ghz components and that the 15 ghz component 7 is the rcs , as is the third 43 ghz component . we calculated the centroid of the first two 43 ghz components , to find an approximate position for the 15 ghz core , and then calculated the p.a . of the 43 ghz rcs from that centroid . the result is shown in figure  [ fig : pa ] . we eliminated one discrepant point at 43 ghz , which was separated by about 20  from nearby 43 ghz points , and one discrepant point at 15 ghz . the correspondence between the two frequencies is generally good , especially after 2005.0 where the agreement is typically within 3 . this further justifies our claim ( paper i ) that the the location of this component is independent of frequency , and that it is a recollimation shock . figure  [ fig : pa ] also contains the p.a . of the 15 ghz ridge line , close to @xmath119 1.0 mas . between 2005.0 and 2010.0 the ridge line p.a . lags the rcs pa , by roughly 0.6 to 1.5 yr . after 2010 the pa of both the rcs and the ridge line stabilizes , and the subsequent variations , with rms amplitude about 3 , may mainly be noise . prior to 2005.0 the variations are faster and more frequent and the lag is erratic . in places there appears to be no lag , but around 2000.0 and again around 2004.0 it is about 0.5 yr . thus it appears that the swinging in pa of the rcs is coupled to the transverse motions of the ridge line . when the rcs is swinging rapidly and strongly , as before 2005 , then so also is the ridge at 1 mas , with an irregular lag in p.a . that sometimes is about a half a year , and at other times is negligible . but when the rcs is swinging more slowly , as after 2005 , then the ridge at 1 mas is also swinging slowly , with a lag of about a year , and after 2010.0 they both are stable , with only small motions that may be dominated by measurement errors .    we suggest that the large transverse waves on the ridge are excited by the swinging in p.a . of the rcs . consider wave a , seen in figure  [ fig : ridge9900 ] . its crest lies near line b and moves downstream at @xmath120 . in 1999.37 the crest is at about @xmath121 mas and at @xmath120 would have been at the rcs ( @xmath122 mas ) around 1998.3 . this is in a data gap at 15 ghz , but at 43 ghz there was a peak in p.a . in mid- or late-1998 . given that in 1999 the time lag between the rcs and the ridge at 1 mas apparently was much less than 1 yr , the association between the peak in the rcs p.a . in 1998 and the crest of wave a is plausible . the fall in p.a . in 1999 and 2000 is seen as the short arrow c in figure  [ fig : ridge9900 ] , and it corresponds to the upstream side of wave a. the downstream side is the advected rise in p.a . of the rcs from mid-1997 to the peak in mid - or late 1998 . the p.a . of the rcs fell from mid-1996 to mid-1997 , and we might expect to find a corresponding crest to the east on wave a , about 1 mas downstream of the main crest to the west . in fact several of the earliest ridge lines in figure  [ fig : ridge9900 ] do show a minor crest to the east at about @xmath123 mas , which is 2 mas or 2 years at 0.92 mas y@xmath64 , downstream of the main crest to the west . a substantial acceleration in the wave speed would be needed for this to match . in any event , we can not speculate usefully on this because it takes place beyond 3 mas , where there is a general bend to the east at all epochs . we conclude that a plausible association can be made between the large swing west then east of the rcs between 1998.0 and 2000.1 , and wave a that is later seen on the ridgeline . a similar connection can be made for wave d , seen in figure  [ fig : ridge0506 ] in 20052006 . it can plausibly be attributed to the large swing of the rcs to the east that began in 2004 and continued into 2005 . this wave does not have a crest as wave a does , but a crude analysis can be made as follows . assume that point @xmath124 on the 2005.71 ridge line is the advected beginning of the wave . with a speed of 1.25 mas y@xmath64 ( table  [ tbl : waves ] ) this means that the swing to the east began around 2003.5 . this date is indicated on the abscissa in figure  [ fig : pa ] . apart from one high point at 2004.1 the p.a . of the rcs falls gradually from 2003.1 until late 2004 , when it must fall abruptly to meet the first point after the data gap in 2005 . this also is seen in figure  [ fig : ridge0506 ] ; the first four epochs have ridge lines that lie together and are straight at @xmath125 out to @xmath126 mas . this is consistent with the rcs p.a . being stationary from roughly 2004.7 to 2005.7 at @xmath127 . this is in a data gap , and this analysis suggests that the rcs p.a . was @xmath128 during most or all of the gap . we conclude that the large swing in p.a . of the rcs from mid-2003 until mid-2005 generated wave d , the largest wave in our data set .    the p.a . of the rcs rose rapidly from 2005.7 to about 2006.5 , but the p.a . of the ridge line rose more slowly , and not as far . from 2007.0 to 2009.0 the p.a . of the inner jet was roughly constant at about -170 , while the p.a . of the rcs slowly dropped to the same value . we do not have a straightforward interpretation of this behavior . we also see in figure  [ fig : ridge_all ] panels ( e ) , ( f ) and ( g ) that following the passage of wave d the jet slowly straightened out . the p.a . of the inner jet @xmath129 propagated as a low - amplitude wave , at roughly the same speed as the large waves , @xmath130 .    as a further complication , during this slow straightening out of the jet we see two more low - amplitude waves . the high - speed wave e ( figure  [ fig : ridge0808 ] ) has no obvious antecedent in the p.a . of the rcs . wave f ( figure  [ fig : ridge0909 ] ) is seen a year after wave e , at the `` usual '' speed of 1.1 mas y@xmath64 . these waves together make a complex set of possibly twisted ridge lines , seen together in figure  [ fig : ridge_all ] panel ( f ) .    in section  [ sec : jet_1013 ] we showed that the waves on the jet subsided in 2010 , and in 2010 - 2012 the ridge line had only a weak variable wiggle . during this time the p.a . of the rcs was essentially constant ; the variations seen in figure  [ fig : pa ] may represent the errors in the measurements , which would be about @xmath131 . these variations in space and time have some regularities , as discussed in section  [ sec : jet_1013 ] , but they do not appear to have a connection to the p.a . of the rcs . in paper i we saw that the component tracks all appear to come from or go through the rcs ( component 7 ) and that they lie in a window centered on @xmath132 @xcite . this now is understood in terms of the waves on the ridge lines , since the components all lie on a ridge . the jet is analagous to a whip with a fixed mean axis being shaken with small amplitudes , in various transverse directions . the whip will occupy a narrow cylinder centered on the axis , and in projection the cylinder becomes our window . in paper i we showed that the magnetic field in the jet of bl lac has a strong transverse component . we assumed that it has a helical form , and that it is likely that the field dominates the dynamics in the jet . this is the condition for the existence of mhd waves that propagate down the jet . we suggested that the moving synchrotron - emitting components are compressions set up by fast and/or slow magnetosonic waves , possibly shocks . now we introduce the third branch of mhd waves in the jet plasma , the alfvn wave , which is a transverse s ( shear ) wave , with the disturbance occurring normal to the propagation direction . in section  [ sec : waves ] we showed that the moving patterns on the jet are transverse waves , and now we suggest that they are alfvn waves . the phase speed of a transverse alfvn wave is given by @xmath133 where @xmath134 is the relativistic scalar alfvn speed , given in equation a6 of paper i , and @xmath135 is the angle between the propagation direction and the magnetic field . the alfvn wave has similar propagation properties ( with respect to the magnetic field direction ) as the slow wave ; i.e. it moves along the field , but not at all normal to it ( @xmath136 ) . note that alfvn waves generally will not produce shocks in an ideal mhd plasma . we now discuss these waves in the jet and present simple models that allow us to estimate the pitch angle @xmath10 of the helix , which we define as the angle between the axis of the helix and the direction of the magnetic field when projected onto that axis .    a simple relation exists for the relativistic phase speeds of the three mhd waves : @xmath137 where @xmath138 is the sound speed ( relative to the speed of light ) , and @xmath139 are the fast , slow , and transverse mhd wave speeds . equation  [ relation ] may be readily verified from equation  [ va ] combined with equations a1 and a2 of paper i. with this result , the three equations for the phase speeds , together with the definitions of the cusp and magnetosonic speeds in equations a3 and a4 in paper i , can be solved for the magnetosonic and alfvn speeds : @xmath140 @xmath141 finally , the propagation angle to the magnetic field @xmath135 can be found from equations  [ va ] and [ alfven ] .    in dealing with this system of equations we are helped with constraints on the mhd wave speeds : @xmath142 , also @xmath143 for an adiabatic sound wave in a relativistic gas . in addition , we adopt a constraint from the one - sidedness of bl lac , @xmath144 , where @xmath145 is the lorentz factor of the beam in the frame of the galaxy ; this gives a jet / counterjet intensity ratio of about @xmath146 for @xmath147 and a spectral index of @xmath148 @xcite . we assume that the three waves travel downstream in the beam frame and parallel to the jet axis . therefore , the propagation angle of all three waves is the pitch angle of the helix itself : @xmath149 . we do not , in fact , measure the wave speeds themselves but rather their apparent speeds in the frame of the galaxy . to relate these to their speeds in the beam frame we first use equation  [ vapp ] and then the relativistic subtraction formula @xmath150 where the superscripts define the coordinate frame . we now have 5 input quantities to the calculation : @xmath151 , @xmath35 and @xmath145 , and with them we can calculate @xmath152 , @xmath153 , @xmath154 , @xmath10 , and the magnetosonic mach number defined as @xmath155 , where @xmath156 is the magnitude of the spatial component of the four - velocity and @xmath157 is the lorentz factor .    to illustrate the relationships among the various waves we show in figure  [ fig : banana ] ( the `` banana diagram '' ) the results for the specific configuration @xmath158 . these values correspond to the fastest superluminal component in bl lac ( paper i ) and to the apparent speeds of the transverse waves noted in section  [ sec : waves ] above . the diagram contains quantities defined in the frame of the beam : sound speed and alfvn speed , @xmath10 the pitch angle of the helix , and @xmath159 the magnetosonic mach number . the diagram is bounded at the left and bottom by @xmath160 and @xmath161 . at the top , for @xmath162 , the boundary traces the curve @xmath163 , but for @xmath164 ( in this case ) , this curve sometimes ventures into a region where there are no solutions for @xmath10 . this region can be eliminated from the banana by continuing the curve for @xmath165 with one that satisfies the criterion @xmath166 at constant @xmath167 , as we have done here . inside the banana our conditions for magnetic dominance @xmath168 and @xmath169 are satisfied everywhere except in a thin quasi - horizontal region at top right , and in a thin quasi - vertical region at left . at the cusp at right @xmath170 , indicating a purely toroidal field and no propagating alfvn waves , regardless of the value of @xmath152 . the banana diagram is set on the plane defined by the lorentz factor of the beam and the apparent ( superluminal ) speed of the slow mhd wave , both measured in the galaxy frame . the location of the banana on this plane is set by the specific set of input parameters as on the top left . figure  [ fig : banana ] shows that knowing the apparent speeds of the three mhd waves and the angle @xmath35 of the jet to the line of sight is not enough to completely determine the jet properties . we must either determine one more quantity or make an assumption about the jet system . we will make two different assumptions for the sound speed , each yielding a simple model . the cases are first , a cold jet , in which the plasma sound speed is negligible ; and the other assumes that @xmath171 .    * model ( a ) : cold plasma *    in paper i we investigated a model of the jet in which an observed slowly moving component with @xmath172 is due to a slow magnetosonic wave whose speed , relative to the jet plasma , is negligible : @xmath173 . this means that the plasma is cold and @xmath174 ( eqn . 4 ) . in this case the apparent slow component speed is the beam speed itself . with this speed for the beam , we then assumed that a fast component was due to a fast magnetosonic wave , and , from the observed apparent speed , we were able to deduce its speed on the jet . this model can be placed in figure  [ fig : banana ] . the model uses @xmath175 , and is located at the dot marked `` * * a * * '' on the boundary of the diagram at @xmath176 . with @xmath177 and @xmath178 , the fast pattern speed is three times greater than the speed of the beam , when the speeds are measured by their lorentz factors . because we now also have a measurement of the apparent transverse alfvn wave propagation speed ( @xmath179 , a typical value from table 1 ) , we can extend this model to include computation of the total alfvn speed @xmath154 , the magnetosonic speed @xmath153 , and the magnetic field pitch angle @xmath10 . with @xmath180 negligible , in the galaxy frame we again have @xmath181 , @xmath182 , and now @xmath183 . then , using equation  [ beta_sub ] , these become in the frame of the beam @xmath184 , @xmath185 , and @xmath186 , yielding @xmath187  a moderate helical magnetic field . since @xmath188 when @xmath189 , we also can calculate the magnetosonic mach number defined in equation  [ beta_ms ] . this yields @xmath190 and qualifies this model as a trans - magnetosonic jet . * model ( b ) : hot plasma *    the plasma hardly can be cold as in model ( a ) because the source is a powerful synchrotron emitter and the electron temperature is probably of order 100 mev ; the electron component of the plasma therefore is probably relativistic . on the other hand , the sound speed may or may not be near @xmath191 , depending on how heavily the plasma is contaminated with heavy , non - relativistic ions . for lack of further information , we choose @xmath192 . but , as seen in figure  [ fig : banana ] , we still need another parameter to establish the solution . we took @xmath193 for model ( a ) to match a slow superluminal component , but if we do that now with @xmath192 it yields @xmath194 and @xmath195 . this value for @xmath196 is less than that typically found in radio beaming studies where @xmath197 @xcite . furthermore , the pitch angle @xmath198 is less than the one estimated from polarization analyses , @xmath165 . ( homan , private communication ) .    to reconcile these values we drop the assumption that the superluminal component with @xmath199 is a slow mhd wave propagating downstream ; it might for example be a reverse mhd shock or wave traveling upstream in the beam frame and seen moving slowly downstream in the galaxy frame . instead , we choose @xmath5 , because it yields acceptable values for @xmath196 and @xmath10 , and matches the speed of a number of superluminal components . the final solution , seen at point @xmath88 in figure  [ fig : banana ] , contains three quantities that are chosen to match observations , @xmath58 , @xmath200 , @xmath201 , and two quantities picked because they are plausible and give reasonable results , @xmath5 , and @xmath202 . the derived quantities are @xmath203 , @xmath204 , @xmath205 , @xmath206 , @xmath207 , @xmath208 , and @xmath209 . the slow and transverse mhd waves are non - relativistic ( @xmath210 ) in the frame of the beam , but the fast wave is mildly relativistic , with @xmath211 . note that relativistic addition to produce the observed speed is non - linear ; @xmath211 plus @xmath212 gives @xmath213 . this is discussed in paper  i.    with @xmath203 and @xmath214 we now calculate the doppler factor @xmath215 , which agrees closely with values in the literature @xcite . also , @xmath216 agrees with estimates from polarization analyses . thus we see that the model for bl lac with alfvn waves on a helical magnetic field is able to explain the moving transverse patterns on the jet of bl lac . it implies a modest lorentz factor for the actual plasma flow ( @xmath217 ) and explains the faster propagation of the components and the transverse disturbances as mhd sonic and alfvn waves , respectively . they are generated primarily at the site of the recollimation shock and propagate downstream on the helical field , each with a speed in the galaxy frame that is the relativistic sum of the wave speed in the beam frame and the plasma flow speed in the galaxy frame . the new model does , however , have a disadvantage . the magnetosonic mach number of 4.7 is rather high and in conflict with the original discussion in paper i on the generation of a collimation shock like c7 : the super - magnetosonic flow emanating from the black hole region should transition to a trans - magnetosonic flow ( @xmath218 ) after it passes through c7 . however , we see in figure  [ fig : banana ] that models with both low mach number and high pitch angle are mutually exclusive : trans - magnetosonic models have @xmath219 , and models with @xmath220 have magnetosonic mach numbers of @xmath221 or more . in order to obtain a more sophisticated model that is compatible both with generation of a post - collimation - shock flow and the polarization observations , it is likely that one or more of the rather restrictive model assumptions in this paper will have to be relaxed . further insight into the propagation of an alfvn wave on a jet can be gained by examining the group velocity , which has only one value , @xmath222 , and is always directed along the magnetic field @xcite . an isolated wave packet will spiral down the jet along the helical magnetic field . a uniform disturbance across the jet will produce a ripple that moves along all the field lines ; i.e. across the jet . the net result is a jump or bend that propagates downstream with speed proportional to the cosine of the pitch angle . this has a close analogy to a transverse mechanical wave on a coiled spring , or slinky . in both cases there is longitudinal tension , provided for the jet by the magnetic field . figure  [ fig : polars ] shows relativistic phase polar diagrams for the two models discussed above and identified in figure  [ fig : banana ] . the diagrams show mhd wave phase speeds in 3-dimensional velocity space with the origin of each at the center of the diagram . each diagram was computed using the relativistic equations a1-a6 in paper i. all surfaces are axisymmetric about the horizontal magnetic axis . in each panel the dotted , solid , and broken lines show respectively the speed - of - light sphere ( unity in all directions ) , the two compressional mhd wave surfaces ( fast [ @xmath223 and slow [ @xmath167 ] ) , and the transverse alfvn wave surface ( @xmath225 ) . unlike the speed of light , the speeds of the mhd waves depend on the polar angle @xmath135 between the propagation and field directions . all three mhd modes are labeled in the left half of the diagrams . the arrows labeled in the right half of the diagrams show the three characteristic wave speeds : sound ( @xmath152 ) , alfvn ( @xmath154 ) , and magnetosonic ( @xmath153 ) , which values are realized along the field for the slow and alfvn modes and normal to the field for the fast mode . as mentioned earlier , the slow and alfvn waves can propagate skew to the field , but not normal to it . some of the relationships among the three types of mhd waves can be seen in figure  [ fig : polars]b . the outer solid loop traces the fast magnetosonic mode , whose speed is a maximum @xmath153 at @xmath226 , and is the same as that of the alfvn wave ( dashed loop ) , @xmath154 , when @xmath227 , provided @xmath168 , where @xmath152 is the sound speed in the plasma . the propagation speed of the alfvn wave is proportional to @xmath228 and this also is approximately true for the slow magnetosonic wave , the inner loop . so far we have been discussing phase polar diagrams in a uniform magnetic field , and now address how this applies to a plasma jet with a helical field . figure  [ fig : helicalfield ] shows a schematic diagram of a helical field jet with the properties of model ( b ) discussed above and in figures  [ fig : banana ] and [ fig : polars ] . the helical field will have a pitch angle of @xmath229 , so the polar diagram in figure  [ fig : polars ] will be rotated by that amount . the propagation direction of the mhd waves points downstream in our model , allowing us to read off the values of their propagation speeds from the polar diagram : @xmath230 , @xmath231 , and @xmath232 . if the helical field and plasma properties are uniform along the jet , the results will be the same everywhere , producing mhd waves with uniform velocities . b is rotated by the angle @xmath233 to align the dashed arrow with the helical field direction on the near side of the jet ( left ) . the propagation speeds of the three mhd waves along the jet axis then can be directly read off the polar diagram ( @xmath230 , @xmath234 , @xmath232 ) . for a uniform helical field one obtains the same results at any point ( e.g. , on the far side of the jet at right . ) [ fig : helicalfield],scaledwidth=50.0% ]    however , there will be a longitudinal current that will cause the field strength and pitch angle to be functions of the radial coordinate @xmath235 . ( see the cut - away view of a plasma rope in figure 6.14 of @xcite , for a simple view of the radial variations in b and @xmath10 . ) but , there should be a cylindrical shell around the axis , covering a modest range of @xmath235 , in which the synchrotron emissivity into the direction of the observer is maximized . we assume that this shell is the dominant region and that the field strength and pitch angle there are the effective values that control the dynamics . ( see @xcite for a discussion of this point . thus , if this is the case , then our dynamical analysis of the waves and a polarization analysis of the emission should result in similar magnetic pitch angle estimates for the magnetic field . a preliminary polarization analysis ( using methods similar to those in @xcite and to be discussed in more detail in the next paper in this series ) produces estimates of at least 60- 70for the pitch angle . this is in agreement with our result for the hot model ( b ) , which gives @xmath236 . the jet of bl lac is highly variable and displays transverse patterns that propagate superluminally downstream on the ridge line . they are not ballistic , like water from a hose , but are constrained , like waves on a whip . the magnetic field is well - ordered with a strong transverse component that we assume is the toroidal part of a helical field . in @xcite we assumed that the helical field provided support for fast- and slow - mode mhd waves whose compressions we see as the superluminal components . we here assume that the moving transverse patterns are alfvn waves propagating on the longitudinal component of the magnetic field . the full set of ridge lines is shown in figure  [ fig : ridge_all ] , and we show six examples of the alfvn waves in figures 7 , 9 , 10 , 11 , and 12 . a movie ( figure  [ fig : movie ] ) provides assistance in studying the motions . the transverse wave activity died down in 2010 and the jet settled to a fixed position angle ( p.a . ) , with a mild wiggle . this wiggle was not stationary , but appeared to oscillate transversely , with amplitude about 0.4 mas . this mild wiggle persisted through the remaining data period , up to 2013.0 . although the transverse propagating alfvn waves were greatly reduced in 2010 - 2013 , the superluminal components , which we identified in paper i as mhd sonic waves , continued roughly as before . figure 2 in paper i shows that during this period they continued with about the same frequency and speed as earlier . furthermore , during the latter half of this period , from about 2011.4 to 2013.0 , bl lac was exceptionally active at shorter wavelengths @xcite , from 1 mm through gamma - rays . this general behavior can fit into our model . we have magnetosonic waves responsible for the superluminal components , and alfvn waves responsible for the moving transverse patterns . these are independent mhd modes , and can be separately excited . we suspect , however , that the increase in short - wavelength activity during the same period as the reduction in alfvn waves ( 20102013 ) is not a coincidence . the velocity of the transverse waves was established by finding characteristic points on the ridge lines where the slope changes , as well as by visual inspection of the delayed superposition of the ridge lines . three of the apparent velocities are near @xmath237 , and one is much faster , with @xmath238 . with @xmath58 and @xmath239 the speeds in the galaxy frame are approximately @xmath240 and in the beam frame @xmath241 . an alfvn wave displaces the jet in the transverse direction , and the observed motion can be converted into a transverse speed . for wave d , the largest wave we observed , the transverse speed , in the galaxy frame , is @xmath242 . this is a rough estimate but safely non - relativistic , and consistent with our assumption that the waves have low amplitude . the timing and direction of some of the the waves are correlated with the p.a . of the recollimation shock ( rcs ) , which swings over 25  in an irregular fashion . it appears that the waves are excited by the swinging of the rcs . this is analogous to exciting a wave on a whip by shaking it . in paper i ( figure 3 ) we saw that the ridge lines occupy a cylinder about 0.7 mas wide and 3 mas long , or 3 ly wide and 120 ly long when a deprojection factor of 10 is used . ( see also @xcite figure 1 . ) we now understand that this cylinder is formed by the transverse waves , whose axes generally are close to the source axis at @xmath243 . the width is set by the amplitude of the largest waves while the length is set by the general bend of the source to the se . we briefly describe the alfvn waves , and provide a method for calculating physical quantities in the jet in terms of the measured wave speeds . we investigate two simple models of the system ; in the first the plasma is cold and the sound speed @xmath161 . this gives results for the lorentz factor and the pitch angle that are in moderate disagreement with results from observations . the second model uses a hot plasma with @xmath244 , and assumes that the slow magnetosonic wave has apparent speed @xmath5 . this yields @xmath245 , pitch angle @xmath236 , alfvn speed @xmath246 , and magnetosonic mach number @xmath247 . this describes a plasma in which the helical magnetic field is strong with a dominant toroidal component .    in our model the lorentz factor for the beam is approximately 4.5 and is smaller than the observed apparent speed of most of the transverse waves as well as the fast superluminal components discussed in paper i. another way to say this is that in most cases the pattern speed is greater than the beam speed . this comes about because the pattern traces a wave traveling downstream on the beam . we conclude that the rapid movements of the transverse patterns in the jet of bl lac can be described as alfvn waves excited at the rcs and propagating downstream on the longitudinal component of a helical magnetic field . the jet can be described as a relativistic , rapidly shaken whip . we suggest that other similar sources be investigated with these ideas in mind . we thank the referee for comments that have improved the manuscript , m. perucho for reading the manuscript and offering helpful suggestions , and the mojave team for comments on the manuscript , and for years of work in producing the data base that makes this work possible . tga acknowledges support by dfg project number os 177/2 - 1 . th was partly supported by the jenny and antti wihuri foundation and by the academy of finland project number 267324 ; ts was partly supported by the academy of finland project 274477 . yyk is partly supported by the russian foundation for basic research ( project 13 - 02 - 12103 ) , research program ofn-17 of the division of physics , russian academy of sciences , and the dynasty foundation . abp was supported by the `` non - stationary processes in the universe '' program of the presidium of the russian academy of sciences . the vlba is a facility of the national radio astronomy observatory , a facility of the national science foundation that is operated under cooperative agreement with associated universities , inc . the mojave program is supported under nasa - fermi grant nnx12a087 g . this study makes use of 43 ghz vlba data from the vlba - bu blazar monitoring program ( vlba - bu - blazar ; http://www.bu.edu/blazars/vlbaproject.html ) , funded by nasa through the fermi guest investigator program . part of this research was carried out at the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration . this research has made use of nasa s astrophysics data system .
Abstract:

Answer: we study the kinematics of ridge lines on the pc - scale jet of the active galactic nucleus bl lac . we show that the ridge lines display transverse patterns that move superluminally downstream , and that the moving patterns are analogous to waves on a whip . their apparent speeds @xmath0 ( units of @xmath1 ) range from 3.9 to 13.5 , corresponding to @xmath2 in the galaxy frame . we show that the magnetic field in the jet is well - ordered with a strong transverse component , and assume that it is helical and that the transverse patterns are alfvn waves propagating downstream on the longitudinal component of the magnetic field . the wave - induced transverse speed of the jet is non - relativistic ( @xmath3 ) . in 2010 the wave activity subsided and the jet then displayed a mild wiggle that had a complex oscillatory behaviour . the alfvn waves appear to be excited by changes in the position angle of the recollimation shock , in analogy to exciting a wave on a whip by shaking the handle . a simple model of the system with plasma sound speed @xmath4 and apparent speed of a slow mhd wave @xmath5 yields lorentz factor of the beam @xmath6 , pitch angle of the helix ( in the beam frame ) @xmath7 , alfvn speed @xmath8 , and magnetosonic mach number @xmath9 . this describes a plasma in which the magnetic field is dominant and in a rather tight helix , and alfvn waves are responsible for the moving transverse patterns .