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: the chromosphere is in a highly dynamic state , varying on timescales of minutes or less . much of the dynamics in the magnetized regions associated with the magnetic network and plage are dominated by short - lived , jet - like features . a whole range of names has been applied to describe these chromospheric features . at the quiet sun limb , they are traditionally called spicules , where they are observed in  as thin , elongated features that develop speeds of 10 - 30 km  s@xmath1  and reach heights of on average 5 - 9 mm during their lifetimes of 3 to 15 minutes @xcite . their widths ( 0.2 - 1 mm ) and dynamics have , until recently , been very close to the resolution limits of observations . as a result of this limitation , the superposition inherent in limb observations , and the difficulty associated with interpreting h@xmath0 data , the properties of spicules have not been very well constrained , which has led to a multitude of theoretical models ( for a review , see , * ? ? ? the same issues have plagued a solid identification of a counterpart to spicules on the disk . observations on the disk in h@xmath0 reveal dark mottles in quiet sun regions . these dark mottles seem similar to spicules , and usually appear in close association with small flux concentrations ( and bright points ) in the magnetic network that outlines the supergranular cells . there has been a long - standing discussion on whether spicules and mottles are the same phenomena . this discussion has focused on the different velocity distributions of mottles and spicules to conclude that mottles are not the disk counterpart of spicules @xcite . however , several groups have used more recent data to argue for a close relationship between mottles and spicules ( see , e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? dynamic , jet - like features can also be seen in and around active region plage regions . these features are shorter ( 1 - 4 mm ) and shorter - lived ( 3 - 6 minutes ) and appear to form a subset of what have traditionally been called active region fibrils . there are also fibrils that do not show jet - like behavior ( fig . [ fig_context ] ) . these are apparently low - lying , heavily inclined and more static loop - like structures that connect plage regions with opposite polarity magnetic flux . in the remainder of this paper , we will not discuss these low - lying fibrils , but focus on the motions of the shorter dynamic fibrils ( hereafter abbreviated as dfs ) . only recently have dfs become the focus of detailed study , with several papers studying their dynamics , influence on the transition region and their oscillatory properties @xcite . in the past , several authors have reported periodicities and oscillations in spicule observations , usually with dominant periods of order 5 minutes . such periodicities in radial velocity , half - width and line intensity have been suggested based on older observations at the limb , but their interpretation has been difficult due to the complicated line formation of h@xmath0 , line - of - sight superposition and limited spatial resolution @xcite . more recently , @xcite have discovered significant oscillatory power in dfs in active regions using high - resolution data from the swedish 1  m solar telescope @xcite . not all dfs show oscillations or recurrence , but a significant number do . most of the power in especially more inclined dfs resides in periods between 4 and 6 minutes @xcite . similar periods have been found by @xcite in mottles in enhanced network or weak plage . @xcite and @xcite developed a new numerical model to explain the cause of these ( quasi)periodic mottles or dfs . they compared their model to data from the transition region and coronal explorer ( trace * ? ? ? * ) and suggested that inclined magnetic fields allow leakage of normally evanescent photospheric p - modes ( which have dominant periods of 5 minutes ) into the atmosphere where they form shocks that drive dfs . a recent paper by @xcite uses extremely high resolution observations obtained at the sst and significantly more advanced numerical modelling to show that dynamic fibrils are indeed caused by magneto - acoustic shock waves .    in this paper , we directly expand on the work by @xcite . here we perform a more comprehensive analysis of the properties and motion of dfs using the same extremely high spatial and temporal resolution sst dataset of @xcite . in addition , we use advanced 2/3d numerical radiative mhd simulations that include the convection zone , photosphere , chromosphere and lower corona to interpret all of the measurements and correlations found in the sst data analysis , many of which have not been described in @xcite . in   [ sec : instrumentation ] we discuss the details of the observations and instrumentation , including a description of a novel , advanced optical setup that allowed a diffraction limited 88 minute long timeseries in  linecenter . we describe the typical temporal evolution and measure properties of dfs in  [ section : dfdescription ] . this section also deals with intriguing regional differences and correlations in df properties that seem to be related to the oscillatory behavior of dfs . to interpret these findings , we use a numerical model that is described in  [ section : discussion ] . this section also summarizes our observational findings and contains a more detailed description of the new model ( introduced by * ? ? ? * ) that explains our observations . our results strongly suggest that dfs naturally form as the result of chromospheric shock waves driven by convective flows and ( global ) oscillations in the photosphere as suggested by @xcite .    in   [ section : qs ] we discuss in more detail how this model and our observations of active region dfs relate to equivalent quiet sun phenomena . the dynamics of quiet sun mottles have been studied extensively by several groups recently @xcite . while the lack of high spatial and temporal resolution data has made it difficult to properly interpret the findings , this recent work seems to have resolved a long - standing issue @xcite regarding the paths mottles follow during their ascent and descent . @xcite and @xcite use  data to establish that mottles follow parabolic paths . in fact , observations by @xcite indicate that spicules at the limb similarly trace parabolic paths . mottles initially appear in the blue wing of extending in the process to lengths of some 2 to 8 mm @xcite . after reaching maximum extent , they fade from the blue wing and appear in the red wing of  to recede generally along the same path . throughout most of their lifetime they appear as dark , elongated features in  linecenter . the measurements of the parabolic paths traced by mottles have been difficult to interpret . the projected decelerations of the top of the mottle are generally observed to be far less than solar gravity . attempts have been made to interpret such low decelerations as projection effects superposed on ballistic flights . however , the projection angles necessary are typically not compatible with measured doppler shifts @xcite , ruling out a purely ballistic interpretation . in addition , a purely ballistic flight would imply large initial velocities that have not been observed @xcite . these observational findings have been difficult to reconcile with the plethora of models for spicules or mottles @xcite . in   [ section : qs ] we describe the similarities between mottles and dfs in more detail , and discuss the applicability of our model to quiet sun mottles . our results are summarized in  [ section : summary ] . the observations were obtained with the swedish 1-m solar telesope on la palma . the solar optical universal polarimeter ( soup , * ? ? ? * ) was used to obtain narrow - band images in the linecenter of  ( filter fwhm 12.8  pm ) . to compensate for seeing deformations in real - time , the observations were aided by the use of the sst adaptive optics system ( ao , * ? ? ? * ) consisting of a tip - tilt mirror , a bimorph mirror with 37 actuators , a correlation tracker , and a shack  hartmann wavefront sensor . as a further measure to reduce seeing aberrations , the data were post - processed with the multi - object multi - frame blind deconvolution image restoration method ( momfbd , * ? ? ? * ) . this method is complementary to the use of ao since it reduces residual aberrations not corrected by the ao and it corrects the areas in the fov that are not used for the ao wavefront sensing . the science cameras consisted of 3 new , high speed cameras of type sarnoff cam1m100 , with a ccd readout time of only 10  ms and an effective exposure time of 15  ms , resulting in a frame rate of approximately 37 frames - per - second ( fps ) . the setup used was specially built for momfbd processing of the data and consisted of one camera behind the soup filter and 2 `` wideband '' cameras in front of the soup filter but after the soup pre - filter ( a 0.8  nm fwhm interference filter centered at 656.3  nm ) . the light for the wideband cameras was split off the main beam by a 90/10 beamsplitter in front of the soup filter and split again evenly for the 2 wideband cameras , which were set up as a phase diversity pair with one camera at 12  mm defocus on a 45  m focal length . a schematic drawing of the red beam optical set - up is shown in fig . [ fig_optical_setup ] .    to ensure simultaneous exposure , a specially designed chopper with a useful aperture of 22@xmath222  mm and a variable duty cycle of @xmath3 - @xmath4 was used as a shutter . to guarantee correct labeling of all camera frames , a unique exposure number was generated by the shutter electronics and stored on the camera pcs . this scheme guarantees that simultaneously exposed frames can be uniquely identified , a feature that is required for momfbd processing but difficult to achieve using the internal clock of the observing pcs at the frame rates used . the large data rate produced by the cameras ( 72  mb@xmath5s@xmath1 per camera ) was reduced by a factor of 3 using real - time data compression , sent from the camera computers to disk - writing computers through ordinary gigabit network connections , and then passed through fiberchannel links to the raid controllers of an atabeast system . more details of the data acquisition , as well as first results from this unique dataset are described by @xcite . the observations were obtained on 04-oct-2005 and consist of a time series of about 78 minute duration , centered on one of the larger spots in the young active region ar10813 at s7,e37 . the observing angle @xmath6 between the line - of - sight and the local vertical is @xmath7 ( @xmath8 ) . the field - of - view ( fov ) was about 64@xmath264  at a pixel scale of 0063 . this corresponds to a fov of 46,400 x 46,400 km , with a pixel scale of 45.7 km . the diffraction limit of the sst at 6563   is 0.165 , i.e. , 120 km . the images from all cameras were jointly processed with the momfbd image restoration method , yielding 2 restored images for each restoration ( 1 wide band and 1 ) . since the common field - of - view of the cameras within an momfbd set is separately aligned to sub - pixel accuracy , the restored images have near - perfect alignment . sets from 12 exposures were used per restoration ( i.e. , 36 images in total ) resulting in a time series with 3  fps . for this study , we selected every third image , giving a time series with a fixed cadence of 1  s between images . the time series were corrected for diurnal field rotation , aligned and de - stretched . the ( local ) offsets were determined on the wide band images and subsequently applied to the  images . the field of view of the time series contains two small sunspots of noaa active region ( ar ) 10813 . the larger of the spots is split into two parts by a lightbridge as seen in the left panel of fig . [ fig_context ] . the spots are surrounded by areas of plage , which includes some pores to the north of the larger spot . the chromosphere , as imaged in h@xmath0 linecenter ( right panel of fig .  [ fig_context ] ) , shows many almost horizontal fibrils that connect the spots to surrounding opposite polarity pores and plage . many of these horizontal fibrils are associated with the superpenumbra of the larger spot . these fibrils are longer ( up to 20 - 40 ) and more static than the dynamic fibrils that are the focus of this paper . most of the dfs are associated with the plage in the region outlined by x=35 - 60and y=20 - 45 . the chromosphere in this region is dominated by relatively short jet - like features that appear to rise and fall within a matter of minutes . during their ascent and descent , individual dfs seem to follow the same , roughly straight , trajectory . at any given time , scores of dfs are visible in this region of interest . dfs often seem to occur in semi - coherent patches which contain several dfs at any given time . the location and orientation of these patches seems to be governed by the topology of the underlying magnetic plage region . within such a patch , it is sometimes difficult to uniquely identify individual dfs because of the semi - coherent behavior over a few arcseconds : some of these coherent fronts show significant substructure . the region which is dominated by dfs appears somewhat brighter in h@xmath0 linecenter than surrounding regions . this may be because hot coronal loops are connected to this region . such a correlation has previously been described by @xcite who find a correlation between the brightness in h@xmath0 linecenter and that of upper transition region moss emission in fe ix / x 171   ( using trace data ) . bright moss emission occurs at the footpoint regions of hot ( @xmath9 mk ) coronal loops . a full disk magnetogram taken by soho / mdi , centered in the middle of the time series , at 09:39 ut on 4-oct-2005 , shows that the two spots and plage in the region of interest share the same polarity ( see contours of fig . [ fig_context ] ) . most of the associated magnetic flux seems to connect to opposite polarity to the north and west of the spots , for example the pores to the north and the plage outside the field of view to the north . this can be deduced from the potential field extrapolation shown in the left panel of fig . [ fig_potext ] . the field extrapolation seems to correlate reasonably well with features outlined by h@xmath0 in some regions , including our region of interest . however , the extrapolation deviates from the h@xmath0 observations in several regions of the active region to the north and east of the spots , as well as close to the pores to the north of the large sunspot . some df properties vary significantly with position in the field of view . several regions were defined : region 1 is at the edge of a plage region where the magnetic field seems more inclined ( see discussion below ) . region 2 covers a dense plage region ( adjacent to region 1 ) where the fields seem more vertical . the average magnetic field strength , as derived from a full - disk mdi magnetogram ( with 2pixels ) , is of order 50 - 100 mx @xmath10 for region 1 , and 150 - 200 mx @xmath10 for region 2 . the other regions contain fewer dfs and their magnetic topology is not as well defined as in regions 1 and 2 . the right panel of fig . [ fig_potext ] illustrates the location of the two main regions ( blue is region 1 , red is region 2 ) . through visual inspection of a time lapse movie of the linecenter images , we found many fibrils in the plage region in proximity to the two spots near the center of the image . we analyzed dfs that occurred throughout the 78 minute long time series . to analyze the properties of the dfs , we manually choose the general direction of a df ( or set of dfs ) . a guideline is drawn in this direction , and a subset of the data , parallel to this guideline , is extracted for further analysis . this procedure results in detailed movies of individual dfs , as well as xt-plots that show the evolution of the extent of the df along the guideline as a function of time . a total of 257 dfs were analyzed in this fashion .    in figs . [ fig_spiclife_1 ] , [ fig_spiclife_2 ] , and [ fig_spiclife_3 ] we present the evolution of three typical dfs . figure  [ fig_spiclife_1 ] shows a df in region 2 . the six panels show that the df rises rapidly to a maximum length in roughly 130  s and recedes along the same path some 260  s after first becoming visible . it is darker than its surroundings , especially at maximum extension , and seems to follow a straight path in the direction that is shared by all surrounding features . the df has some internal structure , most clearly visible at @xmath11  s and @xmath12  s. the maximum extent is relatively modest at some 2 . despite the substructure , the df is thin with a width of some @xmath13 . it is also well isolated even though a number of other dfs are visible in the immediate vicinity . the length of the df as a function of time is well described by a parabola , as shown in the xt-plot in the top panel of fig . [ fig_spicxt ] . this is the case for the vast majority of the dfs that do not suffer from superposition and where the df guideline was chosen with a correct angle , i.e. where the df motion runs along the guideline . all 257 dfs used in the analysis below meet these conditions . each individual df is therefore fitted with a parabola and the deceleration , maximum velocity ( on ascent or descent ) , maximum extent and duration were derived . to determine the beginning and end points of the parabola , we use the following method . the steep parabolic ascent and descent of a typical df ( and the very high cadence of our data ) allows for a robust determination of the beginning and end time of each df . the intersection ( in the xt-plot ) between the parabolic path and a horizontal line at the beginning time or end time is used to determine the location of the `` root '' of the df at , respectively , the beginning and end of its lifetime . we then use these points in the `` xt''-plot to fit a parabola to the path followed by the df . the prevalence of parabolic paths is illustrated in all three panels of fig . [ fig_spicxt ] , where a wide variety of dfs are evident . the df studied in this case ( top panel ) follows its fitted parabola very well with a constant deceleration given by 216  m  s@xmath14 , about two thirds of solar gravity ( 274  m  s@xmath14 ) , a maximum velocity of 27  km  s@xmath1 , and a maximum extent of some 1800  km . the properties for this df and all others discussed in this paper are measured along the projected path . an example of a much wider df in region 2 is shown in fig . [ fig_spiclife_2 ] . this df has extensive internal structure and rises over a front of about @xmath15 , a width that is largely maintained throughout the descending phase of its relatively short lifetime of 188  s. while rising , the df width is not constant along its axis but the df is wider at the bottom . it is typical that such wide dfs do not rise and fall as a rigid body but rather show phase and amplitude variations in velocity at various positions away from the axis . it is also worth noting that this df clearly becomes darker as it rises and brightens somewhat as it falls back along the same path it rose along . the xt-plot in the middle panel of fig . [ fig_spicxt ] reveals that this df also follows a path that is very close to a parabola with a deceleration of 286  m  s@xmath14 , a maximum velocity of 27  km  s@xmath1 , and a maximum extent of 1300  km . a longer lived df from region 1 with a lifetime of some 290  s is shown in fig . [ fig_spiclife_3 ] . this df is located in a region where  images indicate that highly inclined magnetic field lines connect to the region of opposite polarity plage and pores to the west , in the top right corner of fig . [ fig_context ] . the df is wide , roughly 1 , and relatively long . it is also superposed on another df that is evident in the bottom panel of fig . [ fig_spicxt ] . detailed examination of the movie accompanying fig . [ fig_spiclife_3 ] shows the separate evolution of both dfs . the superposition makes it difficult to determine whether these dfs are formed at the exact same location . residual effects of image distortion from atmospheric seeing can cause the jagged appearance of these xt profiles . despite such distortions and the effects of superposition , a good fit for this df can still be obtained with a deceleration of 200  m  s@xmath14 , a maximum velocity of 28  km  s@xmath1 , and a maximum extent of some 2000  km . the properties of the parabolic fits to the xt-plots of 257 dfs were calculated . visual inspection of df time lapse movies confirmed the proper alignment of the dfs and quality of the parabolic fits . the upper left panel of fig . [ fig_hist ] shows that the dfs suffer a deceleration of on average 146  m  s@xmath14  with a significant spread ranging from 40 to 320  m  s@xmath14 and a standard deviation of 56  m  s@xmath14 . this is significantly less than solar gravity . note that the median of all df properties is illustrated and described in the caption of fig . [ fig_hist ] , whereas we enumerate average values in the text . also note that these deceleration values are twice as high as the values reported in @xcite . the latter were mistakenly divided by two . the maximum df lengths are on average @xmath16  km , ranging from 400 to @xmath17  km with a standard deviation of 620  km , as shown in the upper right panel of fig . [ fig_hist ] . the lengths found for dfs are smaller than the widely reported values of @xmath18 to @xmath19  km for spicules at the limb @xcite , but fit well with the values reported in @xcite . on average , the maximum velocity ( on either ascent or descent ) is 18  km  s@xmath1 ranging from 8 to 35  km  s@xmath1  with a standard deviation of 6  km  s@xmath1 , as shown in the middle left panel of fig . [ fig_hist ] . the maximum velocity is never observed to be lower than 8 km  s@xmath1 . df durations are shown in the middle right panel of fig . [ fig_hist ] . the durations vary from 120 to 650  s , with an average of 290  s , and a standard deviation of 85  s. the median duration is lower at 250  s , because of the significant tail of longer duration dfs . in addition to the parabolic fits , we have also measured df widths . to determine the width of a df , we calculate xt-plots ( with x along the direction of the guideline ) for a range of locations that are each offset by one pixel in the direction perpendicular to the guideline . we then find the central axis of the df through visual inspection of the xt-plots : the location with the best - defined xt-plot is chosen as the central axis . starting at this location we find the edges of the df in the direction perpendicular to the df axis . we search for the two positions of maximum intensity gradient furthest from the axis where the df is still visible in the xt-plots . this is done for two different locations along the df , at 50% and at 80% of the maximum extent for a specific time . the width is measured in this fashion at three different times in the df lifetime , at 25% , 50% ( maximum extent ) , and 75% of the df duration . although there is a large spread , we find that the widths of dfs do not vary much with height nor time . the lower left panel of fig . [ fig_hist ] includes all the widths measured and shows a mean of 340  km with a standard deviation of 160  km . while the widest df was measured to be 1100  km , most dfs have widths between the difraction limit of 120  km and 380  km , with a significant tail of wider dfs extending to 700  km . based on the shape of the histogram it seems likely that we are not resolving all dfs . however , a significant fraction is clearly resolved . visual inspection shows that the tops of most dfs are very sharply delineated ( see , e.g. , fig . [ fig_spicxt ] ) . furthermore , inspection of xt-plots shows that dfs often are outlined by a bright ring , as shown in fig . [ fig_bright_ring ] . to determine how sharp the transition at the top of the df is , we calculated the distance along the df axis between the location of maximum intensity gradient and the location of maximum intensity in the bright ring above the df . we refer to this measure as the thickness of the bright ring . it is unclear what causes the bright ring . perhaps it is related to the steep transition from chromospheric to coronal temperatures ? detailed radiative transfer calculations will be necessary to determine what causes the bright ring . regardless of the uncertainty in interpretation , we measured the thickness of the bright ring for all 257 dfs and found an average thickness of 240  km with a standard deviation of 120  km . over half are less than 200  km and many are unresolved , as shown in the lower right panel of fig . [ fig_hist ] .      despite a large spread in properties we find clear differences between df statistics in the two regions defined in the right panel of fig . [ fig_potext ] . on average , dfs in region 1 show lower decelerations , greater lengths , slightly higher maximum velocities , and longer durations than in region 2 ( fig . [ fig_hist ] ) . we also found a slight tendency for dfs to be wider in region 2 ( lower left panel of fig . [ fig_hist ] ) , but found no difference in transition region thickness between the two regions . visual inspection shows that region 2 does show more dfs that are wide with significant substructure , such as the df illustrated in fig . [ fig_spiclife_2 ] .    to further illustrate the difference between region 1 and 2 , maps of the location of dfs , color coded to show various properties , are shown in fig . [ fig_regionalmaps ] . as can be seen , there are clear regional differences in deceleration , duration , and length . there is a prevalence of smaller decelerations in region 1 with typical values around 100  m  s@xmath14 , versus roughly 200  m  s@xmath14  in region 2 . even more clearly , the durations in region 1 are on average about 5  minutes , whereas we find 3  minutes in region 2 . similarly , df lengths in region 1 are typically around @xmath20  km , whereas shorter dfs with lengths around @xmath21  km dominate region 2 . the spatial pattern of maximum velocities is not as clear with a large spread of velocities in both regions .    in principle some of these differences could be caused by projection effects due to the angle between the df axis and the line of sight vector . a df that propagates along a direction that is parallel to the direction of the line of sight would not be observed as a moving feature in these  movies . to measure the absolute extent , velocity and deceleration along the path of the df , the line of sight needs to be perpendicular to the axis of the df . note that the df durations are independent of the line of sight as they do not suffer from these projection effects .    to exclude the possibility that the regional differences in deceleration , velocity and length are caused by projection effects , we try to estimate the absolute direction of propagation of each df . assuming that dfs propagate along the magnetic field , their direction can be estimated from the potential field extrapolation shown in the left panel of fig . [ fig_potext ] . to estimate the effect of projection , we calculate the angle between the line of sight vector and the direction of the magnetic field at a reference height of @xmath22  km above the photosphere . this height was chosen to be very roughly the sum of the top height of the quiescent chromosphere and the typical df length . the angle between extrapolated field lines and the line of sight vector was used to correct the properties of each df for projection . histograms for these corrected values are shown in fig . [ fig_hist_los ] . not surprisingly , all corrected values are larger than the observed values . in addition , the corrected values show the same spatial differences between region 1 and 2 as the uncorrected values . to estimate how reliable the potential field extrapolation is , we measured the mismatch between the direction of the df axis and the direction of the magnetic field as projected onto the surface , and found it to be of order 20  30 degrees . a mismatch of this magnitude can lead to significant errors in calculating the projection angle if the df axis is nearly parallel to the projection of the line of sight vector onto the surface . this means that the corrected values should be considered with caution , especially since the calculations are also based on various assumptions whose validity is unclear . the poor resolution of the mdi magnetogram , the possible non - potentiality of the active region , and the uncertainty in the actual df height with respect to the chosen reference height , all could lead to significant errors in the estimated projection angle . for example , the potential field extrapolation clearly fails in the top part of region 1 where nearly horizontal fields are predicted . such horizontal fields are imcompatible with the visual appearance of  linecenter images . in addition , since horizontal fields at this location are almost parallel to the line of sight vector , they lead to very large and obviously erroneous corrections of the local df parameters .    despite these reservations , it is very difficult to conceive of a magnetic topology with an orientation to the line of sight that would remove the regional differences for all derived parameters . in addition , the regional differences in duration ( which are independent of projection effects ) as well as the correlations found between various parameters in the following subsection , strongly suggest that regional differences are not caused by projection effects . scatterplots of various df properties reveal intriguing correlations . the most striking correlation , also presented in @xcite , is that between the deceleration and the maximum velocity of the dfs , in which a linear relationship is evident ( upper left panel of fig . [ fig_scattergood ] ) . the larger the deceleration the df suffers from , the higher the maximum ( initial or final ) velocity it shows . note that there are , again , significant differences between region 1 and 2 : the dfs in region 2 typically have larger decelerations for a given maximum velocity . a linear fit to the observed relationship between the deceleration and velocity shows that both regions have almost identical slopes . another interesting correlation is that between the maximum length and the duration of a df ( upper right panel of fig . [ fig_scattergood ] ) . the longer dfs tend to have longer lifetimes . as expected from the histograms in the upper and middle right panels of fig . [ fig_hist ] , the dfs in region 1 gravitate towards the upper right of this scatterplot . we also find that the maximum velocity and maximum length of dfs is well correlated ( lower right panel of fig . [ fig_scattergood ] ) : dfs with higher maximum velocity tend to be longer . however , there is a difference in the slope in regions 1 and 2 . a linear fit for region 1 reveals a slope of 87 km per km  s@xmath1 , and 60 km per km  s@xmath1  for region 2 . dfs in region 1 are longer for the same maximum velocity . the deceleration of dfs shows a somewhat less clear correlation with the df duration ( lower left panel of fig . [ fig_scattergood ] ) . the longest - lived dfs typically suffer from the lowest deceleration . this correlation shows a large spread , and is not quite linear for the total population of dfs . it seems that each region may have its own linear relationship . dfs in region 1 show a more steep dependence of duration on deceleration , whereas those in region 2 show a smaller range of durations , for a range of decelerations that is similar in extent to that of region 1 .    based on the correlations described in the above , one would perhaps also expect a correlation between the deceleration and maximum length , or between the maximum velocity and duration . however , such correlations are not clear from the observations , as demonstrated in fig . [ fig_scatterbad ] . the large scatter and different slopes that are evident in the plots that do show correlations may be partially responsible for the lack of clear correlation between deceleration and maximum length or the maximum velocity and duration . however , it is also clear that there are very significant regional differences in fig . [ fig_scatterbad ] , with each region occupying a different area in the scatterplots , and in some cases showing different slopes . it seems probable that such regional differences also are responsible for the lack of clear correlations . when we take the line - of - sight correction into account , none of these correlations , or the lack of correlation , change qualitatively . evidently , the correlation between deceleration and velocity , or that between maximum velocity and length can not change significantly because of projection effects , since both parameters would be corrected by the same factor . the fact that the other correlations , between duration on the one hand , and either deceleration , length or velocity , do not change qualitatively indicates that projection effects do not dominate our dataset .    the same arguments also apply to any kind of selection effects that are caused by the fact that dfs with very large line - of - sight velocities could disappear from or be invisible in our relatively narrow ( 128 m )  linecenter bandpass . since dfs are much wider spectral features ( with widths of order 1000 m , see , e.g. * ? ? ? * ) than the width of the bandpass , this means that only features with line - of - sight velocities of order 20 or more km / s would be invisible or only partially visible in the linecenter timeseries . while in a few cases we can see some fading of the df during its descent , it is still possible to trace the top of the df to its end point . given these arguments , the passband effect can only change the derived correlations marginally .    note that all of the measurements on which the correlations are based , depend on the identification of a single df rising above the chromospheric background . the latter can be quite noisy and sometimes less well defined . filtergrams taken in the wings of h@xmath0 might be helpful for a more exact determination of the root location of each df , but such filtergrams were unavailable in this dataset . this means that the measured length and maximum velocity could in fact be somewhat higher than reported here . however , given the rapid rise and fall of dfs , the measured duration is quite insensitive to this effect . since we are measuring the acceleration by fitting parabolas to the height profile , the acceleration is also insensitive to the position of the chromospheric background . overall , this implies that the relative uncertainty in determining the roots of dfs are unlikely to qualitatively change the correlations reported here . the interpretation of these correlations , or the lack thereof , will be discussed at length in  [ section : discussion ] . the regional differences in df durations evident in the middle right panel of fig . [ fig_hist ] , with 3 minute lifetimes dominating region 2 and lifetimes around 5 minutes more dominant in region 1 , suggest a role for chromospheric waves or oscillations in the formation mechanism of dfs . oscillations in the chromosphere are typically not steady harmonic waves , but wave trains of finite duration . for this reason we perform a wavelet analysis of the  data , since it allows for localization in time of periodic signals . the mother wavelet we use is the complex valued morlet wavelet ( @xmath23 ) which consists of a plane wave modulated by a gaussian , an appropriate shape for wave trains . to estimate the statistical significance of the wavelet power spectra ( 95 % confidence interval ) , we compare them to theoretical spectra for white noise , following @xcite and @xcite . taking into account the cone of influence @xcite , we calculate for each superpixel the number of wavepackets in the  timeseries for which the wavelet power is significant . this is done for a range of wave periods from 150 to 600 seconds . a superpixel is defined as the average of 4x4 original 0.063  pixels . note that the original data is oversampled by a factor of almost 3 with respect to the diffraction limit ( 120 km ) of the sst at the wavelength of . the upper left and upper right panels of fig . [ fig_osc_panels ] show , for each location of the lower left panel , the number of significant wavepackets ( during the  timeseries ) for waves with a period of 180 s and 300 s respectively . the lower right panel shows which wave period dominates , i.e. , has the most wavepackets , at each location of the field of view . from this figure , we find that many locations show significant oscillatory power throughout the timeseries . the dominant period of this power varies across the field of view . clearly visible are the dominant 3 minute oscillations in the sunspots around x=27  and y=37  and at x=45  and y=15 . in addition , the dense plage region which we have designated as region 2 , shows a preponderance of 3 minute oscillations . in contrast , region 1 , adjacent to the dense plage region , is dominated by oscillations with long periods around and above 5 minutes . in this region , the magnetic field is more inclined than in the dense plage region . the longer , less dynamic and more horizontal fibrils , such as those in the superpenumbra of both spots , seem to be dominated by oscillations with periods longer than 10 minutes . it is interesting to note that the dominant period in region 2 corresponds exactly to the df duration we observe in region 2 . apparently , dfs typically live for about 3 minutes in this region , which is also the dominant periodicity of the oscillations . likewise , region 1 is dominated by 5 minute oscillations , and contains dfs with lifetimes of order 5 minutes . since the oscillation measurements are based on the same data , it seems clear that some dfs are periodic or recurring . this was also described by @xcite and suggests a relationship between chromospheric waves and dfs , which will be further explored in  [ section : discussion ] . most dfs follow parabolic paths with a symmetrical ascending and descending phase . these paths are characterized by a large initial velocity , usually of order 15 - 20 km  s@xmath1  ( middle left panel of fig .  [ fig_hist ] ) , that decreases linearly with time until the df has retreated completely . at the end of its life , a df s downward velocity is roughly identical to the initial upward velocity . the linear decrease with time of the velocity , i.e. , the deceleration , is usually between 120 and 280 m  s@xmath14  ( upper left panel of fig . [ fig_hist ] ) . there is a clear linear relationship between the deceleration and the maximum ( down- or upward ) velocity of the df ( upper left panel of fig . [ fig_scattergood ] ) . dfs reach maximum lengths of order 1,000 to 2,000 km ( upper right panel of fig . [ fig_hist ] ) during their lifetimes of 3 to 8 minutes ( middle right panel of fig . [ fig_hist ] ) . the duration and maximum length are well correlated ( upper right panel of fig . [ fig_scattergood ] ) , which indicates that uncertainty about the line - of - sight correction may not play an important role in the regional differences we find . we resolve more than half of the dfs we observe . the average thickness is 240 km ( lower left panel of fig . [ fig_hist ] ) , with a range between the diffraction limit ( 120 km ) and 1000 km . some dfs occur over a wider spatial area , whereas others are very thin and stay thin throughout their lifetime . many of the wider dfs have variable widths as a function of time . these give the visual appearance of a wave front with varying phase speeds or amplitude in the direction perpendicular to the df axis . there is often a brightening above the df , which can be clearly seen in xt-plots ( fig . [ fig_bright_ring ] ) . this brightening might be a sign of increased temperatures at the top of the df . the thickness of this transition is not resolved and is less than 120 km in more than half of the dfs studied ( lower right panel of fig . [ fig_hist ] ) . dfs have significantly different properties in two regions within our field of view . region 2 contains dense plage , i.e. , a collection of strong magnetic field concentrations in which the magnetic field is generally more vertical . region 1 is located at the edge of this dense plage region , where the magnetic field is more inclined , as it connects to relatively nearby opposite polarity plage . the dfs in the dense plage region are shorter ( @xmath24 km ) , have higher decelerations ( @xmath25 m  s@xmath14 ) , slightly lower velocities ( @xmath26 km  s@xmath1 ) and shorter durations of about 3 minutes ( figs . [ fig_hist ] and [ fig_regionalmaps ] ) . in contrast , the dfs in the region where the field is more inclined are longer ( @xmath27 km ) , have lower decelerations ( @xmath28 m  s@xmath14 ) , slightly higher velocities ( @xmath29 km  s@xmath1 ) and longer durations of about 5 minutes ( figs . [ fig_hist ] and [ fig_regionalmaps ] ) . in addition , df lifetimes correlate well with the dominant periodicities observed in both regions : 3 minute oscillations dominate  in the dense plage region , whereas 5 minute oscillations dominate in the region with more inclined field ( fig . [ fig_osc_panels ] ) . our measurements of dfs and oscillatory power strongly suggest that chromospheric waves drive dfs upward . these chromospheric waves are generated by convective flows and ( global ) oscillations in the photosphere and convection zone . the magnetic field concentrations provide a channel for these disturbances to propagate up through the chromosphere . during their path upwards , they shock and drive chromospheric plasma upwards , forming dfs in the process . this mechanism was originally suggested in rudimentary form by @xcite . some models have further studied the effects of the convective part of the photospheric driver on the chromosphere in the so - called rebound - shock model @xcite . more recently , @xcite included both the convective driver and global oscillations in varying magnetic topologies to directly describe their ( lower resolution ) observations of dynamic fibrils . the latter model s predictions agree well with the current findings . this scenario is confirmed by analysis of advanced numerical simulations ( see @xcite and next subsection ) in which dynamic fibril - like features with properties identical to those in our data occur . from the simulations , it is clear that the observed parabolic paths indicate that dfs are shock driven . when a shock hits the transition region ( i.e , the top of a df ) , the transition region is catapulted upwards at a velocity that exceeds the sound speed in the chromosphere ( @xmath30 km  s@xmath1 ) . shock waves generally have velocity profiles in the form of n- or sawtooth shapes ( e.g. * ? ? ? thus , a plasma parcel passing through a shock wave will first experience a sudden kick in velocity - the shock itself  and thereafter a gradual linear deceleration as the shock recedes . assuming a shock wave with amplitude or shock strength of @xmath31 and linear deceleration @xmath32 @xmath33 it is trivially verified that a plasma parcel passing through the shock will describe a parabolic path ; integrating gives a parcel position that moves as @xmath34 df formation from chromospheric shocks is compatible with the fact that we do not find any dfs with initial ( maximum ) velocities of less than 10 km  s@xmath1 . waves with smaller amplitudes will not form shocks and hence dfs . the observed lower limit to the maximum velocities provides an estimate of the sound speed in the upper chromosphere ( 10 km  s@xmath1 ) . its value is compatible with chromospheric temperatures from solar models . more generally , the observed maximum velocities provide a robust estimate of the mach number of chromospheric shocks . our observations indicate that shocks with mach numbers of 1 through at least 3.5 are ubiquitous in the magnetized chromosphere , and drive dfs . analysis of numerical simulations @xcite further shows that several different effects influence the shock formation and evolution so that the observed correlations ( e.g. , between deceleration and velocity ) and regional differences follow naturally . the oscillatory power spectrum of the photospheric driver is dominated by power at 5 minutes . the chromosphere acts as a filter to the input from below , filtering out waves with periods that are longer than the local acoustic cutoff period . in general , the acoustic cutoff period depends on the inclination of the magnetic field lines to the vertical @xcite . under conditions of vertical magnetic field ( e.g. , region 2 ) , this means that the chromosphere is dominated by oscillations and waves with periods at the acoustic cutoff of 3 minutes . as a result , the dfs formed in the dense plage of region 2 are driven by shocks with periods of 3 minutes and thus have lifetimes of order 3 minutes . in contrast , the inclined field , such as that found in region 1 , increases the acoustic cutoff period . this allows the photospheric waves with periods of 5 minutes to leak into the chromosphere , develop into shocks and form dfs with similar lifetimes . since most of the oscillatory power in the photosphere occurs for periods of 5 minutes , the driver in inclined field regions is generally stronger than in regions with more vertical field . this effect is enhanced by the fact that , in dense plage regions , the amplitude of convective flows and global oscillations is generally reduced at the photospheric level . dfs that form along inclined field lines also experience only the component of gravity along the magnetic field , so that decelerations are generally lower in the inclined field region .    the observed correlation between deceleration and maximum velocity can be understood in terms of shock wave physics and the fact that dfs are driven by single shocks . for a given wave period , n-shaped shock waves will show a steeper decline with time in velocity for a greater shock strength . in other words , dfs with higher maximum velocity , i.e. , those that are driven by stronger shocks , will show greater deceleration , for a given df lifetime . it is the combination of these effects that explains the observed regional differences and correlations . region 2 contains more vertical field , so that only 3 minute power drives the chromospheric waves and shocks . given the photospheric power spectrum and the reduced power in dense plage regions , this leads to lower amplitudes , i.e. , lower maximum velocity . in addition , the short lifetimes and more vertical field lead to higher decelerations , and shorter lengths . conversely , region 1 contains more inclined field , allowing the full photospheric peak power at 5 minutes to leak into the chromosphere to drive dfs . as a result , the velocity amplitude is greater , and the longer lifetimes and less vertical field lead to lower decelerations and longer lengths . these findings are confirmed when comparing the observations with recent advanced numerical radiative mhd simulations that include the convection zone , photosphere , chromosphere , transition region and corona @xcite . the equations of mhd are solved using an extended version of the numerical code described in @xcite and in more detail by nordlund & galsgaard at http://www.astro.ku.dk/@xmath35kg . the code functions by using a sixth order accurate scheme to determine partial spatial derivatives . the equations are stepped forward in time using the explicit 3rd order predictor - corrector procedure by @xcite , modified for variable time steps . in order to supress numerical noise , high - order artificial diffusion is added both in the forms of a viscosity and in the form of a magnetic diffusivity . thermal conduction along the magnetic field in the corona and transition region is treated by operator splitting and using a multi - grid solver . the radiative flux divergence from the photosphere and lower chromosphere is obtained using the method of @xcite based on group mean opacities as modified by @xcite to account for scattering . for the upper chromosphere , transition region and corona we assume effectively thin radiative losses . for the middle and upper chromosphere we parametrize the optically thick radiative losses in strong lines and continua from hydrogen and singly ionized calcium by using escape probabilities calculated in 1d with the code of @xcite . further description of the 3d code may be found in @xcite . the 2d simulations described here are run on a grid of @xmath36 points spanning @xmath37  mm@xmath38 . in height , these models cover a region from 2  mm below to 10  mm above the photosphere . the grid has uneven spacing in the vertical , @xmath39 , direction with a grid size of 30  km in the photosphere and chromosphere becoming gradually larger in the corona . we have also run 2d calculations on coarser grids of @xmath40 and @xmath41 points , as well as 3d runs of @xmath42 points spanning @xmath43  mm@xmath44 , with results similar to those reported here .    at the start of a given simulation , a potential magnetic field is added to a model of the solar atmosphere . the resultant model is then allowed to evolve in time until a quasi - steady state is achieved . this takes on the order of a half hour solar time . the models are convectively unstable due the radiative losses in the photosphere . this results in a temperature structure similar to that derived in semi - empirical solar models . the average temperature at the bottom boundary is maintained by setting the entropy of the fluid entering the computational domain . the bottom boundary , based on characteristic extrapolation , is otherwise open , allowing fluid to enter and leave the computational domain as required . the magnetic field at the lower boundary is advected with the fluid . to prevent coronal cooling the upper temperature boundary is set to evolve towards a temperature of @xmath45  k on a timescale of 10  100  s. ( in the 3d models coronal temperatures are maintained by self consistent magnetic dissipation . ) a snapshot of the temperature structure in the lower part of the model including the convection zone , photosphere , and chromosphere is shown in fig . [ sim_snapshot ] . waves generated by convection propagate through the chromosphere , steepen and form shocks roughly 1  mm above the photosphere . in regions above the level where plasma @xmath46 the propagation is modified by the magnetic field and the waves are subjected to reflection , refraction and mode coupling as described by @xcite . of interest to us here are the shocks that propagate along field lines that enter the coronal plasma as indicated in fig  [ sim_snapshot ] . these shocks lift and accelerate the transition region as they propagate into the corona . plotting the motion of the transition region along one such field line as a function of time , e.g. a xt plot , reveals that the position of the transition region describes a parabola ( fig . [ sim_xt ] , compare with fig . [ fig_spicxt ] ) . in the left panel of fig . [ sim_xt ] the logarithm of the temperature is shown , the right panel shows the velocity . it is worth noting that except for a brief moment in the early ascending phase the simulated df shows both upflowing and downflowing velocities along its length at a given time . during descent , the df has only downflowing plasma . note also that near its base the simulated df shows downflows for a major portion of its lifetime . we have measured the transition region motion along several field lines in several ` df - like ' events and fitted the motion with parabolas . the resulting correlation between measured maximum velocities and decelerations ( fig . [ sim_dec - maxv ] ) show the same trend as found in the observed dfs ( upper left panel of fig . [ fig_scattergood ] ) : large decelerations are correlated with large maximum velocities . the range in maximal velocities and decelerations is approximately the same as that found in the observed dfs , and the slope in the correlation is about the same . inspection of fig . [ sim_xt ] also indicates that the simulated dfs have the same range of heights and durations that the observed dfs display . higher resolution 3d simulations with varying magnetic topologies will be necessary to explain the range of df widths , the internal df substructure , and what determines the sometimes sheet - like structure of some dfs . such simulations will also help to fully explain the range of periods and lifetimes that we observe .      how do these findings relate to quiet sun phenomena such as mottles and spicules at the limb ? there are many similaties between mottles and dfs . both phenomena appear as highly dynamic , dark features in the wings and core of . they are associated with magnetic flux concentrations that can be observed as bright points in the blue wing of . more importantly , like dfs , mottles also follow parabolic paths . in an excellent analysis of the dynamics of inclined mottles ( which they refer to as disk spicules ) @xcite find that :    * mottles seem to be generated several hundred kilometers above the photosphere , * mottles undergo real mass motion of order 10 - 30  km  s@xmath1 , extending up to 2 - 10 mm during their lifetime of 2 - 15 minutes , * most mottles are observed to ascend and descend , * mottles in weak plage , or the interior of network cells are more vertical and shorter than the more inclined quiet sun mottles , * the appearance of mottles is varied , they can be thick , thin , tapered or not , curved or straight , * some mottles show downward motion close to their base during the ascending phase , * the decelerations observed in mottles are too small to be consistent with a purely ballistic flight ( at solar gravity ) and the observed doppler shifts , * the largest doppler signal appears at the beginning of the ascending phase , and at the end of receding phase , * the observed velocity profiles are compatible with impulsive acceleration and constant deceleration afterward . * the maximum upward velocity is usually identical in amplitude to the maximum downward velocity , * there is a linear relationship between the duration and the maximum length of mottles , * there is a correlation between the average velocity and maximum length , * nearby mottles are often found to be driven by a common disturbance . these findings agree strikingly well with our observations and modelling of dfs , with only a few discrepancies . the lifetimes and maximum lengths of the mottles @xcite analyzed are typically longer than those of the dfs we describe . we believe this is in part a selection effect , as we chose to focus only on highly dynamic and shorter features in our dataset . we note that our dataset also contains a large number of highly inclined fibrils that are longer lived and less dynamic than the dfs ( see , e.g. , fig . [ fig_osc_panels ] ) . in addition , preliminary analysis of high - quality sst data of quiet sun mottles indicates that the most inclined mottles are on average longer , but that many of the mottles actually show lengths ( 1,000 - 5,000 km ) that are similar to those of dfs . until further detailed analysis , it can however not be excluded that dfs associated with plage are on average somewhat shorter in length and lifetime than quiet sun mottles . preliminary analysis of our numerical simulations suggests that these differences could be related to large scale differences in magnetic topology . further work is necessary to establish the exact relationship between mottles and dfs . however , it should be noted that in addition to these striking similarities with mottles , recent observations of limb spicules @xcite show clear parabolic paths with decelerations and maximum velocities of the same order as those reported for mottles @xcite and those we find for dfs . older observations of limb spicules have also shown evidence for correlations between velocities and maximum lengths @xcite , and between duration and height @xcite . these correlations are similar to the ones we find for dfs . in addition , @xcite has reported on collective behavior of spicules or mottles over several thousand kilometers .    as discussed also by @xcite , these similarities strongly suggest that chromospheric shock waves cause significant excursions of the upper chromosphere in both active region and quiet sun . such a scenario was proposed by @xcite using sst and trace data . further detailed analysis of our numerical simulations can help shed light on several unresolved issues , such as the role of reconnection in forming quiet sun jets , or whether the observed height of limb spicules can be explained by the mechanism proposed here . @xcite suggest that it is possible that limb spicules consist of two populations : jets that are driven by shocks ( as described in this paper ) , and jets caused by reconnection . the latter jets could form a subset that on average is taller than the shock driven jets , and perhaps be a part of a continuous spectrum of reconnection jets that includes surges , macrospicules and   upflow events @xcite . these taller jets may be the features that have uv counterparts ( e.g , observed by sumer ) , which @xcite calls `` spicules '' . the combination of high - resolution chromospheric data and advanced numerical simulations strongly suggests that dfs are formed by upwardly propagating chromospheric oscillations / waves . these waves are generated in the convection zone / photosphere as a result of global p - mode oscillations and convective flows . the magnetic topology acts as a filter , so that only waves of certain periods can propagate into and through the chromosphere , forming shocks along the way , up to regions where h@xmath0 is formed . * dfs follow a parabolic path * dfs in dense plage regions are shorter , slower , undergo larger deceleration and live shorter * dfs in regions adjacent to dense plage regions are longer , faster , undergo less deceleration , and live longer * dfs in dense plage regions show periodicities around 180 s * dfs in regions adjacent to dense plage regions show periodicities around 300 s * the velocity and deceleration of all dfs is well correlated : the higher the velocity , the more deceleration * this correlation is reproduced well in numerical simulations where dfs are driven by shocks caused by upward propagating waves in the chromosphere . * dfs in the simulations follow parabolic paths as well . * the simulations reproduce the observed values of maximum velocities , decelerations and lengths well .    from the results describe above and by @xcite , it seems clear that in active regions , most dynamic fibrils are formed by chromospheric shocks driven by convective flows and oscillations in the photosphere . , g.  b. , bjelksj , k. , korhonen , t.  k. , lindberg , b. , & petterson , b. 2003 , in innovative telescopes and instrumentation for solar astrophysics . keil & s.v . , vol . 4853 , 341350
Abstract:

Answer: we present unprecedented high resolution h@xmath0 observations , obtained with the swedish 1-m solar telescope , that , for the first time , spatially and temporally resolve dynamic fibrils in active regions on the sun . these jet - like features are similar to mottles or spicules in quiet sun . we find that most of these fibrils follow almost perfect parabolic paths in their ascent and descent . we measure the properties of the parabolic paths taken by 257 fibrils , and present an overview of the deceleration , maximum velocity , maximum length and duration , as well as their widths and the thickness of a bright ring that often occurs above dynamic fibrils . we find that the observed deceleration of the projected path is typically only a fraction of solar gravity , and incompatible with a ballistic path at solar gravity . we report on significant differences of fibril properties between those occurring above a dense plage region , and those above a less dense plage region where the magnetic field seems more inclined from the vertical . we compare these findings to advanced numerical 2d radiative mhd simulations , and find that fibrils are most likely formed by chromospheric shock waves that occur when convective flows and global oscillations leak into the chromosphere along the field lines of magnetic flux concentrations . detailed comparison of observed and simulated fibril properties shows striking similarities of the values for deceleration , maximum velocity , maximum length and duration . we compare our results with observations of mottles and find that a similar mechanism is most likely at work in the quiet sun .