Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A novel approach for estimating densities of secretive species from road-survey and spatial-movement data

Accurate estimates of population density are a critical component of effective wildlife conservation and management. However, many snake species are so secretive that their density cannot be determined using traditional methods such as capture–mark–recapture. Thus, the status of most terrestrial snake populations remains completely unknown. We developed a novel simulation-based technique for estimating density of secretive snakes that combined behavioural observations of snake road-crossing behaviour (crossing speed), effort-corrected road-survey data, and simulations of spatial movement patterns derived from radio-telemetry, without relying on mark–recapture. Here, we used radio-telemetry data to parameterise individual-based movement models that estimate the frequency with which individual snakes cross roads and used information on survey vehicle speed and snake crossing speed to determine the probability of detecting a snake, given that it crosses the road transect during a survey. Snake encounter frequencies during systematic road surveys were then interpreted in light of detection probabilities and simulation model results to estimate snake densities and to assess various factors likely to affect abundance estimates. We demonstrated the broad applicability of this approach through a case study of the imperiled southern hognose snake (Heterodon simus) in the North Carolina (USA) Sandhills. We estimated that H. simus occurs at average densities of 0.17 ha –1 in the North Carolina Sandhills and explored the sensitivity of this estimate to assumptions and variation in model parameters. Our novel method allowed us to generate the first abundance estimates for H. simus. We found that H. simus exists at low densities relative to congeners and other mid-sized snake species, raising concern that this species may not only have declined in geographic range, but may also occur at low densities or be declining in their strongholds, such as the North Carolina Sandhills. We present a framework for estimating density of species that have traditionally been considered too secretive to study at the population level. This method will greatly enhance our ability to study and manage a wide variety of snake species and could be applied to other secretive wildlife species that are most frequently encountered during road surveys.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Estimated dose rates to members of the public from external exposure to patients with {sup 131}I thyroid treatment

Purpose: Estimated dose rates that may result from exposure to patients who had been administered iodine-131 ({sup 131}I) as part of medical therapy were calculated. These effective dose rate estimates were compared with simplified assumptions under United States Nuclear Regulatory Commission Regulatory Guide 8.39, which does not consider body tissue attenuation nor time-dependent redistribution and excretion of the administered {sup 131}I. Methods: Dose rates were estimated for members of the public potentially exposed to external irradiation from patients recently treated with {sup 131}I. Tissue attenuation and iodine biokinetics were considered in the patient in a larger comprehensive effort to improve external dose rate estimates. The external dose rate estimates are based on Monte Carlo simulations using the Phantom with Movable Arms and Legs (PIMAL), previously developed by Oak Ridge National Laboratory and the United States Nuclear Regulatory Commission. PIMAL was employed to model the relative positions of the {sup 131}I patient and members of the public in three exposure scenarios: (1) traveling on a bus in a total of six seated or standing permutations, (2) two nursing home cases where a caregiver is seated at 30 cm from the patient’s bedside and a nursing home resident seated 250 cm away from the patient in an adjacent bed, and (3) two hotel cases where the patient and a guest are in adjacent rooms with beds on opposite sides of the common wall, with the patient and guest both in bed and either seated back-to-back or lying head to head. The biokinetic model predictions of the retention and distribution of {sup 131}I in the patient assumed a single voiding of urinary bladder contents that occurred during the trip at 2, 4, or 8 h after {sup 131}I administration for the public transportation cases, continuous first-order voiding for the nursing home cases, and regular periodic voiding at 4, 8, or 12 h after administration for the hotel room cases. Organ specific activities of {sup 131}I in the thyroid, bladder, and combined remaining tissues were calculated as a function of time after administration. Exposures to members of the public were considered for {sup 131}I patients with normal thyroid uptake (peak thyroid uptake of ∼27% of administered {sup 131}I), differentiated thyroid cancer (DTC, 5% uptake), and hyperthyroidism (80% uptake). Results: The scenario with the patient seated behind the member of the public yielded the highest dose rate estimate of seated public transportation exposure cases. The dose rate to the adjacent room guest was highest for the exposure scenario in which the hotel guest and patient are seated by a factor of ∼4 for the normal and differentiated thyroid cancer uptake cases and by a factor of ∼3 for the hyperthyroid case. Conclusions: It was determined that for all modeled cases, the DTC case yielded the lowest external dose rates, whereas the hyperthyroid case yielded the highest dose rates. In estimating external dose to members of the public from patients with {sup 131}I therapy, consideration must be given to (patient- and case-specific) administered {sup 131}I activities and duration of exposure for a more complete estimate. The method implemented here included a detailed calculation model, which provides a means to determine dose rate estimates for a range of scenarios. The method was demonstrated for variations of three scenarios, showing how dose rates are expected to vary with uptake, voiding pattern, and patient location.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Prospects for in vivo estimation of photon linear attenuation coefficients using postprocessing dual-energy CT imaging on a commercial scanner: Comparison of analytic and polyenergetic statistical reconstruction algorithms

Purpose: Accurate patient-specific photon cross-section information is needed to support more accurate model-based dose calculation for low energy photon-emitting modalities in medicine such as brachytherapy and kilovoltage x-ray imaging procedures. A postprocessing dual-energy CT (pDECT) technique for noninvasivein vivo estimation of photon linear attenuation coefficients has been experimentally implemented on a commercial CT scanner and its accuracy assessed in idealized phantom geometries. Methods: Eight test materials of known composition and density were used to compare pDECT-estimated linear attenuation coefficients to NIST reference values over an energy range from 10 keV to 1 MeV. As statistical image reconstruction (SIR) has been shown to reconstruct images with less random and systematic error than conventional filtered backprojection (FBP), the pDECT technique was implemented with both an in-house polyenergetic SIR algorithm, alternating minimization (AM), as well as a conventional FBP reconstruction algorithm. Improvement from increased spectral separation was also investigated by filtering the high-energy beam with an additional 0.5 mm of tin. The law of propagated uncertainty was employed to assess the sensitivity of the pDECT process to errors in reconstructed images. Results: Mean pDECT-estimated linear attenuation coefficients for the eight test materials agreed within 1% of NIST reference values for energies from 1 MeV down to 30 keV, with mean errors rising to between 3% and 6% at 10 keV, indicating that the method is unbiased when measurement and calibration phantom geometries are matched. Reconstruction with FBP and AM algorithms conferred similar mean pDECT accuracy. However, single-voxel pDECT estimates reconstructed on a 1 × 1 × 3 mm{sup 3} grid are shown to be highly sensitive to reconstructed image uncertainty; in some cases pDECT attenuation coefficient estimates exhibited standard deviations on the order of 20% around the mean. Reconstruction with the statistical AM algorithm led to standard deviations roughly 40% to 60% less than FBP reconstruction. Additional tin filtration of the high energy beam exhibits similar pDECT estimation accuracy as the unfiltered beam, even when scanning with only 25% of the dose. Using the law of propagated uncertainty, low Z materials are found to be more sensitive to image reconstruction errors than high Z materials. Furthermore, it is estimated that reconstructed CT image uncertainty must be limited to less than 0.25% to achieve a target linear-attenuation coefficient estimation uncertainty of 3% at 28 keV. Conclusions: That pDECT supports mean linear attenuation coefficient measurement accuracies of 1% of reference values for energies greater than 30 keV is encouraging. However, the sensitivity of the pDECT measurements to noise and systematic errors in reconstructed CT images warrants further investigation in more complex phantom geometries. The investigated statistical reconstruction algorithm, AM, reduced random measurement uncertainty relative to FBP owing to improved noise performance. These early results also support efforts to increase DE spectral separation, which can further reduce the pDECT sensitivity to measurement uncertainty.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1995

Global precipitation estimates based on a technique for combining satellite-based estimates, rain gauge analysis, and NWP model precipitation information

The 'satellite-gauge model' (SGM) technique is described for combining precipitation estimates from microwave satellite data, infrared satellite data, rain gauge analyses, and numerical weather prediction models into improved estimates of global precipitation. Throughout, monthly estimates on a 2.5 degrees x 2.5 degrees lat-long grid are employed. First, a multisatellite product is developed using a combination of low-orbit microwave and geosynchronous-orbit infrared data in the latitude range 40 degrees N - 40 degrees S (the adjusted geosynchronous precipitation index) and low-orbit microwave data alone at higher latitudes. Then the rain gauge analysis is brougth in, weighting each field by its inverse relative error variance to produce a nearly global, observationally based precipitation estimate. To produce a complete global estimate, the numerical model results are used to fill data voids in the combined satellite-gauge estimate. Our sequential approach to combining estimates allows a user to select the multisatellite estimate, the satellite-gauge estimate, or the full SGM estimate (observationally based estimates plus the model information). The primary limitation in the method is imperfections in the estimation of relative error for the individual fields. The SGM results for one year of data (July 1987 to June 1988) show important differences from the individual estimates, including model estimates as well as climatological estimates. In general, the SGM results are drier in the subtropics than the model and climatological results, reflecting the relatively dry microwave estimates that dominate the SGM in oceanic regions.

Huffman, George J.↗
DOE OSTI2019

Estimation of inflow uncertainties in laminar hypersonic double-cone experiments

We propose herein a probabilistic framework for assessing the consistency of an experimental dataset, i.e., whether the stated experimental conditions are consistent with the measurements provided. In case the dataset is inconsistent, our framework allows one to hypothesize and test sources of inconsistencies. This is crucial in model validation efforts. The framework relies on Bayesian inference to estimate experimental settings deemed uncertain, from measurements deemed accurate. The quality of the inferred variables is gauged by its ability to reproduce held-out experimental measurements. We test the correctness of the framework on three double-cone experiments conducted in the CUBRC Inc.'s LENS-I shock tunnel, which have also been numerically simulated successfully. Thereafter, we use the framework to investigate two double-cone experiments (executed in the LENS-XX shock tunnel) which have encountered difficulties when used in model validation exercises. We detect an inconsistency with one of the LENS-XX experiments. In addition, we hypothesize two causes for our inability to simulate LEXS-XX experiments accurately and test them using our framework. We find that there is no single cause that explains all the discrepancies between model predictions and experimental data, but different causes explain different discrepancies, to larger or smaller extent. We end by proposing that uncertainty quantification methods be used more widely to understand experiments and characterize facilities, and we cite three different methods to do so, the third of which we present in this paper.

Thermophysics↗
NASA NTRS1986

Crop acreage estimation using a Landsat-based estimator as an auxiliary variable

The problem of improving upon the ground survey estimates of crop acreages by utilizing Landsat data is addressed. Three estimators, called regression, ratio, and stratified ratio, are studied for bias and variance, and their relative efficiencies are compared. The approach is to formulate analytically the estimation problem that utilizes ground survey data, as collected by the U.S. Department of Agriculture, and Landsat data, which provide complete coverage for an area of interest, and then to conduct simulation studies. It is shown over a wide range of parametric conditions that the regression estimator is the most efficient unless there is a low correlation between the actual and estimated crop acreages in the sampled area segments, in which case the ratio and stratified ratio estimators are better. Furthermore, it is seen that the regression estimator is potentially biased due to estimating the regression coefficient from the training sample segments. Estimation of the variance of the regression estimator is also investigated. Two variance estimators are considered, the large sample variance estimator and an alternative estimator suggested by Cochran. The large sample estimate of variance is found to be biased and inferior to the Cochran estimate for small sample sizes.

Chhikara, R. S.↗
NASA NTRS2003

Estimating Uncertainty in GPCP and TRMM Multi-Satellite Precipitation Estimates

One of the high-priority problems in satellite precipitation estimation is developing algorithms for estimating the errors in precipitation retrievals by individual sensors and subsequent multi-satellite combinations. Classically, we distinguish between "random" and "bias" errors, which do and do not, respectively average to zero over a "big enough" time/space sample. The current operational GPCP and TRMM multi-satellite algorithms are nearly unique in estimating random error for the monthly precipitation estimates from individual sensor systems (including gauge), following Huffman, and then making multi-sensor combinations. No routinely operational global precipitation produces estimates of bias error. Subsequently, a similar scheme has been followed to provide random error estimates for the Multi-satellite Precipitation Analysis (MPA) being computed in real time and after real time for TRMM. The Huffman algorithm for random error is briefly reviewed, including a discussion of the limitations imposed on the algorithm by standard monthly precipitation data sets. Starting from a very simple theoretical treatment of the histogram of precipitation samples in a month, an equation is developed that depends on the estimated average precipitation rate for the month, the number of samples in the month, and two constants. The constants are set separately for each source of precipitation estimate (such as "raingauge") by calibration at selected ground sites. We discuss recent work validating the random error estimates to highlight the successes and limitations of this first-generation approach. We then consider what information is needed from the individual sensor algorithms to facilitate additional accuracy in the estimation of random errors across the time/space span of climate regimes which a global estimation system must handle. In addition, the thorny issue of estimating bias is raised. Finally, the role of error estimates (and the qualitative errors!) in creating combinations of precipitation estimates from different individual sensors is discussed. This issue is particularly important when fine scales in space and time are being considered, say the 0.25 x 0.25-deg 3-hourly estimates in the MPA.

Huffman, G. J.↗
DOE OSTI2018

Effects of exposure estimation errors on estimated exposure-response relations for PM2.5

Highlights: • Recent literature has identified associations between concentrations of fine particulate matter (PM2.5) in ambient air and adverse health effects at estimated exposure concentrations below current air ambient quality standards. • The models and analyses used do not usually quantify uncertainty about exposure estimates. • Realistic errors in exposure estimates can be estimated from available air quality monitoring stations by predicting concentrations at each station from concentrations at neighboring stations. • Exposure estimation errors are large enough to obscure low-dose nonlinearities or thresholds in the true concentration-response (C-R) function, artificially making them appear to be linear at low doses. • Better estimates of C-R functions will require modeling exposure distributions and estimation errors, rather than just estimating average exposure concentrations. Associations between fine particulate matter (PM2.5) exposure concentrations and a wide variety of undesirable outcomes, from autism and auto theft to elderly mortality, suicide, and violent crime, have been widely reported. Influential articles have argued that reducing National Ambient Air Quality Standards for PM2.5 is desirable to reduce these outcomes. Yet, other studies have found that reducing black smoke and other particulate matter by as much as 70% and dozens of micrograms per cubic meter has not detectably affected all-cause mortality rates even after decades, despite strong, statistically significant positive exposure concentration-response (C-R) associations between them. This paper examines whether this disconnect between association and causation might be explained in part by ignored estimation errors in estimated exposure concentrations. We use EPA air quality monitor data from the Los Angeles area of California to examine the shapes of estimated C-R functions for PM2.5 when the true C-R functions are assumed to be step functions with well-defined response thresholds. The estimated C-R functions mistakenly show risk as smoothly increasing with concentrations even well below the response thresholds, thus incorrectly predicting substantial risk reductions from reductions in concentrations that do not affect health risks. We conclude that ignored estimation errors obscure the shapes of true C-R functions, including possible thresholds, possibly leading to unrealistic predictions of the changes in risk caused by changing exposures. Instead of estimating improvements in public health per unit reduction (e.g., per 10 µg/m{sup 3} decrease) in average PM2.5 concentrations, it may be essential to consider how interventions change the distributions of exposure concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Hybrid approach for estimating dynamic states of synchronous generators

Accurate and robust estimation of dynamic states is essential for monitoring and controlling a power grid. To achieve higher accuracy and robustness, this paper proposes a hybrid estimation approach to estimate the dynamic states of synchronous generators using multiple models, classified into a set of local models and a wide-area model. The proposed estimation approach, first, estimates dynamic states using the extended Kalman filter (EKF) based on a local model and a wide-area model. Then, at each time step, probability indexes, which quantify the likelihood of both local and wide-area models, are determined using hypothesis testing based on measurements innovations. Finally, using the weighted fix approach, the states estimated from the local and wide-area models are combined based on their probability indexes. Simulation studies using the two-area four-machine system shows that the proposed approach can improve estimation accuracy and increase robustness against model errors.

Akhlaghi, Shahrokh↗
DOE OSTI2012

Robust Optical Richness Estimation with Reduced Scatter

Reducing the scatter between cluster mass and optical richness is a key goal for cluster cosmology from photometric catalogs. We consider various modifications to the red-sequence matched filter richness estimator of Rozo et al. (2009b), and evaluate their impact on the scatter in X-ray luminosity at fixed richness. Most significantly, we find that deeper luminosity cuts can reduce the recovered scatter, finding that {sigma}{sub ln L{sub X}|{lambda}} = 0.63 {+-} 0.02 for clusters with M{sub 500c} {approx}> 1.6 x 10{sup 14} h{sub 70}{sup -1} M{sub {circle_dot}}. The corresponding scatter in mass at fixed richness is {sigma}{sub ln M|{lambda}} {approx} 0.2-0.3 depending on the richness, comparable to that for total X-ray luminosity. We find that including blue galaxies in the richness estimate increases the scatter, as does weighting galaxies by their optical luminosity. We further demonstrate that our richness estimator is very robust. Specifically, the filter employed when estimating richness can be calibrated directly from the data, without requiring a-priori calibrations of the red-sequence. We also demonstrate that the recovered richness is robust to up to 50% uncertainties in the galaxy background, as well as to the choice of photometric filter employed, so long as the filters span the 4000 {angstrom} break of red-sequence galaxies. Consequently, our richness estimator can be used to compare richness estimates of different clusters, even if they do not share the same photometric data. Appendix A includes 'easy-bake' instructions for implementing our optimal richness estimator, and we are releasing an implementation of the code that works with SDSS data, as well as an augmented maxBCG catalog with the {lambda} richness measured for each cluster.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
DOE OSTI2020

Comparison of equilibrium climate sensitivity estimates from slab ocean, 150-year, and longer simulations

We compare equilibrium climate sensitivity (ECS) estimates from pairs of long (= 800-year) control and abruptly quadrupled CO2 simulations with shorter (150, 300 year) coupled atmosphere-ocean simulations and Slab Ocean Models (SOM). Consistent with previous work, ECS estimates from shorter coupled simulations based on annual averages for years 1-150 underestimate those from SOM (-8% ± 13%) and long (-14% ± 8%) simulations. Analysis of only years 21-150 improved agreement with SOM (-2% ± 14%) and long (-8% ± 10%) estimates. Use of pentadal averages for years 51-150 results in improved agreement with long simulations (-4% ± 11%). While ECS estimates from current generation US models based on SOM and coupled annual averages of years 1-150 range from 2.6°C to 5.3°C, estimates based longer simulations of the same models range from 3.2°C to 7.0°C. Such variations between methods argues for caution in comparison and interpretation of ECS estimates across models.

Dunne, John P.↗
DOE OSTI2014

Using LiF:Mg,Cu,P TLDs to estimate the absorbed dose to water in liquid water around an {sup 192}Ir brachytherapy source

Purpose: The absorbed dose to water is the fundamental reference quantity for brachytherapy treatment planning systems and thermoluminescence dosimeters (TLDs) have been recognized as the most validated detectors for measurement of such a dosimetric descriptor. The detector response in a wide energy spectrum as that of an{sup 192}Ir brachytherapy source as well as the specific measurement medium which surrounds the TLD need to be accounted for when estimating the absorbed dose. This paper develops a methodology based on highly sensitive LiF:Mg,Cu,P TLDs to directly estimate the absorbed dose to water in liquid water around a high dose rate {sup 192}Ir brachytherapy source. Methods: Different experimental designs in liquid water and air were constructed to study the response of LiF:Mg,Cu,P TLDs when irradiated in several standard photon beams of the LNE-LNHB (French national metrology laboratory for ionizing radiation). Measurement strategies and Monte Carlo techniques were developed to calibrate the LiF:Mg,Cu,P detectors in the energy interval characteristic of that found when TLDs are immersed in water around an{sup 192}Ir source. Finally, an experimental system was designed to irradiate TLDs at different angles between 1 and 11 cm away from an {sup 192}Ir source in liquid water. Monte Carlo simulations were performed to correct measured results to provide estimates of the absorbed dose to water in water around the {sup 192}Ir source. Results: The dose response dependence of LiF:Mg,Cu,P TLDs with the linear energy transfer of secondary electrons followed the same variations as those of published results. The calibration strategy which used TLDs in air exposed to a standard N-250 ISO x-ray beam and TLDs in water irradiated with a standard{sup 137}Cs beam provided an estimated mean uncertainty of 2.8% (k = 1) in the TLD calibration coefficient for irradiations by the {sup 192}Ir source in water. The 3D TLD measurements performed in liquid water were obtained with a maximum uncertainty of 11% (k = 1) found at 1 cm from the source. Radial dose values in water were compared against published results of the American Association of Physicists in Medicine and the European Society for Radiotherapy and Oncology and no significant differences (maximum value of 3.1%) were found within uncertainties except for one position at 9 cm (5.8%). At this location the background contribution relative to the TLD signal is relatively small and an unexpected experimental fluctuation in the background estimate may have caused such a large discrepancy. Conclusions: This paper shows that reliable measurements with TLDs in complex energy spectra require a study of the detector dose response with the radiation quality and specific calibration methodologies which model accurately the experimental conditions where the detectors will be used. The authors have developed and studied a method with highly sensitive TLDs and contributed to its validation by comparison with results from the literature. This methodology can be used to provide direct estimates of the absorbed dose rate in water for irradiations with HDR{sup 192}Ir brachytherapy sources.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

National soil organic carbon estimates can improve global estimates

Global inventories of spatial and vertical distribution of soil organic carbon (SOC) stocks are being used in national and global initiatives targeted to mitigate climate change and land degradation impacts. Yet, national level high-resolution estimates of SOC stocks can be useful for improving the accuracy of global SOC inventories. We estimated spatially resolved SOC stocks of surface 0–30 cm and subsurface 30–100 cm layers at a spatial resolution of 30m in tropical Island, Sri Lanka using a legacy harmonized soil database of 122 soil profiles. The national estimates were compared with two global estimates derived from WISE30sec and SoilGrids250m. The tropical Island (land area=64,610 km2) occupying 0.03% of global land area showed a considerable heterogeneity in SOC stocks ranging from 2.0-342.5 Mg ha -1 and 2.7-391.7 Mg ha -1 in the surface and subsurface soil layers, respectively. We found, elevation, precipitation and slope angle as main environmental controllers of the spatial distribution of SOC stocks under tropical climate. Incorporating the pedogenic information (derived from soil series level legacy map, soil orders and suborders) with environmental controllers resulted in better regression models of predicting surface (R 2 =0.61) and subsurface (R 2 =0.81) SOC stocks. Geographically weighted regression kriging derived maps of SOC stocks revealed that 0–100 cm soil layer of the tropical Island stored 500 Tg C contributing for 0.04% of the global SOC stocks. The validation results of our estimates showed low Mean Estimation Error (MEE: surface -1.6 and subsurface -1.6 Mg ha -1 ) and Root Mean Square Error (RMSE: surface 29.5 and subsurface 24.9 Mg ha -1 ) indicating a low bias and satisfactory predictions. The relative improvement of the prediction accuracy of the SOC stocks of our geospatial estimates in the 0–30 cm layer in comparison to SoilGrids250m and WISE30sec data derived SOC stocks were 51.7% and 35.2%, respectively. The SOC stocks predictions of the 30-100 cm soil layer showed even better relative improvement compared to SoilGrids250m (78.4%) and WISE30sec (57.4%) SOC estimates. Compared to estimates of total SOC stocks resulted in this study, WISE30sec data derived SOC stock maps showed 30% over estimation of the C stock in surface 0-30 cm (332 Tg C) and 41% overestimation in 30-100 cm layer (343 Tg C). The over estimation of total SOC stocks by the SoilGrids250 SOC stocks map for the surface 0-30 cm layer was 122% (567 Tg C) and for the 30-100 cm layer it was 209% (750 Tg C). We conclude that the fusion of legacy soil information of SOC stocks with appropriate environmental covariates and pedogenic information derived from legacy area-class soil maps at national level can produce more accurate inventories of spatial and vertical distribution of SOC stocks. As a result, these national inventories have a great potential of upgrading global inventories of SOC stocks.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2014

Modelflow Estimates of Stroke Volume Do Not Correlate With Doppler Ultrasound Estimates During Upright Posture

Orthostatic intolerance affects 60‐80% of astronauts returning from long‐duration missions, representing a significant risk to completing mission‐critical tasks. While likely multifactorial, a reduction in stroke volume (SV) represents one factor contributing to orthostatic intolerance during stand and head up tilt (HUT) tests. Current measures of SV during stand or HUT tests use Doppler ultrasound and require a trained operator and specialized equipment, restricting its use in the field. BeatScope (Finapres Medical Systems BV, The Netherlands) uses a modelflow algorithm to estimate SV from continuous blood pressure waveforms in supine subjects; however, evidence supporting the use of Modelflow to estimate SV in subjects completing stand or HUT tests remain scarce. Furthermore, because the blood pressure device is held extended at heart level during HUT tests, but allowed to rest at the side during stand tests, changes in the finger arterial pressure waveform resulting from arm positioning could alter modelflow estimated SV. The purpose of this project was to compare Doppler ultrasound and BeatScope estimations of SV to determine if BeatScope can be used during stand or HUT tests. Finger photoplethysmography was used to acquire arterial pressure waveforms corrected for hydrostatic finger‐to‐heart height using the Finometer (FM) and Portapres (PP) arterial pressure devices in 10 subjects (5 men and 5 women) during a stand test while simultaneous estimates of SV were collected using Doppler ultrasound. Measures were made after 5 minutes of supine rest and while subjects stood for 5 minutes. Next, SV estimates were reacquired while each arm was independently raised to heart level, a position similar to tilt testing. Supine SV estimates were not significantly different between all three devices (FM: 68+/-20, PP: 71+/-21, US: 73+/-21 ml/beat). Upon standing, the change in SV estimated by FM (‐18+/-8 ml) was not different from PP (‐21+/-12), but both were significantly less than US (‐37+/-16 ml, p<.05). Raising finger BP devices to heart level caused no significant change in SV measured with any of the devices (FM: 1.5+/-19, PP: 1.7+/-26, US: 0.5+/-6), although variability was 3‐6x greater as assessed by both blood pressure devices compared to US. Retrospective analysis of blood pressure data to assess SV in 11 supine subjects revealed significantly different estimates between methods (FM: 95+/-17, US: 75+/-32, p<.05), but the change in SV resulting from HUT was similar between methods (FM: ‐37+/-9, US: ‐40+/-18 ml). However, the correlation coefficient determined from pairs of SV estimated by US and FM was weak (r2=0.03). These data suggest Modelflow cannot be used in lieu of Doppler ultrasound to estimate SV during stand or HUT tests. Further investigation should focus on identifying factors contributing to differences between these measurement techniques in order to make use of a simple method for assessing beat‐by‐beat changes in SV during postural changes, especially during field testing.

Ferguson, Connor R.↗
DOE OSTI2021

Data-Driven Distributed Algorithms for Estimating Eigenvalues and Eigenvectors of Interconnected Dynamical Systems

Here, the paper presents data-driven algorithms to estimate in a distributed manner the eigenvalues, right and left eigenvectors of an unknown linear (or linearized) interconnected dynamic system. In particular, the proposed algorithms do not require the identification of the system model in advance before performing the estimation. As a first step, we consider interconnected dynamical system with distinct eigenvalues. The proposed strategy first estimates the eigenvalues using the well-known Prony method. The right and left eigenvectors are then estimated by solving distributively a set of linear equations. One important feature of the proposed algorithms is that the topology of communication network used to perform the distributed estimation can be chosen arbitrarily, given that it is connected, and is also independent of the structure or sparsity of the system (state) matrix. The proposed distributed algorithms are demonstrated via a numerical example.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Kalman filter data assimilation: Targeting observations and parameter estimation

This paper studies the effect of targeted observations on state and parameter estimates determined with Kalman filter data assimilation (DA) techniques. We first provide an analytical result demonstrating that targeting observations within the Kalman filter for a linear model can significantly reduce state estimation error as opposed to fixed or randomly located observations. We next conduct observing system simulation experiments for a chaotic model of meteorological interest, where we demonstrate that the local ensemble transform Kalman filter (LETKF) with targeted observations based on largest ensemble variance is skillful in providing more accurate state estimates than the LETKF with randomly located observations. Additionally, we find that a hybrid ensemble Kalman filter parameter estimation method accurately updates model parameters within the targeted observation context to further improve state estimation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗