Vertical Velocity Variance Research Articles

Role of large eddies in the breakdown of the Reynolds analogy in an idealized mildly unstable atmospheric surface layer

While the breakdown in similarity between turbulent transport of heat and momentum (or Reynolds analogy) is not disputed in the atmospheric surface layer (ASL) under unstably stratified conditions, the causes of this breakdown are still debated. One reason for the breakdown is differences between how coherent structures transport heat and momentum, and their differing responses to increasing instability. Monin—Obukhov Similarity Theory (MOST), which hypothesizes that only local length‐scales play a role in ASL turbulent transport, implicitly assumes that large‐scale structures are inactive, despite their large energy content. Widely adopted mixing‐length models also rest on this assumption in the ASL. The difficulty of characterizing low‐wavenumber turbulent motions with field measurements motivates the use of high‐resolution Direct Numerical Simulation (DNS), which is free from subgrid‐scale parametrizations and adhoc assumptions near the boundary. Despite the low Reynolds number and idealized geometry of the DNS, DNS‐estimated MOST functions are consistent with ASL field experiments, as are low‐frequency features of the spectra. Parsimonious spectral models for MO stability correction functions for momentum (φm) and heat (φh) are derived, based on idealized vertical velocity variance and buoyancy variance spectra fit to the corresponding DNS spectra. For φm, a spectral model, based only on local length‐scales, matches DNS and field measurements well. In contrast, for φh, the model is substantially biased unless contributions from larger length‐scales are also included. These results are supported by sensitivity analyses based on field measurements that are independent of the DNS. They show that ASL heat transport is not MO‐similar, even under mild stratification, and in the absence of entrainment, non‐stationarity and canopy effects. It further suggests that the breakdown of the Reynolds analogy is at least partially caused by the influence of large eddies on turbulent heat transport.

Influence of Stokes Drift Decay Scale on Langmuir Turbulence

AbstractAccurately scaling Langmuir turbulence (LT) in the ocean surface boundary layer (OSBL) is critical for improving ocean, weather, and climate models. The physical processes by which the structure of LT depends on surface waves’ Stokes drift decay length scale are examined. An idealized model for OSBL turbulent kinetic energy (TKE) provides a conceptual framework with three physical processes: TKE transport, dissipation, and production by the Craik–Leibovich (CL) vortex force (VF) associated with the Stokes drift shear. TKE profiles depend on OSBL depth h, surface roughness length z0, and wavenumber k through the nondimensional parameters kh and kz0. These parameters determine the rate and length scale for the dissipation of TKE produced by the CL-VF. For kz0 ≫ 1, TKE input by the CL-VF is governed by a surface flux with TKE rapidly decaying with depth. Only for kz0 < 1 can TKE penetrate deeper into the OSBL, with the TKE penetration depth controlled by kh. Turbulence-resolving large-eddy simulation results support this conceptual framework and indicate that the dominant Langmuir cell size scales with (kh)−1. Within the depth of dominant Langmuir cells, TKE dissipation is approximately balanced by CL-VF production. Shorter waves contribute less to deeper vertical velocity variance 〈w2〉 because the CL-VF is less effective in generating larger-scale LT. Depth-averaged 〈w2〉 scales with a modified Langmuir number Laϕ = (u*/usϕ)1/2, where u* denotes the water-side surface friction velocity and usϕ is a depth-integrated weighted Stokes drift shear or, equivalently, a spectrally filtered surface Stokes drift.

Identification of Reservoir Thickness and Estimation of Hydrocarbon Reservation used the Pre-Stack Depth Migration (PSDM) in Cikung Area

Information from the geological data shows that the Cikung field has the complex structure which is characterized by carbonate reef. Therefore, it is required in the domain depth migration (PSDM) using interval velocity model is sensitive to variations in vertical and horizontal velocity to identify the thickness of the reservoir and estimation of hydrocarbon reserves. Stages include making the depth structure maps, map making the thickness of the reservoir, and the calculation of hydrocarbon reserves. The results showed that the thickness of the reservoir in Cikung Field, estimated to range from 71 meters to 175 meters with the prospect reservoir is at a depth of 1216 meters to 1247 meters from the surface. Hydrocarbon deposits (oil) is estimated at 1,134 × 106 𝑆𝑇𝐵 or 1,311 × 108 kilo litre by the porosity of 22.6% and 70.7% water saturation.

Full waveform tomography to resolve desert alluvium

This paper presents shear wave velocity distribution developed from full waveform tomography on a seismic dataset recorded in the northwest Las Vegas Valley, Nevada. The setting is a deep and variable alluvial sequence, encompassing gravel, sand, silt and clay, with some carbonate cementation. Seismic data were collected along a ~144-m long array using 72 equally spaced geophones, with an accelerated drop-weight source applied at 6-m intervals. The dataset was interpreted using an advanced 2-D full waveform tomography method. The same dataset was also analyzed using Rayleigh wave dispersion, following the multichannel analysis of surface waves (MASW) method, to develop a single, 1-D velocity profile. Results from the MASW method indicated a mild velocity inversion. The waveform tomography method was able to characterize vertical and lateral velocity variations along the array length to 25m depth. Results showed significantly contrasting high- and low-velocity layers. While depth-averaged shear wave velocities from the two methods were similar, the degree of variability indicated using FWI is supported by strong variation in Standard Penetration Test (SPT) blow counts from a nearby drillhole and by expectations from sediment lithology.

A comparison of vertical velocity variance measurements from wind profiling radars and sonic anemometers

Abstract. Observations of turbulence in the planetary boundary layer are critical for developing and evaluating boundary layer parameterizations in mesoscale numerical weather prediction models. These observations, however, are expensive and rarely profile the entire boundary layer. Using optimized configurations for 449 and 915 MHz wind profiling radars during the eXperimental Planetary boundary layer Instrumentation Assessment (XPIA), improvements have been made to the historical methods of measuring vertical velocity variance through the time series of vertical velocity, as well as the Doppler spectral width. Using six heights of sonic anemometers mounted on a 300 m tower, correlations of up to R2 = 0. 74 are seen in measurements of the large-scale variances from the radar time series and R2 = 0. 79 in measurements of small-scale variance from radar spectral widths. The total variance, measured as the sum of the small and large scales, agrees well with sonic anemometers, with R2 = 0. 79. Correlation is higher in daytime convective boundary layers than nighttime stable conditions when turbulence levels are smaller. With the good agreement with the in situ measurements, highly resolved profiles up to 2 km can be accurately observed from the 449 MHz radar and 1 km from the 915 MHz radar. This optimized configuration will provide unique observations for the verification and improvement to boundary layer parameterizations in mesoscale models.

Elastic internal multiple analysis and attenuation using Marchenko and interferometric methods

Conventional seismic processing aims to create data that contain only primary reflections, whereas real seismic recordings also contain multiples. As such, it is desirable to predict, identify, and attenuate multiples in seismic data. This task is more difficult in elastic (solid) media because mode conversions create families of internal multiples not present in the acoustic case. We have developed a method to predict prestack internal multiples in general elastic media based on the Marchenko method and convolutional interferometry. It can be used to identify multiples directly in prestack data or migrated sections, as well as to attenuate internal multiples by adaptively subtracting them from the original data set. We developed the method on two synthetic data sets, the first composed of horizontal density layers and constant velocities, and the second containing horizontal and vertical density and velocity variations. The full-elastic method is computationally expensive and ideally uses data components that are not usually recorded. We therefore tested an acoustic approximation to the method on the synthetic elastic data from the second model and find that although the spatial resolution of the resulting image is reduced by this approximation, it provides images with relatively fewer artifacts. We conclude that in most cases where cost is a factor and we are willing to sacrifice some resolution, it may be sufficient to apply the acoustic version of this demultiple method.

Applying full-azimuth angle domain imaging to study carbonate reefs at great depths

The main goal of the seismic surveys conducted in the target area in Kazakhstan was to image and detect a major Devonian carbonate barrier reef and to characterize the density and orientation of its main fracture zones. This area is unique in that the carbonate layers occur at great depths, from 5700 to 8500 m, with a complex structure of overburden layers of interbedded clay, sandstone, salt and others (Table 1). These different lithology plays and morphology rocks create strong vertical and lateral velocity variations, resulting in a complex seismic wave phenomenon. In addition, the target carbonate structures contain heterogeneous objects and require high-quality processing of the recorded data. Under such conditions, it is crucial to use full-azimuth, long-offset and dense (high fold) acquisition patterns. Advanced processing sequence tools and high-end depth migrations are required to handle both strong heterogeneity and azimuthal anisotropy effects. Traditional Kirchhoff migrations, even the most accurate ones (wavefront reconstruction, beam), have not been able to provide the required image quality and level of detail required at the target zones (Figure 1). Conventional Kirchhoff migrations generate surface offset-azimuth/offset domain common image gathers (CIG). In this particularly complex area, the correlation between the surface offset-azimuths and the actual, in situ, subsurface slowness-azimuths (azimuth of the incidence/reflected ray pairs at the image points) is relatively poor, leading to significant errors in the estimation of fracture orientation. Additionally, Kirchhoff migrations do not account for multi-pathing (multi-wave path solutions between image points and surface source/receiver locations), which is essential when imaging in areas involving complex wave phenomena. Kirchhoff beam migrations map surface beams with given directions backward to the subsurface. Since they account for the multi-pathing, they provide better results. However, they normally generate the same type of surface offset-azimuth/offset CIGs and therefore cannot be accurately used for azimuthal studies. They also cannot ensure sufficient subsurface illumination, especially at the complex target subsurface regions.

The Effect of Stokes Drift and Transient Rip Currents on the Inner Shelf. Part II: With Stratification

AbstractThis is Part II of a two-part study focused on Stokes drift and transient rip current (TRC) effects on the unstratified (Part I) and stratified (this paper) inner shelf. Part I focuses on funwaveC–Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) coupling and TRC effects on mixing and exchange on an unstratified inner shelf. Here, two simulations (R3 and R4) are performed on a stratified inner shelf and surfzone with typical bathymetry, stratification, and wave conditions. R3 is a COAWST-only simulation (no TRCs), while R4 has funwaveC–COAWST coupling (with TRCs). In R4, TRCs lead to patchy, near-surface cooling, vertical isotherm displacement, and increased water column mixing. For both R3 and R4, the mean Lagrangian circulation has two nearly isolated surfzones and inner-shelf overturning circulation cells, with a stronger, R4, inner-shelf circulation cell. The R4, inner-shelf, vertical velocity variability is 2–3 times stronger than a simulation with TRCs and no stratification. Relative to R3, R4 eddy diffusivity is strongly elevated out to three surfzone widths offshore due to TRCs and TRC-induced density overturns. The R4 inner-shelf stratification is reduced nearshore, and mean isotherms slope more strongly than R3 because of the TRC-enhanced irreversible mixing. At six surfzone widths offshore, both R3 and R4 are in geostrophic balance, explaining the stratified (summertime) observed deviation from Stokes–Coriolis balance. In this region, baroclinic pressure gradients induced by sloping isotherms induce an alongshore geostrophic jet offshore, strongest in R4. In R4, TRCs result in an enhanced (2–10 times) cross-shore exchange velocity across the entire inner shelf, relative to R3. Accurate, stratified, inner-shelf simulations of pollution, larval, or sediment transport must include transient rip currents.

Mechanism of horizontal mass- and salt-exchange between the waters of continental slope and central part of the Black Sea

Estimation of horizontal salt transfer from the central part of the Black Sea to the continental slope is carried out based on numerical modeling. A box hydrodynamic model of the Black Sea based on the POM model is used. Horizontal transport caused by seasonal variability of vertical velocity in the central part of the Black Sea and horizontal turbulent diffusion is considered. In this work we demonstrate that, in the 0- to 30-m layer on average for the year, there is an outflow of salty water under the impact of both factors. In the surface layer, the transport caused by seasonal variability contributes significantly (30%). In the 80- to 150-m layer on average over the year, the waters from the continental slope arrive into the central part of the basin; they are characterized by higher salinity due to the inflow of the Marmara Sea waters across the Bosporus Strait.

Improvement of vertical velocity statistics measured by a Doppler lidar through comparison with sonic anemometer observations

Abstract. Since turbulence measurements from Doppler lidars are being increasingly used within wind energy and boundary-layer meteorology, it is important to assess and improve the accuracy of these observations. While turbulent quantities are measured by Doppler lidars in several different ways, the simplest and most frequently used statistic is vertical velocity variance (w′2) from zenith stares. However, the competing effects of signal noise and resolution volume limitations, which respectively increase and decrease w′2, reduce the accuracy of these measurements. Herein, an established method that utilises the autocovariance of the signal to remove noise is evaluated and its skill in correcting for volume-averaging effects in the calculation of w′2 is also assessed. Additionally, this autocovariance technique is further refined by defining the amount of lag time to use for the most accurate estimates of w′2. Through comparison of observations from two Doppler lidars and sonic anemometers on a 300 m tower, the autocovariance technique is shown to generally improve estimates of w′2. After the autocovariance technique is applied, values of w′2 from the Doppler lidars are generally in close agreement (R2 ≈ 0.95 − 0.98) with those calculated from sonic anemometer measurements.

Closure Schemes for Stably Stratified Atmospheric Flows without Turbulence Cutoff

Abstract Two recently proposed turbulence closure schemes are compared against the conventional Mellor–Yamada (MY) model for stably stratified atmospheric flows. The Energy- and Flux-Budget (EFB) approach solves the budgets of turbulent momentum and heat fluxes and turbulent kinetic and potential energies. The Cospectral Budget (CSB) approach is formulated in wavenumber space and integrated across all turbulent scales to obtain flow variables in physical space. Unlike the MY model, which is subject to a “critical gradient Richardson number,” both EFB and CSB models allow turbulence to exist at any gradient Richardson number and predict a saturation of flux Richardson number () at sufficiently large . The CSB approach further predicts the value of and reveals a unique expression linking the Rotta and von Kármán constants. Hence, all constants in the CSB model are nontunable and stability independent. All models agree that the dimensionless sensible heat flux decays with increasing . However, the decay rate and subsequent cutoff in the MY model appear abrupt. The MY model further exhibits an abrupt cutoff in the turbulent stress normalized by vertical velocity variance, while the CSB and EFB models display increasing trends. The EFB model produces a rapid increase in the ratio of turbulent potential energy and vertical velocity variance as is approached, suggesting a strong self-preservation mechanism. Vertical anisotropy in the turbulent kinetic energy is parameterized in different ways in MY and EFB, but this consideration is not required in CSB. Differences between EFB and CSB model predictions originate from how the vertical anisotropy is specified in the EFB model.

Intensification of Upper-Ocean Submesoscale Turbulence through Charney Baroclinic Instability

AbstractThis study focuses on the description of an oceanic variant of the Charney baroclinic instability, arising from the joint presence of (i) an equatorward buoyancy gradient that extends from the surface into the ocean interior and (ii) reduced subsurface stratification, for example, as produced by wintertime convection or subduction. This study analyzes forced dissipative simulations with and without Charney baroclinic instability (C-BCI). In the former, C-BCI strengthens near-surface frontal activity with important consequences in terms of turbulent statistics: increased variance of vertical vorticity and velocity and increased vertical turbulent fluxes. Energetic consequences are explored. Despite the atypical enhancement of submesoscale activity in the simulation subjected to C-BCI, and contrary to several recent studies, the downscale energy flux at the submesoscale en route to dissipation remains modest in the flow energetic equilibration. In particular, it is modest vis à vis the global energy input to the system, the eddy kinetic energy input through conversion of available potential energy, and the classical inverse cascade of kinetic energy. Linear stability analysis suggests that the southern flank of the Gulf Stream may be conducive to oceanic Charney baroclinic instability in spring, following mode water formation and upper-ocean destratification.

Estimate of Boundary-Layer Depth Over Beijing, China, Using Doppler Lidar Data During SURF-2015

Planetary boundary-layer (PBL) structure was investigated using observations from a Doppler lidar and the 325-m Institute of Atmospheric Physics (IAP) meteorological tower in the centre of Beijing during the summer 2015 Study of Urban-impacts on Rainfall and Fog/haze (SURF-2015) field campaign. Using six fair-weather days of lidar and tower data under clear to cloudy skies, we evaluate the ability of the Doppler lidar to probe the urban boundary-layer structure, and then propose a composite method for estimating the diurnal cycle of the PBL depth using the Doppler lidar. For the convective boundary layer (CBL), a threshold method using vertical velocity variance (sigma _w^2 >0.1,hbox {m}^{2}hbox {s}^{-2}) is used, since it provides more reliable CBL depths than a conventional maximum wind-shear method. The nocturnal boundary-layer (NBL) depth is defined as the height at which sigma _w^2 decreases to 10 % of its near-surface maximum minus a background variance. The PBL depths determined by combining these methods have average values ranging from approx 270 to approx 1500 m for the six days, with the greatest maximum depths associated with clear skies. Release of stored and anthropogenic heat contributes to the maintenance of turbulence until late evening, keeping the NBL near-neutral and deeper at night than would be expected over a natural surface. The NBL typically becomes more shallow with time, but grows in the presence of low-level nocturnal jets. While current results are promising, data over a broader range of conditions are needed to fully develop our PBL-depth algorithms.

Scalar turbulent behavior in the roughness sublayer of an Amazonian forest

Abstract. An important current problem in micrometeorology is the characterization of turbulence in the roughness sublayer (RSL), where most of the measurements above tall forests are made. There, scalar turbulent fluctuations display significant departures from the predictions of Monin–Obukhov similarity theory (MOST). In this work, we analyze turbulence data of virtual temperature, carbon dioxide, and water vapor in the RSL above an Amazonian forest (with a canopy height of 40 m), measured at 39.4 and 81.6 m above the ground under unstable conditions. We found that dimensionless statistics related to the rate of dissipation of turbulence kinetic energy (TKE) and the scalar variance display significant departures from MOST as expected, whereas the vertical velocity variance follows MOST much more closely. Much better agreement between the dimensionless statistics with the Obukhov similarity variable, however, was found for the subset of measurements made at a low zenith angle Z, in the range 0° < |Z| < 20°. We conjecture that this improvement is due to the relationship between sunlight incidence and the “activation–deactivation” of scalar sinks and sources vertically distributed in the forest. Finally, we evaluated the relaxation coefficient of relaxed eddy accumulation: it is also affected by zenith angle, with considerable improvement in the range 0° < |Z| < 20°, and its values fall within the range reported in the literature for the unstable surface layer. In general, our results indicate the possibility of better stability-derived flux estimates for low zenith angle ranges.

Doppler Lidar Observations over a High Altitude Mountainous Site Manora Peak in the Central Himalayan Region

The RAWEX-GVAX field campaign has been carried out from June 2011 to March 2012 over a high altitude site Manora Peak, Nainital (29.4 degrees N; 79.2 degrees E; 1958 m amsl) in the central Himalayas to assess the impacts of absorbing aerosols on atmospheric thermodynamics and clouds. This paper presents the preliminary results of the observations and data analysis of the Doppler Lidar, installed at Nainital. Strong updrafts with vertical winds in the range of similar to 2-4 ms(-1) occurred during the daytime and throughout the season indicating thermally driven convection. On the other hand during nighttime, weak downdrafts persisted during stable conditions. Plan Position Indicator scan of Doppler Lidar showed north-northwesterly winds in the boundary layer. The mixing layer height, derived from the vertical velocity variance, showed diurnal variations, in the range similar to 0.7-1 km above ground level during daytime and very shallow during nighttime.

Variability in ice motion at a land-terminating Greenlandic outlet glacier: the role of channelized and distributed drainage systems

ABSTRACTWe use a combination of field observations and hydrological modelling to examine the mechanisms through which variability in meltwater input affects ice motion at a land-terminating Greenlandic outlet glacier. We find a close agreement between horizontal ice velocity, vertical ice velocity and modelled subglacial water pressure over the course of a melt season. On this basis, we argue that variation in horizontal and vertical ice velocity primarily reflects the displacement of basal ice during periods of cavity expansion and contraction, a process itself driven by fluctuations in basal water pressure originating in subglacial channels. This process is not captured by traditional sliding laws linking water pressure and basal velocity, which may hinder the simulation of realistic diurnal to seasonal variability in ice velocity in coupled models of glacial hydrology and dynamics.

Near-Surface Velocity Model To Enhance PSDM Seismic Imaging of Dukhan Field

This article, written by Special Publications Editor Adam Wilson, contains highlights of paper IPTC 18293, “Near-Surface Velocity Model of Dukhan Field From Multiphysics Survey To Enhance PSDM Seismic Imaging,” by K. Setiyono, Qatar Petroleum; S. Gallo, C. Boulanger, F. Bruere, F. Moreau, and B. Rondeleux, CGG; and J. Snow, Qatar Petroleum, prepared for the 2015 International Petroleum Technology Conference, Doha, Qatar, 7-9 December. The paper has not been peer reviewed. Dukhan field affords an opportunity to assess alternative near-surface velocity-modeling work flows that specifically integrate microgravity (MGR) and vertical electrical sounding (VES). A pilot study of the Dukhan seismic data demonstrates that integration of MGR/VES into prestack depth migration (PSDM) improves the fidelity of seismic images and results in an improved depth match. The resultant MGR/VES PSDM pilot model improved the velocity model between the surface and top Simsima formation on the basis of visible impact on the deeper structure obtained from imaging and better fit between the seismic depth horizons and well measured depths. Introduction The Dukhan field (Fig. 1) along the southwest coast of Qatar is an elongated anticline that has significant hydrocarbon reservoirs in the Jurassic Arab and deeper formations. In 2007, Qatar Petroleum invested in a new full-field, land and shallow-water 3D-seismic survey of the area. Despite considerable and continuous progress in seismic technology, obtaining an accurate, quantitative, well-focused seismic depth image at reservoir level remains challenging. In land seismic, a poor characterization of the near-surface geology is detrimental to the depth imaging of underlying structures. The proper knowledge of the near surface to correct for its perturbations onto the seismic wavefield remains one of the main challenges to be addressed. The near surface influences the kinematics as well as the dynamics of the source-generated seismic waves. Knowledge of the near surface to correct for its induced distortions is critical for a reliable use of seismic images, structures, and amplitude at reservoir level. The near-surface lateral and vertical velocity variations often are not accessible with seismic methods, even to the seismic refraction survey, because their depth of investigation is limited to the top of the shallowest hard layer. In addition, rapid lateral changes and large velocity contrasts are unfavorable to seismic refraction.

Mean-velocity profile of smooth channel flow explained by a cospectral budget model with wall-blockage

A series of recent studies has shown that a model of the turbulent vertical velocity variance spectrum (Fvv) combined with a simplified cospectral budget can reproduce many macroscopic flow properties of turbulent wall-bounded flows, including various features of the mean-velocity profile (MVP), i.e., the “law of the wall”. While the approach reasonably models the MVP’s logarithmic layer, the buffer layer displays insufficient curvature compared to measurements. The assumptions are re-examined here using a direct numerical simulation (DNS) dataset at moderate Reynolds number that includes all the requisite spectral and co-spectral information. Starting with several hypotheses for the cause of the “missing” curvature in the buffer layer, it is shown that the curvature deficit is mainly due to mismatches between (i) the modelled and DNS-observed pressure-strain terms in the cospectral budget and (ii) the DNS-observed Fvv and the idealized form used in previous models. By replacing the current parameterization for the pressure-strain term with an expansive version that directly accounts for wall-blocking effects, the modelled and DNS reported pressure-strain profiles match each other in the buffer and logarithmic layers. Forcing the new model with DNS-reported Fvv rather than the idealized form previously used reproduces the missing buffer layer curvature to high fidelity thereby confirming the “spectral link” between Fvv and the MVP across the full profile. A broad implication of this work is that much of the macroscopic properties of the flow (such as the MVP) may be derived from the energy distribution in turbulent eddies (i.e., Fvv) representing the microstate of the flow, provided the link between them accounts for wall-blocking.

Observed spatiotemporal variability of boundary-layer turbulence over flat, heterogeneous terrain

Abstract. In the spring of 2013, extensive measurements with multiple Doppler lidar systems were performed. The instruments were arranged in a triangle with edge lengths of about 3 km in a moderately flat, agriculturally used terrain in northwestern Germany. For 6 mostly cloud-free convective days, vertical velocity variance profiles were calculated. Weighted-averaged surface fluxes proved to be more appropriate than data from individual sites for scaling the variance profiles; but even then, the scatter of profiles was mostly larger than the statistical error. The scatter could not be explained by mean wind speed or stability, whereas time periods with significantly increased variance contained broader thermals. Periods with an elevated maximum of the variance profiles could also be related to broad thermals. Moreover, statistically significant spatial differences of variance were found. They were not influenced by the existing surface heterogeneity. Instead, thermals were preserved between two sites when the travel time was shorter than the large-eddy turnover time. At the same time, no thermals passed for more than 2 h at a third site that was located perpendicular to the mean wind direction in relation to the first two sites. Organized structures of turbulence with subsidence prevailing in the surroundings of thermals can thus partly explain significant spatial variance differences existing for several hours. Therefore, the representativeness of individual variance profiles derived from measurements at a single site cannot be assumed.

Measurements of Flow Distortion within the IRGASON Integrated Sonic Anemometer and CO$$_2$$/H$$_2$$O Gas Analyzer

Wind-tunnel and field measurements are analyzed to investigate flow distortion within the IRGASON integrated sonic anemometer and CO\(_2\)/H\(_2\)O gas analyzer as a function of wind speed, wind direction and attack angle. The wind-tunnel measurements are complimentary to the field measurements, and the dependence of the wind-tunnel mean-wind-component flow-distortion errors on wind direction agrees well with that of the field measurements. The field measurements exhibit significant overestimation of the crosswind variance and underestimation of the momentum flux with respect to an adjacent CSAT3 sonic, as well as a transfer of turbulent kinetic energy from the streamwise wind component to the cross-stream wind components. In contrast, we find attenuation of only a few percent in the vertical velocity variance and the vertical flux of sonic temperature. The attenuation of the fluxes appears to be caused to a large extent by decorrelation between the horizontal and vertical-velocity components and between the vertical velocity and sonic temperature. Additional flow distortion due to transducer shadowing reduces to some extent the overestimation, but also increases the underestimation of the IRGASON turbulence statistics.