Search NASA⌕ Search

SEARCH · Search NASA

Results for “Random Uncertainty”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

A Satellite Mortality Study to Support Space Systems Lifetime Prediction

Estimating the operational lifetime of satellites and spacecraft is a complex process. Operational lifetime can differ from mission design lifetime for a variety of reasons. Unexpected mortality can occur due to human errors in design and fabrication, to human errors in launch and operations, to random anomalies of hardware and software or even satellite function degradation or technology change, leading to unrealized economic or mission return. This study focuses on data collection of public information using, for the first time, a large, publically available dataset, and preliminary analysis of satellite lifetimes, both operational lifetime and design lifetime. The objective of this study is the illustration of the relationship of design life to actual lifetime for some representative classes of satellites and spacecraft. First, a Weibull and Exponential lifetime analysis comparison is performed on the ratio of mission operating lifetime to design life, accounting for terminated and ongoing missions. Next a Kaplan-Meier survivor function, standard practice for clinical trials analysis, is estimated from operating lifetime. Bootstrap resampling is used to provide uncertainty estimates of selected survival probabilities. This study highlights the need for more detailed databases and engineering reliability models of satellite lifetime that include satellite systems and subsystems, operations procedures and environmental characteristics to support the design of complex, multi-generation, long-lived space systems in Earth orbit.

Fox, George↗
NASA NTRS2021

Spatial Variation of Fine Particulate Matter Levels in Nairobi Before and During the COVID-19 Curfew: Implications for Environmental Justice

The temporary decrease of fine particulate matter (PM(sub 2.5)) concentrations in many parts of the world due to the COVID-19 lockdown spurred discussions on urban air pollution and health. However there has been little focus on sub-Saharan Africa, as few African cities have air quality monitors and if they do, these data are often not publicly available. Spatial differentials of changes in PM(sub 2.5) concentrations as a result of COVID also remain largely unstudied. To address this gap, we use a serendipitous mobile air quality monitoring deployment of eight Sensirion SPS 30 sensors on motorbikes in the city of Nairobi starting on 16 March 2020, before a COVID-19 curfew was imposed on 25 March and continuing until 5 May 2020. We developed a random-forest model to estimate PM(sub 2.5) surfaces for the entire city of Nairobi before and during the COVID-19 curfew. The highest PM(sub 2.5) concentrations during both periods were observed in the poor neighborhoods of Kariobangi, Mathare, Umoja, and Dandora, located to the east of the city center. Changes in PM(sub 2.5) were heterogeneous over space. PM(sub 2.5) concentrations increased during the curfew in rapidly urbanizing, the lower-middle-class neighborhoods of Kahawa, Kasarani, and Ruaraka, likely because residents switched from LPG to biomass fuels due to loss of income. Our results indicate that COVID-19 and policies to address it may have exacerbated existing air pollution inequalities in the city of Nairobi. The quantitative results are preliminary, due to sampling limitations and measurement uncertainties, as the available data came exclusively from low-cost sensors. This research serves to highlight that spatial data that is essential for understanding structural inequalities reflected in uneven air pollution burdens and differential impacts of events like the COVID pandemic. With the help of carefully deployed low cost sensors with improved spatial sampling and at least one reference-quality monitor for calibration, we can collect data that is critical for developing targeted interventions that address environmental injustice in the African context.

Environmental justice↗
NASA NTRS2016

Modeling Forest Biomass and Growth: Coupling Long-Term Inventory and Lidar Data

Combining spatially-explicit long-term forest inventory and remotely sensed information from Light Detection and Ranging (LiDAR) datasets through statistical models can be a powerful tool for predicting and mapping above-ground biomass (AGB) at a range of geographic scales. We present and examine a novel modeling approach to improve prediction of AGB and estimate AGB growth using LiDAR data. The proposed model accommodates temporal misalignment between field measurements and remotely sensed data-a problem pervasive in such settings-by including multiple time-indexed measurements at plot locations to estimate AGB growth. We pursue a Bayesian modeling framework that allows for appropriately complex parameter associations and uncertainty propagation through to prediction. Specifically, we identify a space-varying coefficients model to predict and map AGB and its associated growth simultaneously. The proposed model is assessed using LiDAR data acquired from NASA Goddard's LiDAR, Hyper-spectral & Thermal imager and field inventory data from the Penobscot Experimental Forest in Bradley, Maine. The proposed model outperformed the time-invariant counterpart models in predictive performance as indicated by a substantial reduction in root mean squared error. The proposed model adequately accounts for temporal misalignment through the estimation of forest AGB growth and accommodates residual spatial dependence. Results from this analysis suggest that future AGB models informed using remotely sensed data, such as LiDAR, may be improved by adapting traditional modeling frameworks to account for temporal misalignment and spatial dependence using random effects.

Babcock, Chad↗
NASA NTRS2023

Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures

In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.

Haokui Xu↗
NASA NTRS1987

Design, calibration and error analysis of instrumentation for heat transfer measurements in internal combustion engines

The paper reports the methodology and uncertainty analyses of instrumentation for heat transfer measurements in internal combustion engines. Results are presented for determining the local wall heat flux in an internal combustion engine (using a surface thermocouple-type heat flux gage) and the apparent flame-temperature and soot volume fraction path length product in a diesel engine (using two-color pyrometry). It is shown that a surface thermocouple heat transfer gage suitably constructed and calibrated will have an accuracy of 5 to 10 percent. It is also shown that, when applying two-color pyrometry to measure the apparent flame temperature and soot volume fraction-path length, it is important to choose at least one of the two wavelengths to lie in the range of 1.3 to 2.3 micrometers. Carefully calibrated two-color pyrometer can ensure that random errors in the apparent flame temperature and in the soot volume fraction path length will remain small (within about 1 percent and 10-percent, respectively).

Ferguson, C. R.↗
NASA NTRS2000

Adjoints and Low-rank Covariance Representation

Quantitative measures of the uncertainty of Earth System estimates can be as important as the estimates themselves. Second moments of estimation errors are described by the covariance matrix, whose direct calculation is impractical when the number of degrees of freedom of the system state is large. Ensemble and reduced-state approaches to prediction and data assimilation replace full estimation error covariance matrices by low-rank approximations. The appropriateness of such approximations depends on the spectrum of the full error covariance matrix, whose calculation is also often impractical. Here we examine the situation where the error covariance is a linear transformation of a forcing error covariance. We use operator norms and adjoints to relate the appropriateness of low-rank representations to the conditioning of this transformation. The analysis is used to investigate low-rank representations of the steady-state response to random forcing of an idealized discrete-time dynamical system.

Tippett, Michael K.↗
NASA NTRS1977

The reality of comet groups and pairs

Although the common genetic origin of the Kreutz family of sun-grazing comets has generally been accepted, there remains uncertainty with regard to genetic identity among other groups of comets whose orbital elements are nearly alike. Porter (1952) has listed a number of such groups, and Opik (1971) has made a statistical study of the orbits of 472 comets with aphelion distances beyond Saturn. Opik lists 97 groups that show similarities among their three angular elements, calculates an overall probability of some 10 to the -39th power that these similarities could have occurred by chance, and thus concludes that 60% or more of such comets fall into genetic groups containing from two to seven members. This paper explores the statistical reality of Opik's groups utilizing the Monte Carlo method of statistics as well as ordinary probability theory. The conclusion is reached that except for a few pairs, the similarity among orbital elements within the groups is no greater than random expectation.

Whipple, F. L.↗
NASA NTRS2008

Supersonic/Hypersonic Laminar Heating Correlations for Rectangular and Impact-Induced Open and Closed Cavities

Impact and debris damage to the Space Shuttle Orbiter Thermal Protection System tiles is a random phenomenon, occurring at random locations on the vehicle surface, resulting in random geometrical shapes that are exposed to a definable range of surface flow conditions. In response to the 2003 Final Report of the Columbia Accident Investigation Board, wind tunnel aeroheating experiments approximating a wide range of possible damage scenarios covering both open and closed cavity flow conditions were systematically tested in hypersonic ground based facilities. These data were analyzed and engineering assessment tools for damage-induced fully-laminar heating were developed and exercised on orbit. These tools provide bounding approximations for the damaged-surface heating environment. This paper presents a further analysis of the baseline, zero-pressure-gradient, idealized, rectangular-geometry cavity heating data, yielding new laminar correlations for the floor-averaged heating, peak cavity endwall heating, and the downstream decay rate. Correlation parameters are derived in terms of cavity geometry and local flow conditions. Prediction Limit Uncertainty values are provided at the 95%, 99% and 99.9% levels of significance. Non-baseline conditions, including non-rectangular geometries and flows with known pressure gradients, are used to assess the range of applicability of the new correlations. All data variations fall within the 99% Prediction Limit Uncertainty bounds. Importantly, both open-flow and closed-flow cavity heating are combined into a single-curve parameterization of the heating predictions, and provide a concise mathematical model of the laminar cavity heating flow field with known uncertainty.

Everhart, Joel L.↗
NASA NTRS2016

Estimating the Rate of Occurrence of Renal Stones in Astronauts

Changes in urine chemistry, during and post flight, potentially increases the risk of renal stones in astronauts. Although much is known about the effects of space flight on urine chemistry, no inflight incidence of renal stones in US astronauts exists and the question "How much does this risk change with space flight?" remains difficult to accurately quantify. In this discussion, we tackle this question utilizing a combination of deterministic and probabilistic modeling that implements the physics behind free stone growth and agglomeration, speciation of urine chemistry and published observations of population renal stone incidences to estimate changes in the rate of renal stone presentation. The modeling process utilizes a Population Balance Equation based model developed in the companion IWS abstract by Kassemi et al. (2016) to evaluate the maximum growth and agglomeration potential from a specified set of urine chemistry values. Changes in renal stone occurrence rates are obtained from this model in a probabilistic simulation that interrogates the range of possible urine chemistries using Monte Carlo techniques. Subsequently, each randomly sampled urine chemistry undergoes speciation analysis using the well-established Joint Expert Speciation System (JESS) code to calculate critical values, such as ionic strength and relative supersaturation. The Kassemi model utilizes this information to predict the mean and maximum stone size. We close the assessment loop by using a transfer function that estimates the rate of stone formation from combining the relative supersaturation and both the mean and maximum free stone growth sizes. The transfer function is established by a simulation analysis which combines population stone formation rates and Poisson regression. Training this transfer function requires using the output of the aforementioned assessment steps with inputs from known non-stone-former and known stone-former urine chemistries. Established in a Monte Carlo system, the entire renal stone analysis model produces a probability distribution of the stone formation rate and an expected uncertainty in the estimate. The utility of this analysis will be demonstrated by showing the change in renal stone occurrence predicted by this method using urine chemistry distributions published in Whitson et al. 2009. A comparison to the model predictions to previous assessments of renal stone risk will be used to illustrate initial validation of the model.

gravitational physiology↗
NASA NTRS2003

Structural Reliability Using Probability Density Estimation Methods Within NESSUS

A reliability analysis studies a mathematical model of a physical system taking into account uncertainties of design variables and common results are estimations of a response density, which also implies estimations of its parameters. Some common density parameters include the mean value, the standard deviation, and specific percentile(s) of the response, which are measures of central tendency, variation, and probability regions, respectively. Reliability analyses are important since the results can lead to different designs by calculating the probability of observing safe responses in each of the proposed designs. All of this is done at the expense of added computational time as compared to a single deterministic analysis which will result in one value of the response out of many that make up the density of the response. Sampling methods, such as monte carlo (MC) and latin hypercube sampling (LHS), can be used to perform reliability analyses and can compute nonlinear response density parameters even if the response is dependent on many random variables. Hence, both methods are very robust; however, they are computationally expensive to use in the estimation of the response density parameters. Both methods are 2 of 13 stochastic methods that are contained within the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) program. NESSUS is a probabilistic finite element analysis (FEA) program that was developed through funding from NASA Glenn Research Center (GRC). It has the additional capability of being linked to other analysis programs; therefore, probabilistic fluid dynamics, fracture mechanics, and heat transfer are only a few of what is possible with this software. The LHS method is the newest addition to the stochastic methods within NESSUS. Part of this work was to enhance NESSUS with the LHS method. The new LHS module is complete, has been successfully integrated with NESSUS, and been used to study four different test cases that have been proposed by the Society of Automotive Engineers (SAE). The test cases compare different probabilistic methods within NESSUS because it is important that a user can have confidence that estimates of stochastic parameters of a response will be within an acceptable error limit. For each response, the mean, standard deviation, and 0.99 percentile, are repeatedly estimated which allows confidence statements to be made for each parameter estimated, and for each method. Thus, the ability of several stochastic methods to efficiently and accurately estimate density parameters is compared using four valid test cases. While all of the reliability methods used performed quite well, for the new LHS module within NESSUS it was found that it had a lower estimation error than MC when they were used to estimate the mean, standard deviation, and 0.99 percentile of the four different stochastic responses. Also, LHS required a smaller amount of calculations to obtain low error answers with a high amount of confidence than MC. It can therefore be stated that NESSUS is an important reliability tool that has a variety of sound probabilistic methods a user can employ and the newest LHS module is a valuable new enhancement of the program.

Chamis, Chrisos C.↗
NASA NTRS2000

Retrieving the Vertical Structure of the Effective Aerosol Complex Index of Refraction from a Combination of Aerosol in Situ and Remote Sensing Measurements During TARFOX

The largest uncertainty in estimates of the effects of atmospheric aerosols on climate stems from uncertainties in the determination of their microphysical properties, including the aerosol complex index of refraction, which in turn determines their optical properties. A novel technique is used to estimate the aerosol complex index of refraction in distinct vertical layers from a combination of aerosol in situ size distribution and remote sensing measurements during the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX). In particular, aerosol backscatter measurements using the NASA Langley LASE (Lidar Atmospheric Sensing Experiment) instrument and in situ aerosol size distribution data are utilized to derive vertical profiles of the 'effective' aerosol complex index of refraction at 815 nm (i.e., the refractive index that would provide the same backscatter signal in a forward calculation on the basis of the measured in situ particle size distributions for homogeneous, spherical aerosols). A sensitivity study shows that this method yields small errors in the retrieved aerosol refractive indices, provided the errors in the lidar derived aerosol backscatter are less than 30% and random in nature. Absolute errors in the estimated aerosol refractive indices are generally less than 0.04 for the real part and can be as much as 0.042 for the imaginary part in the case of a 30% error in the lidar-derived aerosol backscatter. The measurements of aerosol optical depth from the NASA Ames Airborne Tracking Sunphotometer (AATS-6) are successfully incorporated into the new technique and help constrain the retrieved aerosol refractive indices. An application of the technique to two TARFOX case studies yields the occurrence of vertical layers of distinct aerosol refractive indices. Values of the estimated complex aerosol refractive index range from 1.33 to 1.45 for the real part and 0.001 to 0.008 for the imaginary part. The methodology devised in this study provides, for the first time a complete set of vertically resolved aerosol size distribution and refractive index data, yielding the vertical distribution of aerosol optical properties required for the determination of aersol-induced radiative flux changes

Redemann, J.↗
NASA NTRS2000

Retrieving the Vertical Structure of the Effective Aerosol Complex Index of Refraction from a Combination of Aerosol in Situ and Remote Sensing Measurements During TARFOX

The largest uncertainty in estimates of the effects of atmospheric aerosols on climate stems from uncertainties in the determination of their microphysical properties, including the aerosol complex index of refraction, which in turn determines their optical properties. A novel technique is used to estimate the aerosol complex index of refraction in distinct vertical layers from a combination of aerosol in situ size distribution and remote sensing measurements during the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX). In particular, aerosol backscatter measurements using the NASA Langley LASE (Lidar Atmospheric Sensing Experiment) instrument and in situ aerosol size distribution data are utilized to derive vertical profiles of the "effective" aerosol complex index of refraction at 815 nm (i.e., the refractive index that would provide the same backscatter signal in a forward calculation on the basis of the measured in situ particle size distributions for homogeneous, spherical aerosols). A sensitivity study shows that this method yields small errors in the retrieved aerosol refractive indices, provided the errors in the lidar-derived aerosol backscatter are less than 30% and random in nature. Absolute errors in the estimated aerosol refractive indices are generally less than 0.04 for the real part and can be as much as 0.042 for the imaginary part in the case of a 30% error in the lidar-derived aerosol backscatter. The measurements of aerosol optical depth from the NASA Ames Airborne Tracking Sunphotometer (AATS-6) are successfully incorporated into the new technique and help constrain the retrieved aerosol refractive indices. An application of the technique to two TARFOX case studies yields the occurrence of vertical layers of distinct aerosol refractive indices. Values of the estimated complex aerosol refractive index range from 1.33 to 1.45 for the real part and 0.001 to 0.008 for the imaginary part. The methodology devised in this study provides, for the first time, a complete set of vertically resolved aerosol size distribution and refractive index data. yielding the vertical distribution of aerosol optical properties required for the determination of aerosol-induced radiative flux changes.

Redemann, J.↗
NASA NTRS1981

The current status of REH theory

A response is made to the evaluation of Fitch (1980) of REH (random evolutionary hits) theory for the evolutionary divergence of proteins and nucleic acids. Correct calculations for the beta hemoglobin mRNAs of the human, mouse and rabbit in the absence and presence of selective constraints are summarized, and it is shown that the alternative evolutionary analysis of Fitch underestimates the total fixed mutations. It is further shown that the model used by Fitch to test for the completeness of the count of total base substitutions is in fact a variant of REH theory. Considerations of the variance inherent in evolutionary estimations are also presented which show the REH model to produce no more variance than other evolutionary models. In the reply, it is argued that, despite the objections raised, REH theory applied to proteins gives inaccurate estimates of total gene substitutions. It is further contended that REH theory developed for nucleic sequences suffers from problems relating to the frequency of nucleotide substitutions, the identity of the codons accepting silent and amino acid-changing substitutions, and estimate uncertainties.

Holmquist, R.↗
NASA NTRS2014

Assimilation of MODIS Snow Cover Through the Data Assimilation Research Testbed and the Community Land Model Version 4

To improve snowpack estimates in Community Land Model version 4 (CLM4), the Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) was assimilated into the Community Land Model version 4 (CLM4) via the Data Assimilation Research Testbed (DART). The interface between CLM4 and DART is a flexible, extensible approach to land surface data assimilation. This data assimilation system has a large ensemble (80-member) atmospheric forcing that facilitates ensemble-based land data assimilation. We use 40 randomly chosen forcing members to drive 40 CLM members as a compromise between computational cost and the data assimilation performance. The localization distance, a parameter in DART, was tuned to optimize the data assimilation performance at the global scale. Snow water equivalent (SWE) and snow depth are adjusted via the ensemble adjustment Kalman filter, particularly in regions with large SCF variability. The root-mean-square error of the forecast SCF against MODIS SCF is largely reduced. In DJF (December-January-February), the discrepancy between MODIS and CLM4 is broadly ameliorated in the lower-middle latitudes (2345N). Only minimal modifications are made in the higher-middle (4566N) and high latitudes, part of which is due to the agreement between model and observation when snow cover is nearly 100. In some regions it also reveals that CLM4-modeled snow cover lacks heterogeneous features compared to MODIS. In MAM (March-April-May), adjustments to snowmove poleward mainly due to the northward movement of the snowline (i.e., where largest SCF uncertainty is and SCF assimilation has the greatest impact). The effectiveness of data assimilation also varies with vegetation types, with mixed performance over forest regions and consistently good performance over grass, which can partly be explained by the linearity of the relationship between SCF and SWE in the model ensembles. The updated snow depth was compared to the Canadian Meteorological Center (CMC) data. Differences between CMC and CLM4 are generally reduced in densely monitored regions.

data assimilation↗
NASA NTRS2017

Effects of Acute Exposures to Carbon Dioxide upon Cognitive Functions

Carbon dioxide (CO2) originates from human metabolism and typically remains about 10-fold higher in concentration on the International Space Station (ISS) than at the earth's surface. There have been recurring complaints by crew members of episodes of "mental viscosity" adversely affecting their performance, and there is evidence from the ISS that associates CO2 levels with reports of headaches by crewmembers. Consequently, flight rules have been employed to control CO2 below 3 mm Hg, which is well below the existing Spacecraft Maximum Allowable Concentration (SMAC) of 10 mm Hg for 24-hour exposures, and 5.3 mm Hg for exposures of 7 to 180 days. Headaches, while sometime debilitating themselves, are also symptoms that can provide evidence that physiological defense mechanisms have been breached, and there is evidence that CO2 has effects at levels below the threshold for headaches. This concern appears to be substantiated in reports that CO2 at concentrations below 2 mm Hg substantially reduced some cognitive functions that are associated with the ability to make complex decisions in conditions that are characterized by volatility, uncertainty, complexity, ambiguity, and delayed feedback. These are conditions that could be encountered by crews in off-nominal situations or during the first missions beyond low earth orbit. Therefore, we set out to determine if decision-making under volatile, uncertain, confusing and ambiguous circumstances, where feedback is delayed or absent, is correlated with low levels of CO2 during acute exposures (several hours) in crew-like subjects and to determine if additional cognitive domains are sensitive to concentrations of CO2 at, or below, current ISS levels by using a test battery that is currently available onboard ISS. We enrolled 22 volunteers (8 females, 14 males) between the ages of 30-55 (38.8 +/- 7.0) years whose training and professional experience reflect that of the astronaut corps. Subjects were divided among 4 study groups. Exposures occurred in the 20-foot environmental chamber in Bldg 7 at JSC. Each group was exposed to each of four concentrations of CO2 (600, 1200, 2500, and 5000 ppm), and the exposure order was balanced across the groups which were randomly assigned to 4 possible complete exposure orders. Study participants and investigators were blinded to the exposure order. Each exposure lasted 4 hours and occurred on the same day each week for 4 consecutive weeks. Testing included the Strategic Management Simulation (SMS) methodology and the Cognition battery of psychometric measures that are being utilized aboard the ISS. Subjects participated in a familiarization session one week prior to the start of exposure sessions. Each exposure session followed the schedule noted in Figure 1 below. On the morning of each exposure, subjects completed the Cognition battery in a conference room at indoor ambient CO2 levels. Subjects then entered the exposure chamber and were acclimated to the environment for 15 minutes before completing another round of the Cognition battery. After a 20 minute rest, subject completed the approx.80 minute SMS. Subjects were allowed another 20 minute rest period before completing the final Cognition test battery in the exposure chamber. Subjects then returned to the conference room where they completed one final round of Cognition. Data from both testing methodologies are currently being analyzed and will be presented. Study investigators have not yet been unblinded to the coded exposure concentrations, but preliminary results suggest differences.

Ryder, V. E.↗
NASA NTRS2019

Modelling above-ground biomass stock over Norway using national forest inventory data with ArcticDEM and Sentinel-2 data

Boreal forests constitute a large portion of the global forest area, yet they are undersampled through field surveys, and only a few remotely sensed data sources provide structural information wall-to-wall throughout the boreal domain. ArcticDEM is a collection of high-resolution (2 m) space-borne stereogrammetric digital surface models (DSM) covering the entire land area north of 60° of latitude. The free-availability of ArcticDEM data offers new possibilities for aboveground biomass mapping (AGB) across boreal forests, and thus it is necessary to evaluate the potential for these data to map AGB over alternative open-data sources (i.e., Sentinel-2). This study was performed over the entire land area of Norway north of 60° of latitude, and the Norwegian national forest inventory (NFI) was used as a source of field data composed of accurately geolocated field plots (n=7710) systematically distributed across the study area. Separate random forest models were fitted using NFI data, and corresponding remotely sensed data consisting of either: i) a canopy height model (ArcticCHM) obtained by subtracting a high-quality digital terrain model (DTM) from the ArcticDEM DSM height values, ii) Sentinel-2 (S2), or iii) a combination of the two (ArcticCHM+S2). Furthermore, we assessed the effect of the forest- and terrain-specific factors on the models’ predictive accuracy. The best model (,i.e., ArcticCHM+S2) explained nearly 60% of the variance of the training set, which translated in the largest accuracy in terms of root mean square error (RMSE=41.4 t/ha). This result highlights the synergy between 3D and multispectral data in AGB modelling. Furthermore, this study showed that despite the importance of ArcticCHM variables, the S2 model performed slightly better than ArcticCHM model. This finding highlights some of the limitations of ArcticDEM, which, despite the unprecedented spatial resolution, is highly heterogeneous due to the blending of multiple acquisitions across different years and seasons. We found that both forest- and terrain-specific characteristics affected the uncertainty of the ArcticCHM+S2 model and concluded that the combined use of ArcticCHM and Sentinel-2 represents a viable solution for AGB mapping across boreal forests. The synergy between the two data sources allowed for a reduction of the saturation effects typical of multispectral data while ensuring the spatial consistency in the output predictions due to the removal of artifacts and data voids present in ArcticCHM data. While the main contribution of this study is to provide the first evidence of the best-case-scenario (i.e., availability of accurate terrain models) that ArcticDEM data can provide for large-scale AGB modelling, it remains critically important for other studies to investigate how ArcticDEM may be used in areas where no DTMs are available as is the case for large portions of the boreal zone.

space-borne imagery↗
NASA NTRS2020

Modeling Whitecaps on Global Scale

Whitecaps play an important role in the surface-atmosphere interactions across the ocean. They are directly linked to the energy dissipation rate during wave breaking and transfer of heat, momentum, and gas/aerosol exchange from the wind to the sea surface. Although the first models of W were dependent only on wind speeds, a large number of diverse models based on wind and sea state which include wave height, wave age, friction velocity, and stability effects have been proposed since then. However, it is recognized that most of the proposed W models have strong systematic (e.g., zonal bias) and random errors when compared against observations. This is partly due to the differences in environmental conditions, measurement techniques, and geographical locations among these studies. But, some of these biases are linked to the inability of the proposed models to capture the variability in W in certain wind/wave regimes. Despite the knowledge of existing biases, W residual relationships from the models with wind and wave fields remain highly uncertain, with residual trends varying between the published studies. Here, we take advantage of the availability of relatively dense observations of W from WindSat microwave satellite retrievals in combination with the University of Miami wave model which was recently incorporated within the NASA GMAO/GEOS system (GEOS-UMWM). We use Windsat W retrievals to assess and constrain the previously published W models and understand the relationships of residuals from models in different wind/wave regimes. We link these unexplained residual variations to additional factors such as swell index, drag coefficient etc and add information to the existing whitecap models. Since Windsat retrievals cover wide range of environmental conditions, it helps to reduce the uncertainties associated with differences in measurement techniques. Regression of wind-wave fields against all Windsat data points (CTL) results in larger residuals for lower wave age and W is overestimated upto ~4% for wave age < 10 and underestimated by upto ~2% as wave age increases. We attest to this bias by considering two approaches. One is to perform regression separately for different stages of wave development such as developing sea, fully developed, and wind sea regimes thereby understanding the sensitivity of regression coefficients to sea state (EXP1). Another is to derive coefficients of W models in EXP1 as a function of additional wind/wave factors such as swell index, drag coefficient, and mean squared slope, deriving more nonlinear W models (EXP2). EXP2 provides reduction in Root Mean Squared Error (RMSE) by 0.1-0.3%. Sea surface drag has a stronger relationship with regression coefficients compared to swell index.These additional factors provide improved parameterizations in different wind and wave age regimes, with smaller unexplained/residual variations in W that has been a major concern in the W community.

Raman, Aishwarya↗
NASA NTRS2010

Probabilistic Assessment of Cancer Risk from Solar Particle Events

For long duration missions outside of the protection of the Earth s magnetic field, space radiation presents significant health risks including cancer mortality. Space radiation consists of solar particle events (SPEs), comprised largely of medium energy protons (less than several hundred MeV); and galactic cosmic ray (GCR), which include high energy protons and heavy ions. While the frequency distribution of SPEs depends strongly upon the phase within the solar activity cycle, the individual SPE occurrences themselves are random in nature. We estimated the probability of SPE occurrence using a non-homogeneous Poisson model to fit the historical database of proton measurements. Distributions of particle fluences of SPEs for a specified mission period were simulated ranging from its 5 th to 95th percentile to assess the cancer risk distribution. Spectral variability of SPEs was also examined, because the detailed energy spectra of protons are important especially at high energy levels for assessing the cancer risk associated with energetic particles for large events. We estimated the overall cumulative probability of GCR environment for a specified mission period using a solar modulation model for the temporal characterization of the GCR environment represented by the deceleration potential (^). Probabilistic assessment of cancer fatal risk was calculated for various periods of lunar and Mars missions. This probabilistic approach to risk assessment from space radiation is in support of mission design and operational planning for future manned space exploration missions. In future work, this probabilistic approach to the space radiation will be combined with a probabilistic approach to the radiobiological factors that contribute to the uncertainties in projecting cancer risks.

Kim, Myung-Hee Y.↗