Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data bias”

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 307 records · Page 17

Comparisons of Satellite and Airborne Altimetry with Ground-Based Data from the Interior of the Antarctic Ice Sheet

A series of traverses has been conducted for validation of the NASA Ice, Cloud, and land Elevation Satellite 2 on the flat interior of the Antarctic ice sheet. GNSS data collected on three separate 88S Traverses intersect 20% of the ICESat-2 reference ground tracks and have precisions of better than ±7 cm and biases of less than ~4 cm. Data from these traverses were used to assess heights from ICESat-2, CryoSat-2, and Airborne Topographic Mapper (ATM). ICESat-2 heights have better than ±3.3 cm bias and better than ±7.2 cm precision. ATM heights have better than 9.3 cm bias and better than ±9.6 cm precision. CryoSat-2 heights have -38.9 cm of bias and ± 47.3 cm precision. These best-case results are from the flat ice-sheet interior but provide a characterization of the quality of satellite and airborne altimetry.

K. M. Brunt↗

Sensitivity of the Tropical Atmospheric Energy Balance to ENSO-Related SST Changes: Comparison of Climate Model Simulations to Observed Responses

This paper focuses on how fresh water and radiative fluxes over the tropical oceans change during ENSO warm and cold events and how these changes affect the tropical energy balance. At present, ENSO remains the most prominent known mode of natural variability at interannual time scales. While this natural perturbation to climate is quite distinct from possible anthropogenic changes in climate, adjustments in the tropical water and energy budgets during ENSO may give insight into feedback processes involving water vapor and cloud feedbacks. Although great advances have been made in understanding this phenomenon and realizing prediction skill over the past decade, our ability to document the coupled water and energy changes observationally and to represent them in climate models seems far from settled (Soden, 2000 J Climate). In a companion paper we have presented observational analyses, based principally on space-based measurements which document systematic changes in rainfall, evaporation, and surface and top-of-atmosphere (TOA) radiative fluxes. Here we analyze several contemporary climate models run with observed SSTs over recent decades and compare SST-induced changes in radiation, precipitation, evaporation, and energy transport to observational results. Among these are the NASA / NCAR Finite Volume Model, the NCAR Community Climate Model, the NCEP Global Spectral Model, and the NASA NSIPP Model. Key disagreements between model and observational results noted in the recent literature are shown to be due predominantly to observational shortcomings. A reexamination of the Langley 8-Year Surface Radiation Budget data reveals errors in the SST surface longwave emission due to biased SSTs. Subsequent correction allows use of this data set along with ERBE TOA fluxes to infer net atmospheric radiative heating. Further analysis of recent rainfall algorithms provides new estimates for precipitation variability in line with interannual evaporation changes inferred from the da Silva, Young, Levitus COADS analysis. The overall results from our analysis suggest an increase (decrease) of the hydrologic cycle during ENSO warm (cold) events at the rate of about 5 W/sq m per K of SST change. Model results agree reasonably well with this estimate of sensitivity. This rate is slightly less than that which would be expected for constant relative humidity over the tropical oceans. There remain, however, significant quantitative uncertainties in cloud forcing changes in the models as compared to observations. These differences are examined in relationship to model convection and cloud parameterizations Analysis of the possible sampling and measurement errors compared to systematic model errors is also presented.

Robertson, Franklin R.↗

Impact of High Resolution SST Data on Regional Weather Forecasts

Past studies have shown that the use of coarse resolution SST products such as from the real-time global (RTG) SST analysis[1] or other coarse resolution once-a-day products do not properly portray the diurnal variability of fluxes of heat and moisture from the ocean that drive the formation of low level clouds and precipitation over the ocean. For example, the use of high resolution MODIS SST composite [2] to initialize the Advanced Research Weather Research and Forecast (WRF) (ARW) [3] has been shown to improve the prediction of sensible weather parameters in coastal regions [4][5}. In an extend study, [6] compared the MODIS SST composite product to the RTG SST analysis and evaluated forecast differences for a 6 month period from March through August 2007 over the Florida coastal regions. In a comparison to buoy data, they found that that the MODIS SST composites reduced the bias and standard deviation over that of the RTG data. These improvements led to significant changes in the initial and forecasted heat fluxes and the resulting surface temperature fields, wind patterns, and cloud distributions. They also showed that the MODIS composite SST product, produced for the Terra and Aqua satellite overpass times, captured a component of the diurnal cycle in SSTs not represented in the RTG or other one-a-day SST analyses. Failure to properly incorporate these effects in the WRF initialization cycle led to temperature biases in the resulting short term forecasts. The forecast impact was limited in some situations however, due to composite product inaccuracies brought about by data latency during periods of long-term cloud cover. This paper focuses on the forecast impact of an enhanced MODIS/AMSR-E composite SST product designed to reduce inaccuracies due data latency in the MODIS only composite product.

Jedlovec, Gary J.↗

On the appropriateness of applying chi-square distribution based confidence intervals to spectral estimates of helicopter flyover data

The validity of applying chi-square based confidence intervals to far-field acoustic flyover spectral estimates was investigated. Simulated data, using a Kendall series and experimental acoustic data from the NASA/McDonnell Douglas 500E acoustics test, were analyzed. Statistical significance tests to determine the equality of distributions of the simulated and experimental data relative to theoretical chi-square distributions were performed. Bias and uncertainty errors associated with the spectral estimates were easily identified from the data sets. A model relating the uncertainty and bias errors to the estimates resulted, which aided in determining the appropriateness of the chi-square distribution based confidence intervals. Such confidence intervals were appropriate for nontonally associated frequencies of the experimental data but were inappropriate for tonally associated estimate distributions. The appropriateness at the tonally associated frequencies was indicated by the presence of bias error and noncomformity of the distributions to the theoretical chi-square distribution. A technique for determining appropriate confidence intervals at the tonally associated frequencies was suggested.

Rutledge, Charles K.↗

New Methods for Retrieval of Chlorophyll Red Fluorescence from Hyperspectral Satellite Instruments: Simulations and Application to GOME-2 and SCIAMACHY

Global satellite measurements of solar-induced fluorescence (SIF) from chlorophyll over land and ocean have proven useful for a number of different applications related to physiology, phenology, and productivity of plants and phytoplankton. Terrestrial chlorophyll fluorescence is emitted throughout the red and far-red spectrum, producing two broad peaks near 683 and 736nm. From ocean surfaces, phytoplankton fluorescence emissions are entirely from the red region (683nm peak). Studies using satellite-derived SIF over land have focused almost exclusively on measurements in the far red (wavelengths greater than 712nm), since those are the most easily obtained with existing instrumentation. Here, we examine new ways to use existing hyperspectral satellite data sets to retrieve red SIF (wavelengths less than 712nm) over both land and ocean. Red SIF is thought to provide complementary information to that from the far red for terrestrial vegetation. The satellite instruments that we use were designed to make atmospheric trace-gas measurements and are therefore not optimal for observing SIF; they have coarse spatial resolution and only moderate spectral resolution (0.5nm). Nevertheless, these instruments, the Global Ozone Monitoring Instrument 2 (GOME-2) and the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY), offer a unique opportunity to compare red and far-red terrestrial SIF at regional spatial scales. Terrestrial SIF has been estimated with ground-, aircraft-, or satellite-based instruments by measuring the filling-in of atmospheric andor solar absorption spectral features by SIF. Our approach makes use of the oxygen (O2) gamma band that is not affected by SIF. The SIF-free O2 gamma band helps to estimate absorption within the spectrally variable O2 B band, which is filled in by red SIF. SIF also fills in the spectrally stable solar Fraunhofer lines (SFLs) at wavelengths both inside and just outside the O2 B band, which further helps to estimate red SIF emission. Our approach is then an extension of previous approaches applied to satellite data that utilized only the filling-in of SFLs by red SIF. We conducted retrievals of red SIF using an extensive database of simulated radiances covering a wide range of conditions. Our new algorithm produces good agreement between the simulated truth and retrievals and shows the potential of the O2 bands for noise reduction in red SIF retrievals as compared with approaches that rely solely on SFL filling. Biases seen with existing satellite data, most likely due to instrumental artifacts that vary in time, space, and with instrument, must be addressed in order to obtain reasonable results. Our 8-year record of red SIF observations over land with the GOME-2 allows for the first time reliable global mapping of monthly anomalies. These anomalies are shown to have similar spatiotemporal structure as those in the far red, particularly for drought-prone regions. There is a somewhat larger percentage response in the red as compared with the far red for these areas that are drought sensitive. We also demonstrate that good-quality ocean fluorescence line height retrievals can be achieved with GOME-2, SCIAMACHY, and similar instruments by utilizing the full complement of radiance measurements that span the red SIF emission feature.

spatiotemporal structure↗

Skill Assessment of a Spectral Ocean-Atmosphere Radiative Model

Ocean phytoplankton, detrital material, and water absorb and scatter light spectrally. The Ocean- Atmosphere Spectral Irradiance Model (OASIM) is intended to provide surface irradiance over the oceans with sufficient spectral resolution to support ocean ecology, biogeochemistry, and heat exchange investigations, and of sufficient duration to support inter-annual and decadal investigations. OASIM total surface irradiance (integrated 200 nm to 4 microns) was compared to in situ data and three publicly available global data products at monthly 1-degree resolution. OASIM spectrally-integrated surface irradiance had root mean square (RMS) difference= 20.1 W/sq m (about 11%), bias=1.6 W/sq m (about 0.8%), regression slope= 1.01 and correlation coefficient= 0.89, when compared to 2322 in situ observations. OASIM had the lowest bias of any of the global data products evaluated (ISCCP-FD, NCEP, and ISLSCP 11), and the best slope (nearest to unity). It had the second best RMS, and the third best correlation coefficient. OASIM total surface irradiance compared well with ISCCP-FD (RMS= 20.7 W/sq m; bias=-11.4 W/sq m, r=0.98) and ISLSCP II (RMS =25.2 W/sq m; bias= -13.8 W/sq m; r=0.97), but less well with NCEP (RMS =43.0 W/sq m ;bias=-22.6 W/sq m; x=0.91). Comparisons of OASIM photosynthetically available radiation (PAR) with PAR derived from SeaWiFS showed low bias (-1.8 mol photons /sq m/d, or about 5%), RMS (4.25 mol photons /sq m/d ' or about 12%), near unity slope (1.03) and high correlation coefficient (0.97). Coupled with previous estimates of clear sky spectral irradiance in OASIM (6.6% RMS at 1 nm resolution), these results suggest that OASIM provides reasonable estimates of surface broadband and spectral irradiance in the oceans, and can support studies on ocean ecosystems, carbon cycling, and heat exchange.

Gregg, Watson, W.↗

Sea Surface Salinity Distribution in the Southern Ocean as Observed from Space

Large-scale spatial and temporal variabilities of sea surface salinity (SSS) in the Southern Ocean from 2011 to 2017 were studied using products derived from microwave sensors on board Aquarius, Soil Moisture and Ocean Salinity (SMOS), and Soil Moisture Active and Passive (SMAP) satellites. Four products, three from Aquarius and one from SMOS, were evaluated and shown to be generally consistent within 0.3 to 0.6 psu and agree favorably with in situ measurements. However, although the Aquarius products show consistent seasonality of SSS with high values of 34.45 psu in October and low values of 33.40 psu in May, theSMOS and SMAP products lack such seasonal variations. This may be caused by larger uncertainties in the SMOS and SMAP data due in part to the lack of concurrent scatterometer measurements that is used to correct for roughness effects. The four products provide similar spatial distributions of SSS with root mean square difference from 0.25 to 0.58 psu. Differences among Aquarius products are mainly due to varying salinity retrieval algorithms, smoothing, and masking of sea ice, while the SMOS product showed the highest SSS deviation that is likely due to the bias adjustment done on the data set. Our analyses show that SSS in the Southern Ocean region has significant meridional variations with the lowest SSS near the ice edge and highest at lower latitudes. The SSS is also lowest in summer indicating the predominant influence of sea ice and glacial melt, but it stays low near ice edges even during the growth season. Salinity together with temperature is the key parameters that control the state and circulation of the World's oceans. The large‐scale distribution of the sea surface salinity of the Southern Ocean is quantified for the first time using Aquarius, Soil Moisture and Ocean Salinity, and Soil Moisture Active and Passive satellite microwave data. Seasonal and interannual changes in the distribution as affected by the sea ice cover, surface temperature, and precipitation are evaluated. Comparative studies were done, using four different products that are publicly available, in conjunction with in situ observations to gain insights into the true nature of the distribution and how consistently the sea surface salinity is depicted by the different products. There are general consistencies in the products, and discrepancies are attributed to different algorithms, smoothing techniques, and sea ice masking. Aquarius data are shown to have higher accuracy than Soil Moisture and Ocean Salinity and Soil Moisture Active and Passive products in part because of having concurrent scatterometer that provides accurate correction to roughness effects.

Garcia-Eidell, Cynthia↗

A noise source identification technique using an inverse Helmholtz integral equation method

A technique is developed which utilizes numerical models and field pressure information to characterize acoustic fields and identify acoustic sources. The numerical models are based on boundary element numerical procedures. Either pressure, velocity, or passive boundary conditions, in the form of impedance boundary conditions, may be imposed on the numerical model. Alternatively, if no boundary information is known, a boundary condition can be left unspecified. Field pressure data may be specified to overdetermine the numerical problem. The problem is solved numerically for the complete sound field from which the acoustic sources may be determined. The model can then be used to idenfify acoustic intensity paths in the field. The solution can be modified and the model used to evaluate design alternatives. In this investigation the method is tested analytically and verified. In addition, the sensitivity of the method to random and bias error in the input data is demonstrated.

Gardner, B. K.↗

Identification of Sea Surface Temperature and Sea Surface Salinity Fronts along the California Coast: Application Using Saildrone and Satellite Derived Products

Coastal upwelling regions are one of the most dynamic areas of the world’s oceans. The California and Baja California Coasts are impacted by both coastal upwelling and the California Current, leading to frontal activity that is captured by gradients in both Sea Surface Temperature (SST) and Sea Surface Salinity (SSS). Satellite data are a great source of spatial data to study fronts. However, biases near coastal areas and coarse resolutions can impair its usefulness in upwelling areas. In this work gradients in SST from NASA Multi-Scale Ultra-High Resolution (MUR) and in two SSS products derived from the Soil Moisture Active Passive (SMAP) NASA mission are compared directly with gradients derived from the Saildrone uncrewed vehicles to validate the gradients as well as to assess their ability to detect known frontal features. The three remotely sensed data sets (MURSST/JPL, SMAP SSS/RSS, SMAP SSS) were co-located with the Saildrone data prior to the calculation of the gradients. Wavelet analysis is used to determine how well the satellite derived SST and SSS products are reproducing the Saildrone derived gradients. Overall results indicate the remote sensing products are reproducing features of known areas of coastal upwelling. Differences between the SST and SSS gradients are mainly associated with the limitations of the microwave derived SSS coverage near land and its reduced spatial resolution. The results are promising for using remote sensing data sets to monitor frontal structure along the California Coast and the application to long term changes in coastal upwelling and dynamics.

Sea Surface Salinity↗

Bias Reduction in Short Records of Satellite Soil Moisture

Although surface soil moisture data from different sources (satellite retrievals, ground measurements, and land model integrations of observed meteorological forcing data) have been shown to contain consistent and useful information in their seasonal cycle and anomaly signals, they typically exhibit very different mean values and variability. These biases pose a severe obstacle to exploiting the useful information contained in satellite retrievals through data assimilation. A simple method of bias removal is to match the cumulative distribution functions (cdf) of the satellite and model data. However, accurate cdf estimation typically requires a long record of satellite data. We demonstrate here that by wing spatial sampling with a 2 degree moving window we can obtain local statistics based on a one-year satellite record that are a good approximation to those that would be derived from a much longer time series. This result should increase the usefulness of relatively short satellite data records.

Reichle, Rolf H.↗

Ionospheric Signatures in Radio Occultation Data

We can extend robustly the radio occultation data record by 6 years (+60%) by developing a singlefrequency processing method for GPS/MET data. We will produce a calibrated data set with profile-byprofile data characterization to determine robust upper bounds on ionospheric bias. Part of an effort to produce a calibrated RO data set addressing other key error sources such as upper boundary initialization. Planned: AIRS-GPS water vapor cross validation (water vapor climatology and trends).

radio occultation↗

Single-Receiver GPS Phase Bias Resolution

Existing software has been modified to yield the benefits of integer fixed double-differenced GPS-phased ambiguities when processing data from a single GPS receiver with no access to any other GPS receiver data. When the double-differenced combination of phase biases can be fixed reliably, a significant improvement in solution accuracy is obtained. This innovation uses a large global set of GPS receivers (40 to 80 receivers) to solve for the GPS satellite orbits and clocks (along with any other parameters). In this process, integer ambiguities are fixed and information on the ambiguity constraints is saved. For each GPS transmitter/receiver pair, the process saves the arc start and stop times, the wide-lane average value for the arc, the standard deviation of the wide lane, and the dual-frequency phase bias after bias fixing for the arc. The second step of the process uses the orbit and clock information, the bias information from the global solution, and only data from the single receiver to resolve double-differenced phase combinations. It is called "resolved" instead of "fixed" because constraints are introduced into the problem with a finite data weight to better account for possible errors. A receiver in orbit has much shorter continuous passes of data than a receiver fixed to the Earth. The method has parameters to account for this. In particular, differences in drifting wide-lane values must be handled differently. The first step of the process is automated, using two JPL software sets, Longarc and Gipsy-Oasis. The resulting orbit/clock and bias information files are posted on anonymous ftp for use by any licensed Gipsy-Oasis user. The second step is implemented in the Gipsy-Oasis executable, gd2p.pl, which automates the entire process, including fetching the information from anonymous ftp

Bertiger, William I.↗

Real-time Sub-cm Differential Orbit Determination of two Low-Earth Orbiters with GPS Bias Fixing

An effective technique for real-time differential orbit determination with GPS bias fixing is formulated. With this technique, only real-time GPS orbits and clocks are needed (available from the NASA Global Differential GPS System with 10-20 cm accuracy). The onboard, realtime orbital states of user satellites (few meters in accuracy) are used for orbit initialization and integration. An extended Kalman filter is constructed for the estimation of the differential orbit between the two satellites as well as a reference orbit, together with their associating dynamics parameters. Due to close proximity of the two satellites and of similar body shapes, the differential dynamics are highly common and can be tightly constrained which, in turn, strengthens the orbit estimation. Without explicit differencing of GPS data, double-differenced phase biases are formed by a transformation matrix. Integer-valued fixing of these biases are then performed which greatly strengthens the orbit estimation. A 9-day demonstration between GRACE orbits with baselines of approx.200 km indicates that approx.80% of the double-differenced phase biases can successfully be fixed and the differential orbit can be determined to approx.7 mm as compared to the results of onboard K-band ranging.

global positioning system (GPS)↗

Combinations of Earth Orientation Measurements: SPACE94, COMB94, and POLE94

A Kalman filter has been used to combine independent measurements of the Earth's orientation taken by the space-geodetic observing techniques of lunar laser ranging, satellite laser ranging, very long baseline interferometry, and the Global Positioning System. Prior to their combination, the data series were adjusted to have the same bias and rate, the stated uncertainties of the measurements were adjusted, and data points considered to be outliers were deleted. The resulting combination, SPACE94, consists of smoothed, interpolated polar motion and UT1-UTC values spanning October 6.0, 1976, to January 27.0, 1995, at 1-day intervals.

Kalman↗

Section drag coefficients from pressure probe transverses of a wing wake at low speeds

An in-flight wing wake section drag investigation was conducted using traversing pitot and static probes. The primary objective was to develop measurement techniques and improve the accuracy of in-flight wing profile drag measurements for low values of dynamic pressure and Reynolds number. Data were obtained on a sailplane for speeds from about 40 knots to 125 knots at chord Reynolds numbers between 1,000,000 and 3,000,000. Tests were conducted with zero flap deflection, deflected flaps, and various degrees of surface roughness, and for smooth and rough atmospheric conditions. Several techniques were used to increase data reliability and to minimize certain bias errors. A discussion of the effects of a total pressure probe in a pressure gradient, and the effects of discrete turbulence levels, on the data presented and other experimental results is also included.

Montoya, L. C.↗

Creating A Consistent Historical NASA POWER Solar Radiation Dataset to Support Renewable Energy, Building Energy Efficiency and Agro-Climatology Decisions

Prediction of Worldwide Energy Resources (POWER) project provides irradiance dataset to support renewable energy, building energy efficiency and agricultural needs. These datasets are derived from Global Energy and Water Cycle Experiment Surface Radiation Budget (GEWEX SRB) and Clouds and the Earth’s Radiant Energy System (CERES SYN1Deg). A systematic bias has been reported between these two datasets for the years with overlapping observations. For obtaining a consistent climate data record spanning the entire time record of observations, it is crucial to understand and remove the bias in the irradiance dataset. Inconsistency in solar radiation data can lead to inaccurate conclusions about solar energy potential and obscure real trends in solar radiation patterns that would impact energy availability assessments. In this study, we adapt quantile mapping approach to remove the systematic bias and to improve reliability of shortwave and longwave irradiance data. We present a validation of the bias corrected data against ground truth. For each 1° latitude and 1° longitude grid box across the globe, we match the CDFs of the reference dataset (CERES SYN1Deg) to that of the SRB dataset, thereby, adjusting the irradiance values to match the empirical distribution of two different measurements. The performance of quantile mapping is evaluated by using the metrics such as Mean Absolute Deviation (MAD). The results indicate that the quantile mapping significantly improves the accuracy and reliability of solar irradiance dataset especially for the weather conditions associated with high cloud cover and extreme irradiance values. The initial range of MAD for the studied sites for daily data was 4 to 19 Wm-2. After correction these reduced to 3 to 7 Wm-2. The findings from this study have important implications for solar energy system design, agricultural planning, and climate modeling community. Reducing the inconsistency and biases in solar irradiance dataset can enable better planning and operation of solar energy systems, leading to increased efficiency and cost-effectiveness. Additionally, this work also contributes to the statistical post-processing techniques in the renewable energy domain and highlights the potential of historical and near-real-time NASA POWER dataset as a valuable resource for solar energy research and applications.

POWER↗

Chandra ACIS Observations of Jovian X-Ray Emission

On November 25 and 26, 1999, the Chandra X-ray spacecraft conducted a set of four 19,000 sec observations of Jupiter. The ACIS-S instrument configuration was used for its good low energy efficiency and spatial resolution. An anomalous response was obtained which was subsequently attributed to strong jovian infrared radiation penetrating the detector and piling up spurious events across the entire X-ray range. However, the pre-observation establishment of an offsetting bias field has allowed the recovery of data from that portion of Jupiter's disc which remained within the elevated portion of the bias field during the observation. This ranges from fewer than 3000 sec to the entire observing time for about 10% of the planet. Auroral emission is seen near both poles in each observation. The northern aurora ia overall more intense than the southern, consistent with prior Einstein and ROSAT Observatory results. The southern aurora shows more modulation with Jupiter's rotation than the northern. Spatial resolution has been improved by at least a factor of two over prior measurements but convincing evidence of structure has not been seen. Lower latitude emission, first observed by ROSAT, is confirmed with flux levels averaging more than a factor of five below peak auroral values. Pronounced variation in the observed emission has occurred over the observing period. The spectral response extends from 0.24 keV, below which noise dominates, to about 1.2 keV. For all four observations the spectrum is clearly enhanced between 0.45 and 0.85 keV. This is apparently unequivocal evidence that Jupiter's X-ray emission is the result of oxygen and perhaps sulfur ions precipitating into the planet's atmosphere, where they undergo charge exchange interactions. The identification of specific transitions lines in the spectrum is among the ongoing efforts. A bremsstrahlung component has not yet been identified.

Garmire, Gordon↗

HZETRN Radiation Transport Validation Using Balloon-Based Experimental Data

The deterministic radiation transport code HZETRN (High charge (Z) and Energy TRaNsport) was developed by NASA to study the effects of cosmic radiation on astronauts and instrumentation shielded by various materials. This work presents an analysis of computed differential flux from HZETRN compared with measurement data from three balloon-based experiments over a range of atmospheric depths, particle types, and energies. Model uncertainties were quantified using an interval-based validation metric that takes into account measurement uncertainty both in the flux and the energy at which it was measured. Average uncertainty metrics were computed for the entire dataset as well as subsets of the measurements (by experiment, particle type, energy, etc.) to reveal any specific trends of systematic over- or under-prediction by HZETRN. The distribution of individual model uncertainties was also investigated to study the range and dispersion of errors beyond just single scalar and interval metrics. The differential fluxes from HZETRN were generally well-correlated with balloon-based measurements; the median relative model difference across the entire dataset was determined to be 30%. The distribution of model uncertainties, however, revealed that the range of errors was relatively broad, with approximately 30% of the uncertainties exceeding ± 40%. The distribution also indicated that HZETRN systematically under-predicts the measurement dataset as a whole, with approximately 80% of the relative uncertainties having negative values. Instances of systematic bias for subsets of the data were also observed, including a significant underestimation of alpha particles and protons for energies below 2.5 GeV/u. Muons were found to be systematically over-predicted at atmospheric depths deeper than 50 g/cm(sup 2) but under-predicted for shallower depths. Furthermore, a systematic under-prediction of alpha particles and protons was observed below the geomagnetic cutoff, suggesting that improvements to the light ion production cross sections in HZETRN should be investigated.

James E Warner↗