Search NASA⌕ Search

SEARCH · Search NASA

Results for “scatterplots”

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.

Evaluation of passive microwave precipitation algorithms in wintertime midlatitude situations

The second intercomparison project of the Global Precipitation Climatology Project (GPCP) examined the estimation of midlatitude, cool-season precipitation. As part of that effort, the authors report here on the results of two microwave techniques, the Goddard scattering algorithm and the physical retrieval algorithm of Kummerow. Results from the estimation of instantaneous rain rate for five overpasses of the Special Sensor Microwave/Imager (SSM/I) are presented in a case study mode to illustrate both the strong and weak points of each technique. These five cases represent a sampling of the various types of precipitating systems observed. Results for the complete set of 20 swaths chosen by the United Kingdom Meteorological Office (UKMO) are then categorized by scatterplots and statistics of instantaneous radar versus microwave-estimated rain rate, rain no-rain contingency tables, and scatterplots of areal coverage of rainfall. Neither algorithm produced a good statistical correlation with the radar data, yet in general, both did well at determining rainy areas. Two reasons are suggested for the low correlation coefficients between both algorithms and the radar data. Time differences between the SSM/I overpass and the radar observations can occasionally account for some of the differences. The primary reason for the low correlations, however, appears to be the predominance of very light rain in the area of interest during the winter. Both algorithms are in good spatial agreement with the radar when the radar data are restricted to rates above 1 mm/h. When all radar rain rates are included, the radar areal coverage increases by as much as a factor of 10 in some cases. Because the Kummerow algorithm does not handle such low rain rates over land very well, and because the Goddard scattering algorithm uses 1 mm/h as the minimum reliably detectable rain rate, regimes that contain large areas of very light rain present inherent difficulties for these retrieval methods. Therefore, the proliferation of low rain rates observed during the experiment is the main contributor to low correlation coefficients and high root-mean-square differences. Mididentification of cold surface (e.g., snow cover) as precipitation was also a problem in several instances.

Negri, Andrew J.↗

Foreshock Langmuir Waves for Unusually Constant Solar Wind Conditions: Data and Implications for Foreshock Structure

Plasma wave data are compared with ISEE 1's position in the electron foreshock for an interval with unusually constant (but otherwise typical) solar wind magnetic field and plasma characteristics. For this period, temporal variations in the wave characteristics can be confidently separated from sweeping of the spatially varying foreshock back and forth across the spacecraft. The spacecraft's location, particularly the coordinate D(sub f) downstream from the foreshock boundary (often termed DIFF), is calculated by using three shock models and the observed solar wind magnetometer and plasma data. Scatterplots of the wave field versus D(sub f) are used to constrain viable shock models, to investigate the observed scatter in the wave fields at constant D(sub f), and to test the theoretical predictions of linear instability theory. The scatterplots confirm the abrupt onset of the foreshock waves near the upstream boundary, the narrow width in D(sub f) of the region with high fields, and the relatively slow falloff of the fields at large D(sub f), as seen in earlier studies, but with much smaller statistical scatter. The plots also show an offset of the high-field region from the foreshock boundary. It is shown that an adaptive, time-varying shock model with no free parameters, determined by the observed solar wind data and published shock crossings, is viable but that two alternative models are not. Foreshock wave studies can therefore remotely constrain the bow shock's location. The observed scatter in wave field at constant D(sub f) is shown to be real and to correspond to real temporal variations, not to unresolved changes in D(sub f). By comparing the wave data with a linear instability theory based on a published model for the electron beam it is found that the theory can account qualitatively and semiquantitatively for the abrupt onset of the waves near D(sub f) = 0, for the narrow width and offset of the high-field region, and for the decrease in wave intensity with increasing D(sub f). Quantitative differences between observations and theory remain, including large overprediction of the wave fields and the slower than predicted falloff at large D(sub f) of the wave fields. These differences, as well as the unresolved issue of the electron beam speed in the high-field region of the foreshock, are discussed. The intrinsic temporal variability of the wave fields, as well as their overprediction based on homogeneous plasma theory, are indicative of stochastic growth physics, which causes wave growth to be random and varying in sign, rather than secular.

Cairns, Iver H.↗

Interactive Visualization of High-Dimensional Petascale Ocean Data

We describe an application for interactive visualization of 5 petabytes of time-varying multivariate data from a high-resolution global ocean circulation model. The input data are 10311 hourly (ocean time) time steps of various 2D and 3D fields from a 22-billion point 1/48- degree “lat-lon cap” configuration of the MIT General Circulation Model (MITgcm). We map the global horizontal model domain onto our 128-screen (8x16) tiled display wall to produce a canonical tiling with approximately one MITgcm grid point per display pixel, and using this tiling we encode the entire time series for multiple native and computed scalar quantities at a collection of ocean depths. We reduce disk bandwidth requirements by converting the model’s floating point data to 16-bit fixed point values, and compressing those values with a lossless video encoder, which together allow synchronized playback at 24 time steps per second across all 128 displays. The application allows dynamic assignment of any two encoded tiles to any display, and has multiple interfaces for quickly specifying various orderly arrangements of tiles. All subsequent rendering is done on the fly, with run time control of colormaps, transfer functions, histogram equalization, and labeling. The two data streams on each screen can be rendered independently and combined in various ways, including blending, differencing, horizontal/ vertical wipes, and checkerboarding. The two data streams on any screen can optionally be displayed as a scatterplot in their joint attribute space. All scatterplots and map-view plots from the same x/y location and depth are linked so they all show the current brushable selection. Ocean scientists have used the system, and have found previously unidentified features in the data.

Interactive↗

Performance evaluation and geologic utility of LANDSAT-4 TM and MSS data

Radiometric calibration accuracy of TM data, radiometric comparison of A-, B-, and P-format data, and geometric registration accuracy of the TM data at enlarged scales were analyzed. Radiometric analysis of the Wind River Basin, Wyoming scene demonstrates that the TM system can be used to extract image reflectance spectra from ground targets following calibration of the system. It was also demonstrated that: (1) image DN (radiance values) vs. ground reflectance calibration scatterplots yield parameters which can be used to constrain atmospheric models and can determine TM radiometric sensitivity; (2) no significant degradation occurs as a result of radiometric and geometric correction by SCROUNGE processing; and (3) TM data can be enlarged to 1:24000 with no major geometric distortions or misregistration problems to USGS topographic maps.

Paylor, E. D.↗

Preliminary Geologic/spectral Analysis of LANDSAT-4 Thematic Mapper Data, Wind River/bighorn Basin Area, Wyoming

A LIDQA evaluation for geologic applications of a LANDSAT TM scene covering the Wind River/Bighorn Basin area, Wyoming, is examined. This involves a quantitative assessment of data quality including spatial and spectral characteristics. Analysis is concentrated on the 6 visible, near infrared, and short wavelength infrared bands. Preliminary analysis demonstrates that: (1) principal component images derived from the correlation matrix provide the most useful geologic information. To extract surface spectral reflectance, the TM radiance data must be calibrated. Scatterplots demonstrate that TM data can be calibrated and sensor response is essentially linear. Low instrumental offset and gain settings result in spectral data that do not utilize the full dynamic range of the TM system.

Lang, H. R.↗

A study of rain estimation methods from space using dual-wavelength radar measurements at near-nadir incidence over ocean

A question arising from the recent interest in spaceborne weather radar is what methods can be used to estimate precipitation parameters from space. In this paper, dual-wavelength airborne radar data obtained from flights conducted during 1988 and 1989 are used to compare rain rates derived from backscattering and attenuation methods. To help interpret the results the surface reference methods are studied by means of scatterplots of the surface cross sections at the two frequencies under rain and no-rain conditions. Approximate criteria are given on combining attenuation and backscattering methods to increase the effective dynamic range of the radar. The dual-wavelength capability of the radar is also used to examine the vertical structure of the precipitation. Another factor affecting the accuracy of the methods is the drop-size distribution. In the final section of the paper a procedure to estimate the profiled drop-size distribution is applied to the measured radar data.

Meneghini, R.↗

Duration versus brightness of gamma-ray bursts: Comparisons between SIGNE and BATSE

We analyze duration and brightness distributions of both the SIGNE Venera 13 and 14 and Burst and Transient Source Experiment (BATSE) gamma-ray burst databases. Choosing T(sub 50) as a measure of the burst duration and using both 64 and 1024 ms peak count rates, we search for correlations between duration, peak brightness, and the ratio V identically equals C(sub 64)/C(sub 1024), proposed as a measure of variability by Lamb, Graziani, & Smith. The duration histogram for SIGNE shows a long-duration peak that is consistent with BATSE, but does not exhibit the short population, instead appearing flat below 0.6 s; the difference is presumably caused by the failure of SIGNE to detect the short, faint bursts that were observed by BATSE. Estimating the instantaneous brightness by C(sub 64), we find that SIGNE confirms the BATSE result that the long and short bursts have similar maximum instantaneous brightnesses. Scatterplots between duration, brightness, and V are consistent for both databases; we show that SIGNE confirms the BATSE observation that there is a lack of bursts that are both bright over 1024 ms and contain a short, bright spike.

Kargatis, V. E.↗

On the Bimodality of ENSO Cycle Extremes

On the basis of sea surface temperature in the El Nino 3.4 region (5 deg. N.,-5 deg. S., 120-170 deg. W.) during the interval of 1950-1997, Kevin Trenberth previously has identified some 16 El Nino and 10 La Nina, these 26 events representing the extremes of the quasi-periodic El Nino-Southern Oscillation (ENSO) cycle. Runs testing shows that the duration, recurrence period, and sequencing of these extremes vary randomly. Hence, the decade of the 1990's, especially for El Nino, is not significantly different from that of previous decadal epochs, at least, on the basis of the frequency of onsets of ENSO extremes. Additionally, the distribution of duration for both El Nino and La Nina looks strikingly bimodal, each consisting of two preferred modes, about 8- and 16-mo long for El Nino and about 9- and 18-mo long for La Nina, as does the distribution of the recurrence period for El Nino, consisting of two preferred modes about 21- and 50-mo long. Scatterplots of the recurrence period versus duration for El Nino are found to be statistically important, displaying preferential associations that link shorter (longer) duration with shorter (longer) recurrence periods. Because the last onset of El Nino occurred in April 1997 and the event was of longer than average duration, onset of the next anticipated El Nino is not expected until February 2000 or later.

Wilson, Robert M.↗

On The Bimodality of ENSO Cycle Extremes

On the basis of sea surface temperature in the El Nino 3.4 region (5N.-5S., 120-170W.) during the interval of 1950-1997, Kevin Trenberth previously has identified some 16 El Nino and 10 La Nina, these 26 events representing the extremes of the quasi-periodic El Nino-Southern Oscillation (ENSO) cycle. Runs testing shows that the duration, recurrence period, and sequencing of these extremes vary randomly. Hence, the decade of the 1990's, especially for El Nino, is not significantly different from that of previous decadal epochs, at least, on the basis of the frequency of onsets of ENSO extremes. Additionally, the distribution of duration for both El Nino and La Nina looks strikingly bimodal, each consisting of two preferred modes, about 8- and 16-months long for El Nino and about 9- and 18-months long for La Nina, as does the distribution of the recurrence period for El Nino, consisting of two preferred modes about 21- and 50- mo long. Scatterplots of the recurrence period versus duration for El Nino are found to be statistically important, displaying preferential associations that link shorter (longer) duration with shorter (longer) recurrence periods. Because the last onset of El Nino occurred in April 1997 and the event was of longer than average duration, onset of the next anticipated El Nino is not expected until February 2000 or later.

Wilson, Robert M.↗

A New, More Physically Based Algorithm, for Retrieving Aerosol Properties over Land from MODIS

The MOD Imaging Spectrometer (MODIS) has been successfully retrieving aerosol properties, beginning in early 2000 from Terra and from mid 2002 from Aqua. Over land, the retrieval algorithm makes use of three MODIS channels, in the blue, red and infrared wavelengths. As part of the validation exercises, retrieved spectral aerosol optical thickness (AOT) has been compared via scatterplots against spectral AOT measured by the global Aerosol Robotic NETwork (AERONET). On one hand, global and long term validation looks promising, with two-thirds (average plus and minus one standard deviation) of all points falling between published expected error bars. On the other hand, regression of these points shows a positive y-offset and a slope less than 1.0. For individual regions, such as along the U.S. East Coast, the offset and slope are even worse. Here, we introduce an overhaul of the algorithm for retrieving aerosol properties over land. Some well-known weaknesses in the current aerosol retrieval from MODIS include: a) rigid assumptions about the underlying surface reflectance, b) limited aerosol models to choose from, c) simplified (scalar) radiative transfer (RT) calculations used to simulate satellite observations, and d) assumption that aerosol is transparent in the infrared channel. The new algorithm attempts to address all four problems: a) The new algorithm will include surface type information, instead of fixed ratios of the reflectance in the visible channels to the mid-IR reflectance. b) It will include updated aerosol optical properties to reflect the growing aerosol retrieved from eight-plus years of AERONE". operation. c) The effects of polarization will be including using vector RT calculations. d) Most importantly, the new algorithm does not assume that aerosol is transparent in the infrared channel. It will be an inversion of reflectance observed in the three channels (blue, red, and infrared), rather than iterative single channel retrievals. Thus, this new formulation of the MODIS aerosol retrieval over land includes more physically based surface, aerosol and radiative transfer with fewer potentially erroneous assumptions.

Levy, Robert C.↗

On the Relation Between Sunspot Area and Sunspot Number

Often, the relation between monthly or yearly averages of total sunspot area, A, and sunspot number, R, has been described using the formula A = 16.7 R. Such a simple relation, however, is erroneous. The yearly ratio of A/R has varied between 5.3 in 1964 to 19.7 in 1926, having a mean of 13.1 with a standard deviation of 3.5. For 1875-1976 (corresponding to the Royal Greenwich Observatory timeframe), the yearly ratio of A/R has a mean of 14.1 with a standard deviation of 3.2, and it is found to differ significantly from the mean for 1977-2004 (corresponding to the United States Air Force/National Oceanic and Atmospheric Administration Solar Optical Observing Network timeframe), which equals 9.8 with a standard deviation of 2.1. Scatterplots of yearly values of A versus R are highly correlated for both timeframes and they suggest that a value of R = 100 implies A=1,538 +/- 174 during the first timeframe, but only A=1,076 +/- 123 for the second timeframe. Comparison of the yearly ratios adjusted for same day coverage against yearly ratios using Rome Observatory measures for the interval 1958-1998 indicates that sunspot areas during the second timeframe are inherently too low.

Wilson, Robert M.↗

GPM Pre-Launch Algorithm Development for Physically-Based Falling Snow Retrievals

In this work we compare and correlate the long time series (Nov.-March) n~easurements of precipitation rate from the Parsivels and 2DVD to the passive (89, 150, 183+/-1, +/-3, +/-7 GHz) observations of NOAA's AMSU-B radiometer. There are approximately 5-8 AMSU-B overpass views of the CARE site a day. We separate the comparisons into categories of no precipitation, liquid rain and falling snow precipitation. Scatterplots between the Parsivel snowfall rates and AMSU-B brightness temperatures (TBs) did not show an exploitable relationship for retrievals. We further compared and contrasted brightness temperatures to other surface measurements such as temperature and relative humidity with equally unsatisfying results. We found that there are similar TBs (especially at 89 and 150 GHz) for cases with falling snow and for non-precipitating cases. The comparisons indicate that surface emissivity contributions to the satellite observed TB over land can add uncertainty in detecting and estimating falling snow. The newest results show that the cloud icc scattering signal in the AMSU-B data call be detected by computing clear air TBs based on CARE radiosonde data and a rough estimate of surface emissivity. That is the differences in computed TI3 and AMSU-B TB for precipitating and nonprecipitating cases are unique such that the precipitating versus ~lon-precipitating cases can be identified. These results require that the radiosonde releases are within an hour of the AMSU-B data and allow for three surface types: no snow on the ground, less than 5 cm snow on the ground, and greater than 5 cm on the ground (as given by ground station data). Forest fraction and measured emissivities were combined to calculate the surface emissivities. The above work and future work to incorporate knowledge about falling snow retrievals into the framework of the expected GPM Bayesian retrievals will be described during this presentation.

Gail Skofronick-Jackson↗

Titan's Surface from Cassini RADAR SAR and High Resolution Radiometry Data of the First Five Flybys

The first five Titan flybys with Cassini's Synthetic Aperture RADAR (SAR) and radiometer are examined with emphasis on the calibration and interpretation of the high-resolution radiometry data acquired during the SAR mode (SAR-radiometry). Maps of the 2-cm wavelength brightness temperature are obtained coincident with the SAR swath imaging, with spatial resolution approaching 6 km. A preliminary calibration shows that brightness temperature in these maps varies from 64 to 89 K. Surface features and physical properties derived from the SAR-radiometry maps and SAR imaging are strongly correlated; in general, we find that surface features with high radar reflectivity are associated with radiometrically cold regions, while surface features with low radar reflectivity correlate with radiometrically warm regions. We examined scatterplots of the normalized radar cross-section sigma(exp o) versus brightness temperature, finding differing signatures that characterize various terrains and surface features. Implications for the physical and compositional properties of these features are discussed. The results indicate that volume scattering is important in many areas of Titan's surface, particularly Xanadu, while other areas exhibit complex brightness temperature variations consistent with variable slopes or surface material and compositional properties.

Titan↗

Estimated Probability of a Cervical Spine Injury During an ISS Mission

Introduction: The Integrated Medical Model (IMM) utilizes historical data, cohort data, and external simulations as input factors to provide estimates of crew health, resource utilization and mission outcomes. The Cervical Spine Injury Module (CSIM) is an external simulation designed to provide the IMM with parameter estimates for 1) a probability distribution function (PDF) of the incidence rate, 2) the mean incidence rate, and 3) the standard deviation associated with the mean resulting from injury/trauma of the neck. Methods: An injury mechanism based on an idealized low-velocity blunt impact to the superior posterior thorax of an ISS crewmember was used as the simulated mission environment. As a result of this impact, the cervical spine is inertially loaded from the mass of the head producing an extension-flexion motion deforming the soft tissues of the neck. A multibody biomechanical model was developed to estimate the kinematic and dynamic response of the head-neck system from a prescribed acceleration profile. Logistic regression was performed on a dataset containing AIS1 soft tissue neck injuries from rear-end automobile collisions with published Neck Injury Criterion values producing an injury transfer function (ITF). An injury event scenario (IES) was constructed such that crew 1 is moving through a primary or standard translation path transferring large volume equipment impacting stationary crew 2. The incidence rate for this IES was estimated from in-flight data and used to calculate the probability of occurrence. The uncertainty in the model input factors were estimated from representative datasets and expressed in terms of probability distributions. A Monte Carlo Method utilizing simple random sampling was employed to propagate both aleatory and epistemic uncertain factors. Scatterplots and partial correlation coefficients (PCC) were generated to determine input factor sensitivity. CSIM was developed in the SimMechanics/Simulink environment with a Monte Carlo wrapper (MATLAB) used to integrate the components of the module. Results: The probability of generating an AIS1 soft tissue neck injury from the extension/flexion motion induced by a low-velocity blunt impact to the superior posterior thorax was fitted with a lognormal PDF with mean 0.26409, standard deviation 0.11353, standard error of mean 0.00114, and 95% confidence interval [0.26186, 0.26631]. Combining the probability of an AIS1 injury with the probability of IES occurrence was fitted with a Johnson SI PDF with mean 0.02772, standard deviation 0.02012, standard error of mean 0.00020, and 95% confidence interval [0.02733, 0.02812]. The input factor sensitivity analysis in descending order was IES incidence rate, ITF regression coefficient 1, impactor initial velocity, ITF regression coefficient 2, and all others (equipment mass, crew 1 body mass, crew 2 body mass) insignificant. Verification and Validation (V&V): The IMM V&V, based upon NASA STD 7009, was implemented which included an assessment of the data sets used to build CSIM. The documentation maintained includes source code comments and a technical report. The software code and documentation is under Subversion configuration management. Kinematic validation was performed by comparing the biomechanical model output to established corridors.

Brooker, John E.↗

Comparison of Monthly IMERG Precipitation Estimates with PACRAIN Atoll Observations

Satellite-based precipitation estimates provide valuable information where surface observations are not readily available, especially over the large expanses of the ocean where in situ precipitation observations are very sparse. This study compares monthly precipitation estimates from the Integrated Multisatellite Retrievals for GPM (IMERG) with gauge observations from 37 low-lying atolls from the Pacific Rainfall Database for the period June 2000–August 2020. Over the analysis period, IMERG estimates are slightly higher than the atoll observations by 0.67% with a monthly correlation of 0.68. Seasonally, DJF shows excellent agreement with a near-zero bias, while MAM shows IMERG is low by 4.6%, and JJA is high by 1.2%. SON exhibits the worst performance, with IMERG overestimating by 6.5% compared to the atolls. The seasonal correlations are well contained in the range 0.67–0.72, with the exception of SON at 0.62. Furthermore, SON has the highest RMSE at 4.70 mm day^(−1), making it the worst season for all metrics. Scatterplots of IMERG versus atolls show IMERG, on average, is generally low for light precipitation accumulations and high for intense precipitation accumulations, with best agreement at intermediate rates. Seasonal variations exist at light and intermediate rate accumulations, but IMERG consistently overestimates at intense precipitation rates. The differences between IMERG and atolls vary over time but do not exhibit any discernible trend or dependence on atoll population. The PACRAIN atoll gauges are not wind-loss corrected, so application of an appropriate adjustment would increase the precipitation amounts compared to IMERG. These results provide useful insight to users as well as valuable information for future improvements to IMERG.

David T. Bolvin↗

Investigation on the Origins of Comets as Revealed through IR High Resolution Spectroscopy: I. Molecular Abundances

We report and analyze molecular abundances updatedin twenty comets by employingmodern data reduction procedures and molecular models. Using boxplots and scatterplots, we examine how the different molecular species are distributed among the comet population, while by means of pie charts we investigate the relative proportions of these molecular species in each comet. We compare these results with the orbital parametersof the selected targets in order to identify trends related to the dynamical history of each comet. In this way, we identify at least three chemical classes based mainly on relative abundances of CO, CH3OH, CH4, C2H6, HCN and NH3. The combination of relative abundances and orbital parameters is then compared with chemical models of planetary system formation and may offer an alternative approach to investigate the origin and evolution of the material in cometary nuclei.Our study also confirms the need to add more objects to our sample in order to improve the statistics, especially for hyper-volatiles (i.e. CH4and CO) in Jupiter Family comets.

M Lippi↗