Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed estimation”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Estimating IDP Origins Using ACLED Political Violence Events: Validation with Lebanon IDP Flow Data

A key input in modeling population distribution and flows following conflict events is the inclusion of Internally Displaced Person (IDP) flows between administrative units within a country. These flows are critical for capturing population movement and redistribution driven by current events, particularly conflict. In some cases, IDP destination data are available while origin data is incomplete or unavailable. This creates a gap in understanding where displacement is occurring, limiting the ability to model population redistribution accurately. Without origin data, it is not possible to reallocate population flows or accurately represent where displacement is occurring within the country. This report evaluates whether Armed Conflict Location Event Data (ACLED) political violence event data can be used to estimate IDP origin distributions when direct origin data are unavailable. The approach is validated using historical IDP flow data from Lebanon, where both origins and destinations are observed. Results show that ACLED event distributions strongly correspond to observed IDP-origin patterns, particularly when using cumulative 60-day event windows. The method is most reliable for identifying major origin districts and approximating proportional origin shares. However, it is not intended to reconstruct exact individual displacement flows, but rather to provide a probabilistic spatial allocation of displacement origins.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Quantile-dependent expressivity of serum C-reactive protein concentrations in family sets

“Quantile-dependent expressivity” occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g., C-reactive protein, CRP) is high or low relative to its distribution. We have previously shown that the heritabilities ( h 2 ) of coffee and alcohol consumption, postprandial lipemia, lipoproteins, leptin, adiponectin, adiposity, and pulmonary function are quantile-specific. Whether CRP heritability is quantile-specific is currently unknown. Methods Serum CRP concentrations from 2,036 sibships and 6,144 offspring-parent pairs were analyzed from the Framingham Heart Study. Quantile-specific heritability from full-sib (β FS , h 2 ={(1 + 8r spouse β FS ) 0.5 - 1}/(2r spouse )) and offspring-parent regression slopes (β OP , h 2 = 2β OP /(1 + r spouse )) were estimated robustly by quantile regression with nonparametric significance determined from 1,000 bootstrap samples. Results Quantile-specific h 2 (±SE) increased with increasing percentiles of the offspring’s age- and sex-adjusted CRP distribution when estimated from β OP ( P trend = 0.0004): 0.02 ± 0.01 at the 10th, 0.04 ± 0.01 at the 25th, 0.10 ± 0.02 at the 50th, 0.20 ± 0.05 at the 75th, and 0.33 ± 0.10 at the 90th percentile, and when estimated from β FS ( P trend = 0.0008): 0.03±0.01 at the 10th, 0.06 ± 0.02 at the 25th, 0.14 ± 0.03 at the 50th, 0.24 ± 0.05 at the 75th, and 0.53 ± 0.21 at the 90th percentile. Conclusion Heritability of serum CRP concentration is quantile-specific, which may explain or contribute to the inflated CRP differences between CRP (rs1130864, rs1205, rs1800947, rs2794521, rs3091244), FGB (rs1800787), IL-6 (rs1800795, rs1800796), IL6R (rs8192284), TNF- α (rs1800629) and APOE genotypes following CABG surgery, stroke, TIA, curative esophagectomy, intensive periodontal therapy, or acute exercise; during acute coronary syndrome or Staphylococcus aureus bacteremia; or in patients with chronic rheumatoid arthritis, diabetes, peripheral arterial disease, ankylosing spondylitis, obesity or inflammatory bowel disease or who smoke.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Convex relaxation for Fokker–Planck equation

We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics by solving the high-dimensional Fokker–Planck equation in terms of low-order cluster moments or marginals. With this approach, we bypass the exponential complexity of estimating the full high-dimensional distribution and directly solve the simplified partial differential equations for low-order moments/marginals. Moreover, the proposed moment/marginal relaxation is fully convex and can be solved via off-the-shelf solvers. We further propose a time-dependent version of the convex programs to study non-equilibrium dynamics. In a specific setting, we show the proposed method can recover a mean-field-type equilibrium density. Numerical results are provided to demonstrate the performance of the proposed algorithm for high-dimensional systems.

Chen, Yian↗
NASA NTRS2016

Probabilistic Damage Characterization Using the Computationally-Efficient Bayesian Approach

This work presents a computationally-ecient approach for damage determination that quanti es uncertainty in the provided diagnosis. Given strain sensor data that are polluted with measurement errors, Bayesian inference is used to estimate the location, size, and orientation of damage. This approach uses Bayes' Theorem to combine any prior knowledge an analyst may have about the nature of the damage with information provided implicitly by the strain sensor data to form a posterior probability distribution over possible damage states. The unknown damage parameters are then estimated based on samples drawn numerically from this distribution using a Markov Chain Monte Carlo (MCMC) sampling algorithm. Several modi cations are made to the traditional Bayesian inference approach to provide signi cant computational speedup. First, an ecient surrogate model is constructed using sparse grid interpolation to replace a costly nite element model that must otherwise be evaluated for each sample drawn with MCMC. Next, the standard Bayesian posterior distribution is modi ed using a weighted likelihood formulation, which is shown to improve the convergence of the sampling process. Finally, a robust MCMC algorithm, Delayed Rejection Adaptive Metropolis (DRAM), is adopted to sample the probability distribution more eciently. Numerical examples demonstrate that the proposed framework e ectively provides damage estimates with uncertainty quanti cation and can yield orders of magnitude speedup over standard Bayesian approaches.

Warner, James E.↗
NASA NTRS2004

Analysis of Upper Air, Ground and Remote Sensing Data For the ATLAS Field Campaign in San Juan, Puerto Rico

The Atlas San Juan Mission was conducted in February 2004 with the main objectives of observing the Urban Heat Island of San Juan, providing high resolution data of the land use for El Yunque Rain Forest and for calibrating remote sensors. The mission was coordinated with NASA staff members at Marshall, Stennis, Goddard, and Glenn. The Airborne Thermal and Land Applications Sensor (ATLAS) from NASA/Stennis, that operates in the visual and IR bands, was used as the main sensor and was flown over Puerto Rico in a Lear 23 jet plane. To support the data gathering effort by the ATLAS sensor, remote sensing observations and upper air soundings were conducted along with the deployment of a number of ground based weather stations and temperature sensors. This presentation focuses in the analysis of this complementary data for the Atlas San Juan Mission. Upper air data show that during the days of the mission the Caribbean mid and high atmospheres were relatively dry and highly stable reflecting positive surface lifted index, a necessary condition to conduct this suborbital campaign. Surface wind patterns at levels below 850mb were dominated by the easterly trades, while the jet stream at the edge of the troposphere dominated the westerly wind at levels above 500mb. The jet stream remained at high latitudes reducing the possibility of fronts. In consequence, only 8.4 mm of precipitation were reported during the entire mission. Observation of soundings located about 150 km apart reflected minimum variations of the boundary layer across the island for levels below 850 meters and a uniform atmosphere for higher levels. The weather stations and the temperature sensors were placed at strategic locations to observe variations across the urban and rural landscapes. Time series plot of the stations' data show that heavily urbanized commercial areas have higher air temperatures than urban and suburban residential areas, and much higher temperatures than rural areas. Temperature differences [dT(U-R)] were obtained by subtracting the values of several stations from a reference urban station, located in the commercial area of San Juan. These time series show that the UHI peaks during the morning between 10:00am and noon to an average of 4.5 C, a temporal pattern not previously observed in similar studies for continental cities. It is also observed a high variability of the UHI with the precipitation patterns even for short events. These results may be a reflection of a large land use density by low level buildings with an apparent absence of significant heat storage effects in the urban areas, and the importance of the surrounding soil and vegetation moisture in controlling the urban tropical climate. The ATLAS data was used to determine albedo and surface temperature patterns on a 10m scale for the study area. These data were used to calibrate the spatial distribution of the surface temperature when using remote sensing images from MODIS (Moderate Resolution Imaging Spectroradiometer). Surface temperatures were estimated using the land surface temperature product MOD11_L2 distributed by the Land Process Distributed Active Archive Center (LP DAAC). These results show the maximum, minimum and average temperatures in San Juan and in the entire Island at a resolution of 1 km. The information retrieved from MODIS for land surface temperatures reflected similar temporal and spatial variations as the weather stations and ATLAS measurements with a highest absolute offset of about 5 C due to the differences between surface and air temperatures.

Gonzalez, J. E.↗
NASA NTRS2006

Studies of the Ability to Hold the Eye in Eccentric Gaze: Measurements in Normal Subjects with the Head Erect

We studied the ability to hold the eyes in eccentric horizontal or vertical gaze angles in 68 normal humans, age range 19-56. Subjects attempted to sustain visual fixation of a briefly flashed target located 30 in the horizontal plane and 15 in the vertical plane in a dark environment. Conventionally, the ability to hold eccentric gaze is estimated by fitting centripetal eye drifts by exponential curves and calculating the time constant (t(sub c)) of these slow phases of gazeevoked nystagmus. Although the distribution of time-constant measurements (t(sub c)) in our normal subjects was extremely skewed due to occasional test runs that exhibited near-perfect stability (large t(sub c) values), we found that log10(tc) was approximately normally distributed within classes of target direction. Therefore, statistical estimation and inference on the effect of target direction was performed on values of z identical with log10t(sub c). Subjects showed considerable variation in their eyedrift performance over repeated trials; nonetheless, statistically significant differences emerged: values of tc were significantly higher for gaze elicited to targets in the horizontal plane than for the vertical plane (P less than 10(exp -5), suggesting eccentric gazeholding is more stable in the horizontal than in the vertical plane. Furthermore, centrifugal eye drifts were observed in 13.3, 16.0 and 55.6% of cases for horizontal, upgaze and downgaze tests, respectively. Fifth percentile values of the time constant were estimated to be 10.2 sec, 3.3 sec and 3.8 sec for horizontal, upward and downward gaze, respectively. The difference between horizontal and vertical gazeholding may be ascribed to separate components of the velocity position neural integrator for eye movements, and to differences in orbital mechanics. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied in a clinical setting to patients who were exposed to environments that may have modified their central integrators and thus require monitoring. Patients with gaze-evoked nystagmus can be flagged by comparing to the above established normative criteria.

Reschke, Millard F.↗
DOE OSTI2023

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Machine Learning-Based Method to Estimate Transformer Primary-Side Voltages with Limited Customer-Side AMI Measurements

Distribution control applications such as volt/var optimization, network reconfiguration, and distribution automation require accurate knowledge of the distribution system state. The lack of sufficient sensors on the primary side of distribution networks often limits the accuracy of the control decisions by these applications. The deployment of advanced metering infrastructure (AMI) provides utilities an opportunity to translate the AMI data on the secondary onto the primary so that it can be used as pseudo-measurements to augment the limited existing measurements on the primary. This paper develops a machine learning based approach for estimating service transformer primary-side voltages by using limited secondary-side AMI measurement. The machine learning model is developed by using random forest algorithm. The estimated primary-side voltages can be used by utilities as pseudo-measurements for distribution control applications. The detailed secondary model topology, which is an essential input data for many existing algorithms, is not required for the proposed method. The performance of the proposed method is validated by using AMI measurements from the field and an actual distribution feeder model of San Diego Gas & Electric Company.

advanced metering infrastructure↗
DOE OSTI2022

Precise ROCOF estimation algorithm for low inertia power grids

The precise estimation of Rate Of Change Of Frequencies (ROCOFs) in a generation trip event can be helpful on power system inertia estimation, fast system response, and accurate event capturing. However, the ROCOF estimations from the existing Synchronized Measurement Devices (SMDs) are usually simply calculated by the finite difference between two adjacent frequency measurement points. The ineluctable noises, disturbances, and spikes from real-world frequency measurements can bring large dynamics to the ROCOF estimation and thus can result in an inaccurate estimation of the initial ROCOF. This issue becomes more serious when the target power grid has low inertia, and a large amount of distributed energy sources are deployed. To address this issue, a precise ROCOF estimation algorithm is designed based on least square method with flexible window size. In addition, a median filter is also designed and applied on the frequency measurements before using the proposed algorithm. The window size and thresholds in the proposed algorithm are determined with historical event data analysis. The proposed algorithm is deployed in a low cost, flexible, and distribution level universal grid analyzer (UGA) platform. In conclusion, multiple experiments are conducted in both a laboratory and the Hawaiian Islands to verify the effectiveness of the proposed algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2011

Evaluating the Performance of Single and Double Moment Microphysics Schemes During a Synoptic-Scale Snowfall Event

Increases in computing resources have allowed for the utilization of high-resolution weather forecast models capable of resolving cloud microphysical and precipitation processes among varying numbers of hydrometeor categories. Several microphysics schemes are currently available within the Weather Research and Forecasting (WRF) model, ranging from single-moment predictions of precipitation content to double-moment predictions that include a prediction of particle number concentrations. Each scheme incorporates several assumptions related to the size distribution, shape, and fall speed relationships of ice crystals in order to simulate cold-cloud processes and resulting precipitation. Field campaign data offer a means of evaluating the assumptions present within each scheme. The Canadian CloudSat/CALIPSO Validation Project (C3VP) represented collaboration among the CloudSat, CALIPSO, and NASA Global Precipitation Measurement mission communities, to observe cold season precipitation processes relevant to forecast model evaluation and the eventual development of satellite retrievals of cloud properties and precipitation rates. During the C3VP campaign, widespread snowfall occurred on 22 January 2007, sampled by aircraft and surface instrumentation that provided particle size distributions, ice water content, and fall speed estimations along with traditional surface measurements of temperature and precipitation. In this study, four single-moment and two double-moment microphysics schemes were utilized to generate hypothetical WRF forecasts of the event, with C3VP data used in evaluation of their varying assumptions. Schemes that incorporate flexibility in size distribution parameters and density assumptions are shown to be preferable to fixed constants, and that a double-moment representation of the snow category may be beneficial when representing the effects of aggregation. These results may guide forecast centers in optimal configurations of their forecast models for winter weather and identify best practices present within these various schemes.

Molthan, Andrew L.↗
DOE OSTI2023

A Multi-Fidelity Gaussian Process Regression Method for Probabilistic Wind Farm Power Curve Estimation

Accurate estimation of the power curve for wind turbines or wind farms is crucial to ensure their efficient operation and management. However, conventional methods for power curve estimation rely either on expensive and infrequent measurements or on low-quality numerical simulations. Moreover, the majority of previous studies on power curve estimation for wind turbines or wind farms focused on deterministic estimation, which provides a point estimate of the relationship between wind speed and power generation. Nevertheless, the deterministic approach fails to consider the inherent uncertainty associated with wind energy production resulting from varying turbine characteristics. This can lead to inaccurate power generation estimation and suboptimal decisions regarding energy management. In this paper, a kernel density estimation (KDE) based Multi-Fidelity Gaussian Process Regression (MFGPR) model is proposed to fuse theoretical power curve data and the ground true measurements to create a mapping of wind speed and wind power. By conducting a case study on an actual wind farm in China, the efficacy of the proposed MFGPR model was demonstrated in characterizing the variability of wind power. The probabilistic MFGPR model was also able to generate confidence intervals that encompassed the measured power, thereby improving the accuracy and confidence in wind power estimation or wind resource assessment. Overall, the proposed MFGPR model offers a reliable approach to integrate high-fidelity ground measurements and theoretical power curve data, resulting in precise wind resource assessment and power estimation.

Gaussian process regression↗
NASA NTRS2003

Satellite Detection of Precipitation over the North Pacific

Intense extratropical winter cyclones often impact the West Coast of North America with strong winds and heavy precipitation. Several times during a winter season, short-term forecasts (24 - 48 hours) of these storms are seriously deficient with central pressure errors in the 10's of hPa and surface low position errors in the 100's of km. For example, 48-hr sea level pressure errors (forecast - observation) at buoy 46005 off the Oregon coast for the 2001 - 2002 winter season is plotted. In addition, two times the standard deviation (determined from pressure errors from the last four winter seasons) are also shown. It is evident from this figure that large forecast errors (i.e. greater than 10 hPa) occurred about 10 times this past winter at buoy 46005 with three events where the errors were 20 hPa. Beside large forecast errors of sea level pressure, numerical forecasts of precipitation for land falling cyclones can also be flawed. This is due in large part to the lack of accurate precipitation information over the ocean. Therefore, remote sensing techniques are the only viable option for obtaining accurate information on the distribution and intensity of precipitation over the North Pacific. Due to the radiative characteristics of precipitation sized hydrometeors at microwave frequencies, microwave sensors are able to detect precipitation over oceanic regions. Past studies have demonstrated the utility of passive microwave rainrate data for locating intense rainfall in rapidly deepening cyclones, in detecting developing polar mesocyclones and in determining frontal bands. There are currently many sources of microwave rainrate data: the Special Sensor Microwave Imager (SSM/I) (currently flying on three platforms), the Advanced Microwave Sounding Unit (AMSU-B) (currently flying on NOAA-15, NOAA-16, and NOAA-17), and the Tropical Rainfall Measuring Mission Microwave Imager (TMI). Data will soon be available from the Advanced Microwave Radiometer-EOS (AMSR-E) on the Aqua platform. In this paper, we present a new technique for mapping rainrate distributions over the North Pacific utilizing rainrate estimates from several microwave sensors and upper-tropospheric winds derived from geosynchronous satellite IR data. The goal of this work is to develop a way to obtain high temporal and spatial rainfall information over the North Pacific. This information will be used to support the verification of model derived precipitation distributions and to support the analysis of in situ measurements of rainfall during the Improvement of Microphysical Parameterization through Observational Verification Experiment (IMPROVE) field campaigns.

Smith, Jeremy↗
DOE OSTI2024

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗
NASA NTRS1987

Assessment Of Advanced Concentrator Photovoltaic Modules

Computer program developed to estimate probabilistically total and component costs to end users of concentrating photovoltaic arrays in 1990's. Based on modified version of SIMRAND program, Simulation of Research and Development projects, -to estimate total cost of system for large number of module designs. Module design that yields minimum cost for total system then chosen as preferred design. To obtain distribution of expected module costs, simulation performed. Component estimates made by SIMRAND include efficiency of module and costs of cell, assembly of cells, lenses, and modules. Written in FORTRAN 77.

Borden, C. S.↗
DOE OSTI2022

Twist-three cross-sections in deeply virtual Compton scattering

We study the deeply virtual Compton scattering process with both twist-two and twist-three Compton form factors and present our cross-sections formulas with all polarization configurations. While the twist-three contributions are generally assumed to be negligible in the literature due to the kinematical suppression, we compare them with the twist-two ones at typical JLab 6 GeV and 12 GeV kinematics as well as EIC kinematics and show their kinematical suppression explicitly, justifying the leading-twist approximation made in the literature. In addition, we also estimate the twist-three Compton form factors using Wandzura-Wilczek relations and inputs of twist-two generalized parton distributions based on a reggeized spectator model. With those estimated Compton form factors, we analyze the kinematical behavior of twist-two and twist-three cross-sections in a wide range of kinematics, and discuss the optimal regions for separating the leading-twist effects from the higher-twist ones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

A new algorithm for microwave delay estimation from water vapor radiometer data

A new algorithm has been developed for the estimation of tropospheric microwave path delays from water vapor radiometer (WVR) data, which does not require site and weather dependent empirical parameters to produce high accuracy. Instead of taking the conventional linear approach, the new algorithm first uses the observables with an emission model to determine an approximate form of the vertical water vapor distribution which is then explicitly integrated to estimate wet path delays, in a second step. The intrinsic accuracy of this algorithm has been examined for two channel WVR data using path delays and stimulated observables computed from archived radiosonde data. It is found that annual RMS errors for a wide range of sites are in the range from 1.3 mm to 2.3 mm, in the absence of clouds. This is comparable to the best overall accuracy obtainable from conventional linear algorithms, which must be tailored to site and weather conditions using large radiosonde data bases. The new algorithm's accuracy and flexibility are indications that it may be a good candidate for almost all WVR data interpretation.

Robinson, S. E.↗
DOE OSTI2020

Data from: Comparing traditional and Bayesian approaches to ecological meta-analysis

Despite the wide application of meta-analysis in ecology, some of the traditional methods used for meta-analysis may not perform well given the type of data characteristic of ecological meta-analyses. We reviewed published meta-analyses on the ecological impacts of global climate change, evaluating the number of replicates used in the primary studies (ni) and the number of studies or records (k) that were aggregated to calculate a mean effect size. We used the results of the review in a simulation experiment to assess the performance of conventional frequentist and Bayesian meta-analysis methods for estimating a mean effect size and its uncertainty interval. Our literature review showed that ni and k were highly variable, distributions were right-skewed, and were generally small (median ni =5, median k=44). Our simulations show that the choice of method for calculating uncertainty intervals was critical for obtaining appropriate coverage (close to the nominal value of 0.95). When k was low (<40), 95% coverage was achieved by a confidence interval based on the t-distribution that uses an adjusted standard error (the Hartung-Knapp-Sidik-Jonkman, HKSJ), or by a Bayesian credible interval, whereas bootstrap or z-distribution confidence intervals had lower coverage. Despite the importance of the method to calculate the uncertainty interval, 39% of the meta-analyses reviewed did not report the method used, and of the 61% that did, 94% used a potentially problematic method, which may be a consequence of software defaults. In general, for a simple random-effects meta-analysis, the performance of the best frequentist and Bayesian methods were similar for the same combinations of factors (k and mean replication), though the Bayesian approaches had higher than nominal (>95%) coverage for the mean effect when k was very low (k<15). Our literature review suggests that many meta-analyses that used z-distribution or bootstrapping confidence intervals may have over-estimated the statistical significance of their results when the number of studies was low; more appropriate methods need to be adopted in ecological meta-analyses.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

A new technique for predicting geosynchronous satellite collision probability

A new technique has been developed to predict the probability of an expired geosynchronous satellite colliding with an active satellite. This new technique employs deterministic methods for modeling the motion of satellites and applies statistical techniques to estimate the collision probability. The collision probability is used to estimate the expected time between collisions based on realistic distributions of expired and active satellites. The primary advantage of this new technique is that realistic distributions can be used in the prediction process instead of uniform distributions as has been used in previous techniques. The expected time between collisions based on a current NORAD database is estimated to be in the hundreds of years.

Mccormick, B.↗