Comparing CO2 Storage and Advection Conditions at Night at Different Carboeuroflux Sites

Anemometer and CO₂ concentration data from temporary campaigns performed at six CARBOEUROFLUX forest sites were used to estimate the importance of non-turbulent fluxes in nighttime conditions. While storage was observed to be significant only during periods of both low turbulence and low advection, the advective fluxes strongly influence the nocturnal CO₂ balance, with the exception of almost flat and highly homogeneous sites. On the basis of the main factors determining the onset of advective fluxes, the ‘advection velocity’, which takes net radiation and local topography into account, was introduced as a criterion to characterise the conditions of storage enrichment/depletion. Comparative analyses of the six sites showed several common features of the advective fluxes but also some substantial differences. In particular, all sites where advection occurs show the onset of a boundary layer characterised by a downslope flow, negative vertical velocities and negative vertical CO₂ concentration gradients during nighttime. As a consequence, vertical advection was observed to be positive at all sites, which corresponds to a removal of CO₂ from the ecosystem. The main differences between sites are the distance from the ridge, which influences the boundary-layer depth, and the sign of the mean horizontal CO₂ concentration gradients, which is probably determined by the source/sink distribution. As a consequence, both positive and negative horizontal advective fluxes (corresponding respectively to CO₂ removal from the ecosystem and to CO₂ supply to the ecosystem) were observed. Conclusive results on the importance of non-turbulent components in the mass balance require, however, further experimental investigations at sites with different topographies, slopes, different land covers, which would allow a more comprehensive analysis of the processes underlying the occurrence of advective fluxes. The quantification of these processes would help to better quantify nocturnal CO₂ exchange rates.     

Estimating the components of the sensible heat budget of a tall forest canopy in complex terrain

Ultrasonic wind measurements, sonic temperature and air temperature data at two heights in the advection experiment MORE II were used to establish a complete budget of sensible heat including vertical advection, horizontal advection and horizontal turbulent flux divergence. MORE II took place at the long-term Carbo-Europe IP site in Tharandt, Germany. During the growing period of 2003 three additional towers were established to measure all relevant parameters for an estimation of advective fluxes, primarily of CO₂. Additionally, in relation to other advection experiments, a calculation of the horizontal turbulent flux divergence is proposed and the relation of this flux to atmospheric stability and friction velocity is discussed. In order to obtain a complete budget, different scaling heights for horizontal advection and horizontal turbulent flux divergence are tested. It is shown that neglecting advective fluxes may lead to incorrect results. If advective fluxes are taken into account, the sensible heat budget based upon vertical turbulent flux and storage change only, is reduced by approximately 30%. Additional consideration of horizontal turbulent flux divergence would in turn add 5–10% to this sum (i.e., the sum of vertical turbulent flux plus storage change plus horizontal and vertical advection). In comparison with available energy horizontal advection is important at night whilst horizontal turbulent flux divergence is rather insignificant. Obviously, advective fluxes typically improve poor nighttime energy budget closure and might change ecosystem respiration fluxes considerably.     

Aspects of CO2 advection measurements

Observations of vegetation–atmosphere exchange of carbon dioxide (CO₂) by the eddy covariance (EC) technique are limited by difficult conditions such as nighttime and heterogeneous terrain. Thus, advective flux components are included into the net ecosystem exchange (NEE) budget. However, advection measurements are experimentally challenging and do not always help to solve the night flux problem of the EC technique. This study investigates alternative methods for the observation of horizontal advection, in particular horizontal concentration gradients, as well as different approaches to coordinate rotation and vertical advection. Continuous high-frequency measurements of the horizontal CO₂ concentration field are employed and compared to the often used discontinuous sequential sampling. Significant differences were found in the case of 30-min mean concentration values between the conventional discontinuous sampling approach and the complete observation of the time series by continuous sampling. Estimates of vertical advection rely on accurate estimates of vertical wind velocity ( $\emph{w}$ ). Therefore, different approaches to the planar fit coordinate rotation have been investigated. Sector-wise rotation was able to eliminate directional dependencies of mean $\emph{w}$ . Furthermore, the effect of the data set length used for rotation (window length) was investigated and was found to have significant impact on estimates of vertical advection, with larger window lengths yielding about 50% larger vertical advection. A sequential planar fit with controlled window length is proposed to give reproducible results. The different approaches to the measurement and calculation of horizontal and vertical advection presented are applied to data obtained during the exchange processes in mountainous region experiment at the FLUXNET site Waldstein–Weidenbrunnen (DE-Bay). Estimates of NEE including advection are compared to NEE from turbulent and storage flux alone without advection. NEE including vertical advection with sector-wise planar fit rotation and controlled window length and including horizontal advection from continuous gradient measurements, which were comprehensively bias corrected by a new approach, did compare well with the expected night flux error, with meteorological drivers of the fluxes and with soil chamber measurements. Unrealistically large and noisy values of horizontal advection from the conventional discontinuous sampling approach, which lead to unrealistic values of NEE, could be eliminated by the alternative approaches presented. We therefore suggest the further testing of those approaches at other sites in order to improve the accuracy of advection measurements and, subsequently, estimates of NEE.     

Horizontal and Vertical Co2 Advection In A Sloping Forest

A system measuring the horizontal and vertical advection was devised and installed in a sloping forest at the Vielsalm site, Belgium. The measurements showed that under stable conditions a flow regime established below the canopy: air flowed horizontally along the slope and entrained the air above the canopy vertically. This movement occurs during stable nights characterised by strongly negative net radiation. It creates negative air concentration gradients in both the vertical and horizontal directions. The advection fluxes associated with these movements are opposite and of a similar order of magnitude. This implies that the horizontal advection cannot be ignored in the carbon budget equation at night. Unfortunately, the large variability of, and considerable uncertainty about, advection fluxes does not enable one to produce estimates of the source term from these equations. Advection measurement systems should be improved in order to enable such estimates to be made. Particular attention should be paid to the estimation of the vertical velocity above the canopy and to the vertical profiles of the horizontal velocity and horizontal CO₂ gradient below the canopy.     

Some methodological questions concerning advection measurements: a case study

A dataset from two campaigns conducted at the Vielsalm experimental site in Belgium was used as a basis for discussing some methodological problems and providing intermediate results on estimating CO₂ advection. The analysis focused on the horizontal [CO₂] gradient and on the vertical velocity w, the variables most affected by uncertainty. The sampling error for half-hourly horizontal [CO₂] gradients was estimated to be 1.3 μmol mol⁻¹. Despite this important random error for half-hour estimations of [CO₂], the mean horizontal [CO₂] gradients in advective conditions were shown to be representative at the ecosystem scale and to extend only to the lowest part of a drainage sub-layer, which developed in the trunk space. By contrast, under daytime conditions, this gradient was shown to be more sensitive to local source heterogeneities. The estimation of the short-term averaged vertical velocity ( was the greater source of error when computing advection terms. The traditional correction methods used to obtain are discussed and a (co)sine correction is tested to highlight the instrumental origin of the offset in w. A comparison of measurements by sonic anemometers placed close together above the canopy showed that the uncertainty on was 0.042 m s⁻¹, which is of the same order of magnitude as the velocity itself. In addition, as the drainage sub-layer is limited to the lowest part of the canopy, the representativeness of is questionable. An alternative computation using the divergence of the horizontal wind speed in the trunk space produced a estimation that was four times lower than the single-point measurement. However, this value gives a more realistic estimate of the vertical advection term and improves the CO₂ budget closure at the site.     

Forest-Air Fluxes Of Carbon, Water And Energy Over Non-Flat Terrain

A field study of surface-air exchange of carbon, water, and energy was conducted at a mid-latitude, mixed forest on non-flat terrain to investigate how to best interpret biological signals from the eddy flux data that may be subject to advective influences. It is shown that during periods of Southwest winds (sector with mild topography), the eddy fluxes are well-behaved in terms of energy balance closure, the existence of a constant flux layer, consistency with chamber observations and the expected abiotic controls on the fluxes. Advective influences are evident during periods with wind from a steep (15%) slope to the Northeast of the tower. These influences appear more severe on CO₂ flux, particularly in stable air, than on the energy fluxes. Large positive flux of CO₂ (> 23 mol m^-2 s^-1) occurs frequently at night. The annual sum of the carbon flux is positive, but the issue about whether the forest is a source of atmospheric carbon remains inconclusive.Attempts are made to assess vertical advectionusing the data collected on a single tower. Over the Southwestsector, vertical advection makes a statistically significant but small contribution to the 30-min energy imbalance and CO₂ flux variations. Contributions by horizontal advection may be larger but cannot be verified directly by the current experimental method.     

A global three-dimensional model of the tropospheric sulfur cycle

The tropospheric part of the atmospheric sulfur cycle has been simulated in a global three-dimensional model. The model treats the emission, transport, chemistry, and removal processes for three sulfur components; DMS (dimethyl sulfide), SO₂ and SO₄ ^2– (sulfate). These processes are resolved using an Eulerian transport model, the MOGUNTIA model, with a horizontal resolution of 10° longitude by 10° latitude and with 10 layers in the vertical between the surface and 100 hPa. Advection takes place by climatological monthly mean winds. Transport processes occurring on smaller space and time scales are parameterized as eddy diffusion except for transport in deep convective clouds which is treated separately. The simulations are broadly consistent with observations of concentrations in air and precipitation in and over polluted regions in Europe and North America. Oxidation of DMS by OH radicals together with a global emission of 16 Tg DMS-S yr^–1 from the oceans result in DMS concentrations consistent with observations in the marine boundary layer. The average turn-over times were estimated to be 3, 1.2–1.8, and 3.2–6.1 days for DMS, SO₂, and SO₄ ^2– respectively.     

On the annual and monthly mean diurnal variations of diffuse solar radiation at a meteorological station in west Africa

Summary In addition to global solar radiationE ^g , the hourly diffuse componentE ^d incident on a horizontal surface has been measured from February 1993 to January 1995 at a meteorological station in tropical West Africa. The measured diffuse solar irradiance data was corrected for shadow band effects. The monthly mean diurnal variations of diffuse solar irradiance obtained for identical months in the two years have been compared and found to be generally consistent. The corresponding monthly mean hourly values ofE ^d for identical months in 1993 and 1994 agreed to within 9% while yielding correlation coefficients greater than 0.960. In addition, the monthly mean daily totals ofE ^d for identical months were found to agree mostly to within 6% and showed virtually the same annual variations in both years. The monthly mean daily total values of diffuse solar radiation for most months in the two years ranged between 7.94 MJm^–2d^–1 and 10.50 MJm^–2d^–1. The monthly mean of daily hourly maximum values ofE ^d obtained for identical months in the two years have been discussed in relation to the dominant atmospheric conditions during these months. The results been presented here have been compared with those of some investigators within and outside the Africa region.With 8 Figures     

Measured and modelled contributions to UV exposures by the albedo of surfaces in an urban environment

Summary The increases in the erythemal UV exposures to horizontal planes and to inclined planes over three surfaces that are found in an urban environment (water, concrete and sand) due to the albedo of these surfaces have been estimated. For the cloud free case, the additional daily estimated UV exposures to a horizontal plane have a maximum value of 222 (Jm^–2)_ER, where the index after the unit is there to indicate that it refers to a biologically effective exposure. In comparison, the daily erythemal UV exposures over a year to a horizontal plane ranged from 425 to 8,321 (Jm^–2)_ER. For a vertical receiving plane that is rotating about a vertical axis, the additional erythemal daily UV exposures for the sub-tropical latitude location of this research for the ranges of solar azimuth angles encountered over the days in each season ranged from 16 to 311 (Jm^–2)_ER, 29 to 566 (Jm^–2)_ER and 46 to 905 (Jm^–2)_ER for water, concrete and sand respectively. The estimated error is ±20% and the calculations are based on clear-sky conditions. The additional erythemal UV averaged over each of the seasons was higher for the receiving plane inclined at 45° below the horizontal plane. In a similar fashion, the vertical surface has the higher additional erythemal UV exposures compared to the surfaces inclined at an angle above the horizontal.     

Accuracy of the relaxed eddy-accumulation technique,evaluated using CO2 flux measurements

A system capable of measuring the fluxes of trace gases was developed. It is based on a simpler version of the eddy-accumulation technique (EA), known as the relaxed eddy-accumulation technique (REA). It accumulates air samples associated with updrafts and downdrafts at a constant flow rate in two containers for later analysis of the trace gas mean concentration. The flux integration is based on the durations of updraft and downdraft events, rather than on the vertical wind velocity (W) as is the case for EA and eddy-correlation (EC) techniques. The flux, calculated by the REA technique, is equal to the difference in the mean concentration of the trace gas of interest between the upward and downward moving eddies, multiplied by the standard deviation of the vertical wind velocity and an empirical coefficient. CO₂ fluxes measured for 162 half-hour periods over a soybean field by both EC and REA techniques showed excellent agreement (coefficient of determination,R ²=0.92). The slope (0.985) and the intercept (–0.042 mg m^–2 s^–1) were not significantly different from 1 and 0, respectively, at the 5% level; and the standard error of estimate was 0.074 mg m^–2 s^–1. It is also shown that the empirical coefficient can be calculated from either latent or sensible heat fluxes. A model describing the effect on this empirical coefficient of not sampling aroundW equal to zero is proposed.Centre for Land and Biological Resources Research Contribution No. 92-212.     

Emissions of marine biogenic sulfur to the atmosphere of northern Europe

Measurements of DMS and other reduced sulfur compounds in surface waters have been carried out from a helicopter in the seas surrounding Scandinavia. Average summer time concentrations of DMS ranged from 70 to 150 ngS L^-1. Simultaneous measurements of biological and physical parameters revealed no correlation between DMS and phytoplankton species, species assemblages, total phytoplankton biomass, chlorophyll a, temperature, and salinity. The only exception was a correlation between DMS concentration, Chrysochromulina spp. belonging to the Prymnesiophyceae, and salinity over a narrow range of salinity in the Baltic Sea.The flux of reduced sulfur to the atmosphere in July in this region is estimated to be 120–170 gS m^-2 d^-1 from the Baltic, 240–810 in the Kattegat/Skagerrak, and 120–690 in the North Sea. Annual fluxes are roughly 100 times higher than these daily fluxes. On an annual basis, biogenic sulfur emissions from the coastal seas are negligible (<1%) compared to the anthropogenic emissions in northern Europe. However, during the summer months, the biogenic sulfur emissions from the seas surrounding the Scandinavian peninsula are estimated to be as high as 20–70% of the anthropogenic emissions in Scandinavia. This makes it of interest to incorporate the biogenic emissions in calculations of long-range transport and deposition of sulfur within the region.Other volatile sulfur species, mainly methyl mercaptan, contribute about 10% of the total flux of reduced sulfur. Estimated fluxes of CS₂ to the atmosphere ranged from 1 gS m^-2 d^-1 in the Baltic Sea to 6 gS m^-2 d^-1 in the North Sea. No emissions for H₂S or COS were detected.     

Actinic flux and photolysis frequency comparison computations using the model PHOTOGT

The accurate radiative transfer model GOMETRAN, initially designed to yield radiances at TOA in the wavelength range 240–790 nm, has been extended to allow for the computation of actinic fluxes down to 175 nm and for the calculation of photolysis frequencies in the atmosphere. The capability of the extended model PHOTOGT (PHOTOGOMETRAN) is demonstrated in a number of successful comparison studies both with recent experiments (ground-based, balloonborne, airborne) and model calculations of radiances, actinic fluxes and photolysis frequencies in the stratosphere and troposphere. In an atmospheric case study, the impact of new quantum yield data for the O₃ » O₂+O(¹ d) photodissociation channel on the photolytic production of O(¹ d) atoms in the lower atmosphere has been quantified.     

Divergence of turbulent fluxes in the surface layer: case of a coastal city

This study quantifies the processes that take place in the layer between the mean building height and the measurement level of an energy balance micrometeorological tower located in the dense old core of a coastal European city. The contributions of storage, vertical advection, horizontal advection and radiative divergence for heat are evaluated with the available measurements and with a three-dimensional, high-resolution meteorological simulation that had been evaluated against observations. The study focused on a summer period characterized by sea-breeze flows that affect the city. In this specific configuration, it appears that the horizontal advection is the dominant term. During the afternoon when the sea breeze is well established, correction of the sensible heat flux with horizontal heat advection increases the measured sensible heat flux up to 100 W m⁻². For latent heat flux, the horizontal moisture advection converted to equivalent latent heat flux suggests a decrease of 50 W m⁻². The simulation reproduces well the temporal evolution and magnitude of these terms.     

Correction Of Eddy-Covariance Measurements Incorporating Both Advective Effects And Density Fluxes

Equations are presented to correct eddy-covariancemeasurements for both fluctuations in density andnon-zero mean advection, induced by convergence ordivergence of flow, and spatial source/sinkinhomogeneity, under steady-state and transientconditions. This correction collapses to theWebb–Pearman–Leuning expression ifthe mean vertical velocity is zero, and formally addsthe Webb–Pearman–Leuning expression to the correctionssuggested by Lee for conditions ofnon-zero vertical velocity and source/sink and meanscalar horizontal homogeneity. The equation requiresmeasurement of the mean vertical gradients of thescalar concentration of interest (air temperature,humidity, CO₂) as well as an accurateestimation of the mean vertical velocity, in additionto the vertical eddy covariance of the scalar. Simplemethods for the approximation of sensor tilt andcomplex terrain flow angle are presented, to allowestimation of non-zero mean vertical velocities. Theequations are applied to data from a maize crop and aforest to give examples of when the correction issignificant. In addition, a term for thethermodynamic expansion energy associated with watervapour flux is derived, which implies that the sonictemperature derived sensible heat flux will accuratelyinclude this contribution.     

Measurement of the Sensible Eddy Heat Flux Based on Spatial Averaging of Continuous Ground-Based Observations

Using the standard eddy-covariance (EC) method to quantify mass and energy exchange at a single location usually results in an underestimation of vertical eddy fluxes at the surface. In order to better understand the reasons for this underestimation, an experimental set-up is presented that is based on spatial averaging of air temperature data from a network of ground-based sensors over agricultural land. For eight days during the 34-day observational period in May and June 2007, additional contributions to the sensible heat flux of more than 50Wm⁻² were measured in the lower surface layer by applying the spatial EC method as opposed to the standard temporal EC method. Smaller but still significant additional sensible heat fluxes were detected for four more days. The additional energy is probably transported in organised convective structures resulting in a mean vertical wind velocity unequal to zero at the tower location. The results show that convective transport contributes significantly to the surface energy budget for measurement heights as low as 2–3 m. Since these structures may be quasi-stationary, they can hardly be captured by a single-location measurement. The spatial EC set-up presented here is capable of quantifying contributions to the sensible heat flux from structures up to the scale of our spatial sensor network, which covered an area 3.5 × 3.5 km. For future experiments aiming at closing the energy balance, the spatial EC method should be employed to measure both the sensible and latent heat fluxes. Experimental determination of the horizontal advection of sensible and latent heat should also be considered, since such transport must occur due to convergence and divergence related to convection.     

Nitrous Oxide Flux Estimates for South-Eastern Australia

Nitrous oxide (N₂O) fluxes for south-easternAustralia have been estimated using a combination ofthe in situ N₂O and radon (Rn) measurementsmade at the Cape Grim Baseline Air Pollution Station,in north-west Tasmania. The average N₂O fluxesfrom the south-eastern mainland of Australia and fromTasmania over the nine years of record analysed (1985–1993) have beenfound to be 130 ± 30 kgN km^-2yr^-1 and 160 ± 45 kgN km^-2yr^-1respectively. These fluxes are larger than expectedand a significant dependence of the flux on rainfallis observed, with greater fluxes in the spring (October–December) andduring periods of positive SouthernOscillation Index. A large flux (1,300 ± 500kgN km^-2 yr^-1) from a nearby island (KingIsland) was also estimated from the data record,indicating a strong source, although the small size ofthe island means that it is not a significant sourcefor Australia.     

A Case Study of CO<Subscript>2</Subscript>, CO and Particles Content Evolution in the Suburban Atmospheric Boundary Layer Using a 2-μm Doppler DIAL,a 1-μm Backscatter Lidar and an Array of In-situ Sensors

A network of remote and in-situ sensors was deployed in a Paris suburb in order to evaluate the mesoscale evolution of the daily cycle of CO₂ and related tracers in the atmospheric boundary layer (ABL) and its relation to ABL dynamics and nearby natural and anthropogenic sources and sinks. A 2-μm heterodyne Doppler differential absorption lidar, which combines measurements of, (1) structure of the atmosphere, (2) radial velocity, and (3) CO₂ differential absorption was a particularly unique element of the observational array. We analyse the differences in the diurnal cycle of CO, CO₂, lidar reflectivity (a proxy for aerosol content) and H₂O using the lidar, airborne measurements in the free troposphere and ground-based measurements made at two sites located few kilometres apart. We demonstrate that vertical mixing dominates the early morning drawdown of CO and aerosol content trapped in the former nocturnal layer but not the H₂O and CO₂ mixing ratio variations. Surface fluxes, vertical mixing and advection all contribute to the ABL CO₂ mixing ratio decrease during the morning transition, with the relative importance depending on the rate and timing of ABL rise. We also show evidence that when the ABL is stable, small-scale (0.1-km vertical and 1-km horizontal) gradients of CO₂ and CO are large. The results illustrate the complexity of inferring surface fluxes of CO₂ from atmospheric budgets in the stable boundary layer.     

Local advection processes in the surface layer of the marginal ice zone

Micro-scale turbulent transport processes over the marginal ice zone have been studied by use of a two-dimensional numerical model. It has been found that internal boundary layers (IBLs) of horizontal mean velocity, temperature, and specific humidity reveal a near field and a far field. In the near field, the change in surface roughness dominates the height and growth rate of a velocity IBL. The change in surface heat flux governs the near field of a temperature and humidity IBL. In the far field, approximately x/¦L _*2 ¦ ~ 20, where L _*2 is the downstream Obukhov length, the downstream stratification more and more influences the growth rate of IBLs basically by modifying the eddy viscosity.Above more complex terrain consisting of an ensemble of ice strips and leads, a merging height h ^M develops, below which the horizontal variability of the surface modification is clearly observed; h ^M varies with the length scale L of surface modification approximately in proportion to h ^M /L ~ 1/20 – 1/10, as a rule of thumb. Above the merging height, an enveloping IBL exists, whose growth depends on the ice cover, i.e., on the integral of surface modification, but changes very little with L.Local advection of momentum, heat, and moisture clearly affects the local surface heat fluxes. Sensible and latent heat fluxes are found to show also a near and far field. However, if areally averaged surface fluxes are to be deduced from grid-averaged flow variables, then details of local advection can be neglected to a reasonably good approximation.     

Eddy Covariance Measurements On Mountain Slopes: The Advantage Of Surface-Normal Sensor Orientation Over A Vertical Set-Up

The measurement of scalar fluxes employing the eddy covariance method is a widely used experimental approach,for which the flow distortion due to obstacles (e.g., sensor mounts and mast)is a well-known but not fully solved problem. In order to reduce flow distortion we installed a sonic anemometer in a surface-normal orientationrelative to the terrain slope, and a second instrument in a verticalposition at a horizontal distance of 1.54 m from the first instrumentWe found a significant reduction in the rotation angle necessary for the coordinaterotation procedure in the x-z plane whencomputing 30-minute flux averages with the surface-normal orientation. In 91% of all cases this rotation angleremained within the angle of incidence of ±10° recommended bythe manufacturer. In contrast, only 24% of the measurements taken with the vertically mounted anemometer were obtained at an angle of incidencewithin ±10°, and 3% were outside the ±30° range specified for an acceptable operation.A data quality test based on the variance of vertical windspeed normalized with friction velocity (_w/u_*) revealed problems for application under stable conditions due to large uncertainties in the determination of the Monin–Obukhov stability parameter z/L. An alternative test using the bulk drag coefficient C_D revealed other problems related to the dependence of C_D on z/z₀, the measuring height normalized by the roughness length, which do not appear to be constantin complex terrain. With both tests, a tendency for a slightly improved dataquality was found for the surface normal set-up, which, however, proved statistically insignificant.It is concluded that the surface-normal set-up of a sonic anemometer significantly reduces flow distortion by thesensor head. Although the surface-normal mounting position therefore appears to be the preferred one, with decreased flow distortion and a slightly improved data quality, no significant differences in turbulent quantities were found between the two set-uppositions. Hence, the consequences for short-term measurements of massand energy fluxes with a surface-normal set-up in complex terrain appearto be relevant only if single flux events are to be inspected, while for long-term measurements of integrated fluxes both the surface-normaland vertical installation of the sonic anemometer are adequate,indicating that eddy covariance measurements in complex terrain are lessdelicate than expected.     

Rapid vertical sampling in stably stratified turbulent shear flows

A new method for obtaining instantaneous vertical profiles of two components of velocity and temperature in thermally stratified turbulent shear flows is presented. In this report, the design and construction of the traversing system will be discussed and results to date will be presented. The method is based on rapid vertical sampling whereby probe sensors are moved vertically at a high speed such that the measurement is approximately instantaneous. The system is designed to collect many measurements for the calculation of statistics such as vertical wave number spectra, mean square vertical gradients, and Thorpe scales. Results are presented for vertical profiles of temperature and compared to vertical profiles measured by single-point Eulerian time averages. The quality of the vertical profiles is found to be good over many profiles. Some comparisons are made between vertical measurements and standard single-point Eulerian measurements for three cases of stably stratified turbulent shear flow in which the initial microscale Reynolds number, Re_λ≈30. In case 1, the mean conditions are characterized by a gradient Richardson number, Ri_g=0.015, for which the flow is “unstable”, meaning the spatially evolving turbulent kinetic energy (E_k) grows. In case 2, Ri_g=0.095, for which the evolving turbulent kinetic energy is almost constant. In case 3, the flow is highly stable, where Ri_g=0.25 and E_k decays with spatial evolution. The measurements indicate anisotropy in the small scales for all cases. In particular, it is found that the ratio grows initially to a maximum and then decays with further evolution. Maximum Thorpe displacements are measured and compared to single-point measures of the vertical scales. It is found that vertical length scales derived from single-point measurements, such as the Ozmidov scale, L_O=(ε/N³)^1/2 and the overturn scale, L_t=θ′/(dT/dz), do not represent well the wide range of overturning scales which are actually present in the turbulence.