Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical error”

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,297 records · Page 72

Optimal Scales for Comparing Satellite and Rain-Gauge Rainfall Estimates for Verification Purposes

In spite of all their problems, rain gauges measure rainfall in such a direct way when compared with other methods of estimating rainfall that comparing their totals to satellite estimates remains an essential tool in the validation of satellite products. Some disagreement between averages of satellite data and rain-gauge data is expected because of the very different sampling patterns of the two systems--the satellite provides only occasional snapshots of large areas, whereas rain gauges provide continuous measurements over very small areas. The comparison of the two requires that some quantitative measure be supplied for the amount of disagreement that can be tolerated due to the differences in sampling. As part of an effort to determine the sampling error of satellite averages, a space-time model for rainfall statistics was developed and its parameters fit to radar data from a field experiment conducted near the Inter-Tropical Convergence Zone (ITCZ) in the eastern Atlantic (GATE). Although the model was intended to represent the statistics of relatively large scale fluctuations of rain, it is surprisingly consistent with the very different scales on which rain gauges observe. It can therefore be used to study some of the issues involved with comparing rain-gauge averages to satellite averages. Its implications for the best time and space scales for comparing the two will be discussed.

Bell, Thomas L.↗

Optimal Strategies for Validating TRMM Rainfall Estimates with Rain-Gauge Data

In spite of all their problems, rain gauges measure rainfall in such a direct way when compared with other methods of estimating rainfall that comparing their totals to satellite estimates remains an essential tool in the validation of satellite products. Some disagreement between averages of satellite data and rain-gauge data is expected because of the very different sampling patterns of the two systems--the satellite provides only occasional snapshots of large areas, whereas rain gauges provide continuous measurements over very small areas. A quantitative estimate of the expected level of disagreement due to the differences in sampling is needed in order for comparisons of the two kinds of averages to be informative. As part of an earlier effort to determine the sampling error of satellite averages, a space-time model for rainfall statistics was developed. Sets of parameter values for the model have been obtained for which the model statistics are found to agree quite well with those of radar data taken over the Atlantic during GATE and the western Pacific during TOGA COARE. Although the model was intended to represent the statistics of relatively large scale fluctuations of rain, it is surprisingly consistent with the very different scales on which rain gauges observe. Implications for maximizing the benefit in comparisons of rain gauge averages and TRMM estimates are derived using the model predictions.

Bell, Thomas L.↗

The (55)Fe X-ray Energy Response of Mercury Cadmium Telluride Near-Infrared Detector Arrays

The accurate determination of a detector's fundamental parameters, including read noise, dark current, and QE, relies on a proper measurement of a detector's conversion gain (e- ADU(exp -1)). Charge coupling effects, such as interpixel capacitance, attenuate photon shot noise and result in an overestimation of conversion gain when implementing the photon transfer technique. An approach involving (55)Fe X-rays provides a potentially straightforward measurement of conversion gain by comparing the observed instrumental counts (ADU) to the known charge (e-) liberated by the X-ray. This technique is already preferred within the CCD community, as the pair production energy for silicon is well established. In contrast, to date the pair production energy is unknown for HgCdTe, a material commonly used for near-infrared detectors. In this paper, we derive a preliminary calibration of the (55)Fe X-ray energy response of HgCdTe using 8 HST WFC3 1.7 micrometers flight grade detectors. Our conversion of the X-ray intensities from counts into electrons implements a technique that restores the 'true' gain via classical propagation of errors. For these detectors, our analysis yields preliminary results of good statistical precision: each Ka event generates 1849 +/- 46 electrons, which corresponds to a pair production energy of 3.21+/-f 0.08 eV. We are continuing to assess potential systematic effects to further refine the accuracy of this result.

Fox, Ori↗

Regional Model Nesting Within GFS Daily Forecasts Over West Africa

The study uses the RM3, the regional climate model at the Center for Climate Systems Research of Columbia University and the NASA/Goddard Institute for Space Studies (CCSR/GISS). The paper evaluates 30 48-hour RM3 weather forecasts over West Africa during September 2006 made on a 0.5 grid nested within 1 Global Forecast System (GFS) global forecasts. September 2006 was the Special Observing Period #3 of the African Monsoon Multidisciplinary Analysis (AMMA). Archived GFS initial conditions and lateral boundary conditions for the simulations from the US National Weather Service, National Oceanographic and Atmospheric Administration were interpolated four times daily. Results for precipitation forecasts are validated against Tropical Rainfall Measurement Mission (TRMM) satellite estimates and data from the Famine Early Warning System (FEWS), which includes rain gauge measurements, and forecasts of circulation are compared to reanalysis 2. Performance statistics for the precipitation forecasts include bias, root-mean-square errors and spatial correlation coefficients. The nested regional model forecasts are compared to GFS forecasts to gauge whether nesting provides additional realistic information. They are also compared to RM3 simulations driven by reanalysis 2, representing high potential skill forecasts, to gauge the sensitivity of results to lateral boundary conditions. Nested RM3/GFS forecasts generate excessive moisture advection toward West Africa, which in turn causes prodigious amounts of model precipitation. This problem is corrected by empirical adjustments in the preparation of lateral boundary conditions and initial conditions. The resulting modified simulations improve on the GFS precipitation forecasts, achieving time-space correlations with TRMM of 0.77 on the first day and 0.63 on the second day. One realtime RM3/GFS precipitation forecast made at and posted by the African Centre of Meteorological Application for Development (ACMAD) in Niamey, Niger is shown.

Druyan, Leonard M.↗

Self-Calibration of BICEP1 Three-Year Data and Constraints on Astrophysical Polarization Rotation

Cosmic microwave background (CMB) polarimeters aspire to measure the faint B-mode signature predicted to arise from inflationary gravitational waves. They also have the potential to constrain cosmic birefringence, rotation of the polarization of the CMB arising from parity-violating physics, which would produce nonzero expectation values for the CMB's temperature to B-mode correlation (TB) and E-mode to B-mode correlation (EB) spectra. However, instrumental systematic effects can also cause these TB and EB correlations to be nonzero. In particular, an overall miscalibration of the polarization orientation of the detectors produces TB and EB spectra which are degenerate with isotropic cosmological birefringence, while also introducing a small but predictable bias on the BB spectrum. We find that BICEP1 three-year spectra, which use our standard calibration of detector polarization angles from a dielectric sheet, are consistent with a polarization rotation of alpha = −2.77deg +/- 0.86deg (statistical) +/- 1.3deg (systematic). We have revised the estimate of systematic error on the polarization rotation angle from the two-year analysis by comparing multiple calibration methods. We also account for the (negligible) impact of measured beam systematic effects. We investigate the polarization rotation for the BICEP1 100 GHz and 150 GHz bands separately to investigate theoretical models that produce frequency-dependent cosmic birefringence. We find no evidence in the data supporting either of these models or Faraday rotation of the CMB polarization by the Milky Way galaxy's magnetic field. If we assume that there is no cosmic birefringence, we can use the TB and EB spectra to calibrate detector polarization orientations, thus reducing bias of the cosmological B-mode spectrum from leaked E-modes due to possible polarization orientation miscalibration. After applying this "self-calibration" process, we find that the upper limit on the tensor-to-scalar ratio decreases slightly, from r < 0.70 to r < 0.65 at 95% confidence.

three-year↗

Trend Estimates of AERONET-Observed and Model-Simulated AOTs Between 1993 and 2013

Recently, temporal changes in Aerosol Optical Thickness (AOT) have been investigated based on model simulations, satellite and ground-based observations. Most AOT trend studies used monthly or annual arithmetic means that discard details of the generally right-skewed AOT distributions. Potentially, such results can be biased by extreme values (including outliers). This study additionally uses percentiles (i.e., the lowest 5%, 25%, 50%, 75% and 95% of the monthly cumulative distributions fitted to Aerosol Robotic Network (AERONET)-observed and ECHAM/MESSy Atmospheric Chemistry (EMAC)-model simulated AOTs) that are less affected by outliers caused by measurement error, cloud contamination and occasional extreme aerosol events. Since the limited statistical representativeness of monthly percentiles and means can lead to bias, this study adopts the number of observations as a weighting factor, which improves the statistical robustness of trend estimates. By analyzing the aerosol composition of AERONET-observed and EMAC-simulated AOTs in selected regions of interest, we distinguish the dominant aerosol types and investigate the causes of regional AOT trends. The simulated and observed trends are generally consistent with a high correlation coefficient (R = 0.89) and small bias (slope+/-2(sigma) = 0.75 +/- 0.19). A significant decrease in EMAC-decomposed AOTs by water-soluble compounds and black carbon is found over the USA and the EU due to environmental regulation. In particular, a clear reversal in the AERONET AOT trend percentiles is found over the USA, probably related to the AOT diurnal cycle and the frequency of wildfires. In most of the selected regions of interest, EMAC-simulated trends are mainly attributed to the significant changes of the dominant aerosols; e.g., significant decrease in sea salt and water soluble compounds over Central America, increase in dust over Northern Africa and Middle East, and decrease in black carbon and organic carbon over Australia.

AERONET↗

Optimization of Microphone Locations for Acoustic Liner Impedance Eduction

Two impedance eduction methods are explored for use with data acquired in the NASA Langley Grazing Flow Impedance Tube. The first is an indirect method based on the convected Helmholtz equation, and the second is a direct method based on the Kumaresan and Tufts algorithm. Synthesized no-flow data, with random jitter to represent measurement error, are used to evaluate a number of possible microphone locations. Statistical approaches are used to evaluate the suitability of each set of microphone locations. Given the computational resources required, small sample statistics are employed for the indirect method. Since the direct method is much less computationally intensive, a Monte Carlo approach is employed to gather its statistics. A comparison of results achieved with full and reduced sets of microphone locations is used to determine which sets of microphone locations are acceptable. For the indirect method, each array that includes microphones in all three regions (upstream and downstream hard wall sections, and liner test section) provides acceptable results, even when as few as eight microphones are employed. The best arrays employ microphones well away from the leading and trailing edges of the liner. The direct method is constrained to use microphones opposite the liner. Although a number of arrays are acceptable, the optimum set employs 14 microphones positioned well away from the leading and trailing edges of the liner. The selected sets of microphone locations are also evaluated with data measured for ceramic tubular and perforate-over-honeycomb liners at three flow conditions (Mach 0.0, 0.3, and 0.5). They compare favorably with results attained using all 53 microphone locations. Although different optimum microphone locations are selected for the two impedance eduction methods, there is significant overlap. Thus, the union of these two microphone arrays is preferred, as it supports usage of both methods. This array contains 3 microphones in the upstream hard wall section, 14 microphones opposite the liner, and 3 microphones in the downstream hard wall section.

Jones, M. G.↗

Evaluation Analysis of NASA SMAP L3 and L4 and SPoRT-LIS Soil Moisture Data in the United States

Soil moisture has a critical role in the development, frequency and persistence of climatic and hydrologic extremes such as drought, heat wave and flooding events. In situ soil moisture data are uneven and sparse in time and space. This highlights the need to utilize other soil moisture sources to fill this spatiotemporal gap. The goal of this study is to validate one satellite-based and two model-based soil moisture datasets with in situ data across the United States. Soil moisture information from the Soil Moisture Active Passive (SMAP) enhanced level 3 (L3) (SMAP L3) and modeled level 4 (L4) (SMAP L4) data at 9-km resolution and Short-term Prediction Research and Transition-Land Information System (SPoRT-LIS) modeled at 3-km resolution were selected for evaluation. SPoRT-LIS is a near real-time, high resolution operational land analysis data. Ground-based data were obtained from the North American Soil Moisture Database (NASMD) for 362 stations. Seven statistical indicators including anomaly, Spearman, and Pearson correlation coefficients, the systematic error (Bias), root mean square error (RMSE), unbiased root mean square error (ubRMSE), and normalized standard deviation (SDV) were used to evaluate the satellite- and model-derived soil moisture data. In addition, the triple collocation (TC) error model was used to measure the error among SMAP L4, SPoRT-LIS and ground-based data. This study assesses which satellite or modeled dataset is most appropriate for specific times and locations to use as a surrogate for in situ observations. Temporal and spatial analysis demonstrated that, overall, SMAP L4 performed better than SMAP L3 and SPoRT-LIS. Strong agreement was observed between SMAP L4 and in situ observations (ρ = 0.53, Bias = −0.006) in all seasons and most regions with various land covers, especially in winter and in the central regions of the United States. For croplands, SMAP L4 presented the best agreement with in situ data, analyzing all period (ρ = 0.60) and non-winter period (ρ = 0.61) separately.

Tavakol, Ameneh↗

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗

GeoXO Atmospheric Composition Instrument (AC) Measurements Conditions Assessment

NOAA’s next generation Geostationary Extended Observations (GeoXO) satellite system will advance Earth observations from geostationary orbit. GeoXO will supply vital information supporting the U.S. weather, ocean, and climate operations. The recommended three-satellite constellation includes spacecraft at GOES-Central position, that will carry an Atmospheric Composition Instrument (ACX) among its payloads; the requirement on the ACX is to scan its region of interest in one hour. An assessment of ACX observational conditions is presented, based on view geometry and solar angle range over the year. The availability for observations of each pixel is estimated per day of the year, as function of the geometric air mass factor (AMF), which can be related to signal to noise ratio using a model of the ACX instrument performance, and ultimately to expected column error for specific atmospheric species. In addition, cloud climatology – based on statistics derived from GOES-16 ABI cloud product data – is considered to determine typical fractions of cloud-free ACX observations and typical cloud optical depths, for each pixel as function of time of day and season. The results of such assessments can be used to prepare for GeoXO ACX data and further the understanding and development of the future instrument operations.

Boryana Efremova↗

Bayesian Model Averaging for Ensemble-Based Estimates of Solvation Free Energies

This paper applies the Bayesian Model Averaging (BMA) statistical ensemble technique to estimate small molecule solvation free energies. There is a wide range methods for predicting solvation free energies, ranging from empirical statistical models to {\it ab initio} quantum mechanical approaches. Each of these methods are based on a set of conceptual assumptions that can affect a method's predictive accuracy and transferability. Using an iterative statistical process, we have selected and combined solvation energy estimates using an ensemble of 17 diverse methods from the SAMPL4 blind prediction study to form a single, aggregated solvation energy estimate. The ensemble design process evaluates the statistical information in each individual method as well as the performance of the aggregate estimate obtained from the ensemble as a whole. Methods that possess minimal or redundant information are pruned from the ensemble and the evaluation process repeats until aggregate predictive performance can no longer be improved. We show that this process results in a final aggregate estimate that outperforms all individual methods by reducing estimate errors by as much as 91% to 1.2 kcal/mol accuracy. We also compare our iterative refinement approach to other statistical ensemble approaches and demonstrate that this iterative process reduces estimate errors by as much as 61%. This work provides a new approach for accurate solvation free energy prediction and lays the foundation for future work on aggregate models that can balance computational cost with predictive accuracy.

solvation, Bayesian model averaging, Electrostatic↗

Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies

The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H 0 , the σ 8 – S 8 tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0 σ tension between the Planck CMB estimate of the Hubble constant H 0 and the SH0ES collaboration measurements. After showing the H 0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ω m , and the amplitude or rate of the growth of structure ( σ 8 , f σ 8 ). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H 0 – S 8 tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Practical and reliable error bars for quantum process tomography

Current techniques in quantum process tomography typically return a single point estimate of an unknown process based on a finite albeit large amount of measurement data. Due to statistical fluctuations, however, other processes close to the point estimate can also produce the observed data with near certainty. Unless appropriate error bars can be constructed, the point estimate does not carry any sound operational interpretation. Here, we provide a solution to this problem by constructing a confidence region estimator for quantum processes. Our method enables reliable estimation of essentially any figure of merit for quantum processes on few qubits, including the diamond distance to a specific noise model, the entanglement fidelity, and the worst-case entanglement fidelity, by identifying error regions which contain the true state with high probability. We also provide a software package, QPtomographer, implementing our estimator for the diamond norm and the worst-case entanglement fidelity. We illustrate its usage and performance with several simulated examples. Furthermore, our tools can be used to reliably certify the performance of, e.g., error correction codes, implementations of unitary gates, or more generally any noise process affecting a quantum system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Solar activity prediction

A statistical study of formulas for predicting the sunspot number several years in advance is reported. By using a data lineup with cycle maxima coinciding, and by using multiple and nonlinear predictors, a new formula which gives better error estimates than former formulas derived from the work of McNish and Lincoln is obtained. A statistical analysis is conducted to determine which of several mathematical expressions best describes the relationship between 10.7 cm solar flux and Zurich sunspot numbers. Attention is given to the autocorrelation of the observations, and confidence intervals for the derived relationships are presented. The accuracy of predicting a value of 10.7 cm solar flux from a predicted sunspot number is dicussed.

Slutz, R. J.↗

Examination of the relationship between outgoing infrared window and total longwave fluxes using satellite data

The relationship between narrowband and broadband thermal radiances is explored to determine the accuracy of outgoing longwave radiation derived from narrowband data. Infrared window (10.2-12.2 microns) data from the Geostationary Operational Environment Satellite are correlated with longwave (5.0-50.0 microns) data from the Earth Radiation Budget Experiment. A simple quadratic fit between the narrowband and longwave fluxes results in standard errors of 4.4-5.3 percent for data that are matched closely in time and space. The use of matched regional flux data with temporal differences up to 59 minutes yields standard errors of 4.1-5.4 percent. About 30 percent of the error may be attributed to limb darkening and spatial and temporal differences in the matched fluxes. The relationship shows a statistically significant dependence on the relative humidity of the atmosphere above the radiating surface. Although this dependency accounts for only about 1 percent of the standard error, it reduces the monthly mean regional errors by more than 10 percent. Data taken over land produced a relationship slightly different from data taken over water. The differences appear to be primarily due to daytime heating of the land surface.

Minnis, Patrick↗

Suboptimal schemes for atmospheric data assimilation based on the Kalman filter

This work is directed toward approximating the evolution of forecast error covariances for data assimilation. The performance of different algorithms based on simplification of the standard Kalman filter (KF) is studied. These are suboptimal schemes (SOSs) when compared to the KF, which is optimal for linear problems with known statistics. The SOSs considered here are several versions of optimal interpolation (OI), a scheme for height error variance advection, and a simplified KF in which the full height error covariance is advected. To employ a methodology for exact comparison among these schemes, a linear environment is maintained, in which a beta-plane shallow-water model linearized about a constant zonal flow is chosen for the test-bed dynamics. The results show that constructing dynamically balanced forecast error covariances rather than using conventional geostrophically balanced ones is essential for successful performance of any SOS. A posteriori initialization of SOSs to compensate for model - data imbalance sometimes results in poor performance. Instead, properly constructed dynamically balanced forecast error covariances eliminate the need for initialization. When the SOSs studied here make use of dynamically balanced forecast error covariances, the difference among their performances progresses naturally from conventional OI to the KF. In fact, the results suggest that even modest enhancements of OI, such as including an approximate dynamical equation for height error variances while leaving height error correlation structure homogeneous, go a long way toward achieving the performance of the KF, provided that dynamically balanced cross-covariances are constructed and that model errors are accounted for properly. The results indicate that such enhancements are necessary if unconventional data are to have a positive impact.

Todling, Ricardo↗

SU-F-T-252: An Investigation of Gamma Knife Frame Definition Error When Using a Pre-Planning Workflow

Purpose: To determine causal factors related to high frame definition error when treating GK patients using a pre-planning workflow. Methods: 160 cases were retrospectively reviewed. All patients received treatment using a pre-planning workflow whereby stereotactic coordinates are determined from a CT scan acquired after framing using a fiducial box. The planning software automatically detects the fiducials and compares their location to expected values based on the rigid design of the fiducial system. Any difference is reported as mean and maximum frame definition error. The manufacturer recommends these values be less than 1.0 mm and 1.5 mm. In this study, frame definition error was analyzed in comparison with a variety of factors including which neurosurgeon/oncologist/physicist was involved with the procedure, number of post used during framing (3 or 4), type of lesion, and which CT scanner was utilized for acquisition. An analysis of variance (ANOVA) approach was used to statistically evaluate the data and determine causal factors related to instances of high frame definition error. Results: Two factors were identified as significant: number of post (p=0.0003) and CT scanner (p=0.0001). Further analysis showed that one of the four scanners was significantly different than the others. This diagnostic scanner was identified as an older model with localization lasers not tightly calibrated. The average value for maximum frame definition error using this scanner was 1.48 mm (4 posts) and 1.75 mm (3 posts). For the other scanners this value was 1.13 mm (4 posts) and 1.40 mm (3 posts). Conclusion: In utilizing a pre-planning workflow the choice of CT scanner matters. Any scanner utilized for GK should undergo routine QA at a level appropriate for radiation oncology. In terms of 3 vs 4 post, it is hypothesized that three posts provide less stability during CT acquisition. This will be tested in future work.

60 APPLIED LIFE SCIENCES↗