Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regression models”

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 811 records · Page 45

Statistical Issues in Galaxy Cluster Cosmology

The number and growth of massive galaxy clusters are sensitive probes of cosmological structure formation. Surveys at various wavelengths can detect clusters to high redshift, but the fact that cluster mass is not directly observable complicates matters, requiring us to simultaneously constrain scaling relations of observable signals with mass. The problem can be cast as one of regression, in which the data set is truncated, the (cosmology-dependent) underlying population must be modeled, and strong, complex correlations between measurements often exist. Simulations of cosmological structure formation provide a robust prediction for the number of clusters in the Universe as a function of mass and redshift (the mass function), but they cannot reliably predict the observables used to detect clusters in sky surveys (e.g. X-ray luminosity). Consequently, observers must constrain observable-mass scaling relations using additional data, and use the scaling relation model in conjunction with the mass function to predict the number of clusters as a function of redshift and luminosity.

Galaxy↗

Monitoring Springs in the Mojave Desert Using Landsat Time Series Analysis

The purpose of this study, based on Landsat satellite data was to characterize variations and trends over 30 consecutive years (1985-2016) in perennial vegetation green cover at over 400 confirmed Mojave Desert spring locations. These springs were surveyed between in 2015 and 2016 on lands managed in California by the U.S. Bureau of Land Management (BLM) and on several land trusts within the Barstow, Needles, and Ridgecrest BLM Field Offices. The normalized difference vegetation index (NDVI) from July Landsat images was computed at each spring location and a trend model was first fit to the multi-year NDVI time series using least squares linear regression.Â

Time Series↗

Landslide Likelihood Prediction using Machine Learning Algorithms

The supply of electricity via power plants is criticalto the operation of many critical infrastructure systems in mod-ern society. Natural hazards can disrupt the power supply, causepower outages that can halt economic growth, and impede emer-gency response until power is restored. The proposed work aimsto predict the landslides likelihood in these critical infrastructurelocations in the Northeastern USA using integrated databases ofexplanatory variables and machine learning algorithms. First,data related to landslides are obtained and merged, includingtopographic, soil moisture, and precipitation-related data. Fiveregression algorithms, namely: Random Forest, Extreme Gradi-ent Boosting (XGBoost), K-Nearest Neighbor regression (KNN),Linear Support Vector Regressor (SVR), and Linear regression,are utilized to predict the landslide probability and evaluatedon the dataset. The accuracy of the models is assessed by usingstatistical metrics such as mean absolute error (MAE), meansquared error (MSE), and root mean squared error (RMSE).The study results show that Random Forest outperformed othermodels with the mutual information feature selection method.It achieved an MSE of 0.0011 with mutual information-basedfeature selection and an MSE of 0.00157 without feature selection.KNN regressor outperformed the other models with an MSEof 0.00139 with correlation-based information selection. Theproposed landslide identification model with Random Forestalgorithm shows outstanding robustness and great potential intackling the landslide likelihood prediction by employing MLalgorithms.

Vasundhara Acharya↗

Improved GEO Derived LW Broadband Fluxes for the Edition5 CERES SYN1Deg Data Product

The NASA Clouds and the Earth's Radiant Energy System (CERES) project provides over 24 years of accurately observed top-of-the-atmosphere (TOA) and surface flux data record for climate monitoring and diagnostic studies. The major level-3 data products include SSF1deg, SYN1deg, FluxByCldTyp, and EBAF data. The CERES SYN1deg product provides 1° gridded observed TOA fluxes and computed surface fluxes based on hourly GEO clouds and derived broadband TOA fluxes. The hourly GEO derived fluxes fill in the temporal gaps between the Terra and Aqua CERES observations. In order to maintain a climate quality dataset the GEO cloud and flux retrievals are consistent over the record. The CERES project is preparing for Edition 5 product providing an opportunity to improve the GEO flux retrieval algorithms. To achieve this, the GEO satellite channel radiances are calibrated against the Aqua-MODIS calibration reference. For SW the GEO channel radiances are then converted to broadband radiances using empirical and theoretical models. The SW broadband radiances are then converted to fluxes using the same CERES angular directional models. For LW the GEO channel radiances are directly converted to LW fluxes utilizing empirical models. The resulting GEO derived broadband fluxes are then normalized to the CERES fluxes regionally by regressing the instantaneous coincident flux pairs. The normalized GEO fluxes are then tied to the CERES instrument calibration, which is very stable over time. This study will attempt to improve the GEO derived LW fluxes based on machine learning technique for the future Edition5 SYN1deg data product.

Moguo Sun↗

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.↗

Satellite-Based Spatiotemporal Trends in PM2.5 Concentrations: China 2004-2013

Three decades of rapid economic development is causing severe and widespread PM2.5(particulate matter (is) less than 2.5 ) pollution in China. However, research on the health impacts of PM2.5 exposure has been hindered by limited historical PM2.5 concentration data. We estimated ambient PM2.5 concentrations from 2004 to 2013 in China at 0.1 deg resolution using the most recent satellite data and evaluated model performance with available ground observations. We developed a two-stage spatial statistical model using the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 aerosol optical depth (AOD) and assimilated meteorology, land use data, and PM2.5 concentrations from China's recently established ground monitoring network. An inverse variance weighting (IVW) approach was developed to combine MODIS Dark Target and Deep Blue AOD to optimize data coverage. We evaluated model predicted PM2.5 concentrations from 2004 to early 2014 using ground observations. The overall model cross-validation R(sup 2) and relative prediction error were 0.79 and 35.6%, respectively. Validation beyond the model year (2013) indicated that it accurately predicted PM(sub 2.5) concentrations with little bias at the monthly (R(sup 2) = 0.73), regression slope = 0.91) and seasonal (R(sup 2) = 0.79), regression slope = 0.92) levels. Seasonal variations revealed that winter was the most polluted season and that summer was the cleanest season. Analysis of predicted PM2.5 levels showed a mean annual increase of 1.97 micro-g/cu cm between 2004 and 2007 and a decrease of 0.46 micro-g/cu cm between 2008 and 2013. Our satellite-driven model can provide reliable historical PM2.5 estimates in China at a resolution comparable to those used in epidemiologic studies on the health effects of long-term PM2.5 exposure in North America. This data source can potentially advance research on PM2.5 health effects in China.

satellite remote sensing↗

Crop status evaluations and yield predictions

A model was developed for predicting the day 50 percent of the wheat crop is planted in North Dakota. This model incorporates location as an independent variable. The Julian date when 50 percent of the crop was planted for the nine divisions of North Dakota for seven years was regressed on the 49 variables through the step-down multiple regression procedure. This procedure begins with all of the independent variables and sequentially removes variables that are below a predetermined level of significance after each step. The prediction equation was tested on daily data. The accuracy of the model is considered satisfactory for finding the historic dates on which to initiate yield prediction model. Growth prediction models were also developed for spring wheat.

Haun, J. R.↗

Searching for Quasi-periodic Oscillations in Astrophysical Transients Using Gaussian Processes

Analyses of quasi-periodic oscillations(QPOs)are important to understanding the dynamic behavior in manyastrophysical objects during transient events like gamma-ray bursts, solarflares, magnetarflares, and fast radiobursts. Astrophysicists often search for QPOs with frequency-domain methods such as(Lomb–Scargle)periodograms, which generally assume power-law models plus some excess around the QPO frequency. Time-series data can alternatively be investigated directly in the time domain using Gaussian process(GP)regression.While GP regression is computationally expensive in the general case, the properties of astrophysical data andmodels allow fast likelihood strategies. Heteroscedasticity and nonstationarity in data have been shown to causebias in periodogram-based analyses. GPs can take account of these properties. Using GPs, we model QPOs as astochastic process on top of a deterministicflare shape. Using Bayesian inference, we demonstrate how to infer GPhyperparameters and assign them physical meaning, such as the QPO frequency. We also perform model selectionbetween QPOs and alternative models such as red noise and show that this can be used to reliablyfind QPOs. Thismethod is easily applicable to a variety of different astrophysical data sets. We demonstrate the use of this methodon a range of short transients: a gamma-ray burst, a magnetarflare, a magnetar giantflare, and simulated solarflare data.

Moritz Hubner↗

ENSO Effect on East Asian Tropical Cyclone Landfall via Changes in Tracks and Genesis in a Statistical Model

Improvements on a statistical tropical cyclone (TC) track model in the western North Pacific Ocean are described. The goal of the model is to study the effect of El Nino-Southern Oscillation (ENSO) on East Asian TC landfall. The model is based on the International Best-Track Archive for Climate Stewardship (IBTrACS) database of TC observations for 1945-2007 and employs local regression of TC formation rates and track increments on the Nino-3.4 index and seasonally varying climate parameters. The main improvements are the inclusion of ENSO dependence in the track propagation and accounting for seasonality in both genesis and tracks. A comparison of simulations of the 1945-2007 period with observations concludes that the model updates improve the skill of this model in simulating TCs. Changes in TC genesis and tracks are analyzed separately and cumulatively in simulations of stationary extreme ENSO states. ENSO effects on regional (100-km scale) landfall are attributed to changes in genesis and tracks. The effect of ENSO on genesis is predominantly a shift in genesis location from the southeast in El Nino years to the northwest in La Nina years, resulting in higher landfall rates for the East Asian coast during La Nina. The effect of ENSO on track propagation varies seasonally and spatially. In the peak activity season (July-October), there are significant changes in mean tracks with ENSO. Landfall-rate changes from genesis- and track-ENSO effects in the Philippines cancel out, while coastal segments of Vietnam, China, the Korean Peninsula, and Japan show enhanced La Nina-year increases.

simulation↗

Satellite Mapping of PM2.5 Episodes in the Wintertime San Joaquin Valley: A "Static" Model Using Column Water Vapor

The use of satellite Aerosol Optical Thickness (AOT) from imaging spectrometers has been successful in quantifying and mapping high PM2.5 (particulate matter mass <2.5µm diameter) episodes for pollution abatement and health studies. However, some regions have high PM2.5 but poor estimation success. The challenges in using Aerosol Optical Thickness (AOT) from imaging spectrometers to characterize PM2.5 worldwide was especially evident in the wintertime San Joaquin Valley (SJV). The SJV's attendant difficulties of high-albedo surfaces and very shallow, variable vertical mixing also occur in other significantly polluted regions around the world. We report on more accurate PM2.5 maps for the whole-winter period in the SJV, Nov 14, 2012?Dec 11, 2013. Intensive measurements by including NASA aircraft were made for several weeks in that winter, the DISCOVER-AQ California mission.We found success with a relatively simple method based on calibration and checking with surface monitors and a characterization of vertical mixing, and incorporating specific understandings of the region's climatology. We estimate PM2.5 to within ~7µgm?3 RMSE and with R values of ~0.9, based on remotely sensed MAIAC (Multi-Angle Implementation of Atmospheric Correction) observations, and that certain further work will improve that accuracy. Mapping is at 1km resolution. This allows a time sequence of mapped aerosols at 1km for cloud-free days. We describe our technique as a "static estimation". Estimation procedures like this one, not dependent on well-mapped source strengths or on transport error, should help full source-driven simulations by deconstructing processes. They also provide a rapid method to create a long-term climatology.Essential features of the technique are (a) daily calibration of the AOT to PM2.5 using available surface monitors, and (b) characterization of mixed-layer dilution using column water vapor (CWV, otherwise "precipitable water"). We noted that on multi-day timescales both water vapor and particles share near-surface sources and both fall to very low values with altitude; indeed, both are largely removed by precipitation. The existence of layers of H2O or aerosol not within the mixed layer adds complexity, but mixed-effects statistical regression captures essential proportionality of PM2.5 and the ratio variable (AOT/CWV). Accuracy is much higher than previous statistical models, and can be extended to the whole Aqua-satellite data record. The maps and time-series we show suggest a repeated pattern for large valleys like the SJV ? progressive stabilization of the mixing height after frontal passages: PM2.5 is somewhat more determined by day-by-day changes in mixing than it is by the progressive accumulation of pollutants (revealed as increasing AOT).

Chatfield, Robert B.↗

Comparison of two atmospheric correction models for a vegetated Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) scene

Current atmospheric correction models applied to imaging spectroscopy data include such methods as residual (scene average) and flat field correction, regression method, and the LOWTRAN 7 method. Due to the limitations of using residual and flat field corrections on vegetated scenes, regression and LOWTRAN 7 are compared. Field measured targets taken at the time of the 13 April, 1989 AVIRIS overflight of Jasper Ridge, California (U.S.) were used to formulate the regression atmospheric correction. Assuming the regressed image represents ground truth, results show that the LOWTRAN 7 method with radiosonde data does not compensate as well for atmospheric water vapor as the regression method, but it may be easier to obtain a posteriori information to perform the LOWTRAN 7 atmospheric correction.

Vandenbosch, Jeannette Marie↗

Statistical Technique for Intermediate and Long-Range Estimation of 13-Month Smoothed Solar Flux and Geomagnetic Index

This report documents the Marshall Space Flight Center (MSFC) 13-month smoothed solar flux (F(sub 10.7)) and geomagnetic index (A(sub p)) intermediate (months) and long-range (years) statistical estimation technique, referred to as the MSFC Lagrangian Linear Regression Technique (MLLRT). Estimates of future solar activity are needed as updated input to upper atmosphere density models used for satellite and spacecraft orbital lifetime predictions. An assessment of the MLLRT computer program's products is provided for 5-year periods from the date estimates were made. This was accomplished for a number of past solar cycles.

Niehuss, K. O.↗

Attribution of Recovery in Lower-Stratospheric Ozone

Multiple satellite and ground-based observations provide consistent evidence that the thickness of Earth's protective ozone layer has stopped declining since 1997, close to the time of peak stratospheric halogen loading. Regression analyses with Effective Equivalent Stratospheric Chlorine (EESC) in conjunction with further analyses using more sophisticated photochemical model calculations constrained by satellite data demonstrate that the cessation of ozone depletion between 18-25 km altitude is consistent with a leveling off of stratospheric abundances of chlorine and bromine, due to the Montreal Protocol and its amendments. However, ozone increases in the lowest part of the stratosphere, from the tropopause to 18 km, account for about half of the improvement in total column ozone during the past 9 years at northern hemisphere mid-latitudes. The increase in ozone for altitudes below 18 km is most likely driven by changes in transport, rather than driven by declining chlorine and bromine. Even with this evidence that the Montreal Protocol and its amendments are having the desired, positive effect on ozone above 18 km, total column ozone is recovering faster than expected due to the apparent transport driven changes at lower altitudes. Accurate prediction of future levels of stratospheric ozone will require comprehensive understanding of the factors that drive temporal changes at various altitudes, and partitioning of the recent transport-driven increases between natural variability and changes in atmospheric structure perhaps related to anthropogenic climate change.

Yang, Eun-Su↗

Attribution of Recovery in Lower-stratospheric Ozone

Multiple satellite and ground-based observations provide consistent evidence that the thickness of Earth's protective ozone layer has stopped declining since 1997, close to the time of peak stratospheric halogen loading. Regression analyses with Effective Equivalent Stratospheric Chlorine (EESC) in conjunction with further analyses using more sophisticated photochemical model calculations constrained by satellite data demonstrate that the cessation of ozone depletion between 18-25 km altitude is consistent with a leveling off of stratospheric abundances of chlorine and bromine, due to the Montreal Protocol and its amendments. However, ozone increases in the lowest part of the stratosphere, from the tropopause to 18 km, account for about half of the improvement in total column ozone during the past 9 years at northern hemisphere mid-latitudes. The increase in ozone for altitudes below 18 km is most likely driven by changes in transport, rather than driven by declining chlorine and bromine. Even with this evidence that the Montreal Protocol and its amendments are having the desired, positive effect on ozone above 18 km, total column ozone is recovering faster than expected due to the apparent transport driven changes at lower altitudes. Accurate prediction of future levels of stratospheric ozone will require comprehensive understanding of the factors that drive temporal changes at various altitudes, and partitioning of the recent transport-driven increases between natural variability and changes in atmospheric structure perhaps related to anthropogenic climate change.

Yang, Eun-Su↗

Novel Analog For Muscle Deconditioning

Existing models of muscle deconditioning are cumbersome and expensive (ex: bedrest). We propose a new model utilizing a weighted suit to manipulate strength, power or endurance (function) relative to body weight (BW). Methods: 20 subjects performed 7 occupational astronaut tasks while wearing a suit weighted with 0-120% of BW. Models of the full relationship between muscle function/BW and task completion time were developed using fractional polynomial regression and verified by the addition of pre- and post-flight astronaut performance data using the same tasks. Spline regression was used to identify muscle function thresholds below which task performance was impaired. Results: Thresholds of performance decline were identified for each task. Seated egress & walk (most difficult task) showed thresholds of: leg press (LP) isometric peak force/BW of 18 N/kg, LP power/BW of 18 W/kg, LP work/ BW of 79 J/kg, knee extension (KE) isokinetic/BW of 6 Nm/Kg and KE torque/BW of 1.9 Nm/kg. Conclusions: Laboratory manipulation of strength / BW has promise as an appropriate analog for spaceflight-induced loss of muscle function for predicting occupational task performance and establishing operationally relevant exercise targets.

Ploutz-Snyder, Lori↗

Novel Analog For Muscle Deconditioning

Existing models (such as bed rest) of muscle deconditioning are cumbersome and expensive. We propose a new model utilizing a weighted suit to manipulate strength, power, or endurance (function) relative to body weight (BW). Methods: 20 subjects performed 7 occupational astronaut tasks while wearing a suit weighted with 0-120% of BW. Models of the full relationship between muscle function/BW and task completion time were developed using fractional polynomial regression and verified by the addition of pre-and postflightastronaut performance data for the same tasks. Splineregression was used to identify muscle function thresholds below which task performance was impaired. Results: Thresholds of performance decline were identified for each task. Seated egress & walk (most difficult task) showed thresholds of leg press (LP) isometric peak force/BW of 18 N/kg, LP power/BW of 18 W/kg, LP work/BW of 79 J/kg, isokineticknee extension (KE)/BW of 6 Nm/kg, and KE torque/BW of 1.9 Nm/kg.Conclusions: Laboratory manipulation of relative strength has promise as an appropriate analog for spaceflight-induced loss of muscle function, for predicting occupational task performance and establishing operationally relevant strength thresholds.

Ploutz-Snyder, Lori↗

Modeled soft X-ray solar irradiances

Solar soft X-rays have historically been inaccurately modeled in both relative variations and absolute magnitudes by empirical solar extreme ultraviolet (EUV) irradiance models. This is a result of the use of a limited number of rocket data sets which were primarily associated with the calibration of the AE-E satellite EUV data set. In this work, the EUV91 solar EUV irradiance model has been upgraded to improve the accuracy of the 3.0 to 5.0 nm relative irradiance variations. The absolute magnitude estimate of the flux in this wavelength range has also been revised upwards. The upgrade was accomplished by first digitizing the SOLRAD 11 satellite 4.4 to 6.0 nm measured energy flux data set, then extracting and extrapolating a derived 3.0 to 5.0 nm photon flux from these data, and finally by performing a correlation between these derived data and the daily and 81-day mean 10.7 cm radio flux emission using a multiple linear regression technique. A correlation coefficient of greater than 0.9 was obtained between the dependent and independent data sets. The derived and modeled 3.0 to 5.0 nm flux varies by more than an order of magnitude over a solar cycle, ranging from a flux below 1 x 10(exp 8) to a flux greater than 1 x 10(exp 9) photons per sq cm per sec. Solar rotational (27-day) variations in the flux magnitude are a factor of 2. The derived and modeled irradiance absolute values are an order of magnitude greater than previous values from rocket data sets related to the calibration of the AE-E satellite.

Tobiska, W. Kent↗

Numerical Calculation of the Drag Force Acting on an Insoluble Particle Moving in Front of a Solidifying Interface

The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. Whether a particle will be engulfed or continuously pushed by SLI is essentially determined by the balance of forces acting on the particle. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that not only the steady-state approach is inappropriate to model the interaction process but also that even at steady-state the theoretical solution underestimates the drag force. A regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gave a relationship of the form F(sub d) (exp num) = (square root of 3) pi eta V(sub p)((R(sub p) / d) (exp 10 gamma/3)) as compared to the theoretical solution F(sub d) (exp theor) = 3 X (square root of 2) pi eta V(sub p)((R(sub p) / d) (exp 3/2)) where F(sub d) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma is approximately equal to 0.577), and d is the width of the gap between the particle and the SLI. Model validation as well as the influence of the interface shape on the value of F(sub d) will also be discussed.

Catalina, Adrian V.↗