the tmt project conducted an extensive site monitoring campaign to identify the most suitable site to host the observatory . knowing the vertical distribution of the optical turbulence strength ( ts ) , as expressed by the structure constant of the refractive index @xmath0 , in earth s atmosphere  the turbulence profile ( tp or @xmath1 )  is of great importance for modern astronomical observing techniques .
especially adaptive optics ( ao ) instrumentation in its various forms relies in its inital design on a good assumption of the typical altitude distribution of the ts above the telescope .
ao will become an important observing technique of the future extremely large telescopes ( elts ) and it is therefore essential to obtain measurements of the tps above the elt candidate sites .
such measurements can in principle be obtained in - situ by means of balloon borne turbulence sensors ( @xcite ) .
but a continuous monitoring of the tp using balloons would be very costly .
therefore , site testing programs use remote sensing techniques to obtain measures of the total ts .
the instruments used during such campaigns have to be rugged and should only require small size telescopes .
these requirements limited early site testing programes to measurements of the total seeing only and were first based on visual observations ( @xcite ) , photoelectric ( @xcite ) or photographic techniques ( @xcite ) .
a radical improvement of measuring the total seeing was obtained by the development of the differential image motion monitor ( dimm , @xcite ) .
remote sensing of the actual tp based on measurements of stellar scintillations has been done since some time by means of the scintillation detection and ranging ( scidar ) technique ( @xcite ) .
this technique , however , requires telescopes with apertures of at least 1  m for being effective .
therefore , in recent years several other techniques , like slope detection and ranging ( slodar , @xcite ) or multi aperture scintillation sensor ( mass , @xcite ) have been developed to overcome the requirement for such large apertures and still obtain tps with reasonable altitude coverage and resolution .
the tmt site testing employed combined mass - dimm instruments ( @xcite ) to assess the low resolution tp through the entire atmosphere . in order to measure the tp within the lowest few hundred meters above the sites , acoustic turbulence profilers sodars ( sound detection and ranging , @xcite ) were also employed .    here
we report on the results obtained with these turbulence profilers during the tmt site testing . in section  [ sitesdata ] the candidate sites ,
the deployed instrumentation and data are described .
sections  [ massbasedprofiles ] and [ sodarbasedprofiles ] show and discuss the observed `` typical '' tps . sections  [ seasonalprofiles ] , [ nightprofiles ] and [ meteorelations ] discuss seasonal and nighttime changes and the dependence on the ground wind speed of the ts at the different altitudes .
after a preselection based on satellite data and general topographic conditions ( @xcite ) , tmt deployed site testing stations on five candidate sites : 13n on mauna kea ( hawaii ) , san pedro mrtir or spm ( mexico ) , cerro tolar , cerro armazones and cerro tolonchar ( all in northern chile ) . the general properties of these mountains are summarized in table  [ sites ] .
more detailed information about these sites is given in @xcite and @xcite .
these stations have been operated for at least two years on each mountain , with a considerable time span being simultaneous between all stations ( see table  [ sites ] ) .
all tmt site testing stations were equipped with a custom built 35  cm cassegrain telescope , manufactured by teleskoptechnik halfmann .
these telescopes were mounted on 6.5  m tall towers , resulting in an elevation of the telescopes of approximately 7  m above ground .
all components of the tmt site testing equipment are summarized in @xcite and @xcite . on each telescope a combined mass - dimm device ( @xcite ) was deployed .
the dimm channel of these instruments measures the differential image motion of the star , which then results in the total seeing from the telescope level to the top of the atmosphere ( @xcite ) .
the used dimm system and its precision is described in detail by @xcite , who reported the seeing precision  the comparability between the tmt dimms  to be better than 002 .
the mass ( multi aperture scintillation sensor ) reconstructs , by measuring the spatial structure of the flying shadows , a tp at six altitudes @xmath2 km above the telescope ( @xcite ) .
the mass tps consist of the integrals of @xmath0 over each altitude bin , weighted with functions which peak at the the layer height @xmath3 and go to zero at the neighboring mass layer altitudes ( thus causing some overlap between neighboring layers ) .
therefore , the mass tps are given in @xmath4 , which is equivalent to the seeing of the layer to the power of @xmath5 .
from such a profile the seeing which would be seen by an observer 500  m above the telescope to the top of the atmosphere is computed .
also the isoplanatic angle @xmath6 can directly be obtained from these profiles .
the precision of the strength of the individual mass layers was found by @xcite to be better than @xmath7  m@xmath8 and of the mass seeing to be better than 005 .
all turbulence values in this paper are given at a wavelength of @xmath9 nm and corrected for zenith distance .
from the difference between mass and dimm seeing data , the seeing originating between the site testing telescope and the first mass layer can be computed by @xmath10 we refer to this lowest layer seeing as ground layer seeing ( gl ) . in order to obtain the gl , only mass data were used for which dimm data exist within 60  seconds .
this is usually the case as the dimm channel is used to acquire the target star ( see below ) .
when the gl is so weak that its strength is comparable to the measurement accuracy of mass and dimm , the mass seeing can become larger than the dimm seeing .
we treat these cases by applying the exponent ( @xmath11 ) in eq .
[ glformula ] on the absolute of the difference and reverting the sign of the result in case @xmath12 , same as in @xcite .
even though negative gl values are unphysical , not taking them into account would bias the resulting statistics ; in particular the low seeing percentiles .    the atmospheric turbulence time constant @xmath13 requires the additional knowledge of the vertical wind speed profile @xmath14 at the time of observation , which is usually unknown
however , @xcite found that mass data can in principle be used to infer a @xmath15 value based on the differential exposure time scintillation index .
but this @xmath15 has to be properly calibrated to obtain the real @xmath16 and is generally unreliable for @xmath17ms in the case of the tmt mass systems ( @xcite ) .
a thorough investigation of the @xmath16 measurements obtained by the tmt site testing is given by @xcite .
here , we make use of the uncalibrated @xmath15 data but only for data selection purposes for specific analyses within section  [ massbasedprofiles ] .
mass and dimm measurements are triggered simultaneously and the dimm integrates the stellar light for 36  sec , whereas the mass samples 60  sec . including processing time , this results in a mass and dimm measurement every 7090  sec . as the dimm ccd is used to acquire the target star , typically for more than 95% of the collected mass data a dimm measurement was taken within 60  sec . however , sometimes the dimm or mass computations failed to provide valid results .
this is in part due to data taken under cloudy conditions .
as the robotic system operated the telescope whenever the environmental conditions did allow a safe operation , a high data collection efficiency was assured ( @xcite ) .
nightly operations commenced approximately one hour after sunset and ceased approximately one hour before sunrise .
the tmt site testing project also operated acoustic turbulence profilers sodars ( sound detection and ranging , i.e. , @xcite ) at the candidate sites .
these systems emit sound pulses into the atmosphere and obtain from the backscatter signal , the turbulence and wind speed profiles above the site .
two sodar systems , sfas and xfas ( both manufactured by scintec ) , were operated .
the sfas profiler , measures from 10  m to 200  m above the sites with an altitude resolution of 5  m. xfas profiler range from 40  m to 800  m with a resolution of 20  m. the tmt site testing project operated three sfas and xfas pairs which were rotated among most sites ( no sodar was operated on tolar ) .
these systems obtained turbulence and wind speed profiles every 30  min .
the details on the operations , calibrations and accuracy are given in @xcite .
table  [ sites ] indicates the time spans covered by the data shown in this work .
figure  [ datapts ] shows the amount of simultaneous mass - dimm data collected at each site during each month of a year .
even though data collection was ongoing on several sites until mid 2008 , only the data obtained during the indicated periods were considered in the tmt site testing process and are presented in this work .
the site testing stations at cerro tolar , mauna kea 13n , spm and cerro tolonchar were dismantled in august 2007 , june 2008 , august 2008 and november 2008 , respectively . at some of these sites , mass - dimm data collection
had already been terminated before these dates .
the collected mass observations allow the computation of the characteristic turbulence profile above each of the candidate sites .
there is , however , no unique definition of a _
characteristic _ vertical turbulence profile that is ideally suited for all applications . in the following
we describe two different definitions of representative turbulence profiles . in both methods
we include the gl layer strength in these profiles and therefore present only data when mass and dimm data were available simultaneously .
as described before , this does not significantly bias or limit the data sample .
the first way is to compute the median and mean ts for each mass layer over the entire data set .
the statistics  25%ile , median and 75%ile  of the individual layer tss are shown in tab .
[ profilestable ] . for completeness
that table also shows the mean and rms values for each layer .
these data represent the statistical properties of each individual layer , but these profiles can not be considered as realistic or representative . in particular because the layers will not be entirely independent from each other .
for example , for the lowest mass layers , which are spaced close to each other , one might expect that shearing forces could cause interaction between these layers .    the other way to define
a tp is to compute the ts of the layers under particular seeing , isoplanatic angle and coherence time conditions .
for example , we looked for the typical tp for the @xmath18%ile of the dimm seeing , which is inverse proportional to the fried parameter @xmath19 . to this end
, we first identified @xmath20% of all dimm data with values closest to @xmath19 .
then , we compute the statistics of the mass data which were taken simultaneous with the selected dimm measurements . in this study
we set @xmath21  % , e.g. , @xmath22  % of all data around the 25 , 50 and 75 dimm seeing percentiles were used to obtain the turbulence profile .
the so defined profiles for the median @xmath23 and @xmath16 conditions are shown in the panels b , c and d of figure  [ massprofiles ] . for reference
, tables  [ profiledimmseltable ] , [ profiletheta0seltable ] and [ profiletau0seltable ] in the appendix , also show the mean and medians of the profiles computed around the 25 , 50 and 75 percentiles of respective integrated turbulence parameter .
it must be noted that even though these profiles represent the typical vertical distribution of turbulence during the times when a certain value of an integrated turbulence parameter is observed , computing that parameter from these profiles does not result in the value of the respective percentile of that parameter computed from the overall data sample .
for example in the case of tolar , computing the total seeing from the median profile from mass data selected around the median dimm seeing results in 054 .
but the overall dimm seeing statistics results in a median seeing of 063 ( @xcite ) .
this is because in the first case the statistics is done before the exponent of @xmath11 is applied to the layer ts , whereas the median seeing is obtained after that step .
these turbulence profiles now allow to understand the differences in the integrated turbulence parameters between the sites ( @xcite ) .
the three sites in chile show a very similar vertical distribution of the ts .
compared to the northern hemisphere sites , 13n on mauna kea and spm , the chilean sites show a weaker ground layer ts by approximately 50% . however , the ts above the chilean sites is stronger than the northern sites at altitudes above 4  km . at these altitudes 13n
shows the weakest ts , even less than spm whose 16  km layer is approximately as strong as that of the chilean sites .
this explains why spm shows a mass seeing similarly low as 13n but an isoplanatic angle of the same level as the chilean sites . because strong weight is put on the highest layers in the isoplanatic angle calculation ( @xmath24 )
, this makes the 16  km mass layer dominant for this specific parameter and results in 13n showing the largest isoplanatic angle . for the seeing calculation all altitudes
are weighted equally , explaining the low mass seeing of spm .
the different altitude weighting is also visible in the selected turbulence profiles .
the spread between the profiles selected from the various @xmath25 percentiles is larger than the spread between the @xmath26 selected profiles .
as example , figure  [ t6profiles ] shows the profiles from 13n ( from tables  [ profiledimmseltable ] and [ profiletheta0seltable ] ) .
it can be seen that the difference between the 25%ile and 75% gl ts median values selected from @xmath25 is approximately a factor of 30 larger than in the case of the @xmath26 selected profiles .
again , the strong weight put on the heighest layers ( @xmath24 ) in the @xmath26 calculation , causes that mainly variations of the highest layers will drive the spread of the distribution of the @xmath26 selected tp .
the seeing selected profiles reflect the variability of all layers equally , thus explaining the much wider statistical distribution of the seeing selected tps .
as the chilean sites show a weaker gl than the northern sites , but the northern sites show a weaker mass seeing than the chilean sites , a crossover between the tps above the northern and southern hemisphere sites has to occur within the lowest 500  m. we show now that the sodar based profiles allow to detect this crossover . in order to compare the sodar results directly to the mass  dimm measurements , we selected only mass - dimm observations obtained during the 30  minute sodar data acquisition times . for a thorough discussion on the data acquistion , analysis and calibration of the sodar data we refer to @xcite . figure  [ sodarprofiles ] shows the seeing which an observer would see at various altitudes above the ground up to 200  m. to obtain this profile , the xfas profile was integrated from 200  m to 500  m and added to the median mass seeing value during the time of observation .
this mass seeing is indicated by the markers to the right of the sodar profiles .
this 200  m xfas+mass seeing was then added to the sfas layers , resulting in the profiles shown in figure  [ sodarprofiles ] .
the profiles were first added and the median was calculated afterwards .
the observed dimm seeing taken simultaneously with the sodar observations is indicated by the markers to the left of the profiles , at 7.5  m above ground .
even though the sodar profiles start at 10  m above ground , the 10  m datapoint is an integral from 7.5  m to 12.5  m. and , as the sodars were always mounted approximately 1  m to 3  m below the bases of the 6.5  m tall telescope towers , it can be assumed that the sodars do not miss turbulence compared to the dimm .
the sodars were typically operated during much less time than the mass - dimm on the sites .
these profiles can therefore not be taken as representative on absolute scales for the sites .
nevertheless , these profiles allow to draw some interesting conclusions .
as can be seen from figure  [ sodarprofiles ] , the sodar profiles provide a good match between the mass and dimm observations .
this demonstrates the quality of the data analysis and calibration ( @xcite ) .
the agreement between sfas , xfas , mass and dimm is better than 10% , showing that these three methods , each based on a different physical concept , measures the ts with good accuracy .
the only site where the agreement is worse is spm .
this can be attributed to the increased acoustic noise due to the trees surrounding the site , which affect the sodar measurements ( see @xcite for details ) .
it can also be seen in figure  [ sodarprofiles ] that the profiles obtained at the chilean sites are similar to each other .
their tps show that the bulk of the gl is located below 40  m and 60  m at tolonchar and armazones , respectively . in comparison ,
the northern sites start from a weaker ts at higher altitudes , but below 100  m they increase much stronger than the chilean sites . between approximately 60  m and 90  m , their profiles are crossing the ts of the chilean sites .    from the presented turbulence profiles we conclude ,
that the topography and surface properties of the sites dominate the optical ts over the lowest hundred meters above ground , which are the main drivers for the overall seeing .
the site at spm is surrounded by trees and we believe that the therefore increased surface roughness , heat capacity , and release of humidity , cause the enhanced ts ; this raises the effective height of the ground . as 13n
is located approximately 150  m below the summit ridge of mauna kea , it falls below the wake region when winds are coming from eastern directions , as reported in @xcite . in comparison ,
the chilean sites are located on the actual summit of the candidate mountains .
the terrain shows low surface roughness and is basically free of vegetation .
these sites are are very dry and all of similar soil properties .
these arguments hold specifically for armazones and tolar .
cerro tolonchar is different in the sense that it represents a summit plateau with an extension of approximately 500  m in north  south and 250  m in east  west direction and that some sparse vegetation is present .
but as the location of the tolonchar site testing equipment is on the northern most edge of this plateau and winds are predominantely coming from the north - west ( @xcite ) , the air above that location resembles closely the free air flow .
in recent years , indications for seasonal variability of ts above other observatories were found ( @xcite , @xcite ) . clearly revealing the presence of seasonal variability in the ts at various altitudes will help to shed light on the mechanisms driving optical turbulence and will be essential for validating numerical models which link meteorological modeling to optical ts . to look for seasonal variations of the vertical turbulence distribution throughout the atmosphere in the data we present here
, we constructed what we call a standard year .
for each month of the year , the median ts of each layer was computed from all measurements taken during the corresponding months .
the resulting standard years for the turbulence distribution above the tmt candidate sites are shown in figure  [ txcnyear ] . due to technical problems , bad weather and the daylight seasons
, not all months contain the same amount of data . to judge the significance of features seen in the median ts data
, we again refer to figure  [ datapts ] which shows how many mass - dimm profiles were used for the computation of ts in each month of the standard year .
several features in figure  [ txcnyear ] are coherent at several sites and in various layers to a level of a few @xmath27  m@xmath8 .
this indicates that the resolution of the mass is likely better than a few @xmath27  m@xmath8 . as can be seen , the seasonal behaviour of the mass layers below approximately 4  km is anticorrelated with the mass layers at 8  km .
the 4  km and 16  km layer exhibit less pronounced seasonal variations at some sites , if any .
the 8  km mass layer ( @xmath28 10  km to 12  km a.s.l . ) at the sites falls closest to the altitude domain of the tropopause .
the lower mass layers thus trace the troposphere , whereas the 16  km mass layer ( @xmath28 18  km to 20  km a.s.l . )
rather traces the lower stratosphere . that the parts of the atmosphere below and above this altitude behave differently
is therefore not unexpected .
we will discuss these two altitude regimes separately .
first we discuss the chilean sites . above tolar
all layers up to 4  km show the lowest ts between may and august , raising to their highest values during the summer months . above armazones
this behaviour is less pronounced and only clearly detectable in the 0.5  km and 1  km layer , maybe some weak variation is visible in the 2  km layer as well .
tolonchar shows seasonal trends also only in its 0.5  km and 1  km layer .
the ts of the 2  km and 4  km mass layers at tolonchar show very low ts during april and may . as at tolonchar
only few data are available during the months of may and june , the feature in may represents a single year only and it is questionable whether it can be taken as typical for this site .
still , there appears to be a maximum altitude of approximately 6  km asl below which seasonal variations  weak ts in winter and strong ts during summer  can be detected in the atmosphere above northern chile .    at these low altitudes , local terrain effects will play an imporant role .
like the proximity of tolar to the pacific ocean which acts as heat reservoir ; and tolonchar rising steeply out of its surrounding arid terrain , ranging already into the westerly winds which dominate at higher elevations ( @xcite ) .
but also the large scale atmospheric conditions might be the driver of the ts .
the `` altiplano winter '' phenomenon is a result of the shift of the inner tropical convergence zone southwards over the amazonia region . over south america
this results in increased low altitude winds and convective , monsoon - like , activity ( @xcite , @xcite ) . during summer
the chilean sites are located within the region where the parts of these northerly winds which cross the andes and the anti - cyclonic southerlies meet and turn westwards off the pacific coast .
even though wind speeds in the region of the sites are low , one might expect that , as these flows merge , will result in increased large scale wind shear ; in addition to the convectivity of the notherlies .
the situation at the sites of spm and 13n is somewhat different to what is observed in chile .
the layers up to 2  km , in particular the 0.5  km one , are basically turbulence free until may .
the 0.5  km layer then becomes strong ( up to @xmath7 m@xmath8 ) from june onwards .
this behaviour does not show a smooth seasonal variation , like temperature or humidity measured at these sites .
only the 1  km and 2  km layers above 13n show a rather smooth seasonal variation on a very low level ( few @xmath27  m@xmath29 ) .
these variations are opposite to what is observed in chile : strong ts in winter and weak ts during summer . at spm
the seasonal variation is not reaching into the 4  km mass layer , whereas 13n might show some faint seasonal trend of this layer .
interestingly , the behaviour of the 0.5  km layer does not correspond to the gl .
for example , during jja the 0.5  km layer becomes strongest at spm whereas its gl is weakest during these months .
both sites are not far to the pacific ocean and we expect it to play an important role in the ts . but their local terrain is quite complex ( vegetation at spm and 13n being not on the actual summit ) , not allowing an easy qualitative assessment of the ts behaviour . for the large scale flows it can be said that both sites are located in the flow of the hawaiian high which extends further south during the summer months . during the winter it shifts towards north and the north - westerly flow in the area of hawaii
interacts with the south - easterly wind which is driven by the aleutian low ( @xcite ) .
this might explain the seasonal variation in the 1  km and 2  km layers at 13n .      of the higher altitudes ,
the 8  km layer shows the stronger seasonal variation , whereas the 16  km layer shows essentialy no seasonal trend .
the very strong seasonal variation of the 8  km layer , of up to @xmath30  m@xmath8 , is the main driver for the seasonal variability of the isoplanatic angle at some of the sites ( @xcite ) .    with altitudes of the tmt candidate sites ranging from 2.3 to 4.5  km asl
, the 8  km mass layer traces the domain of the jetstream , whose typical altiude is around 200  mb ( corresponding to @xmath3111  km asl ) .
an influence of the jetstream activity , or more general of the wind speed , on the ts has been discussed since some time in the literature ( @xcite , @xcite , @xcite ) . to see whether the wind speed is indeed one of the main drivers for the ts at these altitudes
, we make use of the vertical wind speed profiles provided by the the ncep / ncar reanalysis project ( @xcite ) .
the used locations of the ncep models are given in table  [ ncepcoord ] .
the ncep wind speed profiles are provided every six hours and at 16 altitudes , ranging from 761  m to 25  km asl . of these , only 13 or 14 levels are above the candidate sites , depending on the sites elevation . to correlate these data with the mass observations , the wind speed profiles have to be convolved with the mass weighting function of the corresponding layer .
this was done after interpolating the wind speed profile linearly onto a grid with 100  m vertical resolution .
again , we computed a standard year for these wind speed profiles , but only from those data for which mass - dimm observations were taken within @xmath323  hours . figure  [ ncepcn2 ] shows the resulting scatter plots for the standard year data of the @xmath33  km layer above the sites
. in particular armazones shows an excellent linear correlation .
tolar and tolonchar do show a similar increase of the ts with increasing wind speed , however , with more scatter than armazones . in the case of tolonchar
this enhanced scatter is in part due to low statistics during certain months ( see above ) .
nevertheless , the slope ( increase of ts with wind speed , marked as parameter @xmath34 in fig .  [ ncepcn2 ] ) is at the same level at all three sites . on the other hand , above spm and
13n the ts remains at very low levels up to @xmath35 m/s and only the higher wind speeds show a trend towards stronger @xmath36m@xmath8 values .
in other words , for the same level of wind speed , the turbulence strength at the northern sites is only about half of what is observed above chile .
also the rate of increase is only half of what is observed at the southern sites .
we tested similar correlations for the other mass layers but these do not show significant correlations .
this confirms the findings by @xcite who showed a direct correlation between the 200  mbar wind speed and the total seeing observed at telescopes on mauna kea and la silla .
@xcite deduced a threshold of 20  m / s at the 200  mbar level for good total seeing conditions .
this is remarkably similar to our observations taking into account that @xcite had to employ correction factors for the various seeing values and could only use the total seeing . @xcite later found from balloon data a correlation of @xmath16 with the 200  mbar wind speed above pachn .
this led to a number of studies of high altitude winds above other existing and potential observatory sites ( e.g. , @xcite , @xcite ) .
it is in principle understandable that high wind speeds are more likely to develop stronger ts , but our observations of the weaker increase of ts above the sites of 13n and spm are at first puzzling and are in contrast to the common threshold found in the @xcite data .
this could indicate that the mass misses to detect the jetstream above spm and 13n .
but as the mass calculates the total seeing from the integral of the reconstructed @xmath4 values and also from the total scintillation , it is possible to test the overall integration scheme of the mass system .
we found that both methods deliver values of the total seeing which agree to within 10% : the armazones median mass seeing from the scintillation calulation resulted in 040 , as compared to 043 from the @xmath0 reconstruction ( for 13n mauna kea the numbers read 030 and 034 , respectively ) .
we believe that this , in combination with the previously shown agreement with other instruments , demonstrates that no turbulence is acutally missed by the mass . on the other hand , the interpolated wind speed will certainly not be very accurate as the original ncep altitude grid is not dense .
we experimented with different interpolation schemes , as well as using the raw wind speed from the closest level to 8  km above site altitude and all resulted in very similar figures .    the observed differences between the correlations would then show that the 200  mb wind speed alone is not a good indicator on a global scale for the median ts at these altitudes and that the jetstream above the northern sites during the times of our observations was somewhat less turbulent than above the chilean mountains .
this might be caused by the lack of a large continental land mass and a high mountain range , as in the case of the andes mountains . in northern
chile , the airflow over the south american continent , which rises steeply out of the sea , can be expected to show enhanced formation of atmospheric waves .
such waves can also be forced by differential heating patterns like the proximity of the atacama desert to the pacific ocean .
atmospheric waves are commonly thought to result in an increased ts .
unfortunately , such waves and their turbulence production can not yet be modeled properly , even though comparisons between models and observations show promising results ( e.g. , @xcite ) .
our results might indicate that it is the actual interaction between wave and jet stream which is driving the ts .
to investigate the typical behaviour of the ts during the night we computed , similar to the standard year , a standard night .
it consists of the median ts of all data collected during each hour after sunset .
figure  [ txcnnight ] shows the results for the sites . to avoid artefacts introduced by the seasonal variations of the times of sunrise and sunset , the graphs in figure  [ txcnnight ] represent the merged behaviour for the hourly medians computed with respect to sunset and sunrise .
we found that the curves computed with respect to sunrise and sunset match best when midnight is set to 5  hours and the standard night for which data are available , is therefore 11  hours long .
the jumps of the curves at 5  hours in figure  [ txcnnight ] are the offset between the curves computed with respect to sunset and sunrise .
again , it can be noted that coherent structures are apparent in these graphs on a scale of a few @xmath27  m@xmath8 , consistent with what was found in section  [ seasonalprofiles ] .
also , the median night time evolution of the ts indicates that the atmosphere above approximately 4  km above the sites is basically decoupled from air below .
the high altitude layers remain close to constant throughout the night , whereas the lower layers show a decreasing ts with time .
only 13n and tolar exhibit almost no change during the course of the night in the lower layers . at spm
the 0.5  km layers at this site shows a decrease of its ts until one two hours before sunrise .
the ts in the layers up to 2  km above armazones decrease during the first half of the night and then remain at almost constant values .
finally , the air above tolonchar shows a continuous drop of ts throught the night , which maybe only levels off in the 0.5  km layer .
the observed behaviour appears similar to a cooling process and the differences between the sites on how the ts evolves below the tropopause , are related to the evolution of wind shear and temperature gradients .
we do not have measurements in order to compute the richardson number or any other stability parameter in this altitude range .
nevertheless , we make the assumption that the temperature and wind speed measured at the sites are related to their respective gradients throughout the air above the site .
this assumption is only strictly true below the height of the boundary layer which , based on the sodar profiles presented earlier , we expect to be in the order of 100  m. figure  [ mednightwst ] shows the standard night for the wind speed and temperature , recorded 2  m above ground simultaneous to the tp measurements .
it can be seen that the temperature descreases , similarly at all sites .
in contrast , the nightly evolution of the wind speed shows remarkable differences between the sites .
in particular armazones shows an increase during the second half of the night , whereas tolonchar shows a continuous decrease of the measured wind speed .
similar to tolonchar , but at a lower level , is the behaviour at 13n and tolar .
note that our wind speed measurements are unreliable at spm due to the site being surrounded by trees .
therefore , we will not try to interpret the behaviour at spm .    from these observations
we suspect that cooling of the terrain during the night is important for the decrease of the ts .
this is in particular the case for continental sites such as tolonchar and armazones .
we attribute the constant ts at armazones during the second half of the night to increased wind shear . due to its much higher heat capacity and its better efficiency to carry thermal energy below its surface ( @xcite ) ,
the sea surface temperature of the pacific ocean remains almost constant during night time .
this certainly affects the turbulent fluxes within the air above and close to the ocean .
we therefore suspect the proximity of the sites of tolar and 13n to the pacific ocean as reason for the almost constant ts observed at these sites .
spm being located in the center of the peninsula baja california appears to be rather an inland case .
the wind speed plays a crucial role as turbulence driver by introducing shear in the flow . here
the turbulence measurements of each mass layer are correlated to the wind speed which was measured 2  m above ground by the automated weather stations .
these weather stations were equipped with cup anemometers of type an3 by monitor sensors . in principle
we could have also used the wind measurements taken by the csat-3 sonic anemometers which were installed close to the telescopes at 7  m above ground ( @xcite ) .
the advantage of sonic anemometer measurements is that mechanical anemometers are prone to problems induced by dust or ice in the mechanisms .
however , we found the two systems to agree very well during most of the time .
a particular case is tolonchar , where a strong gradient exists between the 2  m and 7  m level , due to the steep terrain drop towards the prevailing wind direction .
still , we show here only the aws anemometer data , as the sonic anemometers were installed only during the second half of the site testing program and thus do not cover the full time span of mass - dimm measurements .
the panels of figure  [ cn2vsws ] show the dependence of the tss of the individual mass layers on the wind speed .
the median @xmath4 value was computed in 1  m / s wind speed bins for wind speeds up to the operational limit of the site testing telescope .
this limit was typically 15  m / s , only the tolonchar setup was limited to 12  m / s due to the strong gradient mentioned before .
we expect that the measured wind speed at 2  m above ground is indicative of the wind speed within the lowest mass layers . for the higher layers this assumption
can not be expected to hold , as the wind adjusts here to the geostrophic conditions which differ from the situation close to the ground .
nevertheless , some correlation might still be present . to visualize this behaviour figure  [ awsncepcorr ]
shows the correlation coefficient of each ncep model wind speed with the measured wind speed 2  m above the sites .
again , measured wind speeds were simultaneous to mass - dimm observations and the median was calculated for data taken within @xmath323 hours around the ncep model times .
the strongest correlations , which is always positive , are found at altitudes slightly above the sites .
the correlations break down at altitudes below the sites .
as was seen in the overall profiles , as well as in @xcite , the ground layer at all sites shows the strongest ts by up to a factor of a few .
this is not unexpected , as the mass - dimm defined ground layer covers ( most of ) the atmospheric boundary layer , i.e. , the layer in which turbulence production is dominated by the interaction between air and ground .
shearing forces are strong in this layer as the air speed at surface has to reach zero . also the ground might heat or cool the air above it by means of the sensible heat flux , introducing temperature gradients .
as was already shown in @xcite , the gl shows the lowest @xmath37 in the range of wind speeds between approximately 2  m / s and 8  m / s for the chilean sites and approximately 2  m / s and 3  m / s for 13n mauna kea .
spm shows a continuous increase of @xmath37 of the gl with increasing wind speed .
as that site is surrounded by trees the wind speed measurements , especially at low values , are doubtful
. nevertheless , the observed increase of the gl with wind speed appears reasonable also at that site .
the behaviour of the 0.5 , 1 and 2  km layers above the sites seems to be dominated by the orography of the surrounding terrain .
these layers above tolar remain at the same level of ts up to 10  m / s . at armazones the 0.5  km layer rises strongly , whereas its 1 and 2  km layers rise less steep .
these trends are similar to what happens above 13n . above tolonchar
the 0.5  km and 1  km layers remain at very low ts , but the 2  km layer rises up to 5  m / s and then shows a quite strong increase by @xmath38  m@xmath8 towards 7  m / s .
tolonchar is about 2000 m above the level of the salar de atacama ( a flat dry salt lake ) , which is in the prevailing wind direction and within 20  km from the mountain base .
these 20  km are smoothly sloped terrain ranging up 500m-800 m below the tolonchar summit .
therefore , one could assume that the wind reaching the summit effectively resembles an altitude of 1000 - 2000 m above the surface and the 0.5 and 1  km mass layers trace the free air flow between 2 and 3  km above the effective terrain .
tolar and armazones reach out similarly but not as high and steep as tolonchar . in the case of armazones a mountain range which is slightly higher ( 3400  m asl ) is located in the prevailing wind direction .
we speculate that these orographic features cause the onset of waves , driving the ts at these altitudes .
this picture seems reasonable also for the 13n site , as it is below the actual summit . as
our spm wind speed data are not reliable due to the presence of trees , we rather do not want to fit it into this picture .
nevertheless , the increase in ts of the 0.5  km layer appears reasonable .
the 8 and 16  km layers show as well an increase in their ts with wind speed at the sites , in particular tolonchar . from figure  [ awsncepcorr ]
it turns out that tolonchar shows the highest correlation of all sites between the measured and modeled wind speeds .
this indicates that its altitude brings it into the regime of the geostrophic westerly flows which are probably better modelled by ncep and are less dependent on the local orography and therefore the correlation is higher .
the mass and sodar based turbulence profiles collected at the tmt candidate sites have been presented .
the average profiles explain the origin of the observed behaviour of the integrated turbulence parameters . at all sites
the turbulence strength of the lowest 500  m  the ground layer  is found to be the dominating contributor to the overall seeing .
this layer is stronger at 13n and spm as compared to the chilean sites . on the other hand , 13n and
spm show lower turbulence strengths at higher altitudes , resulting in a weaker mass seeing than at the chilean sites .
the sodar measurements , even though these cover much less time than the mass observations , support this finding .
they indicate that the altitude dependence of the overall seeing at 13n and spm is steeper than at the chilean sites between approximately 60  m and 100  m above ground .
this results in a stronger total seeing below this altitude at the northern sites as compared to the chilean sites .
the agreement found between the measurements of mass , dimm and sodar of 10% demonstrates that these techniques do all have a good accuracy to measure @xmath0 .
however , we expect the sodar measurements to be less accurate ( @xcite ) than the mass - dimm measurements .
the isoplanatic angle differences between the sites depends on the distribution of the high altitude turbulence . in this respect ,
spm shows similar high altitude turbulence strengths as the chilean sites and therefore a similar isoplanatic angle , despite its mass seeing being comparable to 13n .
seasonal and nocturnal variations of the median turbulence strength recorded by the mass have also been shown .
these data indicate that the resolution of the mass is in the order of a few @xmath27  m@xmath8 .
clear seasonal variations are observed up to 2  km and at 8  km above the sites .
the 4  km and 16  km layers show only little or no seasonal trends .
the 8  km mass layer turbulence strength was found to correlate well with the wind speeds at these altitudes above the chilean sites . as seasonal turbulence strengths of the 8  km and the lower mass layers
are basically anit correlated , it can be expected that ground layer adaptive optics will be most efficient during the summer months at all sites .
nocturnal variations are typically confined to regions below the 4  km mass layer .
these trends seem to indicate that orographic features and surface energy properties on scales of tens of km are important drivers for the turbulence strength at altitudes below the tropopause .
in particular the close proximity of the pacific ocean to some of the sites is dominating the behaviour of the turbulence strengths at altitudes above approximately 200  m.    astronomical turbulence profilers provide a powerful tool to sense clear sky atmopsheric turbulence .
the collected data are thus not only providing the necessary information for the decision at which site to build tmt and for the design of future astronomical instruments .
these data have also to be analysed in view of their relation to meteorological mechanisms .
the present study could only provide a brief overview of some aspects of what can be extracted from such data .
more detailed studies using climatological and meteorological data , computational models in combination with new observations from turbulence profilers should result in a better understanding of the actual drivers of optical turbulence in earths atmosphere .
the tmt project gratefully acknowledges the support of the tmt partner institutions .
they are the association of canadian universities for research in astronomy ( acura ) , the california institute of technology and the university of california .
this work was supported as well by the gordon and betty moore foundation , the canada foundation for innovation , the ontario ministry of research and innovation , the national research council of canada , the natural sciences and engineering research council of canada , the british columbia knowledge development fund , the association of universities for research in astronomy ( aura ) and the u.s
. national science foundation .
barletti r. , ceppatelli g. , paterno l. , righini a. , speroni n. , 1977 , , _ 54 _ , 649 birkle k. , elssser h. , neckel th .
, schnur g. , schwarze b. , 1976 , , _ 46 _ , 397 butterley t. , wilson r.w . ,
sarazin m. , 2006 , , _ 369 _ , 835 carrasco e. , avila r. , carramina a. , 2005 , , _ 117 _ , 104 crescenti g.h . , 1997 , bull .
78 _ , 651 cuevas o. , et al . , 2009 , `` meteorological study of an event of occurence of optical turbulence at the tmt candidate sites '' , in optical turbulence
meteorology meets astronomy , masciadri e. & sarazin m. ( eds . ) , imperial college press , _ in press _
els s.g . , schck m. , seguel j. , tokovinin a. , kornilov v. , riddle r. , skidmore w. , travouillon t. , vogiatzis k. , blum r. , bustos e. , gregory b. , vasquez j. , walker d. , gillett p. , 2008
, , _ 47 _ , 14 , 2610 garcia - lorenzo b. , fuensalida j.j . , muoz - tuon c. , mendizabal e. , 2005 , , _ 356 _ , 849 irwin j.b . , 1966 , ,
_ 71 _ , 28 jumper g.y .
, murphy e.a . ,
ruggiero f.h .
, roadcap j.r . , ratkokwski a.j . , vernin j. , trinquet h. , 2007 , environ .
fluid mech . , _ 7 _ , 351 kalnay e. , kanamitsu m. , kistler r. , et al . , 1996 , bull .
soc . , _ 77 _ , 437 kornilov v. , tokovinin a. , shatsky n. , voziakova o. , potanin s. , & safonov b. , 2007 , , _ 382 _ , 1268 masciadri e. & garfias t. , 2001 , , _ 366 _ , 708 masciadri e. & egner s. , 2006 , , _ 118 _ , 1604 riddle r. , schck m. , skidmore w. , 2006 , spie , 6267 riddle r. , schck m. , skidmore w. , els s. , travouillon t. , bustos e. , seguel j. , walker d. , gillett p. , 2009
, _ in prep .
_ sarazin m. & roddier f. , 1990 , , _ 227 _ , 294 sarazin m. & tokovinin a. , 2002 , `` the statistics of isoplanatic angle and adaptive optics time constant derived from dimm data '' , in beyond conventional adaptive optics , raggazoni r. , hubin n. , esposito s. ( eds . ) , eso conference and workshop proceedings , vol .
58 , 321 satyamurty p. , nobre c.a . ,
silva dias p.l . , 1998 , `` meteorology of the tropics
south america '' , in ` meteorology of the southern hemisphere ' , koroly d.j & vincent d.g .
( eds . ) , meteorological monographs , vol .
49 , 119 schck m. , els s. , riddle r. , skidmore w. , travouillon t. , 2007 , , _
31 _ , 10 schck m. , els s. , riddle r. , skidmore w. , travouillon t. , blum r. , bustos e. , chanan g. , djorgovski g. , gillett p. , gregory b. , nelson j. , otrola a. , seguel j. , vasquez j. , walker a. , walker d. , wang l. , 2009 , pasp , _ submitted _ skidmore w. , els s. , schck m. , riddle r. , travouillon t. , seguel j. , bustos e. , walker d. , 2009 , pasp , _ in prep . _
skidmore w. , els s. , riddle r. , schck m. , travouillon t. , 2009a , pasp , _ in prep .
_ smith t.m . ,
& reynolds r.w . , 1998 , j. climate , _ 11 _ , 3320 stock j. , 1964 , `` procedures for location of astronomical observatory sites '' , in ` le choix des sites dobservatoires astronomiques ' , rosch j. ( ed . ) , iau symp .
19 stull r.b . , 1988 , `` an introduction to boundary
layer meteorology '' , kluwer academic publishers tokovinin a. , 2002 , , _ 41 _ , 957 tokovinin a. & kornilov v. , 2007 , , _ 381 _ , 1179 travouillon t. , els s. , riddle r. , schck m. , skidmore w. , blum r. , bustos e. , gillett p. , gregory b. , seguel j. , vasquez j. , walker d. , 2009a , , _ in prep .
_ travouillon t. , els s. ,
riddle r. , schck m. , skidmore w. , 2009b , pasp , _ submitted _ walker a. , djorgovski g. , blum r. , schck m. , els s. , riddle r. , skidmore w. , travouillon t. , bustos e. , deyoung d. , nelson j. , seguel j. , vasquez , vogiatzis k. , walker d. , pasp , _ in prep . _ vernin j. & roddier f. , 1973 , josa , _ 63 _ , 270 vincent d.g .
, 1998 , `` meteorology of the tropics  pacific ocean '' , in ` meteorology of the southern hemisphere ' , koroly d.j & vincent d.g .
( eds . ) , meteorological monographs , vol .
49 , 101 vernin j. , 1986 , spie , _ 628 _ , 142 wang l. , schck m. , chanan g. , et al .
2007 , , _ 46 _ , 25 , 6460 zou j. & lau k .-
, 1998 , journal of climate , _ 11 _ , 1020    crrrrrrrrrrr site name & latitude & longitude & elevation & + and telescope # & [ deg n ] & [ deg w ] & [ m a.s.l . ] & mass - dimm & sodar + cerro tolar  t1 & -21.9639 & 70.0997 & 2290 & oct . 2003  apr .
2007 & + cerro armazones
t2 & -24.5800 & 70.1833 & 3064 & nov .
2004  feb .
2008 & sep .
2005  dec .
2005 + cerro tolonchar  t3 & -23.9333 & 67.9750 & 4480 & jan .
2006  feb .
2008 & mar .
2007  feb .
2008 + spm
t4 & 31.0456 & 115.4691 & 2830 & oct .
2004  feb .
2008 & may 2006  jun .
2006 + 13n mauna kea  t6 & 19.8330 & 155.4810 & 4050 & jun .
2005  feb .
2008 & oct .
2005  dec .
2007 +    crrrrrrrrrrr site name & latitude & longitude + and telescope # & [ deg
n ] & [ deg w ] + cerro tolar
t1 & -21.25 & 71.25 + cerro armazones  t2 & -23.75 & 71.25 + cerro tolonchar  t3 & -23.75 & 68.75 + spm
t4 & 31.25 & 116.25 + 13n mauna kea  t6 & 18.75 & 156.25 +
c c c c c c c height & + @xmath39km@xmath40 $ ] & 25%ile & median & 75%ile & & mean & rms +   + 0.0 & 5.26e-14 & 1.11e-13 & 2.02e-13 & & 1.62e-13 & 2.88e-13 + 0.5 & 4.21e-18 & 3.84e-15 & 3.09e-14 & & 4.26e-14 & 1.37e-13 + 1.0 & 3.52e-16 & 3.21e-15 & 1.37e-14 & & 2.74e-14 & 1.02e-13 + 2.0 & 1.78e-16 & 3.88e-15 & 3.49e-14 & & 3.77e-14 & 8.69e-14 + 4.0 & 9.80e-15 & 2.77e-14 & 5.52e-14 & & 4.42e-14 & 5.79e-14 + 8.0 & 4.09e-15 & 2.54e-14 & 7.72e-14 & & 5.05e-14 & 6.21e-14 + 16.0 & 2.29e-14 & 3.51e-14 & 5.49e-14 & & 4.57e-14 & 3.96e-14 +   + 0.0 & 6.04e-14 & 1.19e-13 & 2.15e-13 & & 1.71e-13 & 2.77e-13 + 0.5 & 1.23e-17 & 9.78e-15 & 5.58e-14 & & 6.93e-14 & 1.85e-13 + 1.0 & 6.59e-17 & 1.99e-15 & 1.10e-14 & & 2.38e-14 & 9.17e-14 + 2.0 & 5.27e-16 & 7.16e-15 & 3.23e-14 & & 3.28e-14 & 7.44e-14 + 4.0 & 6.34e-15 & 1.93e-14 & 4.09e-14 & & 3.33e-14 & 4.66e-14 + 8.0 & 8.66e-15 & 3.03e-14 & 7.96e-14 & & 5.63e-14 & 6.87e-14 + 16.0 & 1.96e-14 & 3.04e-14 & 4.84e-14 & & 4.05e-14 & 3.68e-14 +   + 0.0 & 4.72e-14 & 1.04e-13 & 1.69e-13 & & 1.14e-13 & 1.93e-13 + 0.5 & 2.11e-21 & 1.76e-15 & 2.60e-14 & & 3.71e-14 & 1.21e-13 + 1.0 & 2.87e-17 & 1.29e-15 & 8.27e-15 & & 2.38e-14 & 9.18e-14 + 2.0 & 9.02e-16 & 1.29e-14 & 6.21e-14 & & 6.18e-14 & 1.30e-13 + 4.0 & 1.65e-14 & 3.87e-14 & 7.59e-14 & & 6.06e-14 & 7.11e-14 + 8.0 & 6.75e-15 & 3.47e-14 & 1.03e-13 & & 6.49e-14 & 7.41e-14 + 16.0 & 2.24e-14 & 3.62e-14 & 6.12e-14 & & 5.77e-14 & 6.96e-14 +   + 0.0 & 1.75e-13 & 2.76e-13 & 5.09e-13 & & 4.88e-13 & 7.05e-13 + 0.5 & 1.67e-21 & 3.10e-15 & 3.62e-14 & & 6.39e-14 & 2.66e-13 + 1.0 & 1.56e-21 & 5.39e-16 & 5.93e-15 & & 2.22e-14 & 1.09e-13 + 2.0 & 9.82e-19 & 4.30e-16 & 9.41e-15 & & 2.55e-14 & 8.61e-14 + 4.0 & 8.00e-16 & 1.37e-14 & 4.27e-14 & & 3.55e-14 & 5.86e-14 + 8.0 & 2.70e-17 & 1.06e-14 & 5.57e-14 & & 4.15e-14 & 6.57e-14 + 16.0 & 2.25e-14 & 3.60e-14 & 6.10e-14 & & 5.43e-14 & 6.00e-14 +   + 0.0 & 1.38e-13 & 2.41e-13 & 4.29e-13 & & 3.63e-13 & 5.07e-13 + 0.5 & 1.22e-20 & 5.11e-15 & 4.21e-14 & & 7.71e-14 & 3.11e-13 + 1.0 & 1.88e-16 & 2.46e-15 & 1.13e-14 & & 4.03e-14 & 1.83e-13 + 2.0 & 3.26e-17 & 8.98e-16 & 1.18e-14 & & 2.55e-14 & 8.15e-14 + 4.0 & 4.81e-15 & 1.72e-14 & 4.12e-14 & & 3.48e-14 & 5.56e-14 + 8.0 & 4.34e-16 & 6.21e-15 & 3.35e-14 & & 3.07e-14 & 5.58e-14 + 16.0 & 1.27e-14 & 2.03e-14 & 3.48e-14 & & 3.02e-14 & 3.53e-14 +    c c c c c c c c c height & & & + @xmath39km@xmath40 $ ] & 25% @xmath25 & 50% @xmath25 & 75% @xmath25 & & 25% @xmath25 & 50% @xmath25 & 75% @xmath25 +   + 0.0 & 9.41e-14 & 1.31e-13 & 1.81e-13 & & 8.72e-14 & 1.20e-13 & 1.68e-13 + 0.5 & 1.17e-15 & 4.29e-15 & 9.85e-15 & & 1.24e-14 & 2.35e-14 & 4.62e-14 + 1.0 & 2.14e-15 & 3.53e-15 & 6.03e-15 & & 6.20e-15 & 1.24e-14 & 2.80e-14 + 2.0 & 9.79e-16 & 5.26e-15 & 2.08e-14 & & 8.15e-15 & 2.14e-14 & 4.98e-14 + 4.0 & 2.04e-14 & 2.99e-14 & 4.51e-14 & & 2.54e-14 & 3.69e-14 & 6.19e-14 + 8.0 & 2.09e-14 & 3.15e-14 & 4.24e-14 & & 3.76e-14 & 5.27e-14 & 6.73e-14 + 16.0 & 3.32e-14 & 3.66e-14 & 3.96e-14 & & 3.90e-14 & 4.55e-14 & 5.33e-14 +   + 0.0 & 1.02e-13 & 1.37e-13 & 2.05e-13 & & 9.64e-14 & 1.29e-13 & 1.91e-13 + 0.5 & 2.03e-15 & 1.02e-14 & 3.71e-14 & & 1.20e-14 & 2.94e-14 & 7.86e-14 + 1.0 & 9.61e-16 & 2.17e-15 & 5.49e-15 & & 4.29e-15 & 9.47e-15 & 2.62e-14 + 2.0 & 2.90e-15 & 1.04e-14 & 1.99e-14 & & 1.01e-14 & 2.42e-14 & 4.54e-14 + 4.0 & 1.41e-14 & 2.14e-14 & 3.24e-14 & & 1.84e-14 & 2.75e-14 & 4.69e-14 + 8.0 & 2.49e-14 & 3.69e-14 & 4.52e-14 & & 3.89e-14 & 5.76e-14 & 7.48e-14 + 16.0 & 2.75e-14 & 3.09e-14 & 3.53e-14 & & 3.23e-14 & 3.85e-14 & 4.72e-14 +   + 0.0 & 1.03e-13 & 1.24e-13 & 1.41e-13 & & 9.59e-14 & 1.17e-13 & 1.33e-13 + 0.5 & 8.53e-16 & 2.46e-15 & 4.17e-15 & & 1.15e-14 & 2.23e-14 & 4.61e-14 + 1.0 & 1.02e-15 & 1.33e-15 & 2.81e-15 & & 4.19e-15 & 7.76e-15 & 2.24e-14 + 2.0 & 3.17e-15 & 1.30e-14 & 4.32e-14 & & 1.23e-14 & 2.95e-14 & 6.76e-14 + 4.0 & 2.91e-14 & 3.96e-14 & 5.98e-14 & & 3.39e-14 & 4.68e-14 & 7.55e-14 + 8.0 & 2.14e-14 & 3.86e-14 & 7.46e-14 & & 3.65e-14 & 6.01e-14 & 8.67e-14 + 16.0 & 3.07e-14 & 3.66e-14 & 4.61e-14 & & 3.38e-14 & 4.28e-14 & 6.01e-14 +   + 0.0 & 2.06e-13 & 3.09e-13 & 5.55e-13 & & 1.94e-13 & 2.90e-13 & 5.07e-13 + 0.5 & 1.07e-15 & 2.35e-15 & 6.73e-15 & & 1.56e-14 & 2.75e-14 & 6.10e-14 + 1.0 & 2.52e-16 & 6.51e-16 & 1.33e-15 & & 4.04e-15 & 9.02e-15 & 2.56e-14 + 2.0 & 1.89e-16 & 3.80e-16 & 1.67e-15 & & 4.55e-15 & 1.03e-14 & 3.71e-14 + 4.0 & 9.77e-15 & 1.62e-14 & 2.80e-14 & & 1.71e-14 & 3.00e-14 & 5.49e-14 + 8.0 & 3.91e-15 & 1.38e-14 & 2.86e-14 & & 1.90e-14 & 3.75e-14 & 6.12e-14 + 16.0 & 3.08e-14 & 3.62e-14 & 4.76e-14 & & 3.55e-14 & 4.47e-14 & 6.60e-14 +   + 0.0 & 1.83e-13 & 2.85e-13 & 4.89e-13 & & 1.67e-13 & 2.56e-13 & 4.28e-13 + 0.5 & 1.45e-15 & 4.89e-15 & 1.22e-14 & & 1.78e-14 & 3.43e-14 & 6.49e-14 + 1.0 & 1.33e-15 & 2.36e-15 & 5.03e-15 & & 5.06e-15 & 1.02e-14 & 3.24e-14 + 2.0 & 2.93e-16 & 9.77e-16 & 3.62e-15 & & 4.75e-15 & 1.33e-14 & 3.21e-14 + 4.0 & 1.25e-14 & 1.95e-14 & 2.58e-14 & & 1.95e-14 & 3.03e-14 & 4.48e-14 + 8.0 & 4.30e-15 & 7.30e-15 & 1.29e-14 & & 1.67e-14 & 2.64e-14 & 4.46e-14 + 16.0 & 1.89e-14 & 2.06e-14 & 2.34e-14 & & 2.40e-14 & 2.72e-14 & 3.48e-14 +    c
c c c c c c c c height & & & + @xmath39km@xmath40 $ ] & 25% @xmath26 & 50% @xmath26 & 75% @xmath26 & & 25% @xmath26 & 50% @xmath26 & 75% @xmath26 +   + 0.0 & 1.32e-13 & 1.20e-13 & 9.28e-14 & & 2.00e-13 & 1.82e-13 & 1.47e-13 + 0.5 & 7.43e-15 & 5.92e-15 & 9.65e-16 & & 4.45e-14 & 5.10e-14 & 4.33e-14 + 1.0 & 3.08e-15 & 3.86e-15 & 3.67e-15 & & 1.82e-14 & 2.52e-14 & 3.28e-14 + 2.0 & 2.17e-15 & 2.93e-15 & 8.35e-15 & & 2.39e-14 & 3.42e-14 & 5.18e-14 + 4.0 & 2.67e-14 & 3.42e-14 & 2.87e-14 & & 3.24e-14 & 4.39e-14 & 4.84e-14 + 8.0 & 7.49e-15 & 2.18e-14 & 6.76e-14 & & 1.53e-14 & 3.14e-14 & 7.11e-14 + 16.0 & 2.68e-14 & 3.83e-14 & 5.47e-14 & & 2.52e-14 & 3.60e-14 & 5.28e-14 +   + 0.0 & 1.32e-13 & 1.23e-13 & 1.06e-13 & & 1.90e-13 & 1.76e-13 & 1.68e-13 + 0.5 & 1.53e-14 & 1.37e-14 & 5.12e-15 & & 6.45e-14 & 8.82e-14 & 7.45e-14 + 1.0 & 2.35e-15 & 2.60e-15 & 2.56e-15 & & 1.73e-14 & 2.25e-14 & 3.08e-14 + 2.0 & 4.86e-15 & 7.31e-15 & 1.17e-14 & & 2.05e-14 & 2.95e-14 & 3.97e-14 + 4.0 & 1.83e-14 & 2.26e-14 & 2.27e-14 & & 2.29e-14 & 3.00e-14 & 3.75e-14 + 8.0 & 1.24e-14 & 2.85e-14 & 7.16e-14 & & 1.75e-14 & 3.46e-14 & 7.57e-14 + 16.0 & 2.25e-14 & 3.30e-14 & 4.86e-14 & & 2.15e-14 & 3.15e-14 & 4.69e-14 +   + 0.0 & 1.17e-13 & 1.12e-13 & 8.79e-14 & & 1.40e-13 & 1.32e-13 & 1.05e-13 + 0.5 & 7.59e-15 & 4.53e-15 & 1.17e-18 & & 4.20e-14 & 4.16e-14 & 3.11e-14 + 1.0 & 1.87e-15 & 2.42e-15 & 5.66e-16 & & 1.25e-14 & 2.15e-14 & 1.96e-14 + 2.0 & 7.00e-15 & 8.38e-15 & 2.16e-14 & & 3.18e-14 & 3.87e-14 & 6.41e-14 + 4.0 & 3.34e-14 & 4.58e-14 & 3.60e-14 & & 3.85e-14 & 5.44e-14 & 5.51e-14 + 8.0 & 1.09e-14 & 2.55e-14 & 9.83e-14 & & 1.79e-14 & 3.55e-14 & 1.02e-13 + 16.0 & 2.56e-14 & 3.98e-14 & 5.88e-14 & & 2.42e-14 & 3.76e-14 & 5.66e-14 +   + 0.0 & 2.39e-13 & 2.61e-13 & 3.22e-13 & & 3.44e-13 & 3.90e-13 & 5.36e-13 + 0.5 & 5.03e-15 & 4.08e-15 & 1.73e-15 & & 3.91e-14 & 5.42e-14 & 6.44e-14 + 1.0 & 7.48e-16 & 9.78e-16 & 5.68e-16 & & 8.99e-15 & 1.51e-14 & 3.12e-14 + 2.0 & 3.09e-16 & 3.93e-16 & 7.34e-16 & & 9.43e-15 & 1.85e-14 & 3.65e-14 + 4.0 & 1.27e-14 & 1.86e-14 & 2.07e-14 & & 1.95e-14 & 3.29e-14 & 4.58e-14 + 8.0 & 9.21e-16 & 1.01e-14 & 4.77e-14 & & 6.73e-15 & 2.00e-14 & 5.60e-14 + 16.0 & 2.51e-14 & 3.89e-14 & 6.20e-14 & & 2.41e-14 & 3.72e-14 & 6.02e-14 +   + 0.0 & 2.32e-13 & 2.44e-13 & 2.42e-13 & & 3.33e-13 & 3.69e-13 & 3.59e-13 + 0.5 & 7.29e-15 & 1.14e-14 & 4.59e-15 & & 4.10e-14 & 7.89e-14 & 1.12e-13 + 1.0 & 1.79e-15 & 2.68e-15 & 3.46e-15 & & 1.62e-14 & 2.76e-14 & 5.44e-14 + 2.0 & 5.80e-16 & 8.99e-16 & 1.10e-15 & & 1.02e-14 & 2.06e-14 & 3.38e-14 + 4.0 & 1.30e-14 & 2.12e-14 & 2.90e-14 & & 1.79e-14 & 2.95e-14 & 4.25e-14 + 8.0 & 1.48e-15 & 5.93e-15 & 2.39e-14 & & 4.89e-15 & 1.32e-14 & 3.51e-14 + 16.0 & 1.50e-14 & 2.29e-14 & 3.70e-14 & & 1.43e-14 & 2.15e-14 & 3.47e-14 +    c c c c c c c c c height & & & + @xmath39km@xmath40 $ ] & 25% @xmath16 & 50% @xmath16 & 75% @xmath16 & & 25% @xmath16 & 50% @xmath16 & 75% @xmath16 +   + 0.0 & 1.41e-13 & 1.11e-13 & 7.97e-14 & & 2.02e-13 & 1.76e-13 & 1.29e-13 + 0.5 & 9.42e-15 & 4.13e-15 & 7.54e-16 & & 5.24e-14 & 4.72e-14 & 3.42e-14 + 1.0 & 4.12e-15 & 3.02e-15 & 3.07e-15 & & 3.26e-14 & 2.32e-14 & 2.42e-14 + 2.0 & 5.87e-15 & 2.69e-15 & 4.23e-15 & & 3.86e-14 & 3.50e-14 & 4.38e-14 + 4.0 & 2.94e-14 & 2.81e-14 & 2.86e-14 & & 3.87e-14 & 4.00e-14 & 4.88e-14 + 8.0 & 8.51e-15 & 2.39e-14 & 6.83e-14 & & 1.98e-14 & 3.63e-14 & 7.49e-14 + 16.0 & 2.80e-14 & 3.64e-14 & 4.83e-14 & & 3.17e-14 & 4.01e-14 & 5.46e-14 +   + 0.0 & 1.31e-13 & 1.12e-13 & 1.01e-13 & & 1.87e-13 & 1.58e-13 & 1.58e-13 + 0.5 & 1.36e-14 & 1.25e-14 & 6.33e-15 & & 6.86e-14 & 7.97e-14 & 7.56e-14 + 1.0 & 1.99e-15 & 2.35e-15 & 2.35e-15 & & 1.76e-14 & 2.35e-14 & 2.69e-14 + 2.0 & 4.75e-15 & 7.67e-15 & 1.12e-14 & & 1.98e-14 & 2.89e-14 & 3.93e-14 + 4.0 & 1.59e-14 & 2.12e-14 & 2.49e-14 & & 2.21e-14 & 2.91e-14 & 3.64e-14 + 8.0 & 1.13e-14 & 2.84e-14 & 6.51e-14 & & 1.83e-14 & 3.73e-14 & 7.31e-14 + 16.0 & 2.36e-14 & 3.06e-14 & 3.95e-14 & & 2.67e-14 & 3.50e-14 & 4.60e-14 +   + 0.0 & 1.20e-13 & 1.01e-13 & 6.26e-14 & & 1.43e-13 & 1.23e-13 & 7.65e-14 + 0.5 & 4.47e-15 & 1.49e-15 & 7.96e-18 & & 4.05e-14 & 3.79e-14 & 3.52e-14 + 1.0 & 1.53e-15 & 1.60e-15 & 8.31e-16 & & 1.57e-14 & 1.91e-14 & 2.66e-14 + 2.0 & 6.94e-15 & 8.55e-15 & 2.46e-14 & & 3.28e-14 & 4.06e-14 & 6.38e-14 + 4.0 & 2.98e-14 & 3.63e-14 & 5.92e-14 & & 3.70e-14 & 4.54e-14 & 7.11e-14 + 8.0 & 1.57e-14 & 3.91e-14 & 1.10e-13 & & 2.82e-14 & 5.07e-14 & 1.12e-13 + 16.0 & 3.14e-14 & 3.71e-14 & 5.35e-14 & & 3.55e-14 & 4.29e-14 & 6.85e-14 +   + 0.0 & 2.29e-13 & 2.95e-13 & 3.87e-13 & & 3.24e-13 & 4.26e-13 & 6.38e-13 + 0.5 & 5.56e-15 & 3.70e-15 & 1.40e-15 & & 4.61e-14 & 5.30e-14 & 7.65e-14 + 1.0 & 8.20e-16 & 8.35e-16 & 2.40e-16 & & 1.30e-14 & 1.82e-14 & 3.75e-14 + 2.0 & 3.71e-16 & 3.99e-16 & 4.32e-16 & & 1.47e-14 & 2.52e-14 & 3.23e-14 + 4.0 & 1.15e-14 & 1.67e-14 & 1.91e-14 & & 2.46e-14 & 3.46e-14 & 4.22e-14 + 8.0 & 1.32e-15 & 8.79e-15 & 5.51e-14 & & 1.14e-14 & 2.45e-14 & 6.61e-14 + 16.0 & 2.77e-14 & 3.63e-14 & 5.53e-14 & & 3.23e-14 & 4.23e-14 & 6.58e-14 +   + 0.0 & 2.46e-13 & 2.50e-13 & 2.21e-13 & & 3.70e-13 & 3.95e-13 & 3.74e-13 + 0.5 & 1.04e-14 & 1.06e-14 & 2.51e-15 & & 5.74e-14 & 8.16e-14 & 1.55e-13 + 1.0 & 2.50e-15 & 2.88e-15 & 3.99e-15 & & 2.68e-14 & 4.07e-14 & 7.53e-14 + 2.0 & 8.59e-16 & 8.83e-16 & 2.31e-15 & & 1.97e-14 & 2.31e-14 & 4.04e-14 + 4.0 & 1.52e-14 & 2.20e-14 & 3.55e-14 & & 2.30e-14 & 3.18e-14 & 5.17e-14 + 8.0 & 2.84e-15 & 9.31e-15 & 4.18e-14 & & 9.47e-15 & 2.17e-14 & 5.48e-14 + 16.0 & 1.75e-14 & 2.35e-14 & 3.68e-14 & & 2.03e-14 & 2.77e-14 & 4.22e-14 +