the interaction of extragalactic jets with their environment leads , under rather general circumstances , to the stratification of the beam of the jet in the direction normal to its velocity .
the morphology of fr i sources ( e.g. , m87 , @xcite ; 3c 31 , @xcite ; 3c 296 , @xcite ) has been explained in terms of a jet whose dynamics is dominated by a boundary shear layer . in such layers the emissivity in radio and optical maps peaks and , in some cases ( e.g. , mrk 501 ; @xcite ) , the spectrum is rather flat , suggesting that the acceleration of non - thermal particles takes place right at the boundary region .
further evidence in favor of such stratification is provided by radio polarization measurements , which indicate that , towards the jet boundary , the magnetic field is highly ordered and parallel to the jet axis @xcite .
a radial stratification of jets in fr  ii radio galaxies has also been observed @xcite .
the edge brightening found in 3c353 has been interpreted by @xcite as due to the doppler hiding of the emission of the central spine of the jet , suggesting that most of the observed radiation comes from the jet boundary layer . on the other hand ,
the _ rails of low polarization _ found by @xcite close the the jet boundary indicate that the magnetic field in the layer is either axial or toroidal but not radial . even at parsec scale ,
the fr  ii radio jet of 1055 + 018 exhibits a transverse structure consistent with a _ spine - shear boundary layer _ jet morphology @xcite .
also sources which are in the borderline dividing fr  i and fr  ii sources show evidences of radial jet stratification .
a couple of examples of such sources are 3c  15 and cen  a. in the first case the jet is generally narrower in the optical than in the radio @xcite , and simple spine - sheath model may account for the polarization angle differences seen in the optical and radio data . for cen  a , the recent x - ray data @xcite
can also be properly explained by a stratified jet model with a radially decreasing velocity field .
the aforementioned observational evidence along with some others ( e.g. , @xcite ) , suggests that the velocity in some jet boundary layers is smaller ( but still relativistic ) than the velocity of the beam itself . on the other hand ,
laboratory experiments show that , almost unavoidably , turbulent interaction layers may develop as jets propagate into a viscous medium ( e.g. , @xcite ) .
these laboratory shear layers display a radial velocity profile roughly monotonic in which the velocity of the jet core smoothly decreases until it vanishes at the external medium .    from a theoretical point of view
, the reynolds number of the jet flow may reach values of @xmath0 and , hence , it seems unquestionable that such flows will quickly develop a turbulent boundary layer which spreads into the flow and leads to the entrainment and acceleration of ambient gas @xcite .
however , the nature and the amount of viscosity in relativistic jets is still largely unknown . most probably , the effective viscosity in the lateral interaction between the jet an the external medium is a mixture of the turbulent eddy viscosity , magnetic viscosity ( based on the finite size of the larmor radius ; see , e.g. , @xcite ) and _ cosmic - ray _ viscosity ( i.e. , the viscosity originated due to the diffusion of the momentum carried by energetic particles ; see , e.g. , @xcite ) .
whichever the origin of the viscosity is , the boundary layer itself results from the nonlinear growth of kelvin - helmholtz ( kh ) instabilities on the jet surface . the fact that relativistic jets are prone to these instabilities was demonstrated theoretically in the pioneer work of @xcite .
later , many other papers have explored the stability of relativistic jets with respect to kh instabilities in the vortex - sheet approximation ( e.g. , @xcite ) . @xcite and @xcite presented the first attempts to study ( in the linear regime ) the influence of a finite thickness shear layer in classical and relativistic jets , respectively .
very recently , @xcite has also considered the stability properties ( in the linear regime ) of spine - sheath magnetized , relativistic jets in which the magnetic field is parallel to the flow velocity .
the development of kh modes in sheared relativistic ( slab ) flows from the linear ( analytic ) regime to the non - linear regime by means of relativistic hydrodynamic simulations has been carried forward by @xcite .
alternative approaches to the linear analysis of the stability of stratified relativistic jets and sheared relativistic jets have been performed by @xcite , @xcite @xcite and @xcite , respectively .
further support in favor of radial jet stratification as a consequence of the interaction between a jet and its environment is provided by three - dimensional hydrodynamic simulations of relativistic , large scale jets @xcite and two - dimensional simulations of relativistic magnetized jets @xcite .
these simulations show that radially stratified jets are developed out of initially uniform beams .
these numerical works display also a rather smooth transition layer with decreasing transverse velocity towards the jet boundary and high specific internal energy .
it must be noted that the viscosity in such simulations is of numerical origin , and mimics only qualitatively actual viscous flows .
in addition to the existence of _ environmental _ reasons to produce a radial flow stratification in astrophysical jets , this effect can also be a natural consequence of the jet formation process .
jet launching from accretion disks may directly lead to a certain mass flux profile and/or magnetic flux profile within the jet .
this possibility has been verified by means of axisymmetric mhd simulations both in the classical @xcite , as well as in the general relativistic @xcite regime .
@xcite show that , under physically plausible conditions ( see sect .  [
sec : rmhd_boost ] ) , the lateral structure of the interaction layer of a relativistic flow can be richer than previously thought .
this is also the case in relativistic magnetized jets @xcite .
due to a purely relativistic phenomenon ( the conversion of specific enthalpy into bulk lorentz factor across a flow discontinuity ) , the jet velocity may increase right at the contact discontinuity ( cd ) separating the jet and the external medium . beyond the cd the velocity decreases steeply to zero in the radial direction .
the growth of the velocity or , equivalently , of the lorentz factor of the beam is associated to the development of a rarefaction wave that emerges form the cd where the density and the pressure decrease with respect to their corresponding values at the jet spine .
the profiles exhibited by the hydrodynamic variables in these layers are non - monotonic , in contrast with the smooth monotonic profiles that have been typically discussed in the literature .
thereby , we will refer to them as _ anomalous _ or _
ar _ shear layers .    as mentioned above , jet shear layers , represent natural sites for particle acceleration , providing high - energy cosmic rays and influencing the dynamics of relativistic jets in extragalactic radio sources by forming cosmic - ray cocoons @xcite .
the efficiency of the acceleration process in these turbulent shear layers depends on the particle mean free path and on the velocity structure .
@xcite ( 2002 ; so02 hereafter ) performed a through study of both the acceleration processes acting at jet boundaries and the resulting observational effects .
they former authors considered the possibility that turbulent standard boundary layers of kiloparsec scale jets may substantially contribute to the radiative output of the jet . in the present work
we will closely follow so02 working hypothesis , however we will replace the monotonic shear layer kinematic structure they assumed by an anomalous one as suggested by @xcite . the aim being to show whether ar - shear layers imprint any distinctive feature in the radiation produced by the jet as compared with standard boundary layers .
we show in this paper that , indeed , the radiative output of jets bound by anomalous shear layers significantly differs from that of ordinary boundary layers . thus observations may confirm or rule out the existence of such anomalous shear layers in kiloparsec scale jets and , consistently , this
may be used to constrain hardly known physical parameters in these regions .
the plan of the paper is to first show [   [ sec : rmhd_boost ] ] the influence of magnetic fields in the profiles of the physical variables in anomalous boundary layers . in this way
we extend the results of @xcite and @xcite to the adequate parameter range . in sect .
[ sec : model ] we summarize the basic model developed by so02 and adapt it to account for the kinematic differences when applied to anomalous shear layers whose profiles will be discussed in sect .
[ sec : rmhd_boost ] .
section  [ sec : sed ] provides the spectral energy distribution obtained for the different models of kinematic shear layers shown in the previous section .
we discuss our results in sect .
[ sec : discussion ] and sum them up in sect .
[ sec : summary ] .
the dynamics of a relativistic magnetohydrodynamic jet in an external medium can be treated as the motion of two fluids , one of which ( the jet ) is much hotter and at higher ( or equal ) pressure than the other one ( the external medium ) , and is moving with a large tangential velocity with respect to the cold , slowly moving fluid .
furthermore , the jet can be magnetized while the external medium magnetic field is much smaller than that of the jet .
@xcite found that if the specific enthalpy of the jet is sufficiently large , a net conversion of internal to kinetic energy can be produced through a purely relativistic channel .
the reason for its exclusively relativistic character is that , the evolution of contact or tangential discontinuities ( like , e.g. , the one that separates laterally the jet from the external medium ) , is governed by a genuine coupling between the specific enthalpy and the lorentz factor at both sides of the discontinuity which does not exist in newtonian ( magneto-)hydrodynamics . such mechanism yields a substantial hydrodynamic boost along the boundary layer between the jet and the external medium . basically , the boost along the jet lateral boundaries is due to the work done by the external medium on an overpressured , relativistic jet .
@xcite extended the validity of the results of @xcite to the case in which the jet is magnetized and conclude that the presence of a jet magnetic field ( either poloidal or toroidal ) enlarges the boost that might be obtained by purely hydrodynamic means .
the magnetohydrodynamic boost is larger for jets carrying toroidal than poloidal magnetic fields .
the present analysis is aimed to show what is the qualitative variation of the lorentz factor and of the magnetic field across a jet boundary layer .
therefore , following @xcite , we model the interaction between the jet and the external medium as a one dimensional riemann problem in cartesian coordinates .
furthermore , we restrict to the case in which the magnetic field is perpendicular to the flow speed in the jet .
as demonstrated by @xcite and @xcite this simple model retains the basic features of the phenomenon .
certainly , a more detailed model would require three dimensional relativistic magnetohydrodynamic simulations .
we solve exactly the riemann problem set by the jump in the conditions between the jet and the external medium assuming that the magnetic field is perpendicular to the velocity field .
for this purpose we use the method devised by @xcite . the two uniform `` left '' ( jet ) and `` right '' ( external medium ) states which set the proposed riemann problems possess different and discontinuous magnetohydrodynamical properties : the rest - mass density @xmath1 , the total pressure @xmath2 ( @xmath3 and @xmath4 being the magnetic and the gas pressure , respectively ) , the components of the velocity normal @xmath5 and tangential @xmath6 to the initial discontinuity , and the magnetization @xmath7 . in the following , we will use for the quantities in the jet and in the external medium the subscripts `` l '' and `` r '' , respectively .
like in the purely hydrodynamic case , depending on the conditions of the left and right states two qualitatively and quantitatively different solutions develop ( fig .
[ fig : sl_structure_magnetic ] ) . on the one hand ,
a pattern of waves can be formed , which we indicate as @xmath8 ( fig .
[ fig : sl_structure_magnetic ] left panel ) , where @xmath9 refers to the shock propagating towards the left sweeping up the jet , @xmath10 to the shock moving towards the right crossing the external medium and @xmath11 to the contact discontinuity between the two . on the other hand ,
the pattern can change for sufficiently large velocities parallel to the jet axis ( tangential ) and in this case the inward - moving shock is replaced by a rarefaction wave thus producing a @xmath12 pattern ( fig .
[ fig : sl_structure_magnetic ] right panel ) .
once again , we point out that the change in the wave pattern does only happen in relativistic ( magneto-)hydrodynamics . only in the relativistic regime the evolution of a riemann problem depends upon the components of the velocity parallel to the initial discontinuity or , equivalently , upon the lorentz factor at both sides of the initial discontinuity .
furthermore , only in the relativistic regime there is a coupling between the momentum and energy fluxes that depends upon the specific enthalpy of initial left and right states .
we notice that the variation of the magnetization in the intermediate left state is less than a few percent when a @xmath12 solution forms ( see fig .
[ fig : sl_structure_magnetic ] right panel ) . in case
a @xmath13 pattern develops , the magnetization in the shocked left intermediate state @xmath14 grows due to the shock compression , but , in any case , @xmath15 .
although @xcite have already explored the influence of a dynamically relevant magnetic field in the jet on the aforementioned boost , they restricted their study to only four cases and to extremely hot jets ( @xmath16 ) , which is appropriate for gamma - ray burst jets ( see , e.g. , @xcite ) .
here we consider a parameter space more adequate for kiloparsec scale jets .
more precisely , we consider _ warm _ or _ cold _ jets with a ratio @xmath17 and a magnetization @xmath18 , which is varied between @xmath19 ( slightly magnetized jet ) and @xmath20 ( poynting flux dominated jet ) is considered in their model mhdb . ] .
the left state will have comparable but larger total pressure than the external medium @xmath21 , will be under dense @xmath22 , with a fixed bulk lorentz factor @xmath23^{-1/2}=10 $ ] .
we choose to fix the external medium ( right state ) to be cold ( @xmath24 ) , non - magnetized and non - moving .
precisely , we take for the external medium @xmath25 , @xmath26 , @xmath27 , where @xmath28 is a normalization constant which allows us to be scale - free .
we point out that with this parametrization we provide physically plausible conditions for both the jet and the external medium . for convenience , hereafter we will assume @xmath29 , unless stated otherwise .
we explore the resulting riemann solution as we vary the component of the jet velocity normal to the shear layer . ] and the jet magnetization @xmath18 .
we point out that models with @xmath30 are actually rather cold , since @xmath31 ( e.g. , for @xmath32 , @xmath33 ) .
according to the most broadly accepted view , even mhd - generated jets may become matter dominated at kiloparsec scale ( e.g. , @xcite ) .
however , there are others who advocate jet models which are basically electromagnetic entities at such large scales ( e.g. , @xcite ) . in view of these two possible extremes
, we have covered the @xmath34-parameter space with jet models which are both matter dominated and poyinting - flux dominated at kiloparsec scale . in fig .
[ fig : betal_dep ] we show the value of the lorentz factor reached at the state left to the contact discontinuity @xmath35 .
this intermediate state is made out of boosted or deboosted jet matter depending respectively on whether a @xmath36 or a @xmath13 solution develops .
figure  [ fig : betal_dep ] illustrates how a magnetohydrodynamic boost , where @xmath37 , exists below a critical value of the normal velocity ( @xmath38 in the case considered here ) whereby a @xmath12 solution develops .
the critical value of @xmath39 is independent of the jet magnetization ( note the crossing of all the solutions at the point where @xmath40 ) .    the magnetohydrodynamic boost is larger if the magnetization of the jet is increased : while for a poorly magnetized jet ( @xmath41 ) the boost may increase the bulk lorentz factor of the layer by @xmath42 , for a strongly magnetized jet ( @xmath43 ) the lorentz factor increase can be @xmath44 .
increasing the magnetization beyond @xmath45 does not produce larger boosts ( the effect saturates ) .
below @xmath46 the results are almost indistinguishable from the case @xmath47 and we have decided not to include more lines in the fig .
[ fig : betal_dep ] for the sake of readability .
the bottom line of this parametric study is that for conditions realizable in kiloparsec scale jets , one may find an increase of the lorentz factor in the transition layer between the jet and the external medium ranging from @xmath48 , while the magnetization in that layer is basically the same as in the jet if a @xmath36 pattern occurs .
furthermore , since only the strength of the magnetohydrodynamic boost , but not its existence , depends on the topology of the field ( see @xcite ) , we point out that it also happens in the presence of randomly oriented magnetic fields .
our model follows very closely that devised by so02 .
we assume that non - thermal particles can be accelerated at the boundary layers of relativistic jets ( e.g. , @xcite ) , thereby producing high - energy cosmic rays . as mentioned in sect .
[ sec : intro ] , the main goal of this paper is to replace the prototype shear layer assumed in so02 by an anomalous shear layer of the type suggested by @xcite and whose basic features have been outlined in sect .
[ sec : rmhd_boost ] . at the most basic level
, this means to replace the kinematic shear layer structure assumed by so02 , namely , a jet sheath where the bulk lorentz factor ( @xmath49 ) decreases monotonically ( fig .
[ fig : sl_structure]a ) , by a boundary layer where the lorentz factor develops a spike in which it is larger than at the jet core ( fig .
[ fig : sl_structure]b ) .
we shall assume that the magnetic field is uniform in the whole shear layer , which is compatible with the results of our previous section and a reasonable assumption in our simple model .
the working hypothesis made by so02 are also valid in our case since , the only difference between monotonic and anomalous layers is simply kinematic ( see sect .  [
sec : kinematic ] ) . thus , following so02 , the maximum electron lorentz factor which can be obtained by the combination of acceleration of particles along the shear layer and the radiative cooling of such particles is @xmath50 where @xmath51 is the magnetic field in microgauss and @xmath52cms@xmath53 is the alfvn speed @xmath54 in units of @xmath55cms@xmath53 . for kiloparsec scale jets , typical values of these two parameters
are @xmath56 and @xmath57 , which yields @xmath58 .
we assume that the resulting spectral energy distribution of the electrons accelerated at the layer consists of a low energy power law @xmath59 , with @xmath60 , finished with a high - energy component modeled as a nearly mono energetic peak at @xmath61 due to the pile - up of accelerated particles caused by losses ( see fig .  1 of so02 ) .
in terms of the dirac delta function , @xmath62 , at energies above the injection energy , @xmath63 , the considered quasy - stationary , averaged spectrum of electrons can be cast in the form @xmath64 where @xmath65 and @xmath66 are normalization parameters whose values ( @xmath67@xmath68 , @xmath69@xmath68 ) result from assuming equipartition between the magnetic field energy density @xmath70 and the energy densities of both electron spectral components : @xmath71    we point out that in so02 , the power - law part of the electron spectrum was also endowed with an exponential cut - off instead of with a heaviside @xmath72 function in spite of their eqs .  6 and 7 . the @xmath72 function was used in their work only to compute analytic estimates , not for the numerical integrations necessary to calculate their spectral energy distributions @xcite .
assuming equipartition between the particle and magnetic field energy in relativistic sheared jets has been proved to be a rather solid theoretical assumption ( see , e.g. , @xcite ) .
however , in some cases the values of @xmath65 and @xmath66 derived from eq .
[ eq : normalization ] and based on equipartition arguments may yield only lower bounds of the true values .
@xcite argue that powerful jets in quasars and fr ii objects can be far from the minimum energy condition , and the field strength very likely exceeds the equipartition value .
on the other hand , upper limits to the inverse compton radiation of the jet in m87 imposed by hess and hegra cerenkov telescopes also indicate that the magnetic field can not be weaker than the equipartition value @xcite .      the relevant radiative processes taking place at the boundary layer of a relativistic , kpc - scale jet are , on the one hand , the synchrotron ( syn ) emission of the spectral family represented by eq .
[ eq : n(gamma ) ] and , on the other hand , their inverse compton ( ic ) cooling due to the interaction with the previously produced synchrotron photons ( synchrotron self compton -ssc- ) or with seed photons of the cosmic microwave background ( external compton -ec- ) .
these radiative processes are the dominant for large scale relativistic ( @xmath73 ) jets at distances from the galactic nucleus @xmath74kpc .
following so02 , we assume that the magnetic field is randomly oriented in the shear layer . ] . additionally , we neglect synchrotron self - absorption and , therefore , we limit the computation of the resulting spectra to frequencies @xmath75hz .    restricted to the case of an isotropic electron distribution @xmath76 the synchrotron emissivity averaged over a randomly oriented magnetic field @xmath77 can be computed from @xmath78 where @xmath79 and @xmath80 is a combination of the modified second order bessel functions @xcite as @xmath81 .
\label{eq : r(x)}\ ] ]    the ic photon emissivity in the source frame is assumed to be uniform .
] , @xmath82 , at a certain photon energy @xmath83 ( in units of the electron rest mass , @xmath84 , @xmath85 ) and scattering direction @xmath86 , is given by @xmath87 where @xmath88 is the seed photon number density as a function of the incident photon energy @xmath89 and the photon direction @xmath90 , @xmath91 is the electron energy distribution , @xmath92 is the angle between the electron and the incident photon directions , @xmath93 is the electron lorentz factor corresponding to a velocity @xmath94 and @xmath95 its direction in the plasma source frame ( e.g. , @xcite ) .
the rest - frame emissivity in cgs units is @xmath96 .    in the calculation of the emissivity due to ic process ( ssc or ec )
, we limit ourselves to the thompson regime , in which the scattering cross section @xmath97@xmath98 is independent of the seed photon energy , which is an adequate approximation in the regime @xmath99 .
following @xcite , the scattered photons will be beamed along the direction of the electron motion ( i.e. @xmath100 ) and will emerge after scattering with an average final energy @xmath101 .
we assume as so02 that the synchrotron radiation and the electron distribution are isotropic and , therefore , @xmath102 and @xmath103 .
then , in eq .  [ eq : n_ic ] one can use @xmath104 in order to obtain the following expression of @xmath105 : @xmath106 using eqs .
[ eq : n(gamma ) ] ,  [ eq : n_ssc ] , the relation @xmath107 and the optically thin synchrotron intensity @xmath108 , where @xmath109 is the emitting region linear size , the integral form of the synchrotron self - compton emissivity @xmath110 is written as @xmath111 where @xmath112 is the lower limit imposed by the restriction @xmath113 , which guaranties that we stay within the thompson regime .
the cmb radiation is strongly anisotropic in the source frame which is moving at a bulk lorentz factor @xmath49 , providing an external source of seed photons which move opposite to the jet direction , i.e. , @xmath114 .
thereby , @xmath115 , where @xmath95 is the direction of the colliding electron in the comoving frame of the shear layer , similarly to what we made for the case of ssc radiation . in the thomson regime
@xmath116 $ ] , since the scattered photon energy is @xmath117 , where @xmath118 is the energy of the incident cmb photon in units of the electron rest - mass .
hence , eq .  [ eq : n_ic ] can be cast as @xmath119\   .
\label{eq : n_ec}\ ] ]    in order to compute the ec emissivity produced by the interaction of cmb photons with the power - law part of the electron distribution , we model the cmb field as a monochromatic radiation with photon energies equal to the average value of a thermal black body spectrum @xcite but blue - shifted by a factor @xmath49 , i.e. , @xmath120 , and a number density @xmath121 .
the value that we adopt for the energy density of the cmb radiation in the rest frame of the cmb is @xmath122erg@xmath68 . with these approximations ,
the ec emissivity associated to the power - law part of the electron distribution can be computed analytically @xmath123^{3/2 }            \exp{\left[-\left(\frac{\epsilon}{\gamma ( 1+\mu )              \langle \epsilon_{\rm cmb}^\ast \rangle \gamma_{\rm eq}^2 } \right)^{1/2 }             \right]}\nu^{-1/2}\:\ : } \nonumber \\   & & \mbox{if }           \gamma_0 ^ 2 \le \displaystyle {            \frac{h\nu}{2.7kt^\ast\gamma(1+\mu)}}\ , .
\label{eq : j_ec_1}\end{aligned}\ ] ]    for the monoenergetic electron peak component , we use the exact black body spectrum transformed to the source frame . in this case , the expression of @xmath124 in the boundary layer comoving frame is ( see eq .  2.58 of @xcite ) @xmath125^{-1}\ , , \label{eq : n_cmb}\ ] ] where @xmath126 is the boltzmann constant
. substituting eq .
[ eq : n_cmb ] and @xmath127 into eq.[eq : n_ec ] , taking into account that @xmath128 and the restriction imposed by the thompson regime ( @xmath129 ) , one obtains @xmath130    the total ec emissivity is the sum of eqs .
[ eq : j_ec_1 ] and [ eq : j_ec_2 ] , i.e. , @xmath131      we model a large scale jet as a uniform cylindrical spine of radius @xmath132 with a lorentz factor @xmath133 followed by a transition layer where the lorentz factor changes linearly with radius until @xmath134 , where the flow velocity either becomes zero ( similar to the so02 model ) or reaches a maximum from which it abruptly decreases to zero ( figs .  [
fig : sl_structure ] , [ fig : scheme ] ) . for completeness
, we also consider a boundary layer where the lorentz factor is uniform .
more explicitly @xmath135 where we take @xmath136 for our reference models .
the subscript @xmath137 labels the case of a monotonic shear layer ( @xmath138 ) , @xmath139 represents the case of an anomalous shear layer ( @xmath140 ) and @xmath141 denotes the case of a jet where the boundary layer is uniform and has the same lorentz factor as the core ( @xmath142 ) .
the kinematic idealization of the case corresponding to an anomalous layer is a rough prototype of the structure displayed by the lorentz factor when a @xmath36 solution develops ( fig .
[ fig : sl_structure_magnetic ] ) .
we neglect the variation of the magnetic field across the boundary layer because it is rather small ( see sect .  [ sec : rmhd_boost ] ) and we assume that the magnetic field is uniform and equal to @xmath143 g .    in order to properly normalize the emissivity properties of the different jet models , we choose that all the shear layers transport the same kinetic power as the monotonically decreasing one with a jet radius @xmath144kpc and a density and pressure equal to the values chosen for the left state of the models shown in sect .
[ sec : rmhd_boost ] .
this means that we adjust the external jet radius of the other two cases to @xmath145kpc and @xmath146kpc . since , typically , for extragalactic jets , @xmath147g@xmath68 @xcite , the fluxes of mass and of magnetic field in the jet axial direction are @xmath148    for each cylindrical shell in the shear layer of radius @xmath149 and moving with a bulk lorentz factor @xmath150 , we compute the local source frame emissivities for the radiative processes considered in sect .  [
sec : model ] ( eqs .
[ eq : j_syn ] , [ eq : j_ssc ] and [ eq : j_ec ] ) . under the assumption
that the radiating electrons are distributed uniformly in the boundary layer ( for a comoving observer ) , the emissivity of the sub - layer at the distance @xmath149 from the jet axis can be computed as @xmath151 where @xmath152^{-1}$ ] is the doppler factor associated to a sub - layer moving with the lorentz factor @xmath153^{-1/2}$ ] at an angle @xmath154 with respect to the line of sight .
the subscript @xmath155 syn , ssc or ec specifies the radiative process .
the observed flux density from the boundary layer volume @xmath156 appropriately modified by the flow - beaming patterns is @xmath157 @xmath158 being the distance to the observer ( we take @xmath159mpc unless stated otherwise ) .
@xmath160 corresponds to a volume element of the shear layer with an observed length @xmath161 , which we assume to be @xmath162kpc .
the integrals involved in the evaluation of eq .
[ eq : s(nu , theta ) ] are performed numerically and , in order to compute the ssc contribution to the flux density , we estimate the ssc emissivity ( eq .
[ eq : j_ssc ] ) by taking @xmath163 .
our main goal is to outline the observational differences between the emission of large - scale jets laterally endowed by different types of shear layer .
the observed spectral energy distribution of the radiation emitted by different kinematic boundary layers in a prototype kpc - scale jet is show in fig .
[ fig : sed_ar - so ] for two different viewing angles @xmath164 and @xmath165 .
we neglect both the absorption of very high energy ( vhe ) gamma - rays during the propagation to the observer and any contribution from a separate non - thermal electron population accelerated at shocks happening at the jet core .
the choice of a small and large viewing angle is motivated by the fact that for small viewing angles , @xmath166 , the larger doppler boosting provided to the emitted radiation by the larger lorentz factor in ar layers than in monotonic layers makes the total flux density of the former case also slightly larger ( fig .
[ fig : sed_ar - so ] left panel ) .
the differences in the total flux density are smaller in the synchrotron dominated part of the spectrum than in the ec dominated one . at larger angles ( fig .
[ fig : sed_ar - so ] right panel ) , the flux density predicted for standard layers is 2 to 5 orders of magnitude larger than for uniform or ar layers due to the doppler hiding of the later two cases .
the differences are very important in the ec dominated part of the spectrum ( i.e. , in the gamma - ray regime ) .
both , anomalous and uniform boundary layers generate a rather similar spectral pattern , particularly at small viewing angles . as we observe the jet at larger values of @xmath167 ( fig .
[ fig : sed_ar - so ] right panel ) the smallest lorentz factor of the uniform case provides a smaller doppler hiding which results into a larger flux density than in case of ar layers .
considering specifically the flux density produced by anomalous layers ( fig .
[ fig : sed_ar - thetas ] ) , we note two main differences with the standard case studied by so02 . first , the separation of the ssc and ec peaks corresponding to ic emission of the power - law distribution of the electrons ( eq .  [ eq : n(gamma ) ] ) is larger in the former case . for jets with ar layers the frequency separation between the ec1 and ssc1 peaks grows with increasing viewing angle from less than one order of magnitude ( @xmath168 ) to more than two orders ( @xmath169 ) . in case of jets endowed by standard layers ( fig .
[ fig : sed_uniform - thetas ] ) , the frequency separation between the ec1 and ssc1 peaks also grows with the viewing angle , but much less than in the previous case : the frequency separation between the two peaks is smaller than a factor 8 at @xmath170 .
if a sufficiently large number of tev photons were detected , this difference could be used to distinguish between jets with standard and anomalous layers .
however , we repeat here the cautionary note of so02 : the spectral energy distributions ( seds ) computed in this paper correspond to the most optimistic scenario , with a highly relativistic jet spine and very efficient acceleration creating electrons with large lorentz factors up to @xmath171 . either smaller values of @xmath133 or @xmath172 would reduce the observed very high energy gamma - ray flux and shift the corresponding peaks to lower energies .
the second most important difference between ar and standard radiating boundary layers is that the exponential decay of the synchrotron spectrum is shifted towards smaller frequencies ( delving deeper into the very soft x - ray regime ) .
this has a big impact in the effective x - ray spectral index as we will discuss below .
as noted by @xcite and confirmed by so02 , extragalactic jets with a monotonic boundary layer may show a jet - to - counterjet radio - to - optical brightness ratio ( @xmath173 ) which is smaller than if the jet were uniform .
we find that such an assertion has to be modified for jets limited by ar - boundary layers . in the latter case , at small viewing angles ( namely , @xmath174 for the considered parametrization ) the jet / counterjet asymmetry is @xmath175 times larger than the one corresponding to a uniform jet with the same kinetic power ( fig .
[ fig : jet2counterjet ] ) . at larger viewing angles ,
the brightness asymmetry turns out to be almost indistinguishable from the case of an equivalent uniform jet ( see in fig .  [
fig : jet2counterjet ] how the solid and dashed lines practically overlap each other for @xmath176 if @xmath136 or for @xmath177 in case @xmath178 ) .
certainly , since the brightness asymmetry is a result of the doppler boosting ( hiding ) of the jet ( counterjet ) , the smaller is the value of the maximum lorentz factor in the anomalous shear layer ( @xmath179 ) , the smaller will be the maximum asymmetry at small @xmath167 .
if the lorentz factor of the jet spine is small , the jet - to - counterjet brightness ratio reduces in absolute value ( see the grey lines corresponding to @xmath180 , in fig .  [
fig : jet2counterjet ] ) .
, for the case of @xmath178 we take @xmath181 .
] however , even for the relatively small value of @xmath178 the jet - to - counterjet brightness ratio of the anomalous shear layer case is larger than the one corresponding to the standard case at @xmath136 independent of the viewing angle ( e.g. , at @xmath182 , @xmath183 } \simeq 4.35 $ ] for the ar model , while @xmath184 } \simeq 4.25 $ ] ) .
if there were a handle on the value of the lorentz factor of the jet , this difference might allow us to discriminate observationally between jets with anomalous and standard shear layers .
however , given the similarity of the jet - to - counterjet brightness ratio between jets with anomalous layers and jets with _
top - hat _ profiles ( uniform ) , it might be practically impossible to disentangle from this unique value whether the jet is shielded by an anomalous layer or by no layer at all .    considering the results obtained for standard boundary layers ,
so02 concluded that the mildly relativistic velocities inferred from the observed brightness asymmetries at tens of kiloparsec scales , may correspond to a slower boundary layer and not necessarily to the faster jet core , which might be highly relativistic at the observed distances .
if anomalous shear layers may happen in actual jets , slower spines are preferred and relatively moderate maximum shear layer lorentz factors might also be invoked to explain the observed jet - to - counterjet asymmetries .    in our models ,
the x - ray radiated power of kpc - scale jets bounded by standard layers is dominated by the high - energy electron bump spectral component , implying that the resulting synchrotron x - ray flux lies above the extrapolated radio - to - optical continuum ( see so02 and figs .
[ fig : sed_ar - so ] , [ fig : sed_uniform - thetas ] ) .
this is also true for jets endowed by anomalous boundary layers ( see fig .
[ fig : sed_ar - thetas ] at @xmath185hz ) . for both types of layers ,
the spectral slope at x - ray frequencies may be rather different from that of the power - law at radio - to - optical frequencies , and the difference depends also on the viewing angle .
this implies that the effective x - ray spectral index @xmath186 computed between @xmath187kev and @xmath188kev as @xmath189 \over \log \left [ \nu^{\ast}_2 \ , / \ ,        \nu^{\ast}_1 \right ] } , \label{eq : alphax}\ ] ] can be significantly distinct .
we point out that @xmath190 is computed in this paper as the value of the flux density calculated using eq .
[ eq : s(nu , theta ) ] instead of so02 who only take @xmath191 , with @xmath192 given by eq .
[ eq : r(x ) ] , and @xmath193 .    for highly relativistic jet cores ( @xmath136 ) , fig .
[ fig : alphax ] ( black lines ) shows that the most noticeable difference between jets with anomalous or uniform shear layers and jets with standard limiting boundaries , is the fact that the effective x - ray spectral index of the jet @xmath194 is larger than the corresponding to the counterjet @xmath195 for jet inclinations @xmath196 . for standard boundary layer jets @xmath197 , @xmath198 holds ( i.e. , the counter jet has a steeper x - ray continuum as compared to the jet spectrum ) .
furthermore , jets with uniform and ar layers show a much larger$ ] , @xmath199 is 6 to 7 times larger than @xmath200 . ]
jet and counter jet effective x - ray spectral index than jets with standard layers if the inclination is @xmath201 . indeed , @xmath194 in the former two types of boundary layers is even larger than the @xmath202 corresponding to jets with monotonic layers for viewing angles @xmath201 .
the large values of the effective x - ray spectral index in jets with uniform and ar layers is due to the fact that the decay of the spectrum after the synchrotron peak crosses the x - ray band and moves towards optical frequencies for inclinations @xmath203 ( fig .
[ fig : sed_ar - thetas ] ) .
consistently , the flux density at 5kev decays abruptly and , thus @xmath204 grows .    in case of jets with more moderate spine velocities ( @xmath180 ; fig .
[ fig : alphax ] , grey lines ) , there is no viewing angle for which @xmath205 , independent on the shear layer model .
however , it remains true that the counterjet x - ray spectral index of jets bearing uniform or ar layers is much larger than the one corresponding to standard boundary layer jets ( there is a factor of 2 to 4 difference between them ) .
taking together the results for models with jet spines flowing at @xmath178 and 10 , it turns out that the range of variability of both @xmath194 and @xmath206 is larger for jets flanked by uniform and ar layers than for jets with standard boundary layers .
interestingly , the large - scale x - ray emission observed from several radio - loud agns exhibits rather different spectral characteristics . for instance , x - ray spectra of the known jets in quasars are very flat ( e.g. , @xmath207 for 3c 207 , @xcite ; @xmath208 for pks 1127 , @xcite ; @xmath209 for 3c 273 and pks 0637 ) .
the low effective spectral index is consistent with small jet inclinations ( c.f .
[ fig : alphax ] , inset ) independent of the shear layer model adopted . on the other hand , jets in radio galaxies tend to display relatively steep x - ray spectra , @xmath210 ( @xcite for 3c 66b , @xcite for b2 0206 and b2 0755 ) .
these values correspond to large inclination angles on fig .
[ fig : alphax ] , and/or lower jet lorentz factors for jets with standard boundary layers .
but , if uniform or ar emitting layers are considered , the viewing angles needed to account for such values of the x - ray spectral index are smaller ( @xmath211 $ ] if @xmath136 , or @xmath212 $ ] if @xmath178 ) .
very recently , @xcite have reported values of the x - ray spectral index @xmath213 at distances of @xmath214kpc from the nucleus in cen a. these large values of @xmath215 are not within reach of jet models flanked by standard layers ( the maximum of @xmath194 is 1.84 at @xmath170 for the model with @xmath216 ) .
nonetheless , jets with uniform or ar boundaries may yield @xmath217 for @xmath218 ( @xmath219 ) if @xmath136 ( @xmath180 ) .
indeed , considering that the inclination of cen a is @xmath220 @xcite , our results suggest that , if a uniform or an ar emitting boundary layer were responsible for the emission of cen a at distances from the nucleus @xmath221kpc , small values of the spine lorentz factor ( @xmath222 ) are preferred to explain the large values of @xmath215 observed at such scales .
although very recently deep _
chandra _ of observations of cen a @xcite have provided us with spectra of some counterjet x - ray features , this is not the usual case for large scale counterjets .
there are , however , some estimates of the lower limit of the x - ray brightness asymmetry , like e.g. , for the 3c 66b ( @xmath223 , @xcite ) , pictor a ( @xmath224 , @xcite ; @xmath225 @xcite ) , or pks 1127 - 145 ( @xmath226 , @xcite ) .
our model of ar or uniform radiating boundary layer tends to favor larger x - ray brightness asymmetries between the jet and the counterjet .
the reason being that the x - ray emission of the jet is dominated by the synchrotron radiation from the monoenergetic , hard ( @xmath227 ) electron component , while the counterjet emission is dominated by the ec radiation of the power - law part of the electron spectrum .
so02 chose the parameters of their radiating boundary layers such that they exhibit clearly the main characteristics of their model .
we do not repeat here their discussion .
instead , we focus on the discussion of the choice of the new parameters of our model .
so02 anticipated that the form of the @xmath150 radial profile and of the related spatial variation of the acceleration efficiency can significantly influence the beaming pattern and the intensity of the boundary layer emission .
we actually have confirmed this point by replacing the monotonically decaying profile of @xmath150 by an idealization of the lorentz factor profile that results from ar boundaries .
this results in a different relative normalization of the ssc and ec components presented in figs .
[ fig : sed_ar - so]-[fig : sed_uniform - thetas ] .
the uniform layer model could be considered as a limiting case of ar layers when @xmath228 .
thus , when looking at the plots of the jet - to - counterjet ratio ( fig .  [ fig : jet2counterjet ] ) or the x - ray effective spectral index ( fig .  [ fig : alphax ] ) , any model including an ar layer with @xmath229 $ ] ( if @xmath136 ; black lines in figs .
[ fig : jet2counterjet ] and [ fig : alphax ] ) or @xmath230 $ ] ( if @xmath178 ; grey lines in figs .
[ fig : jet2counterjet ] and [ fig : alphax ] ) , would display a graph located between the boundaries set by the solid and dashed lines of the corresponding figures .
since the area enclosed by these two line types ( solid and dashed ) is rather small , it turns out that our choice of the value of @xmath231 does not significantly influence the results .
more precisely , if for the considered values of @xmath133 ( 3 or 10 ) we would have picked up any other value of @xmath232 $ ] , neither the value of the jet - to - counterjet brightness ratio , nor @xmath215 , nor even the spectrum would have changed appreciably . since
, keeping all other physical conditions fixed , a decrease of @xmath233 yields smaller values of @xmath234 , the results we obtain shall be qualitatively valid both for magnetized or non - magnetized jets .
the jet - to - counterjet effective x - ray spectral index asymmetry and its angular dependence ( fig .  [ fig : alphax ] ) depends on which of the two spectral components at x - ray frequencies dominates : synchrotron radiation of the monoenergetic , high - energy electron component or the ec radiation from low energy electrons ( @xmath235 ) . in its turn
, the dominance of any of these two processes depends on the exact shape of the electron spectrum at the highest energies , on the viewing angle , on the jet lorentz factor and on @xmath150 .
figure  [ fig : alphax ] serves as an example of the comparative behavior of different boundary layer models ( in which all the parameters are the same except @xmath150 ) .      in this work ,
we compare the observational signatures imprinted by different kinematic models of shear layers of kiloparsec scale jets on the radiation of ultrarelativistic electrons accelerated at such boundary layers .
the ( simple ) radiative model of a jet boundary layer follows that of so02 .
alternative ( more sophisticated ) models to compute the boundary layer radiation from multizone models ( see , e.g. , @xcite ) will be considered elsewhere .
the most important difference among distinct kinematic models is the radial profile of the lorentz factor @xmath150 across the shear layer .
we have considered three different profiles .
the first one is the usually assumed for large scale jets , namely , a monotonically decaying lorentz factor profile ( as in so02 ) .
the second one is motivated by the results of @xcite and @xcite , namely , a non - monotonic profile where the lorentz factor reaches a maximum at the limit between the jet and the external medium ( ar layer ) . for completeness
, the case of a uniform jet with a sharp edge is considered .
the last case can be regarded as a limit of the ar profile when the maximum lorentz factor in the boundary layer equals the lorentz factor of the jet spine . along the way
, we have conveniently developed the results of @xcite and @xcite to include the effects of dynamically important magnetic fields in the parameter range which is adequate for large scale jets . from such development
we have obtained and idealized kinematic model of anomalous boundary layers in magnetized extragalactic jets .    not surprisingly , we find that the jet - to - counter jet brightness ratio at radio frequencies is larger for jets flanked by ar layers than for uniform or standard layers with the same jet core lorentz factor .
this results from the fact that the differences in the radial profile of the velocity field determine the beaming pattern of the boundary layer emission .
thereby , this effect has to be taken into account when the jet bulk lorentz factor is inferred from the jet - to - counterjet brightness asymmetry .
the differences between ar and uniform boundary layer jets are small in the framework of our simple model . in practice
, these small differences might render their observational distinction very difficult .
comparing these two models with standard radiating boundary layers the differences are much larger .
several independent clues can be used to distinguish ( observationally ) among them :    1 .
if a sufficiently large number of tev photons were detected , it would be possible to distinguish between jets with standard and anomalous layers by looking at the separation of the ssc and ec peaks .
the separation between the two peaks is larger for jets with ar or uniform layers . 2 .
a large jet - to - counter jet brightness ratio ( about two orders of magnitude larger than in standard boundary layer jets ) is expected for jets bounded by ar or uniform layers .
3 .   for large viewing angles
the effective x - ray spectral index is much larger for ar or uniform jets than for jets with standard boundary layers .
there can exist an inversion of the jet and counterjet x - ray spectral indices for jets seen at large viewing angles .
however , for most inclination angles and moderate bulk lorentz factors of the jet core , the effective x - ray spectral index of the jet is much smaller than that of the counterjet . indeed ,
if there is a hint on the bulk lorentz factor of the jet ( e.g. , because superluminal proper motions are detected ) , the jet - to - counterjet spectral x - ray index ratio may tell us whether the radiating layer is uniform , standard or ar .
we note that for large lorentz factors of the jet core ( @xmath236 ) , the largest ratio @xmath237 corresponds to jets with uniform layers , whilst the smaller one should be identified with standard boundary layers . for more moderate core lorentz factor ,
the ratio @xmath238 is larger for jets flanked by ar layers than for jets with uniform or standard layers . in spite of these facts ,
the differences in the effective spectral x - ray index between models with uniform and with ar layers are rather small and , they might render their observational distinction very difficult unless a very careful analysis of the observational data is performed .
we kindly acknowledge the work of jerome ferrand during his internship at the departamento de astronoma y astrofsica .
we also thank l. stawarz for very useful discussions and helpful hints .
maa is a ramn y cajal fellow of the spanish ministry of education and science .
he also acknowledges partial support from the spanish ministry of education and science ( aya2004 - 08067-c03-c01 , aya2007 - 67626-c03 - 01 , csd2007 - 00050 ) .
pm is at the university of valencia with a european union marie curie incoming international fellowship ( meif - ct-2005 - 021603 ) .
the authors thankfully acknowledge the computer resources , technical expertise and assistance provided by the barcelona supercomputing center - centro nacional de supercomputacin .
aloy , m.a . , ibez , j.m . ,
mart , j.m . ,
mller , e. 1999 , supl .
, 122 , 151 aloy , m.a . ,
ibez , j.m . ,
mart , j.m . ,
gmez , j.l . , & mller , e. 1999 , , 523 , l125 aloy , m.a . ,
pons , j.a .
, ibez , j.m . 1999 ,
comm . , 120 , 115 aloy , m.a . , gmez , j.l . ,
ibez , j.m . ,
mart , j.m .
, & mller , e. 2000 , , 528 , l85 aloy , m.a . , mller , e. , ibez , j.m . , mart , j.m . , & macfadyen , a. 2000 , , 531 , l119 aloy , m.a . ,
ibez , j.m . , miralles , j.a . , & urpin , v. 2002 , , 396 , 693 aloy , m.a . ,
mart , j.m . ,
gmez , j.l . ,
agudo , i. , mller , e. , & ibez , j.m .
2003 , , 585 , l109 m.  a. , janka h .-
, mller e. , 2005 , a&a , 436 , 273 aloy , m.a . , & rezzolla , l. 2006 , , 640 , l115 m.  a. , obergaulinger , m. , 2007 , revista mexicana de astronoma y astrofsica , 30 , 96 attridge , j.m . ,
roberts , d.h . , & wardle , j.f.c .
1999 , , 518 , l87 baan , w.a . 1980 , , 239 , 433 begelman , m.c . ,
li , z .- y .
1994 , , 426 , 269 birkinshaw , m. 1991 , , 252 , 505 , r. , aloy , m.  a. , janka , h.  -th . , & mueller , e. 2007 , , 463 , 51 blandford , r.g .
2002 , in eso astrophysics symposia .
m. gilfanov , r. sunyaev , & e. churazov .
springer - verlag , p. 381
blandford , r. 2003 , ap&ss , 288 , 155 blumenthal , g.r . , gould , r.j .
1970 , rev .
phys . , 42 , 237 brown , g.  l. & roshko , a. 1974 , j. fluid mech . , 64(4 ) , 775 brunetti , g. , bondi , m. , comastri , a. , setti , g. 2001 , , 381 , 795 cawthorne , t.v .
, wardle , j.f.c . , roberts , d.h . , & gabuzda , d.c .
1993 , , 416 , 519 chiaberge , m. , celotti , a. , capetti , a. , & ghisellini , g. 2000 , , 358 , 104 crusius , a. , & schlickeiser , r. 1986 , , 164 , l16 dermer , c.d .
1995 , , 446 , l63 de young , d.  s. 1993 , , 405 , l13 dulwich , f. , worrall , d.m . ,
birkinshaw , m. , padgett , c.a . , perlman , e.s .
2007 , , 374 , 1216 earl , j.a . ,
jokipii , j.r . ,
morfill , g. 1988 , , 331 , l91 fendt , c. 2006 , , 651 , 272 fendt , c. , ouyed , r. 2004 , , 608 , 378 ferrari , a. 1998 , ara&a , 36 , 539 ferrari , a. , trussoni , e. , & zaninetti , l. 1978 , , 64 , 43 ferrari , a. , trussoni , e. , & zaninetti , l. 1982 , , 198 , 1065 hanasz , m. , & sol , h. 1996 , , 315 , 355 hardcastle , m.j . , alex&er , p. , pooley , g.g . , & riley , j.m .
1997 , , 288 , l1 hardcastle , m.j . , birkinshaw , m. , & worrall , d.m .
2001 , , 326 , 1499 hardcastle , m.j . ,
croston , j.h .
2005 , , 363 , 649 hardcastle , m.j . , kraft , r.p .
, sivakoff , g.r . , et al .
2007 , , 670 , l81 hardee , p.e .
1979 , , 234 , 47 hardee , p.e .
2007 , , 664 , 26 kataoka , j. , & stawarz , l. 2005 , , 622 , 797 kataoka , j. , & stawarz , l. , aharonian , f. , takahara , f. , ostrowski , m. , edwards , p.g .
2006 , , 641 , 158 komissarov , s.s .
1990 , sval , 16 , 284 laing , r.a .
1996 , in _ ` energy transport in radio galaxies & quasars ' , asp conference series vol .
100 _ , eds .
hardee , a.h .
bridle & j.a .
zensus , san francisco leismann , t. , antn , l. , aloy , m. a. , mller , e. , mart , j. m. , miralles , j. a. , ibez , j. m. 2005 , a&a , 436 , 503 mckinney , j.c . , narayan , r. ( 2007 ) , , 375 , 513 melinani , z. , keppens , r. 2007 , , 475 , 785 mimica p. , aloy m - a .
, mller e. , brinkmann w. 2004 , a&a , 418 , 947 mimica p. , aloy m - a . ,
mller e. , brinkmann w. 2005 , a&a , 441 , 103 mimica p. , aloy m - a . ,
mller e. 2007 , a&a , 466 , 93 mizuno , y. , hardee , p. , hartmann , d.h .
, nishikawa , k .-
, & zhang , b. 2008 , , 672 , 72 ostrowski , m. 1990 , , 238 , 435 ostrowski , m. 1998 , , 335 , 134 ostrowski , m. 2000 , , 312 , 579 owen , f.n . ,
hardee , p.e . , & cornwell , t.j .
1989 , , 340 , 698 perlman , e.s . , biretta , j.a . , fang , z. , sparks , w.b . , &
macchetto , f.d .
1999 , , 117 , 2185 , m. , hanasz , m. , mart , j.  m. , & miralles , j.  a. 2007 , phys .
e , 75 , 056312 pudridtz , r.e . , rogers , c.s .
, ouyed , r. 2006 , , 365 , 1131 rieger , f.  m. , & duffy , p. 2004
, , 617 , 155 romero , r. , mart , j.m . , pons , j.a . ,
ibez , j.m .
, miralles , j.a .
( 2005 ) , j. fluid mech .
544 , 323 siemiginowska , a. , bechtold , j. , aldcroft , t.l . ,
elvis , m. , harris , d.a . ,
dobrzycki , a. 2002 , , 570 , 543 skibo , j.b . ,
dermer , c.d . ,
kinzer , r.l .
1994 , apj , 426 , l23 stawarz ,  . , & ostrowski , m. 2002 , , 578 , 763 ( so02 ) stawarz ,  .
, siemiginowska , a. , ostrowski , m. , & sikora , m. 2005 , , 626 , 120 stawarz ,  .
2007 , private communication .
swain , m.r . ,
bridle , a.h . , & baum , s.a .
1998 , , 507 , l29 turland , b.  d. , & scheuer , p.  a.  g. 1976 , , 176 , 421 urpin , v. 2002 , , 385 , 14 urpin , v. 2006 , , 455 , 779 wilson , a.s . ,
young , a.j .
, shopbell , p.l .
2001 , , 547 , 740 worrall , d.m .
, birkinshaw , m. , & hardcastle , m.j .
2001 , , 326 , l7