Search NASA⌕ Search

SEARCH · Search NASA

Results for “variance component analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Full covariance of CMB and lensing reconstruction power spectra

CMB and lensing reconstruction power spectra are powerful probes of cosmology. However, they are correlated, since the CMB power spectra are lensed, and the lensing reconstruction is constructed using CMB multipoles. This paper discusses how we perform a full analysis of the auto- and cross-covariances, including polarization power spectra and minimum-variance lensing estimators, and compare with simulations of idealized future CMB-S4 observations. Covariances sourced by fluctuations in the unlensed CMB and instrumental noise can largely be removed by using a realization-dependent subtraction of lensing reconstruction noise, leaving a relatively simple covariance model that is dominated by lensing-induced terms and well described by a small number of principal components. The correlations between the CMB and lensing power spectra will be detectable at the level of $\sim 5σ$ for a CMB-S4 mission, and neglecting them could underestimate some parameter error bars by several tens of percent. However, we found that the inclusion of external priors or data sets to estimate parameter error bars can make the impact of the correlations almost negligible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Fine structure and time variation of the quiet sun at 1.3 cm

The Hat Creek two-element interferometer has been used to study the quiet sun at 22 GHz. A statistical analysis of output of the interferometer clearly shows the existence of time variations on the quiet sun with time scales of approximately 180 s. The observations suggest that the fine structure on the quiet sun might consist of two components - one which varies with the time scales of approximately 3 min and the other being relatively stable. The average visibility amplitude indicates that the fine structure on the sun has a typical angular size of approximately 6 sec. The observation that the variance and the mean of the visibility amplitudes depend in the same way on the projected baseline suggests that the transient sources have angular size similar to the average size of the fine structure on the quiet sun. Power spectra of the output of the interferometer show no significant periodicity.

Rao, A. P.↗
NASA NTRS2014

Quantitative Estimation of the Impact of European Teleconnections on Interannual Variation of East Asian Winter Temperature and Monsoon

The impact of European teleconnections including the East AtlanticWest Russia (EA-WR), the Scandinavia (SCA), and the East Atlantic (EA) on East Asian winter temperature variability was quantified and compared with the combined effect of the Arctic Oscillation (AO), the Western Pacific (WP), and the El-Nino Southern Oscillation (ENSO), which are originated in the Northern Hemispheric high-latitudes or the Pacific. Three European teleconnections explained 22-25 percent of the total monthly upper-tropospheric height variance over Eurasia. Regression analysis revealed warming by EA-WR and EA and cooling by SCA over mid-latitude East Asia during their positive phase and vice versa. Temperature anomalies were largely explained by the advective temperature change process at the lower troposphere. The average spatial correlation over East Asia (90-180E, 10-80N) for the last 34 winters between observed and reconstructed temperature comprised of AO, WP and ENSO effect (AWE) was approximately 0.55, and adding the European teleconnection components (ESE) to the reconstructed temperature improved the correlation up to approximately 0.64. Lower level atmospheric structure demonstrated that approximately five of the last 34 winters were significantly better explained by ESE than AWE to determine East Asian seasonal winter temperatures. We also compared the impact between EA-WR and AO on the 1) East Asian winter monsoon, 2) cold surge, and 3) the Siberian high. These three were strongly coupled, and their spatial features and interannual variation were somewhat better explained by EA-WR than AO. Results suggest that the EA-WR impact must be treated more importantly than previously thought for a better understanding of East Asian winter temperature and monsoon variability.

East Asian Winter Monsoon↗
DOE OSTI2016

Planck 2015 results: XVI. Isotropy and statistics of the CMB

In this paper, we test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The “Cold Spot” is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Finally, where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A new process sensitivity index to identify important system processes under process model and parametric uncertainty

Hydrological models are always composed of multiple components that represent processes key to intended model applications. When a process can be simulated by multiple conceptual-mathematical models (process models), model uncertainty in representing the process arises. While global sensitivity analysis methods have been widely used for identifying important processes in hydrologic modeling, the existing methods consider only parametric uncertainty but ignore the model uncertainty for process representation. To address this problem, this study develops a new method to probe multimodel process sensitivity by integrating the model averaging methods into the framework of variance-based global sensitivity analysis, given that the model averaging methods quantify both parametric and model uncertainty. A new process sensitivity index is derived as a metric of relative process importance, and the index includes variance in model outputs caused by uncertainty in both process models and model parameters. For demonstration, the new index is used to evaluate the processes of recharge and geology in a synthetic study of groundwater reactive transport modeling. The recharge process is simulated by two models that converting precipitation to recharge, and the geology process is also simulated by two models of different parameterizations of hydraulic conductivity; each process model has its own random parameters. The new process sensitivity index is mathematically general, and can be applied to a wide range of problems in hydrology and beyond.

Dai, Heng↗
DOE OSTI2023

O’Hare Airport Short-Term Ground Transportation Modal Demand Forecast Using Gaussian Processes

Here, the principal objective of this study is to analyze the spatial and temporal variation of ground transportation airport demand and provide demand forecast to inform planning capability and explore alternatives for investments to accommodate airport growth. Because of its good adaptability and strong generalization ability for dealing with high-dimensional input, small-sample, and nonlinear spatial data, Gaussian process (GP) regression is used to provide forecast estimates using data from transportation network company (TNC) trips and urban rail passengers at Chicago's O'Hare International Airport. TNC airport trips differ significantly, with three times more distance, more than twice the travel time, and half of the share requests compared with nonairport trips. This highlights the need for separate demand models. Hourly analysis of the rail service indicates that this is likely heavily used by airport workers, whereas TNC services focus on travelers because of variations in the peak demand hours. Heteroscedastic GP regression is implemented because of differences in trip variance between night and day hours. Estimates are given for weekdays and weekend trips, and the 95% confidence intervals are calculated. The introduction of flight schedule information into the models shows marginal improvements in their performance. However, fitting a GP regression becomes computationally expensive with increased sample size and the introduction of spatial components. Transportation planners and policymakers can use the results and methods implemented in this study to optimize transportation assets and provide long-range simulations of the current and future conditions in the area.

42 ENGINEERING↗
NASA NTRS1983

LANDSAT-D thematic mapper image dimensionality reduction and geometric correction accuracy

When principal component analysis of a subscene of a section of the Sacramento River showed lower correlation among the TM spectral components that were observed for the MSS spectral components, principal component analysis was applied to a LANDSAT 2 MSS subscene of the same area for comparison purposes. Correlation coefficient matrices indicate the pairwise similarity and correlation of the data for the spectral components. The principal components transformation matrix, indicates the weights applied to the original components to generate the transformed components. The first two TM components can be described as visible and near infrared. For the MSS data, the first transformed component is roughly the average of the four original components. The second transformed component is roughly the difference between the visible and infrared components. Tables show that 97.0% of the variance in an MSS image is contained in only two transformed components.

Ford, G. E.↗
DOE OSTI2024

Multilevel Monte Carlo Estimators For Derivative-Free Optimization Under Uncertainty

Optimization is a key tool for scientific and engineering applications; however, in the presence of models affected by uncertainty, the optimization formulation needs to be extended to consider statistics of the quantity of interest. Optimization under uncertainty (OUU) deals with this endeavor and requires uncertainty quantification analyses at several design locations; i.e., its overall computational cost is proportional to the cost of performing a forward uncertainty analysis at each design location. An OUU workflow has two main components: an inner loop strategy for the computation of statistics of the quantity of interest, and an outer loop optimization strategy tasked with finding the optimal design, given a merit function based on the inner loop statistics. Here, in this work, we propose to alleviate the cost of the inner loop uncertainty analysis by leveraging the so-called multilevel Monte Carlo (MLMC) method, which is able to allocate resources over multiple models with varying accuracy and cost. The resource allocation problem in MLMC is formulated by minimizing the computational cost given a target variance for the estimator. We consider MLMC estimators for statistics usually employed in OUU workflows and solve the corresponding allocation problem. For the outer loop, we consider a derivative-free optimization strategy implemented in the SNOWPAC library; our novel strategy is implemented and released in the Dakota software toolkit. We discuss several numerical test cases to showcase the features and performance of our approach with respect to its Monte Carlo single fidelity counterpart.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Improving Self-Driving Labs: Quantifying System-Level Experiment Repeatability and Broadening Instrument-Level Compatibility

Modular Autonomous Research System (MARS) is a self-driving laboratory (SDL) which performs wet-lab science with peptide-lanthanide combinations in an automated and, ultimately, an autonomous manner to aid in soil analysis for domestic lithium mining. Autonomous experimentation involves automated experimentation, experiment planning, and active learning. MARS consists of a 6-axis robotic arm (UR5e) on a linear rail, pipette robots (Opentrons 2), and microplate readers. These components transport, operate on, and collect data with chemical solutions in standard labware. For effective autonomy, MARS must perform system-level labware operations repeatably, plan experiments autonomously, and be portable between research-domains. Repeatability is evaluated by labware placement precision, such that future operations can properly locate labware, as well as the elapsed time, so that low variance mean estimates of experiment duration can inform high-level researcher decision making. Autonomous experiment planning is the next step to decouple experimentation from human management; however, there is a conflict between the ideal system-level experiment goals and the constraints imposed by instruments’ limitations. Sub-domain portability is a long-term goal to extend MARS’ research beyond the chemistry of peptide-lanthanide binding to other sub-domains without having to invest significant overhead to system retrofitting. To address these goals, we manually trained the robotic arm labware placement and modelled statistical failurerate and uncertainty Additionally, we benchmarked the duration and variance of each experiment sub-operation as a heuristic for research decision making. Next, we use a parameterized geometric program (PGP) approach to design experiments that optimize system-level objectives and satisfy instrument-level constraints. Lastly, we proposed a Python framework to maximize MARS’ extensibility to other scientific sub-domains through a JSON-based experiment specification.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Linear minimum variance filters applied to carrier tracking

A new approach is taken to the problem of tracking a fixed amplitude signal with a Brownian-motion phase process. Classically, a first-order phase-lock loop (PLL) is used; here, the problem is treated via estimation of the quadrature signal components. In this space, the state dynamics are linear with white multiplicative noise. Therefore, linear minimum-variance filters, which have a particularly simple mechanization, are suggested. The resulting error dynamics are linear at any signal/noise ratio, unlike the classical PLL. During synchronization, and above threshold, this filter with constant gains degrades by 3 per cent in output rms phase error with respect to the classical loop. However, up to 80 per cent of the maximum possible noise improvement is obtained below threshold, where the classical loop is nonoptimum, as demonstrated by a Monte Carlo analysis. Filter mechanizations are presented for both carrier and baseband operation.

Gustafson, D. E.↗
NASA NTRS1999

Clustering of the Diffuse Infrared Light from the COBE DIRBE Maps: Power Spectrum Analysis and Excess Isotropic Component of Fluctuations - 3

The cosmic infrared background (CIB) radiation is the cosmic repository for energy release throughout the history of the universe. The spatial fluctuations of the CIB resulting from galaxy clustering are expected to be at least a few percent on scales of a degree, depending on the luminosity and clustering history of the early universe. Using the all-sky data from the COBE DIRBE instrument at wavelengths 1.25 - 100 microns we attempt to measure the CIB fluctuations. In the near-IR, foreground emission is dominated by small scale structure due to stars in the Galaxy. There we find a strong correlation between the amplitude of the fluctuations and Galactic latitude after removing bright foreground stars. Using data outside the Galactic plane (absolute value of b > 20 deg) and away from the center (90 deg < l < 270 deg) we extrapolate the amplitude of the fluctuations to cosec absolute value of b = 0. We find a positive intercept of delta.F(sub rms) = 15.5(sup +3.7, sub -7.0), 5.9(sup +1.6, sub -3.7), 2.4(sup +0.5, sub -0.9), 2.0(sup +0.25, sub -0.5) nW/sq m.sr at 1.25, 2.2, 3.5 and 4.9 microns respectively, where the errors are the range of 92% confidence limits. For color subtracted maps between band 1 and 2 we find the isotropic part of the fluctuations at 7.6(sup +1.2, sub -2.4) nW/sq m.sr. Based on detailed numerical and analytic models, this residual is not likely to originate from the Galaxy, our clipping algorithm, or instrumental noise. We demonstrate that the residuals from the fit used in the extrapolation are distributed isotropically and suggest that this extra variance may result from structure in the CIB. We also obtain a positive intercept from a linear combination of maps at 1.25 and 2.2 microns. For 2 deg < theta < 15 deg, a power-spectrum analysis yields limits of (theta/5 deg) x delta.F(sub rms)(theta) < 6, 2.5, 0.8, 0.5 nW/sq m.sr at 1.25, 2.2, 3.5 and 4.9 microns respectively. From 10 - 100 microns, the dominant foregrounds are emission by dust in the Solar system and the Galaxy. There the upper limits on the CIB fluctuations are below 1 nW/sq m.sr and are lowest (< equal 0.5 nW/sq m.sr) at 25 microns.

Kashlinsky, A.↗
DOE OSTI2023

Optimal Frequency-domain Analysis for Spacecraft Time Series: Introducing the Missing-data Multitaper Power Spectrum Estimator

While the Lomb–Scargle periodogram is foundational to astronomy, it has a significant shortcoming: the variance in the estimated power spectrum does not decrease as more data are acquired. Statisticians have a 60 yr history of developing variance-suppressing power spectrum estimators, but most are not used in astronomy because they are formulated for time series with uniform observing cadence and without seasonal or daily gaps. Here we demonstrate how to apply the missing-data multitaper power spectrum estimator to spacecraft data with uniform time intervals between observations but missing data during thruster fires or momentum dumps. The F-test for harmonic components may be applied to multitaper power spectrum estimates to identify statistically significant oscillations that would not rise above a white noise–based false alarm probability. Multitapering improves the dynamic range of the power spectrum estimate and suppresses spectral window artifacts. We show that the multitaper–F-test combination applied to Kepler observations of KIC 6102338 detects differential rotation without requiring iterative sinusoid fitting and subtraction. Significant signals reside at harmonics of both fundamental rotation frequencies and suggest an antisolar rotation profile. Next we use the missing-data multitaper power spectrum estimator to identify the oscillation modes responsible for the complex "scallop-shell" shape of the K2 light curve of EPIC 203354381. We argue that multitaper power spectrum estimators should be used for all time series with regular observing cadence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Godiva IV Simulated Radiation Field Characterization and Variance Reduction

Godiva IV is a system comprised of highly enriched uranium alloyed with molybdenum in the form of fuel plate rings. The reactor, along with its predecessors, was designed with the unique ability to satisfy interests in the super-prompt-critical reactor operation space. Originally, the reactor was part of the Los Alamos Critical Experiments Facility (LACEF) at Technical Area-18 (TA-18). The radiation field around Godiva at this facility was well characterized and understood. As a fast neutron system, the neutron spectrum in and around Godiva was close to a Watt Fission spectrum. The Kiva where Godiva IV was located at LACEF was made of thin, sheet metal walls which did not contribute significantly to the neutron spectrum. Following the transition of LACEF to the National Critical Experiments and Research Center (NCERC) in Nevada, Godiva-IV was moved from TA-18 to the Device Assembly Facility (DAF) at the Nevada National Security Site (NNSS). Part of this move brought renewed interest in radiation field characterization. The new facility introduced significant changes to the environment surrounding Godiva, and preliminary foil irradiation results suggested that the room contribution to the neutron spectrum was significant. Unlike at TA-18, a large thermal neutron signature was added to the fast spectrum from Godiva due to significant room return. A primary goal due to the additional complexity that the room return adds to the Godiva IV radiation emission spectrum was the development of an efficient Monte Carlo N-Particle (MCNP) calculation capable of characterizing the neutron spectrum anywhere in the room around Godiva. A campaign of activation foil irradiations and analysis were completed to support the validation of the MCNP model. The modeling of these foils in MCNP can be easily done with a standard volumetric neutron flux tally. However, given the multitude of locations and reaction rates to be modeled, further steps must be taken to increase the efficiency of these calculations in MCNP. During this study, a benchmark model currently under development for Godiva IV was used. A qualitative assessment of the thermal neutron contributors was performed using spatial neutron distribution plots. Additional detail was added to the model based on the qualitative results showing the thermal spectrum’s large sensitivity to hydrogenous material. Neutron energy spectra was evaluated at discrete locations in the room around Godiva to quantify the relative contribution of various components. It was discovered that the concrete walls are the largest contributor to the thermal signature, with minor contributions from plastic components surrounding Godiva. Following these results, two different variance reduction techniques were implemented to improve the problem efficiency in these calculations. In the first approach, an F5 point detector tally was implemented in the standard Godiva IV criticality problem. The second approach involved a weight-window generator implementation with an F5 point detector tally in a fixed source problem. The weight window implementation reduced the runtime from 42739.55 minutes to 1803.34 minutes (computer time), compared to the F5 KCODE implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

A kernel machine method for detecting higher order interactions in multimodal datasets: Application to schizophrenia

Technological advances are enabling us to collect multimodal datasets at an increasing depth and resolution while with decreasing labors. Understanding complex interactions among multimodal datasets, however, is challenging. Here in this study, we tested the interaction effect of multimodal datasets using a novel method called the kernel machine for detecting higher order interactions among biologically relevant multimodal data. Using a semiparametric method on a reproducing kernel Hilbert space, we formulated the proposed method as a standard mixed-effects linear model and derived a score-based variance component statistic to test higher order interactions between multimodal datasets. The method was evaluated using extensive numerical simulation and real data from the Mind Clinical Imaging Consortium with both schizophrenia patients and healthy controls. Our method identified 13-triplets that included 6 gene-derived SNPs, 10 ROIs, and 6 gene-specific DNA methylations that are correlated with the changes in hippocampal volume, suggesting that these triplets may be important for explaining schizophrenia-related neurodegeneration. The performance of the proposed method is compared with the following methods: test based on only first and first few principal components followed by multiple regression, and full principal component analysis regression, and the sequence kernel association test. With strong evidence (p-value ≤0.000001), the triplet (MAGI2, CRBLCrus1.L, FBXO28) is a significant biomarker for schizophrenia patients. This novel method can be applicable to the study of other disease processes, where multimodal data analysis is a common task.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1987

Analysis and interpretation of satellite fragmentation data

The velocity perturbations of the fragments of a satellite can shed valuable information regarding the nature and intensity of the fragmentation. A feasibility study on calculating the velocity perturbations from existing equations was carried out by analyzing 23 major documented fragmentation events. It was found that whereas the calculated values of the radial components of the velocity change were often unusually high, those in the two other orthogonal directions were mostly reasonable. Since the uncertainties in the radial component necessarily translate into uncertainties in the total velocity change, it is suggested that alternative expressions for the radial component of velocity be sought for the purpose of determining the cause of the fragmentation from the total velocity change. The calculated variances in the velocity perturbations in the two directions orthogonal to the radial vector indicate that they have the smallest values for collision induced breakups and the largest values for low-intensity explosion induced breakups. The corresponding variances for high-intensity explosion induced breakups generally have values intermediate between those of the two extreme categories. A three-dimensional plot of the variances in the two orthogonal velocity perturbations and the plane change angle shows a clear separation between the three major types of breakups. This information is used to reclassify a number of satellite fragmentation events of unknown category.

Tan, Arjun↗
DOE OSTI2023

A Vertically Resolved Analysis of Radiative Feedbacks on Moist Static Energy Variance in Tropical Cyclones

Abstract A vertically resolved moist static energy (MSE) variance budget framework is used to diagnose processes associated with the development of tropical cyclones (TCs) in a general circulation model (GCM) under realistic boundary conditions. Previous studies have shown that interactions between radiation and MSE promote TC development. Here, we examine the vertical contributions of radiation and its interactions with MSE by performing several mechanism-denial experiments in which synoptic-scale radiative interactions are suppressed either in the boundary layer or in the free troposphere. Partly suppressing radiative interactions results in a reduction in global TC frequency. However, the magnitude of reduction and structure of the feedback depend on the intensity and structure of the TCs in these mechanism-denial experiments, indicating that both the magnitude and the vertical location of radiative interactions can impact global TC frequency. Using instantaneous 6-hourly outputs, an explicit computation reveals distinct spatial patterns of the advection term: the vertical component is positive in the mid- to upper troposphere, which reflects an upward transport of MSE by deep convection, whereas the horizontal component is positive in the boundary layer. These results illustrate the impact of the vertical distribution of radiative interactions and vertically varied contribution of the advection term in the development of TCs.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗