Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “uncertainty Quantification”

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 433 records · Page 24

Attribution of Flux Partitioning Variations Between Land Surface Models Over the Continental U.S.

Accurate quantification of the terrestrial evapotranspiration (ET) components of plant transpiration (T), soil evaporation (E) and evaporation of the intercepted water (I) is necessary for improving our understanding of the links between the carbon and water cycles. Recent studies have noted that, among the modeled estimates, large disagreements exist in the relative contributions of T, E and I to the total ET. As these models are often used in data assimilation environments for incorporating and extending ET relevant remote sensing measurements, understanding the sources of inter-model differences in ET components is also necessary for improving the utilization of such remote sensing measurements. This study quantifies the contributions of two key factors explaining inter-model disagreements to the uncertainty in total ET: (1) contribution of the local partitioning and (2) regional distribution of ET. The analysis is conducted by using outputs from a suite of land surface models in the North American Land Data Assimilation System (NLDAS) configuration. For most of these models, transpiration is the dominant component of the ET partition. The results indicate that the uncertainty in local partitioning dominates the inter-model spread in modeled soil evaporation E. The inter-model differences in T are dominated by the uncertainty in the distribution of ET over the Eastern U.S. and the local partitioning uncertainty in the Western U.S. The results also indicate that uncertainty in the T estimates is the primary driver of total ET errors. Over the majority of the U.S., the contribution of the two factors of uncertainty to the overall uncertainty is non-trivial.

evapotranspiration partitioning; transpiration; so↗

NEWS Climatology Project: The State of the Water Cycle at Continental to Global Scales

NASA's Energy and Water Cycle Study (NEWS) program fosters collaborative research towards improved quantification and prediction of water and energy cycle consequences of climate change. In order to measure change, it is first necessary to describe current conditions. The goal of the NEWS Water and Energy Cycle Climatology project is to develop "state of the global water cycle" and "state of the global energy cycle" assessments based on data from modern ground and space based observing systems and data integrating models. The project is a multiinstitutional collaboration with more than 20 active contributors. This presentation will describe results of the first stage of the water budget analysis, whose goal was to characterize the current state of the water cycle on mean monthly, continental scales. We examine our success in closing the water budget within the expected uncertainty range and the effects of forcing budget closure as a method for refining individual flux estimates.

Rodell, Matthew↗

Transfer-AE: A novel autoencoder-based impact detection model for structural digital twin

Accurately detecting the location and intensity of impacts is crucial for ensuring structural safety. Currently, AI-based structural impact detection methods are widely used for their excellent detection accuracy. However, their generalization capability is limited by the scenarios present in the training data. Many complex and dangerous impact scenarios are difficult to conduct real-world experiments on to collect sufficient samples. To capture all impact scenarios and fully leverage the advantages of AI-based detection technologies, advanced methods involve combining real-world structural monitoring data with corresponding numerical models to construct digital twins. These methods continuously refine the created numerical models with limited real-world data and provide diverse impact scenarios through numerical model simulations. However, there are inevitable differences between digital models and physical models that are challenging to correct through mechanical means. This discrepancy in data distribution between the two models significantly hinders the application of digital twin technology in impact/event identification tasks. To address this challenge, this study proposes a novel model based on autoencoders, named Transfer-AE. Transfer-AE encodes the common features of digital twins in the latent space to bridge the uncertainty gap at a macro scale between numerical models and physical models and synchronously fits the magnitude and location of the impact load in the decoder. This enables consistent detection results for the same impact event, whether the sample comes from the numerical model or the physical model. Transfer-AE includes two operating modes: Mode 1 has a fixed computational complexity with stable inference speed, but the training cost and difficulty increase with data distribution. Mode 2's computational complexity increases with data distribution, but it has a fixed training cost and speed. In both cases involving the geodesic dome structure simulating a deep space habitat and the IASC-ASCE benchmark structure, Transfer-AE demonstrated the best performance in impact localization and quantification tasks compared to mainstream domain-adaptive transfer models.

Chengjia Han↗

Effect of Baseline Period on Quantification of Climate Extremes over the United States

Extreme weather events, including heatwaves, heavy precipitation, and drought, have a large impact on society through human health, destruction of infrastructure, ecological change, and economic losses. Reanalyses such as NASA’s Modern Era Retrospective Analysis, version 2 (MERRA-2) are a valuable tool for analyzing past extreme events to determine their underlying causes and how extremes have changed over the past four decades. The detection of extreme events relies on a threshold for precipitation or temperature that is derived using a reference period, which can then be used to determine how extreme an event was or what the return period is. Operational centers typically use a 30-year climatology period that shifts in time every ten years, while the World Meteorological Organization suggests that the maximum amount of data should be included for the detection of extreme events due to their rare occurrence. As global and regional climate continues to change, the interpretation of extreme events is reliant on the baseline period that is used for the underlying thresholds and can be a source of uncertainty for the policy making community. Three baseline periods – 1981-2010, 1991-2020, and 1981-2020 – will be used to compute percentiles of temperature and precipitation across the contiguous United States and will then by employed to determine monthly indices representing heatwaves, cold spells, and extreme precipitation events. A spatial and temporal analysis of the resulting extreme weather indices will demonstrate the appropriateness for each baseline period for the evaluation of extreme events.

Extremes↗

Two Decades Observing Smoke Above Clouds in the South-Eastern Atlantic Ocean: Deep Blue Algorithm Updates and Validation with ORACLES Field Campaign Data

This study presents and evaluates an updated algorithmfor quantification of absorbing aerosols above clouds(AACs) from passive satellite measurements. The focus isbiomass burning in the south-eastern Atlantic Ocean duringthe 2016 and 2017 ObseRvations of Aerosols above CLoudsand their intEractionS (ORACLES) field campaign deployments.The algorithm retrieves the above-cloud aerosoloptical depth (AOD) and underlying liquid cloud opticaldepth and is applied to measurements from the Sea-viewingWide Field-of-view Sensor (SeaWiFS), Moderate ResolutionImaging Spectroradiometer (MODIS), and Visible InfraredImaging Radiometer Suite (VIIRS) from 1997 to 2017. AirborneNASA Ames Spectrometers for Sky-Scanning, Sun-Tracking Atmospheric Research (4STAR) and NASA LangleyHigh Spectral Resolution Lidar 2 (HSRL2) data collectedduring ORACLES provide important validation for spectralAOD for MODIS and VIIRS; as the SeaWiFS missionended in 2010, it cannot be evaluated directly. The 4STARand HSRL2 comparisons are complementary and reveal performancegenerally in line with uncertainty estimates providedby the optimal estimation retrieval framework used. Atpresent the two MODIS-based data records seem the mostreliable, although there are differences between the deployments,which may indicate that the available data are not yetsufficient to provide a robust regional validation. Spatiotemporalpatterns in the data sets are similar, and the time seriesare very strongly correlated with each other (correlationcoefficients from 0.95 to 0.99). Offsets between the satellitedata sets are thought to be chiefly due to differences in absolutecalibration between the sensors. The available validationdata for this type of algorithm are limited to a small numberof field campaigns, and it is strongly recommended that suchairborne measurements continue to be made, both over thesouthern Atlantic Ocean and elsewhere.

Sayer, Andrew M.↗

Thinking Outside of the Blue Marble: Novel Ocean Applications Using the VIIRS Sensor

While planning for future space-borne sensors will increase the quality, quantity, and duration of ocean observations in the years to come, efforts to extend the limits of sensors currently in orbit can help shed light on future scientific gains as well as associated uncertainties. Here, we present several applications that are unique to the polar orbiting Visual Infrared Imaging Radiometer Suite (VIIRS), each of which challenge the threshold capabilities of the sensor and provide lessons for future missions. For instance, while moderate resolution polar orbiters typically have a one day revisit time, we are able to obtain multiple looks of the same area by focusing on the extreme zenith angles where orbital views overlap, and pair these observations with those from other sensors to create pseudo-geostationary data sets. Or, by exploiting high spatial resolution (imaging) channels and analyzing patterns of synoptic covariance across the visible spectrum, we can obtain higher spatial resolution bio-optical products. Alternatively, non-traditional products can illuminate important biological interactions in the ocean, such as the use of the Day-Night-Band to provide some quantification of phototactic behavior of marine life along light polluted beaches, as well as track the location of marine fishing vessel fleets along ocean fronts. In this talk, we explore ways to take full advantage of the capabilities of existing sensors in order to maximize insights for future missions.

ocean color↗

Exploring the Production of NOx by Lightning and Its Impact on Tropospheric Ozone

Our quantitative understanding of free tropospheric (FT) chemistry is quite poor. State-of-the-art regional air quality models (e.g., US EPA's CMAQ) perform very poorly in simulating FT chemistry, with Uniform ozone around 70 ppb throughout the FT in summer, while ozonesonde data show much higher levels of ozone and much spatial-temporal structure. Such models completely neglect lightning-NOx (LNOx) emissions (the most significant source of NOx in the FT), and also contain large uncertainties in the specifications of intercontinental transport, stratosphere-troposphere exchange (STE) and PBLFT exchange (PFTE). Global air chemistry models include LNOx, but in very crude fashion, with the frequency and distribution of lightning being based on modeled cloud parameters (hence large uncertainty), lightning energetics being assumed to be constant for all flashes (literature value, while in reality there is at least a two-orders of magnitude variability from flash-to-flash), and the production of NOx in the surrounding heated air, per Joule of heating, being assumed to be constant also (literature value, while in fact it is a non-linear function of the dissipated heat and local air density, p). This situation is commonly blamed on paucity of pertinent observational data, but for the USA, there is now a wealth of surface- and satellite-based data of lightning available to permit much improved observation-based estimation of LNOx emissions. In the FT, such NOx has a long residence time, and also the ozone production efficiency from NOx there is considerably higher than in the PBL. It is, therefore, of critical importance in FT chemistry. This paper will describe the approach and data products of an ongoing NSSTC project aimed at a much-improved quantification of not only LNOx production on the scale of continental USA based on local and regional lightning observations, but also of intercontinental transport, STE and PFTE, all in upgraded simulations of tropospheric transport and chemistry. In our approach for LNOx, (a) we utilize continuous observed lightning information from the NLDN ground network and from satellite imagers (OTD and LIS) to quantify lightning frequency and distribution at the spatial-temporal scales of models such as CMAQ; (b) we develop new methodologies to quantify flash-specific lightning energy dissipation as heat (epsilon) using data from the research-grade lightning measurement facility at NASA-KSC, and to parameterize epsilon based on regional lightning monitoring data (ground- and satellite-based); and, (c) we develop a new parameterization of NOx production as a function of epsilon and rho. Based on such observation-based information, we are working to develop a gridded, episodic LNOx emissions inventory for the USA for use in models like CMAQ. We are also developing approaches for global(MOZART)- regional(CMAQ) chemistry coupling to improve intercontinental transport and STE. Finally, we are developing new methodologies for assimilation of satellite-observed (GOES) clouds into meteorological modeling (MM5), to improve PFTE and to optimize co-location of cloud convection and observed lightning. We will incorporate these improvements in CMAQ simulations over the USA to better understand FT processes and chemistry, and its impact on ground-level ozone.

Gillani, Noor↗

The NEWS Water Cycle Climatology

NASA's Energy and Water Cycle Study (NEWS) program fosters collaborative research towards improved quantification and prediction of water and energy cycle consequences of climate change. In order to measure change, it is first necessary to describe current conditions. The goal of the first phase of the NEWS Water and Energy Cycle Climatology project was to develop "state of the global water cycle" and "state of the global energy cycle" assessments based on data from modern ground and space based observing systems and data integrating models. The project was a multi-institutional collaboration with more than 20 active contributors. This presentation will describe the results of the water cycle component of the first phase of the project, which include seasonal (monthly) climatologies of water fluxes over land, ocean, and atmosphere at continental and ocean basin scales. The requirement of closure of the water budget (i.e., mass conservation) at various scales was exploited to constrain the flux estimates via an optimization approach that will also be described. Further, error assessments were included with the input datasets, and we examine these in relation to inferred uncertainty in the optimized flux estimates in order to gauge our current ability to close the water budget within an expected uncertainty range.

Rodell, Matthew↗

Estimating Mangrove Aboveground Biomass from Airborne LiDAR Data: a Case Study from the Zambezi River Delta

Mangroves are ecologically and economically important forested wetlands with the highest carbon (C) density of all terrestrial ecosystems. Because of their exceptionally large C stocks and importance as a coastal buffer, their protection and restoration has been proposed as an effective mitigation strategy for climate change. The inclusion of mangroves in mitigation strategies requires the quantification of C stocks (both above and belowground) and changes to accurately calculate emissions and sequestration. A growing number of countries are becoming interested in using mitigation initiatives, such as REDD+ (reducing emissions from deforestation and forest degradation), in these unique coastal forests. However, it is not yet clear how methods to measure C traditionally used for other ecosystems can be modified to estimate biomass in mangroves with the precision and accuracy needed for these initiatives. Airborne Lidar (ALS) data has often been proposed as the most accurate way for larger scale assessments but the application of ALS for coastal wetlands is scarce, primarily due to a lack of contemporaneous ALS and field measurements. Here, we evaluated the variability in field and Lidar-based estimates of aboveground biomass (AGB) through the combination of different local and regional allometric models and standardized height metrics that are comparable across spatial resolutions and sensor types, the end result being a simplified approach for accurately estimating mangrove AGB at large scales and determining the uncertainty by combining multiple allometric models.We then quantified wall-to-wall AGB stocks of a tall mangrove forest in the Zambezi Delta, Mozambique. Our results indicate that the Lidar H100 height metric correlates well with AGB estimates, with R(exp 2) between 0.80 and 0.88 and RMSE of 33% or less. When comparing Lidar H100 AGB derived from three allometric models, mean AGB values range from 192 Mg ha(exp −1) up to 252 Mg ha(exp −1).We suggest the best model to predict AGB was based on the East Africa specific allometry and a power-based regression that used Lidar H100 as the height input with an R(exp 2) of 0.85 and an RMSE of 122 Mg ha(exp −1) or 33%. The total AGB of the Lidar inventoried mangrove area (6654 ha) was 1 350 902 Mg with a mean AGB of 203 Mg ha(exp −1) +/-166 Mg ha(exp −1). Because the allometry suggested here was developed using standardized height metrics, it is recommended that the models can generate AGB estimates using other remote sensing instruments that are more readily accessible over other mangrove ecosystems on a large scale, and as part of future carbon monitoring efforts in mangroves.

Fatoyinbo, Temilola↗

Dynamic remapping decisions in multi-phase parallel computations

The effectiveness of any given mapping of workload to processors in a parallel system is dependent on the stochastic behavior of the workload. Program behavior is often characterized by a sequence of phases, with phase changes occurring unpredictably. During a phase, the behavior is fairly stable, but may become quite different during the next phase. Thus a workload assignment generated for one phase may hinder performance during the next phase. We consider the problem of deciding whether to remap a paralled computation in the face of uncertainty in remapping's utility. Fundamentally, it is necessary to balance the expected remapping performance gain against the delay cost of remapping. This paper treats this problem formally by constructing a probabilistic model of a computation with at most two phases. We use stochastic dynamic programming to show that the remapping decision policy which minimizes the expected running time of the computation has an extremely simple structure: the optimal decision at any step is followed by comparing the probability of remapping gain against a threshold. This theoretical result stresses the importance of detecting a phase change, and assessing the possibility of gain from remapping. We also empirically study the sensitivity of optimal performance to imprecise decision threshold. Under a wide range of model parameter values, we find nearly optimal performance if remapping is chosen simply when the gain probability is high. These results strongly suggest that except in extreme cases, the remapping decision problem is essentially that of dynamically determining whether gain can be achieved by remapping after a phase change; precise quantification of the decision model parameters is not necessary.

Nicol, D. M.↗

Quantification of Aerosol Type, and Sources of Aerosols Over the Indo- Gangetic Plain

Differences and similarities in aerosol characteristics, for the first time, over two environmentally distinct locations in Indo-Gangetic plain (IGP) e Kanpur (KPR) (urban location) and Gandhi College (GC) (rural site) are examined. Aerosol optical depths (AODs) exhibit pronounced seasonal variability with higher values during winter and premonsoon. Aerosol fine mode fraction (FMF) and Ångstrom exponent (a) are higher over GC than KPR indicating relatively higher fine mode aerosol concentration over GC. Higher FMF over GC is attributed to local biomass burning activities. Analysis of AOD spectra revealed that aerosol size distribution is dominated by wide range of fine mode fractions or mixture of modes during winter and postmonsoon, while during premonsoon and monsoon coarse mode aerosols are more abundant. Single scattering albedo (SSA) is lower over GC than KPR. SSA spectra reveals the abundance of fine mode (coarse mode) absorbing (scattering) aerosols during winter and postmonsoon (premonsoon and monsoon). Spectral SSA features reveal that OC contribution to enhanced absorption is negligible. Analysis shows that absorbing aerosols can be classified as Mostly Black Carbon (BC), and Mixed BC and Dust over IGP. Mixed BC and dust is always higher over KPR, while Mostly BC is higher over GC throughout the year. The amount of long range transported dust exhibits a gradient between KPR (higher) and GC (lower). Results on seasonally varying aerosol types, and absorbing aerosol types and their gradients over an aerosol hotspot are important to tune models and to reduce the uncertainty in radiative and climate impact of aerosols.

Organic Carbon↗

On the Statistical Analysis of X-ray Polarization Measurements

In many polarimetry applications, including observations in the X-ray band, the measurement of a polarization signal can be reduced to the detection and quantification of a deviation from uniformity of a distribution of measured angles of the form alpha plus beta cosine (exp 2)(phi - phi(sub 0) (0 (is) less than phi is less than pi). We explore the statistics of such polarization measurements using both Monte Carlo simulations as well as analytic calculations based on the appropriate probability distributions. We derive relations for the number of counts required to reach a given detection level (parameterized by beta the "number of sigma's" of the measurement) appropriate for measuring the modulation amplitude alpha by itself (single interesting parameter case) or jointly with the position angle phi (two interesting parameters case). We show that for the former case when the intrinsic amplitude is equal to the well known minimum detectable polarization (MDP) it is, on average, detected at the 3sigma level. For the latter case, when one requires a joint measurement at the same confidence level, then more counts are needed, by a factor of approximately equal to 2.2, than that required to achieve the MDP level. We find that the position angle uncertainty at 1sigma confidence is well described by the relation sigma(sub pi) equals 28.5(degrees) divided by beta.

Strohmayer, T. E.↗

PACE Technical Report Series, Volume 7: Ocean Color Instrument (OCI) Concept Design Studies

Extending OCI hyperspectral radiance measurements in the ultraviolet to 320 nm on the blue spectrograph enables quantitation of atmospheric total column ozone (O3) for use in ocean color atmospheric correction algorithms. The strong absorption by atmospheric ozone below 340 nm enables the quantification of total column ozone. Other applications are possible but were not investigated due to their exploratory nature and lower priority.The first step in the atmospheric correction processing, which converts top-of-the-atmosphere radiances to water-leaving radiances, is removal of the absorbance by atmospheric trace gases such as water vapor, oxygen, ozone and nitrogen dioxide. Details of the atmospheric correction process currently used by the Ocean Biology Processing Group (OBPG) and will be employed for PACE with appropriate modifications, are described by Mobley et al. [2016]. Atmospheric ozone absorbs within the visible to near-infrared spectrum between ~450 nm and 800nm and most appreciably between 530 nm and 650 nm, a spectral region critical for maintaining NASA's chlorophyll-a climate data record and for PACE algorithms planned to characterize phytoplankton community composition and other ocean color products.While satellite-based observations will likely be available during PACE's mission lifetime, the difference in acquisition time with PACE, the coarseness in their spatial resolution, and differences in viewing geometries will introduce significant levels of uncertainties in PACE ocean color data products.

Cetinic, Ivona↗

Detection and Retrieval of Multi-Layered Cloud Properties Using Satellite Data

Four techniques for detecting multilayered clouds and retrieving the cloud properties using satellite data are explored to help address the need for better quantification of cloud vertical structure. A new technique was developed using multispectral imager data with secondary imager products (infrared brightness temperature differences, BTD). The other methods examined here use atmospheric sounding data (CO2-slicing, CO2), BTD, or microwave data. The CO2 and BTD methods are limited to optically thin cirrus over low clouds, while the MWR methods are limited to ocean areas only. This paper explores the use of the BTD and CO2 methods as applied to Moderate Resolution Imaging Spectroradiometer (MODIS) and Advanced Microwave Scanning Radiometer EOS (AMSR-E) data taken from the Aqua satellite over ocean surfaces. Cloud properties derived from MODIS data for the Clouds and the Earth's Radiant Energy System (CERES) Project are used to classify cloud phase and optical properties. The preliminary results focus on a MODIS image taken off the Uruguayan coast. The combined MW visible infrared (MVI) method is assumed to be the reference for detecting multilayered ice-over-water clouds. The BTD and CO2 techniques accurately match the MVI classifications in only 51 and 41% of the cases, respectively. Much additional study is need to determine the uncertainties in the MVI method and to analyze many more overlapped cloud scenes.

Minnis, Patrick↗

Top-down Estimate of Dust Emissions Through Integration of MODIS and MISR Aerosol Retrievals With the Geos-chem Adjoint Model

Predicting the influences of dust on atmospheric composition, climate, and human health requires accurate knowledge of dust emissions, but large uncertainties persist in quantifying mineral sources. This study presents a new method for combined use of satellite-measured radiances and inverse modeling to spatially constrain the amount and location of dust emissions. The technique is illustrated with a case study in May 2008; the dust emissions in Taklimakan and Gobi deserts are spatially optimized using the GEOSChem chemical transport model and its adjoint constrained by aerosol optical depth (AOD) that are derived over the downwind dark-surface region in China from MODIS (Moderate Resolution Imaging Spectroradiometer) reflectance with the aerosol single scattering properties consistent with GEOS-chem. The adjoint inverse modeling yields an overall 51% decrease in prior dust emissions estimated by GEOS-Chem over the Taklimakan-Gobi area, with more significant reductions south of the Gobi Desert. The model simulation with optimized dust emissions shows much better agreement with independent observations from MISR (Multi-angle Imaging SpectroRadiometer) AOD and MODIS Deep Blue AOD over the dust source region and surface PM10 concentrations. The technique of this study can be applied to global multi-sensor remote sensing data for constraining dust emissions at various temporal and spatial scales, and hence improving the quantification of dust effects on climate, air quality, and human health.

MODIS aerosol retrievals↗

A Quantitative Study of Oxygen as a Metabolic Regulator

An acute reduction in oxygen (O2) delivery to a tissue is generally associated with a decrease in phosphocreatine, increases in ADP, NADH/NAD, and inorganic phosphate, increased rates of glycolysis and lactate production, and reduced rates of pyruvate and fatty acid oxidation. However, given the complexity of the human bioenergetic system and its components, it is difficult to determine quantitatively how cellular metabolic processes interact to maintain ATP homeostasis during stress (e.g., hypoxia, ischemia, and exercise). Of special interest is the determination of mechanisms relating tissue oxygenation to observed metabolic responses at the tissue, organ, and whole body levels and the quantification of how changes in tissue O2 availability affect the pathways of ATP synthesis and the metabolites that control these pathways. In this study, we extend a previously developed mathematical model of human bioenergetics to provide a physicochemical framework that permits quantitative understanding of O2 as a metabolic regulator. Specifically, the enhancement permits studying the effects of variations in tissue oxygenation and in parameters controlling the rate of cellular respiration on glycolysis, lactate production, and pyruvate oxidation. The whole body is described as a bioenergetic system consisting of metabolically distinct tissue/organ subsystems that exchange materials with the blood. In order to study the dynamic response of each subsystem to stimuli, we solve the ordinary differential equations describing the temporal evolution of metabolite levels, given the initial concentrations. The solver used in the present study is the packaged code LSODE, as implemented in the NASA Lewis kinetics and sensitivity analysis code, LSENS. A major advantage of LSENS is the efficient procedures supporting systematic sensitivity analysis, which provides the basic methods for studying parameter sensitivities (i.e., changes in model behavior due to parameter variation). Sensitivity analysis establishes relationships between model predictions and problem parameters (i.e., initial concentrations, rate coefficients, etc). It helps determine the effects of uncertainties or changes in these input parameters on the predictions, which ultimately are compared with experimental observations in order to validate the model. Sensitivity analysis can identify parameters that must be determined accurately because of their large effect on the model predictions and parameters that need not be known with great precision because they have little or no effect on the solution. This capability may prove to be important in optimizing the design of experiments, thereby reducing the use of animals. This approach can be applied to study the metabolic effects of reduced oxygen delivery to cardiac muscle due to local myocardial ischemia and the effects of acute hypoxia on brain metabolism. Other important applications of sensitivity analysis include identification of quantitatively relevant pathways and biochemical species within an overall mechanism, when examining the effects of a genetic anomaly or pathological state on energetic system components and whole system behavior.

Radhakrishnan, Krishnan↗

Airborne Lidar-Based Estimates of Tropical Forest Structure in Complex Terrain: Opportunities and Trade-Offs for REDD+

Background: Carbon stocks and fluxes in tropical forests remain large sources of uncertainty in the global carbon budget. Airborne lidar remote sensing is a powerful tool for estimating aboveground biomass, provided that lidar measurements penetrate dense forest vegetation to generate accurate estimates of surface topography and canopy heights. Tropical forest areas with complex topography present a challenge for lidar remote sensing. Results: We compared digital terrain models (DTM) derived from airborne lidar data from a mountainous region of the Atlantic Forest in Brazil to 35 ground control points measured with survey grade GNSS receivers. The terrain model generated from full-density (approx. 20 returns/sq m) data was highly accurate (mean signed error of 0.19 +/-0.97 m), while those derived from reduced-density datasets (8/sq m, 4/sq m, 2/sq m and 1/sq m) were increasingly less accurate. Canopy heights calculated from reduced-density lidar data declined as data density decreased due to the inability to accurately model the terrain surface. For lidar return densities below 4/sq m, the bias in height estimates translated into errors of 80-125 Mg/ha in predicted aboveground biomass. Conclusions: Given the growing emphasis on the use of airborne lidar for forest management, carbon monitoring, and conservation efforts, the results of this study highlight the importance of careful survey planning and consistent sampling for accurate quantification of aboveground biomass stocks and dynamics. Approaches that rely primarily on canopy height to estimate aboveground biomass are sensitive to DTM errors from variability in lidar sampling density.

Airborne lidar↗

A Probabilistic Asteroid Impact Risk Model

Asteroid threat assessment requires the quantification of both the impact likelihood and resulting consequence across the range of possible events. This paper presents a probabilistic asteroid impact risk (PAIR) assessment model developed for this purpose. The model incorporates published impact frequency rates with state-of-the-art consequence assessment tools, applied within a Monte Carlo framework that generates sets of impact scenarios from uncertain parameter distributions. Explicit treatment of atmospheric entry is included to produce energy deposition rates that account for the effects of thermal ablation and object fragmentation. These energy deposition rates are used to model the resulting ground damage, and affected populations are computed for the sampled impact locations. The results for each scenario are aggregated into a distribution of potential outcomes that reflect the range of uncertain impact parameters, population densities, and strike probabilities. As an illustration of the utility of the PAIR model, the results are used to address the question of what minimum size asteroid constitutes a threat to the population. To answer this question, complete distributions of results are combined with a hypothetical risk tolerance posture to provide the minimum size, given sets of initial assumptions. Model outputs demonstrate how such questions can be answered and provide a means for interpreting the effect that input assumptions and uncertainty can have on final risk-based decisions. Model results can be used to prioritize investments to gain knowledge in critical areas or, conversely, to identify areas where additional data has little effect on the metrics of interest.

Asteroid↗