Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

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 793 records · Page 44

Conical-Domain Model for Estimating GPS Ionospheric Delays

The conical-domain model is a computational model, now undergoing development, for estimating ionospheric delays of Global Positioning System (GPS) signals. Relative to the standard ionospheric delay model described below, the conical-domain model offers improved accuracy. In the absence of selective availability, the ionosphere is the largest source of error for single-frequency users of GPS. Because ionospheric signal delays contribute to errors in GPS position and time measurements, satellite-based augmentation systems (SBASs) have been designed to estimate these delays and broadcast corrections. Several national and international SBASs are currently in various stages of development to enhance the integrity and accuracy of GPS measurements for airline navigation. In the Wide Area Augmentation System (WAAS) of the United States, slant ionospheric delay errors and confidence bounds are derived from estimates of vertical ionospheric delay modeled on a grid at regularly spaced intervals of latitude and longitude. The estimate of vertical delay at each ionospheric grid point (IGP) is calculated from a planar fit of neighboring slant delay measurements, projected to vertical using a standard, thin-shell model of the ionosphere. Interpolation on the WAAS grid enables estimation of the vertical delay at the ionospheric pierce point (IPP) corresponding to any arbitrary measurement of a user. (The IPP of a given user s measurement is the point where the GPS signal ray path intersects a reference ionospheric height.) The product of the interpolated value and the user s thin-shell obliquity factor provides an estimate of the user s ionospheric slant delay. Two types of error that restrict the accuracy of the thin-shell model are absent in the conical domain model: (1) error due to the implicit assumption that the electron density is independent of the azimuthal angle at the IPP and (2) error arising from the slant-to-vertical conversion. At low latitudes or at mid-latitudes under disturbed conditions, the accuracy of SBAS systems based upon the thin-shell model suffers due to the presence of complex ionospheric structure, high delay values, and large electron density gradients. Interpolation on the vertical delay grid serves as an additional source of delay error. The conical-domain model permits direct computation of the user s slant delay estimate without the intervening use of a vertical delay grid. The key is to restrict each fit of GPS measurements to a spatial domain encompassing signals from only one satellite. The conical domain model is so named because each fit involves a group of GPS receivers that all receive signals from the same GPS satellite (see figure); the receiver and satellite positions define a cone, the satellite position being the vertex. A user within a given cone evaluates the delay to the satellite directly, using (1) the IPP coordinates of the line of sight to the satellite and (2) broadcast fit parameters associated with the cone. The conical-domain model partly resembles the thin-shell model in that both models reduce an inherently four-dimensional problem to two dimensions. However, unlike the thin-shell model, the conical domain model does not involve any potentially erroneous simplifying assumptions about the structure of the ionosphere. In the conical domain model, the initially four-dimensional problem becomes truly two-dimensional in the sense that once a satellite location has been specified, any signal path emanating from a satellite can be identified by only two coordinates; for example, the IPP coordinates. As a consequence, a user s slant-delay estimate converges to the correct value in the limit that the receivers converge to the user s location (or, equivalently, in the limit that the measurement IPPs converge to the user s IPP).

Sparks, Lawrence↗

Non-Targeted Effects and the Dose Response for Heavy Ion Tumorigenesis

BACKGROUND: There is no human epidemiology data available to estimate the heavy ion cancer risks experienced by astronauts in space. Studies of tumor induction in mice are a necessary step to estimate risks to astronauts. Previous experimental data can be better utilized to model dose response for heavy ion tumorigenesis and plan future low dose studies. DOSE RESPONSE MODELS: The Harderian Gland data of Alpen et al.[1-3] was re-analyzed [4] using non-linear least square regression. The data set measured the induction of Harderian gland tumors in mice by high-energy protons, helium, neon, iron, niobium and lanthanum with LET s ranging from 0.4 to 950 keV/micron. We were able to strengthen the individual ion models by combining data for all ions into a model that relates both radiation dose and LET for the ion to tumor prevalence. We compared models based on Targeted Effects (TE) to one motivated by Non-targeted Effects (NTE) that included a bystander term that increased tumor induction at low doses non-linearly. When comparing fitted models to the experimental data, we considered the adjusted R2, the Akaike Information Criteria (AIC), and the Bayesian Information Criteria (BIC) to test for Goodness of fit.In the adjusted R2test, the model with the highest R2values provides a better fit to the available data. In the AIC and BIC tests, the model with the smaller values of the summary value provides the better fit. The non-linear NTE models fit the combined data better than the TE models that are linear at low doses. We evaluated the differences in the relative biological effectiveness (RBE) and found the NTE model provides a higher RBE at low dose compared to the TE model. POWER ANALYSIS: The final NTE model estimates were used to simulate example data to consider the design of new experiments to detect NTE at low dose for validation. Power and sample sizes were calculated for a variety of radiation qualities including some not considered in the Harderian Gland data set and with different background tumor incidences. We considered different experimental designs with varying number of doses and varying low doses dependant on the LET of the radiation. The optimal design to detect a NTE for an individual ion had 4 doses equally spaced below a maximal dose where bending due to cell sterilization was < 2%. For example at 100 keV/micron we would irradiate at 0.03 Gy, 0.065 Gy, 0.13 Gy, and 0.26 Gy and require 850 mice including a control dose for a sensitivity to detect NTE with 80% power. Sample sizes could be improved by combining ions similar to the methods used with the Harderian Gland data.

Chappelli, Lori J.↗

Mission X: Train Like an Astronaut Pilot Study

Mission X: Train Like an Astronaut is an international educational challenge focusing on fitness and nutrition as we encourage students to "train like an astronaut." Teams of students (aged 8-12) learn principles of healthy eating and exercise, compete for points by finishing training modules, and get excited about their future as "fit explorers." The 18 core exercises (targeting strength, endurance, coordination, balance, spatial awareness, and more) involve the same types of skills that astronauts learn in their training and use in spaceflight. This first-of-its-kind cooperative outreach program has allowed 14 space agencies and various partner institutions to work together to address quality health/fitness education, challenge students to be more physically active, increase awareness of the importance of lifelong health and fitness, teach students how fitness plays a vital role in human performance for exploration, and inspire and motivate students to pursue careers in STEM fields. The project was initiated in 2009 in response to a request by the International Space Life Sciences Working Group. USA, Netherlands, Italy, France, Germany, Austria, Colombia, Spain, and United Kingdom hosted teams for the pilot this past spring, and Japan held a modified version of the challenge. Several more agencies provided input into the preparations. Competing on 131 teams, more than 3700 students from 40 cities worldwide participated in the first round of Mission X. OUTCOMES AND BEST PRACTICES Members of the Mission X core team will highlight the outcomes of this international educational outreach pilot project, show video highlights of the challenge, provide the working group s initial assessment of the project and discuss the future potential of the effort. The team will also discuss ideas and best practices for international partnership in education outreach efforts from various agency perspectives and experiences

Lloyd, Charles W.↗

Earth as an Extrasolar Planet: Earth Model Validation Using EPOXI Earth Observations

The EPOXI Discovery Mission of Opportunity reused the Deep Impact flyby spacecraft to obtain spatially and temporally resolved visible photometric and moderate resolution near-infrared (NIR) spectroscopic observations of Earth. These remote observations provide a rigorous validation of whole disk Earth model simulations used to better under- stand remotely detectable extrasolar planet characteristics. We have used these data to upgrade, correct, and validate the NASA Astrobiology Institute s Virtual Planetary Laboratory three-dimensional line-by-line, multiple-scattering spectral Earth model (Tinetti et al., 2006a,b). This comprehensive model now includes specular reflectance from the ocean and explicitly includes atmospheric effects such as Rayleigh scattering, gas absorption, and temperature structure. We have used this model to generate spatially and temporally resolved synthetic spectra and images of Earth for the dates of EPOXI observation. Model parameters were varied to yield an optimum fit to the data. We found that a minimum spatial resolution of approx.100 pixels on the visible disk, and four categories of water clouds, which were defined using observed cloud positions and optical thicknesses, were needed to yield acceptable fits. The validated model provides a simultaneous fit to the Earth s lightcurve, absolute brightness, and spectral data, with a root-mean-square error of typically less than 3% for the multiwavelength lightcurves, and residuals of approx.10% for the absolute brightness throughout the visible and NIR spectral range. We extend our validation into the mid-infrared by comparing the model to high spectral resolution observations of Earth from the Atmospheric Infrared Sounder, obtaining a fit with residuals of approx.7%, and temperature errors of less than 1K in the atmospheric window. For the purpose of understanding the observable characteristics of the distant Earth at arbitrary viewing geometry and observing cadence, our validated forward model can be used to simulate Earth s time dependent brightness and spectral properties for wavelengths from the far ultraviolet to the far infrared.brightness

Robinson, Tyler D.↗

Analysis of ALTAIR 1998 Meteor Radar Data

We describe a new analysis of a set of 32 UHF meteor radar traces recorded with the 422 MHz ALTAIR radar facility in November 1998. Emphasis is on the velocity measurements, and on inferences that can be drawn from them regarding the meteor masses and mass densities. We find that the velocity vs altitude data can be fitted as quadratic functions of the path integrals of the atmospheric densities vs distance, and deceleration rates derived from those fits all show the expected behavior of increasing with decreasing altitude. We also describe a computer model of the coupled processes of collisional heating, radiative cooling, evaporative cooling and ablation, and deceleration - for meteors composed of defined mixtures of mineral constituents. For each of the cases in the data set we ran the model starting with the measured initial velocity and trajectory inclination, and with various trial values of the quantity mPs 2 (the initial mass times the mass density squared), and then compared the computed deceleration vs altitude curves vs the measured ones. In this way we arrived at the best-fit values of the mPs 2 for each of the measured meteor traces. Then further, assuming various trial values of the density Ps, we compared the computed mass vs altitude curves with similar curves for the same set of meteors determined previously from the measured radar cross sections and an electrostatic scattering model. In this way we arrived at estimates of the best-fit mass densities Ps for each of the cases. Keywords meteor ALTAIR radar analysis 1 Introduction This paper describes a new analysis of a set of 422 MHz meteor scatter radar data recorded with the ALTAIR High-Power-Large-Aperture radar facility at Kwajalein Atoll on 18 November 1998. The exceptional accuracy/precision of the ALTAIR tracking data allow us to determine quite accurate meteor trajectories, velocities and deceleration rates. The measurements and velocity/deceleration data analysis are described in Sections II and III. The main point of this paper is to use these deceleration rate data, together with results from a computer model, to determine values of the quantities mPs 2 (the meteor mass times its material density squared); and further, by combining these m s 2 values with meteor mass estimates for the same set of meteors determined separately from measured radar scattering

Zinn, J.↗

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗

Fermi/GBM Observations of SGRJ0501 + 4516 Bursts

We present our temporal and spectral analyses of 29 bursts from SGRJ0501+4516, detected with the Gamma-ray Burst Monitor onboard the Fermi Gamma-ray Space Telescope during the 13 days of the source activation in 2008 (August 22 to September 3). We find that the T(sub 90) durations of the bursts can be fit with a log-normal distribution with a mean value of approx. 123 ms. We also estimate for the first time event durations of Soft Gamma Repeater (SGR) bursts in photon space (i.e., using their deconvolved spectra) and find that these are very similar to the T(sub 90)s estimated in count space (following a log-normal distribution with a mean value of approx. 124 ms). We fit the time-integrated spectra for each burst and the time-resolved spectra of the five brightest bursts with several models. We find that a single power law with an exponential cutoff model fits all 29 bursts well, while 18 of the events can also be fit with two black body functions. We expand on the physical interpretation of these two models and we compare their parameters and discuss their evolution. We show that the time-integrated and time-resolved spectra reveal that E(sub peak) decreases with energy flux (and fluence) to a minimum of approx. 30 keV at F = 8.7 x 10(exp -6)erg/sq cm/s, increasing steadily afterwards. Two more sources exhibit a similar trend: SGRs J1550 - 5418 and 1806 - 20. The isotropic luminosity, L(sub iso), corresponding to these flux values is roughly similar for all sources (0.4 - l.5 x 10(exp 40) erg/s.

Lin, Lin↗

Algorithm for Compressing Time-Series Data

An algorithm based on Chebyshev polynomials effects lossy compression of time-series data or other one-dimensional data streams (e.g., spectral data) that are arranged in blocks for sequential transmission. The algorithm was developed for use in transmitting data from spacecraft scientific instruments to Earth stations. In spite of its lossy nature, the algorithm preserves the information needed for scientific analysis. The algorithm is computationally simple, yet compresses data streams by factors much greater than two. The algorithm is not restricted to spacecraft or scientific uses: it is applicable to time-series data in general. The algorithm can also be applied to general multidimensional data that have been converted to time-series data, a typical example being image data acquired by raster scanning. However, unlike most prior image-data-compression algorithms, this algorithm neither depends on nor exploits the two-dimensional spatial correlations that are generally present in images. In order to understand the essence of this compression algorithm, it is necessary to understand that the net effect of this algorithm and the associated decompression algorithm is to approximate the original stream of data as a sequence of finite series of Chebyshev polynomials. For the purpose of this algorithm, a block of data or interval of time for which a Chebyshev polynomial series is fitted to the original data is denoted a fitting interval. Chebyshev approximation has two properties that make it particularly effective for compressing serial data streams with minimal loss of scientific information: The errors associated with a Chebyshev approximation are nearly uniformly distributed over the fitting interval (this is known in the art as the "equal error property"); and the maximum deviations of the fitted Chebyshev polynomial from the original data have the smallest possible values (this is known in the art as the "min-max property").

Hawkins, S. Edward, III↗

Analyses of Cometary Silicate Crystals: DDA Spectral Modeling of Forsterite

Comets are the Solar System's deep freezers of gases, ices, and particulates that were present in the outer protoplanetary disk. Where comet nuclei accreted was so cold that CO ice (approximately 50K) and other supervolatile ices like ethane (C2H2) were preserved. However, comets also accreted high temperature minerals: silicate crystals that either condensed (greater than or equal to 1400 K) or that were annealed from amorphous (glassy) silicates (greater than 850-1000 K). By their rarity in the interstellar medium, cometary crystalline silicates are thought to be grains that formed in the inner disk and were then radially transported out to the cold and ice-rich regimes near Neptune. The questions that comets can potentially address are: How fast, how far, and over what duration were crystals that formed in the inner disk transported out to the comet-forming region(s)? In comets, the mass fractions of silicates that are crystalline, f_cryst, translate to benchmarks for protoplanetary disk radial transport models. The infamous comet Hale-Bopp has crystalline fractions of over 55%. The values for cometary crystalline mass fractions, however, are derived assuming that the mineralogy assessed for the submicron to micron-sized portion of the size distribution represents the compositional makeup of all larger grains in the coma. Models for fitting cometary SEDs make this assumption because models can only fit the observed features with submicron to micron-sized discrete crystals. On the other hand, larger (0.1-100 micrometer radii) porous grains composed of amorphous silicates and amorphous carbon can be easily computed with mixed medium theory wherein vacuum mixed into a spherical particle mimics a porous aggregate. If crystalline silicates are mixed in, the models completely fail to match the observations. Moreover, models for a size distribution of discrete crystalline forsterite grains commonly employs the CDE computational method for ellipsoidal platelets (c:a:b=8.14x8.14xl in shape with geometrical factors of x:y:z=1:1:10, Fabian et al. 2001; Harker et al. 2007). Alternatively, models for forsterite employ statistical methods like the Distribution of Hollow Spheres (Min et al. 2008; Oliveira et al. 2011) or Gaussian Random Spheres (GRS) or RGF (Gielen et al. 200S). Pancakes, hollow spheres, or GRS shapes similar to wheat sheaf crystal habit (e.g., Volten et al. 2001; Veihelmann et al. 2006), however, do not have the sharp edges, flat faces, and vertices seen in images of cometary crystals in interplanetary dust particles (IDPs) or in Stardust samples. Cometary forsterite crystals often have equant or tabular crystal habit (J. Bradley). To simulate cometary crystals, we have computed absorption efficiencies of forsterite using the Discrete Dipole Approximation (DDA) DDSCAT code on NAS supercomputers. We compute thermal models that employ a size distribution of discrete irregularly shaped forsterite crystals (nonspherical shapes with faces and vertices) to explore how crystal shape affects the shape and wavelength positions of the forsterite spectral features and to explore whether cometary crystal shapes support either condensation or annealing scenarios (Lindsay et al. 2012a, b). We find forsterite crystal shapes that best-fit comet Hale-Bopp are tetrahedron, bricks or brick platelets, essentially equant or tabular (Lindsay et al. 2012a,b), commensurate with high temperature condensation experiments (Kobatake et al. 2008). We also have computed porous aggregates with crystal monomers and find that the crystal resonances are amplified. i.e., the crystalline fraction is lower in the aggregate than is derived by fitting a linear mix of spectral features from discrete subcomponents, and the crystal resonances 'appear' to be from larger crystals (Wooden et al. 2012). These results may indicate that the crystalline mass fraction in comets with comae dominated by aggregates may be lower than deduced by popular methods that only emoy ensembles of discrete crystals.

Wooden, Diane↗

Assessment of Mars Atmospheric Temperature Retrievals from the Thermal Emission Spectrometer Radiances

Motivated by the needs of Mars data assimilation. particularly quantification of measurement errors and generation of averaging kernels. we have evaluated atmospheric temperature retrievals from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) radiances. Multiple sets of retrievals have been considered in this study; (1) retrievals available from the Planetary Data System (PDS), (2) retrievals based on variants of the retrieval algorithm used to generate the PDS retrievals, and (3) retrievals produced using the Mars 1-Dimensional Retrieval (M1R) algorithm based on the Optimal Spectral Sampling (OSS ) forward model. The retrieved temperature profiles are compared to the MGS Radio Science (RS) temperature profiles. For the samples tested, the M1R temperature profiles can be made to agree within 2 K with the RS temperature profiles, but only after tuning the prior and error statistics. Use of a global prior that does not take into account the seasonal dependence leads errors of up 6 K. In polar samples. errors relative to the RS temperature profiles are even larger. In these samples, the PDS temperature profiles also exhibit a poor fit with RS temperatures. This fit is worse than reported in previous studies, indicating that the lack of fit is due to a bias correction to TES radiances implemented after 2004. To explain the differences between the PDS and Ml R temperatures, the algorithms are compared directly, with the OSS forward model inserted into the PDS algorithm. Factors such as the filtering parameter, the use of linear versus nonlinear constrained inversion, and the choice of the forward model, are found to contribute heavily to the differences in the temperature profiles retrieved in the polar regions, resulting in uncertainties of up to 6 K. Even outside the poles, changes in the a priori statistics result in different profile shapes which all fit the radiances within the specified error. The importance of the a priori statistics prevents reliable global retrievals based a single a priori and strongly implies that a robust science analysis must instead rely on retrievals employing localized a priori information, for example from an ensemble based data assimilation system such as the Local Ensemble Transform Kalman Filter (LETKF).

Hoffman, Matthew J.↗

Communicating the Threat of a Tropical Cyclone to the Eastern Range

The 45th Weather Squadron (45 WS) has developed a tool to help visualize the Wind Speed Probability product from the National Hurricane Center (NHC) and to help communicate that information to space launch customers and decision makers at the 45th Space Wing (45 SW) and Kennedy Space Center (KSC) located in east central Florida. This paper reviews previous work and presents the new visualization tool, including initial feedback as well as the pros and cons. The NHC began issuing their Wind Speed Probability product for tropical cyclones publicly in 2006. The 45 WS uses this product to provide a threat assessment to 45 SW and KSC leadership for risk evaluations with an approaching tropical cyclone. Although the wind speed probabilities convey the uncertainty of a tropical cyclone well, communicating this information to customers is a challenge. The 45 WS continually strives to provide the wind speed probability information to customers in a context which clearly communicates the threat of a tropical cyclone. First, an intern from the Florida Institute of Technology (FIT) Atmospheric Sciences department, sponsored by Scitor Corporation, independently evaluated the NHC wind speed probability product. This work was later extended into a M.S. thesis at FIT, partially funded by Scitor Corporation and KSC. A second thesis at FIT further extended the evaluation partially funded by KSC. Using this analysis, the 45 WS categorized the probabilities into five probability interpretation categories: Very Low, Low, Moderate, High, and Very High. These probability interpretation categories convert the forecast probability and forecast interval into easily understood categories that are consistent across all ranges of probabilities and forecast intervals. As a follow-on project, KSC funded a summer intern to evaluate the human factors of the probability interpretation categories, which ultimately refined some of the thresholds. The 45 WS created a visualization tool to express the timing and risk for multiple locations in a single graphic. Preliminary results on an on-going project by FIT will be included in this paper. This project is developing a new method of assigning the probability interpretation categories and updating the evaluation of the performance of the NHC wind speed probability analysis.

Winters, Katherine A.↗

Comparison of the Variations of Sunspot Number, Number of Sunspot Groups, and Sunspot Area, 1875 -2013

Examined are the yearly variations and ratios of sunspot number, the number of sunspot groups, and the total corrected sunspot area for the interval 1875-2013. While yearly sunspot number independently correlates strongly (r = 0.98) with the yearly number of sunspot groups (y = -2 + 11.99x) and the total corrected sunspot area (y = 5 + 0.059x), the strongest correlation (Ry12 = 0.99) is the one based on the bivariate fit of sunspot number against the combined variations of the number of sunspot groups and sunspot area (y = 1 + 5.88x1 + 0.031x2, where y refers to sunspot number, x1 refers to the number of sunspot groups, and x2 refers to the sunspot area). While all cycle minima based on the bivariate fit are concurrent with the observed minimum in sunspot number, cycle maxima are sometimes found to differ. For sunspot cycles 12, 19, 20, and 23, cycle maximum is inferred to have occurred in 1884, 1958, 1970, and 2002, respectively, rather than in 1883, 1957, 1968, and 2000, based on the observed sunspot number. Also, cycle 19's maximum amplitude based on observed sunspot number seems too high in comparison to that found using the bivariate fit. During the 139-year interval 1875-2013, the difference between the observed and predicted sunspot number based on the bivariate fit is <1 standard error of estimate (se) (<6.4) for 111 years, between 1 and <2 se (6.4 to <12.8) for 28 years, and ≥2 se (≥12.8) for only 4 years, these years being 1957 (16.6), 1978 (-15.8), 1980 (23), and 1982 (-16.3). For sunspot cycle 24, the difference between observed and predicted values has been only -0.7 and 3.2 (≤0.5 se).

Wilson, Robert M.↗

A New Global Regression Analysis Method for the Prediction of Wind Tunnel Model Weight Corrections

A new global regression analysis method is discussed that predicts wind tunnel model weight corrections for strain-gage balance loads during a wind tunnel test. The method determines corrections by combining "wind-on" model attitude measurements with least squares estimates of the model weight and center of gravity coordinates that are obtained from "wind-off" data points. The method treats the least squares fit of the model weight separate from the fit of the center of gravity coordinates. Therefore, it performs two fits of "wind- off" data points and uses the least squares estimator of the model weight as an input for the fit of the center of gravity coordinates. Explicit equations for the least squares estimators of the weight and center of gravity coordinates are derived that simplify the implementation of the method in the data system software of a wind tunnel. In addition, recommendations for sets of "wind-off" data points are made that take typical model support system constraints into account. Explicit equations of the confidence intervals on the model weight and center of gravity coordinates and two different error analyses of the model weight prediction are also discussed in the appendices of the paper.

strain-gage balance↗

Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results

We present the final nine-year maps and basic results from the Wilkinson Microwave Anisotropy Probe (WMAP) mission. The full nine-year analysis of the time-ordered data provides updated characterizations and calibrations of the experiment. We also provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate cosmic microwave background (CMB) anisotropy from foreground emission, and both types of signals are analyzed in detail.We provide new point source catalogs as well as new diffuse and point source foreground masks. An updated template-removal process is used for cosmological analysis; new foreground fits are performed, and new foreground reduced are presented.We nowimplement an optimal C(exp -1)1 weighting to compute the temperature angular power spectrum. The WMAP mission has resulted in a highly constrained Lambda-CDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that big bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (N(sub eff) = 3.84 +/- 0.40). The model fit also implies that the age of the universe is (sub 0) = 13.772 +/- 0.059 Gyr, and the fit Hubble constant is H(sub 0) = 69.32 +/- 0.80 km/s/ Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity, reported earlier by the WMAP team, now has high statistical significance (n(sub s) = 0.9608+/-0.0080); and the universe is close to flat/Euclidean (Omega = −0.0027+0.0039/−0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter Lambda-CDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further.With no significant anomalies and an adequate goodness of fit, the inflationary flat Lambda-CDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

dark matter↗

Ultrasonic Measurement and Monitoring of Loads in Bolts used in Structural Joints

The paper is an overview of work by the author in measuring and monitoring loads in bolts using an ultrasonic extensometer. A number of cases of bolted joints are covered. These include, a clamped joint with clearance fit between the bolt and hole, a clamped joint with the bolt in an interference fit with the hole, a flanged joint which allows the flange and bolt to bend; and a shear joint in a clevis and tang configuration. These applications were initially developed for measuring and monitoring preload in the NASA Space Shuttle Orbiter critical joints but are also applicable for monitoring loads in other critical bolted joints of structures such as transportation bridges and other aerospace structures. The paper explains how to set‐up a model to estimate the load factor and accuracy for the ultrasonic preload application in a clamped joint with clearance fit. The ultrasonic preload application for clamped joint with bolt in an interference fit can also be used to measure diametrical interference between the bolt shank and hole; and interference pressure on the bolt shank. Model and experimental data are given to demonstrate use of ultrasonic measurements in a shear joint. A bolt in a flanged joint experiences both tensile and bending loads. This application involves measurement of bending and tensile preload in a bolt. The ultrasonic beam bends due to the bending load on the bolt. A numerical technique to compute the trace of ultrasonic ray is presented.

Koshti, Ajay↗

Quantification of ln-Flight Physical Changes: Anthropometry and Neutral Body Posture

Currently, NASA does not have sufficient in-flight anthropometric data to assess the impact of changes in body shape and size. For developing future planetary and reduced-gravity suits, NASA needs to quantify the impacts of microgravity on anthropometry and body posture to ensure optimal crew performance, fit, and comfort. To obtain data on these changes, circumference, length, height, breadth, and depth for body segments (chest, waist, bicep, thigh, calf) from astronauts for preflight, in-flight, and post-flight conditions needs to be collected. Once these data have been collected, pre-flight, in-flight, and post-flight anthropometric values will be compared, yielding microgravity factors. The neutral body posture (NBP) will also be measured, to determine body posture (joint angle) changes between subjects throughout the duration of a mission. Data collection, starting with Increments 37/38, is still in progress but has been completed for 6 out of 9 subjects. NASA suit engineers and NASA's Extravehicular Activity (EVA) Project Office have identified that suit fit in microgravity could become an issue. It has been noted that crewmembers often need to adjust their suit sizing once they are in orbit. This adjustment could be due to microgravity effects on anthropometry and postural changes, and is necessary to ensure optimal crew performance, fit, and comfort in space. To date, the only data collected to determine the effects of microgravity on physical human changes were collected during Skylab 4, the Apollo-Soyuz Test Project (ASTP), Space Shuttle mission STS-57, and a recent HRP study on seated height changes due to spinal elongation (Spinal Elongation, Master Task List [MTL] #221). The Skylab 4, ASTP, and the STS-57 studies found that, according to photographs, a distinct NBP exists. The still photographs showed a distinguishable posture with the arms raised and the shoulders abducted; in addition, the knees are flexed, with noticeable hip flexion, and the foot is plantar flexed [1,2]. This combination is considered to be the standard set of body joint angles for an NBP in microgravity. A recent simulated microgravity NBP study [3] showed individual variability and inconsistencies in defining NBP. This variation may be influenced by spinal growth, and other potential anthropometry factors such as spinal curvature, age, and gender. Data on the variation of this posture data is required for all kinds of space device designs (such as suits, habitat, and mobility aids). The method proposed in this study considers the dynamic nature of body movement and will use a measurement technique to continually monitor posture and develop a probability likelihood of the neutral posture and how the NBP postures are affected by anthropometry. Additionally, Skylab studies found that crewmembers experienced a stature growth of up to 3 percent. The data included 3 crewmembers and showed that a biphasic stature growth occurs once the crewmember enters into weightlessness. However, the HRP Spinal Elongation study showed that crewmembers could experience about a 6 percent growth in seated height and a 3 percent stature growth, when exposed to microgravity. The results of that study prove that not all anthropometric measurements have the same microgravity percent growth factor. For EVA and suit engineers to properly update the sizing protocol for microgravity, they need additional anthropometric data from space missions. Hence, this study is aimed at gathering additional in-flight anthropometric measurements, such as length, depth, breadth, and circumference, to determine the changes to body shape and size caused by microgravity effects. It is anticipated that by recording the potential changes to body shape and size, NASA will develop a better suit sizing protocol for the International Space Station and other space missions. In essence, this study will help NASA quantify the impacts of microgravity on anthropometry to ensure optimal crew performance, fit, and comfort. This study will use standard anthropometry data collection techniques, 3D laser scanning, digital still photography, and video data, and perform photogrammetric analyses to determine the changes that occur to the body shape and size, and to NBP, when the human body is exposed to a microgravity environment.

Young, K. S.↗

Mission X: Train Like an Astronaut Challenge

The Mission X: Train Like an Astronaut Challenge was developed in 2011 to encourage proper exercise and nutrition at an early age by teaching young people to live and eat like space explorers. The strong correlation between an unhealthy childhood diet and adolescent fitness, and the onset of chronic diseases as an adult is the catalyst for Mission X. Mission X is dedicated to assisting people on a global scale to live healthier lifestyles and learn about human space exploration. The Mission X: Train Like an Astronaut 2015 (MX15) International Challenge hosted almost 40,000 children on 800 teams, 28 countries affiliated with 12 space agencies. The MX15 website included 17 languages. MX15, the fifth annual international fitness challenges sponsored by the NASA Human Research Program worked with the European Space Agency and other space agencies from around the world. In comparison to MX14, MX15 expanded to include four additional new countries, increased the number of students by approximately 68% and the number of teams by 29%. Chile' and South Korea participated in the new fall Astro Charlie Walk Around the Earth Challenge. Pre-challenge training materials were made more readily available from the website. South Korea completed a prospective assessment of the usability of the MX content for improving health and fitness in 212 preschool children and their families. Mission X is fortunate to have the support of the NASA, ESA and JAXA astronaut corps. In MX15, they participated in the opening and closing events as well as while on-board the International Space Station. Italian Astronaut Samantha Cristoretti participated as the MX15 Astronaut Ambassador for health and fitness providing the opening video and other videos from ISS. United Kingdom Astronaut Tim Peake and US Astronaut Kate Rubins have agreed to be the MX Ambassadors for 2016 and 2017 respectively. The MX15 International Working Group Face-to-Face meeting and Closing Event were held at the Agenzia Spaziale Italiana (ASI) in Rome, Italy. A record number of twenty-eight countries participated. Austria and Norway have offered to host the 2016 and 2017 working group meetings. MX16 planning began with the working group meetings and areas of improvement will include another second early challenge to accommodate countries in the Southern Hemisphere, recommended changes to the MX website, development of a more defined approach to metrics, a change to the format of future MX International Working Group meetings, and proposed new activities to be developed by the MX International Educator Working Group. We look forward to welcoming many new participants in 2016!

Lloyd, Charles W.↗

Vertical Structure and Optical Properties of Titans Aerosols from Radiance Measurements Made Inside and Outside the Atmosphere

Prompted by the detection of stratospheric cloud layers by Cassini's Composite Infrared Spectrometer (CIRS; see Anderson, C.M., Samuelson, R.E. [2011]. Icarus 212, 762-778), we have re-examined the observations made by the Descent Imager/Spectral Radiometer (DISR) in the atmosphere of Titan together with two constraints from measurements made outside the atmosphere. No evidence of thin layers (<1 km) in the DISR image data sets is seen beyond the three previously reported layers at 21 km, 11 km, and 7 km by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M.G. [2009]. Icarus 199, 442-448). On the other hand, there is evidence of a thicker layer centered at about 55 km. A rise in radiance gradients in the Downward-Looking Visible Spectrometer (DLVS) data below 55 km indicates an increase in the volume extinction coefficient near this altitude. To fit the geometric albedo measured from outside the atmosphere the decrease in the single scattering albedo of Titan's aerosols at high altitudes, noted in earlier studies of DISR data, must continue to much higher altitudes. The altitude of Titan's limb as a function of wavelength requires that the scale height of the aerosols decrease with altitude from the 65 km value seen in the DISR observations below 140 km to the 45 km value at higher altitudes. We compared the variation of radiance with nadir angle observed in the DISR images to improve our aerosol model. Our new aerosol model fits the altitude and wavelength variations of the observations at small and intermediate nadir angles but not for large nadir angles, indicating an effect that is not reproduced by our radiative transfer model. The volume extinction profiles are modeled by continuous functions except near the enhancement level near 55 km altitude. The wavelength dependence of the extinction optical depth is similar to earlier results at wavelengths from 500 to 700 nm, but is smaller at shorter wavelengths and larger toward longer wavelengths. A Hapke-like model is used for the ground reflectivity, and the variation of the Hapke single scattering albedo with wavelength is given. Fits to the visible spectrometers looking upward and downward are achieved except in the methane bands longward of 720 nm. This is possibly due to uncertainties in extrapolation of laboratory measurements from 1 km-am paths to much longer paths at lower pressures. It could also be due to changes in the single scattering phase functions at low altitudes, which strongly affect the path length through methane that the photons travel. We demonstrate the effects on the model fits by varying each model parameter individually in order to illustrate the sensitivity of our determination of each model parameter.

Titan↗