Search NASA⌕ Search

SEARCH · Search NASA

Results for “trend estimation”

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 307 records · Page 17

Applications of OMI Tropospheric NO(Sub 2) Data: Air Quality Trends in the US and Lightning NO(x) Source Strength

This talk will be presented in two parts: 1) an analysis of tropospheric column NO2 trends in the eastern half of the United States over the period 2005 to 2009 and 2) estimation of lightning NO(x) production rates based on OMI observations and lightning flash rate data. The air quality trends in the eastern US will be determined for specific subregions using tropospheric column NO2 data from OMI for 2005 through 2008 and from GOME-2 for 2007 through 2009. This period is characterized by significant NO(x) emission reductions at power plants within most of this region. The air quality trends will be compared with those estimated from continuous emission monitoring data from the power plants compiled by the US Environmental Protection Agency. OMI NO2 data have also been used to estimate lightning NOx production per flash in selected storms near Costa Rica and Panama during the 2007 NASA TC4 field campaign and over the continental US, Gulf of Mexico, and western Atlantic during the summers of 2005 and 2006. The lightning signal is extracted from the OMI data through a custom retrieval in which an NO2 profile representative of convective outflow is used in the airmass factor calculation and the background NO2 column is subtracted from the tropospheric column. When combined with NO(x)/NO2 ratios from the NASA GMT model and observed flash rates, the resulting estimates of NO(x) production per flash are comparable to those estimated obtained from analyses of aircraft data and cloud-resolving modeling.

Pickering, Kenneth↗

150 Shades of Green: Using the Full Spectrum of Remote Sensing Reflectance to Elucidate Color Shifts in the Ocean

This article proposes a simple and intuitive classification system by which to define full spectral remote sensing reflectance (Rrs(λ)) data with a quantitative output that enables a more manageable handling of spectral information for aquatic science applications. The weighted harmonic mean of the Rrs(λ) wavelengths outputs an Apparent Visible Wavelength (in units of nanometers), representing a one-dimensional geophysical metric of color that is inherently correlated to spectral shape. This dimensionality reduction of spectral information combined with the output along a continuum of wavelength values offers a robust and user-friendly means to describe and analyze spectral Rrs(λ) in terms of spatial and temporal trends and variability. The uncertainty in the algorithm's estimation of spectral shape is demonstrated on a global scale, in addition to the utility of the algorithm to discern spectral-spatial-temporal trends in the ocean, on a per-pixel basis for the entire 22 year continuous ocean color (SeaWiFS and MODIS-Aqua) time-series. This technique can be applied to datasets of varying multi- and hyper-spectral resolutions, providing continuity between heritage and future satellite sensors, and further enabling an effective means of elucidating similarities or differences in complex spectral signatures within the constraints of two dimensions. This straightforward means of conceptualizing multi-dimensional variability can help maximize the potential of the spectral information embedded in remote sensing data.

ocean color↗

Impact of Satellite Viewing-Swath Width on Global and Regional Aerosol Optical Thickness Statistics and Trends

We use the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite aerosol optical thickness (AOT) product to assess the impact of reduced swath width on global and regional AOT statistics and trends. Alongtrack and across-track sampling strategies are employed, in which the full MODIS data set is sub-sampled with various narrow-swath (approximately 400-800 km) and single pixel width (approximately 10 km) configurations. Although view-angle artifacts in the MODIS AOT retrieval confound direct comparisons between averages derived from different sub-samples, careful analysis shows that with many portions of the Earth essentially unobserved, spatial sampling introduces uncertainty in the derived seasonal-regional mean AOT. These AOT spatial sampling artifacts comprise up to 60%of the full-swath AOT value under moderate aerosol loading, and can be as large as 0.1 in some regions under high aerosol loading. Compared to full-swath observations, narrower swath and single pixel width sampling exhibits a reduced ability to detect AOT trends with statistical significance. On the other hand, estimates of the global, annual mean AOT do not vary significantly from the full-swath values as spatial sampling is reduced. Aggregation of the MODIS data at coarse grid scales (10 deg) shows consistency in the aerosol trends across sampling strategies, with increased statistical confidence, but quantitative errors in the derived trends are found even for the full-swath data when compared to high spatial resolution (0.5 deg) aggregations. Using results of a model-derived aerosol reanalysis, we find consistency in our conclusions about a seasonal-regional spatial sampling artifact in AOT Furthermore, the model shows that reduced spatial sampling can amount to uncertainty in computed shortwave top-ofatmosphere aerosol radiative forcing of 2-3 W m(sup−2). These artifacts are lower bounds, as possibly other unconsidered sampling strategies would perform less well. These results suggest that future aerosol satellite missions having significantly less than full-swath viewing are unlikely to sample the true AOT distribution well enough to obtain the statistics needed to reduce uncertainty in aerosol direct forcing of climate.

remote sensing↗

Calibration of TOPEX/POSEIDON at Platform Harvest

We present estimates for the mean bias of the TOPEX/POSEIDON NASA altimeter (ALT) and the Centre National d'Etudes Spatiales altimeter (SSALT) using in-situ data gathered at Platform Harvest during the first 36 cycles of the mission. Data for 21 overflights of the ALT and six overflights of the SSALT have been analyzed. The analysis includes an independent assessment of in-situ measurements of sea level, the radial component of the orbit, wet tropospheric path delay, and ionospheric path delay. (The sign convention used is such that, to correct the geophysical data record values for sea level, add the bias algebraically. Unless otherwise stated, the uncertainty in a given parameter is depicted by +/- sigma(sub x), where sigma(sub x) is the sample standard deviation of x about the mean.) Tide gauges at Harvest provide estimates of sea level with an uncertainty of +/- 1.5 cm. The uncertainty in the radial component of the orbit is estimated to be +/- 1.3 cm. In-situ measurements of tropopsheric path delay at Harvest compare to within +/- 1.3 cm of the TOPEX/POSEIDON microwave radiometer, and in-situ measurements of the ionospheric path delay compare to within -0.4 +/- 0.7 cm of the dual-frequency ALT and 1.1 +/- 0.6 cm of Doppler orbitography and radiopositioning integrated by satellite. We obtain mean bias estimates of -14.5 +/- 2.9 cm for the ALT and +0.9 +/- 3.1 cm for the SSALT (where the uncertainties are based on the standard deviation of the estimated mean (sigma(sub bar x/y), which is derived from sample statistics and estimates for errors that cannot be observed). These results are consistent with independent estimates for the relative bias between the two altimeters. A linear regression applied to the complete set of data shows that there is a discernable secular trend in the time series for the ALT bias estimates. A preliminary analysis of data obtained through cycle 48 suggests that the apparent secular drift may be the result of a poorly sampled annual signal.

Christensen, E. J.↗

Main bearing response in a waked 15-MW floating wind turbine in below-rated conditions

Increased wind turbine size raises unknowns related to structural flexibility. Moreover, moving to deeper waters, component reliability becomes more critical. This work investigates main bearing response dependence on drivetrain flexibility and wake impingement in a two-turbine wind farm. A 15-MW floating direct-drive turbine is considered. Large eddy simulations (LES) are employed to model neutral, stable and unstable atmospheric conditions at below-rated mean wind speed, while the engineering codes OpenFAST and FAST.Farm simulate turbine and wake behavior. Results indicate significant sensitivities in fatigue estimates to lateral distance between the upstream and downstream turbine. The trends are most substantial in stable conditions, where the waked downwind main bearing sees twice the fatigue damage estimates of the upstream turbine for one position and 50% for another. Main bearing fatigue sensitivity to drivetrain flexibility is minor, while properly including generator rotor inertia loads is important for the axial forces of the locating (axially fixed) bearing, especially in stable conditions.

17 WIND ENERGY↗

Contrasting long-term temperature trends reveal minor changes in projected potential evapotranspiration in the US Midwest

Abstract Warming generally leads to increased evaporative demand, altering the amount of water needed for growing crops. For the Midwest, some studies have suggested that reaching yield targets by 2050 will not be possible without additional precipitation or large expansion of irrigation. Here, we show that this claim is not supported by the historical summer climate trends, which indicate that the warming of daily average temperatures is largely driven by increases in minimum temperatures, while maximum temperatures have decreased. This has translated into a net decrease in vapor pressure deficit (VPD) and potential evapotranspiration (PET). With the increasing rainfall, this suggests that crop water deficits have likely become less frequent in the region despite the warming climate. By projecting these trends into 2050 and ancillary use of a crop model, we estimate minor changes in PET that would have minimal effects on corn yields (<6%) under persistence of these trends.

54 ENVIRONMENTAL SCIENCES↗

Challenge Problem 1: Preliminary Model Development and Assessment of Flexible Heat Transfer Modeling Approaches

This report presents the modeling progress of a university consortium studying flexible modeling for turbulent heat transfer. In the context of advanced nuclear reactor design, where challenges include non-unity Prandtl fluids, free and mixed convection, and transitional flow, computational fluid dynamics methods are underutilized due to low confidence in modeling approaches and high computational cost. This report evaluates the performance of engineering CFD models in diabatic, buoyant turbulent flow. It finds that all turbulence models including several variants of the k-ε and k-ω models struggle to predict accurate turbulent momentum and heat transfer in such flows. The Nusselt numbers have been compared between the models and the DNS data, where calculations have been performed for each case and trends have shown a good agreement between DNS estimated Nusselt numbers and available correlations as well as experimental data. A novel DNS correlation for high Pr cases as they are transitioning from mixed convection to natural convection has been proposed. While Nusselt number errors relative to DNS range from 20% to 50%, the models capture similar trends to DNS with respect to Nusselt suppression and amplification under varying levels of buoyancy effect. This report also highlights model form error as a significant contributor to CFD predictions and proposes a framework for quantifying model error and improving confidence in CFD calculations. To improve the predictive capability of engineering CFD models, data-driven approaches for turbulence models are investigated. Theoretical frameworks based on the invariant tensor / vector basis neural networks for prediction of Reynolds stresses and turbulent heat fluxes are employed. The models are developed using direct numerical simulations data for forced convection flows of different fluids in vertical planar channel domain. The framework is implemented in spectral element solvers Nek5000 / nekRS and has shown a potential for future development and consideration of mixed convection flows.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Initial Sensorimotor and Cardiovascular Data Acquired from Soyuz Landings: Establishing a Functional Performance Recovery Time Constant

INTRODUCTION Testing of crew responses following long-duration flights has not been previously possible until a minimum of more than 24 hours after landing. As a result, it has not been possible to determine the trend of the early recovery process, nor has it been possible to accurately assess the full impact of the decrements associated with long-duration flight. To overcome these limitations, both the Russian and U.S. programs have implemented joint testing at the Soyuz landing site. This International Space Station research effort has been identified as the functional Field Test, and represents data collect on NASA, Russian, European Space Agency, and Japanese Aerospace Exploration Agency crews. RESEARCH The primary goal of this research is to determine functional abilities associated with long-duration space flight crews beginning as soon after landing as possible on the day of landing (typically within 1 to 1.5 hours). This goal has both sensorimotor and cardiovascular elements. To date, a total of 15 subjects have participated in a 'pilot' version of the full 'field test'. The full version of the 'field test' will assess functional sensorimotor measurements included hand/eye coordination, standing from a seated position (sit-to-stand), walking normally without falling, measurement of dynamic visual acuity, discriminating different forces generated with the hands (both strength and ability to judge just noticeable differences of force), standing from a prone position, coordinated walking involving tandem heel-to-toe placement (tested with eyes both closed and open), walking normally while avoiding obstacles of differing heights, and determining postural ataxia while standing (measurement of quiet stance). Sensorimotor performance has been obtained using video records, and data from body worn inertial sensors. The cardiovascular portion of the investigation has measured blood pressure and heart rate during a timed stand test in conjunction with postural ataxia testing (quiet stance sway) as well as cardiovascular responses during sensorimotor testing on all of the above measures. We have also collected motion sickness data associated with each of the postflight tests. When possible rudimentary cerebellar assessment was undertaken. In addition to the immediate post-landing collection of data, postflight data has been acquired twice more within 24 hours after landing and measurements continue until sensorimotor and cardiovascular responses have returned to preflight normative values (approximately 60 days postflight). SUMMARY The level of functional deficit observed in the crew tested to date is more severe than expected, clearly triggered by the acquisition of gravity loads immediately after landing when the demands for crew intervention in response to emergency operations will be greatest. Measureable performance parameters such as ability to perform a seat egress, recover from a fall or the ability to see clearly when walking, and related physiologic data (orthostatic responses) are required to provide an evidence base for characterizing programmatic risks and the degree of variability among crewmembers for exploration missions where the crew will be unassisted after landing. Overall, these early functional and related physiologic measurements will allow the estimation of nonlinear sensorimotor and cardiovascular recovery trends that have not been previously captured.

Reschke, M. F.↗

When Will MISR Detect Rising High Clouds?

It is predicted by both theory and models that high-altitude clouds will occur higher in the atmosphere as a result of climate warming. This produces a positive longwave feedback and has a substantial impact on the Earth's response to warming. This effect is well established by theory, but is poorly constrained by observations, and there is large spread in the feedback strength between climate models. For this study, we use the NASA Multi-angle Imaging SpectroRadiometer (MISR) to examine changes in Cloud-Top-Height (CTH). MISR uses a stereo-imaging technique to determine CTH. This approach is geometric in nature and insensitive to instrument calibration and therefore is well suited for trend analysis and studies of variability on long time scales. In this article we show that the current MISR record does have an increase in CTH for high-altitude cloud over Southern Hemisphere (SH) oceans but not over Tropical or the Northern Hemisphere (NH) oceans. We use climate model simulations to estimate when MISR might be expected to detect trends in CTH, that include the NH. The analysis suggests that according to the models used in this study MISR should detect changes over the SH ocean earlier than the NH, and if the model predictions are correct should be capable of detecting a trend over the Tropics and NH very soon (3–10 years). This result highlights the potential value of a follow-on mission to MISR, which no longer maintains a fixed equator crossing time and is unlikely to be making observations for another 10 years.

54 ENVIRONMENTAL SCIENCES↗

Arctic Sea Ice Parameters from AMSR-E Data using Two Techniques, and Comparisons with Sea Ice from SSM

We use two algorithms to process AMSR-E data in order to determine algorithm dependence, if any, on the estimates of sea ice concentration, ice extent and area, and trends and to evaluate how AMSR-E data compare with historical SSM/I data. The monthly ice concentrations derived from the two algorithms from AMSR-E data (the AMSR-E Bootstrap Algorithm, or ABA, and the enhanced NASA Team algorithm, or NT2) differ on average by about 1 to 3%, with data from the consolidated ice region being generally comparable for ABA and NT2 retrievals while data in the marginal ice zones and thin ice regions show higher values when the NT2 algorithm is used. The ice extents and areas derived separately from AMSR-E using these two algorithms are, however, in good agreement, with the differences (ABA-NT2) being about 6.6 x 10(exp 4) square kilometers on average for ice extents and -6.6 x 10(exp 4) square kilometers for ice area which are small compared to mean seasonal values of 10.5 x 10(exp 6) and 9.8 x 10(exp 6) for ice extent and area: respectively. Likewise, extents and areas derived from the same algorithm but from AMSR-E and SSM/I data are consistent but differ by about -24.4 x 10(exp 4) square kilometers and -13.9 x 10(exp 4) square kilometers, respectively. The discrepancies are larger with the estimates of extents than area mainly because of differences in channel selection and sensor resolutions. Trends in extent during the AMSR-E era were also estimated and results from all three data sets are shown to be in good agreement (within errors).

Comiso, Josefino C.↗

Regional Characteristics of Observable Foreshocks

Abstract Measures of foreshock occurrence are systematically examined using earthquake catalogs for eight regions (Italy, southern California, northern California, Costa Rica, Onshore Japan, Alaska, Turkey, and Greece) after imposing a magnitude ≥3.0 completeness level. Foreshocks are identified using three approaches: a magnitude-dependent space + fixed-time windowing method, a nearest-neighbor clustering method, and a modified magnitude-dependent space + variable-time windowing method. The method with fixed-time windows systematically yields higher counts of foreshocks than the other two clustering methods. We find similar counts of foreshocks across the three methods when the magnitude aperture is equalized by including only earthquakes in the magnitude range M*−2≤ M< M*, in which M* is the mainshock magnitude. For most of the catalogs (excluding Italy and southern California), the measured b-values of the foreshocks of all region-specific mainshocks are lower by 0.1–0.2 than b-values of respective aftershocks. Allowing for variable-time windows results in relatively high probabilities of having at least one foreshock in Italy (∼43%–56%), compared to other regional catalogs. Foreshock probabilities decrease to 14%–41% for regions such as Turkey, Greece, and Costa Rica. Similar trends are found when requiring at least five foreshocks in a sequence to be considered. Estimates of foreshock probabilities for each mainshock are method dependent; however, consistent regional trends exist regardless of method, with regions such as Italy and southern California producing more observable foreshocks than Turkey and Greece. Some regions with relatively high background seismicity have comparatively low probabilities of detectable foreshock activity when using methods that account for variable background, possibly due to depletion of near-failure fault conditions by background activity.

Geochemistry & Geophysics↗

Understanding how chemical structure affects ignition-delay-time $\phi$-sensitivity

$\phi$-sensitivity is the change in ignition delay time (IDT) with respect to the fuel-to-air equivalence ratio ($\phi$). High $\phi$-sensitivity is a desirable fuel property for applications in advanced compression ignition and multi-mode engine designs. Understanding how $\phi$-sensitivity depends on chemical structure is essential for selecting promising biofuels from the ever-growing list of proposed candidates. Here, we investigate the effect of chemical structure on $\phi$-sensitivity with experiment, simulation, and theory. Experimental Advanced Fuel Ignition Delay Analyzer (AFIDA) measurements for 2,4-dimethylpentane and diisopropyl ether provide evidence that branching and functional groups strongly impact $\phi$-sensitivity. Further insights into this dependence are obtained with 0-D kinetic simulations with existing mechanisms for n-pentane, diethyl ether, 3-pentanone, n-heptane, 2-methylhexane, 2,4-dimethylpentane, and 2,2,3-trimethylbutane. Quantum mechanical (QM) G4 calculations of low-temperature reactions help explain the observed experimental and simulation trends. Specifically, these QM calculations provide theoretical estimates of the ketohydroperoxide (KHP) dissociation rates, the HO 2 formation rates from peroxy radical (ROO), and the “cross-over” temperatures, i.e., the temperature at which ROO dissociation is favored compared to hydroperoxyl radical (QOOH) formation. Each of these reaction rates is compared to the n-alkane reference point to determine the impact of branching and different functional groups. Although kinetic mechanisms typically assume that KHP dissociation rates are invariant of chemical environment, our QM results suggest that this rate can span a range of roughly two orders of magnitude. We also discuss the importance of including the peroxy-hydroperoxy (OO-OOH) hydrogen transfer reaction for branched ethers. Finally, the insights gained assist in proposing a highly $\phi$-sensitive compound, namely, isopropyl propyl ether.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying the impacts of micro- and mild- hybrid vehicle technologies on fleetwide fuel economy and electrification

Micro- and mild-hybridization (jointly labeled as M-HEV) is gaining popularity as a cost-effective technology for fuel economy improvement, but whether and how M-HEV may compete against less efficient conventional internal combustion engine vehicles (ICEV), more efficient full hybrid electric vehicles (HEV), and plug-in electric vehicles (PEV) is not well understood. As such, this study aims at evaluating the impact of the market adoption of M-HEV on the average fuel economy of the new vehicle fleet and on the sales share of PEVs. The study reviews recent sales trends and market forecasts, and uses published estimates of manufacturing cost and fuel economy of M-HEV with an existing discrete choice model (Market Acceptance of Advanced Automotive Technologies or MA3T) to project the market penetration and impacts of M-HEV under different scenarios of M-HEV choice positions, designed to enhance conclusion robustness. It is found that among engine-based powertrain choices, micro-HEV appears to be the most cost-effective, followed by ICEVs, mild-HEV and finally full HEVs. M-HEV technologies are likely to improve fleetwide average fuel economy without significant adverse effects on sales of plug-in electric vehicles, and are likely to remain highly competitive outside PEVs through 2050.

33 ADVANCED PROPULSION SYSTEMS↗

Fossil fuel CO 2 emissions over metropolitan areas from space: A multi-model analysis of OCO-2 data over Lahore, Pakistan

Urban areas, where gathering more than 55% of the global population, alone contributed to more than 70% of anthropogenic fossil fuel carbon dioxide (CO 2ff ) emissions. Accurate quantification of CO 2ff emissions from urban areas is of great importance to the formulation of global warming mitigation policies to achieve carbon neutrality by 2050. Satellite-based inversion techniques are unique among “top-down” approaches, potentially allowing us to track CO 2ff emission changes over cities globally. However, its accuracy is still limited by incomplete background information, cloud blockages, aerosol contaminations, and uncertainties in models and emission inventories used as prior. To evaluate the current potential of space-based quantification techniques, we present the first attempt to monitor long-term changes in CO 2ff emissions based on the OCO-2 satellite measurements of column-averaged dry-air mole fractions of CO 2 (X CO2 ) over a fast-growing Asian metropolitan area: Lahore, Pakistan. We first examined the OCO-2 data availability at global scale. About 17% of OCO-2 soundings are marked as high-quality soundings by quality flags over the global 70 most populated cities over the period 2014-2019. Cloud blockage and aerosol contamination are the two main causes of data loss. As an attempt to recover additional soundings, we evaluated the effectiveness of OCO-2 quality flags at the city level by comparing three flux quantification methods (WRF-Chem, X-STILT, and flux cross-sectional integration method), all based on the Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC) product. The satellite/bottom-up emissions (OCO-2/ODIAC) ratios of the high-quality tracks better converged across the three methods compared to the all-data tracks with reduced uncertainties in emissions. Thus, OCO-2 quality flags are useful filters of low-quality OCO-2 retrievals at local scales, although originally designed for global-scale studies. All three methods consistently suggested that the ratio medians are greater than 1, which implies that the ODIAC slightly underestimated the CO2ff emissions over Lahore. Additionally, our estimation of the a posteriori CO2ff emission trend was about 734 kt C/year (i.e., an annual 6.7% increase). 10,000 Monte Carlo simulations of the Mann-Kendall upward trend test showed that less than 10% prior uncertainty for 8 tracks (or less than 20% prior uncertainty for 25 tracks) is required to achieve a greater-than-50% trend significant possibility at a 95% confidence level. It implies that the trend is driven by the prior and not due to the assimilation of OCO-2 retrievals. The key to improving the role of satellite data in CO 2 emission trend detection lies in collecting more frequent high-quality tracks near metropolitan areas to achieve significant constraints from X CO2 retrievals.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Synthetic method of analogues for emerging infectious disease forecasting

The Method of Analogues (MOA) has gained popularity in the past decade for infectious disease forecasting due to its non-parametric nature. In MOA, the local behavior observed in a time series is matched to the local behaviors of several historical time series. The known values that directly follow the historical time series that best match the observed time series are used to calculate a forecast. This non-parametric approach leverages historical trends to produce forecasts without extensive parameterization, making it highly adaptable. However, MOA is limited in scenarios where historical data is sparse. This limitation was particularly evident during the early stages of the COVID-19 pandemic, where the emerging global epidemic had little-to-no historical data. In this work, we propose a new method inspired by MOA, called the Synthetic Method of Analogues (sMOA). sMOA replaces historical disease data with a library of synthetic data that describe a broad range of possible disease trends. This model circumvents the need to estimate explicit parameter values by instead matching segments of ongoing time series data to a comprehensive library of synthetically generated segments of time series data. We demonstrate that sMOA has competitive performance with state-of-the-art infectious disease forecasting models, out-performing 78% of models from the COVID-19 Forecasting Hub in terms of averaged Mean Absolute Error and 76% of models from the COVID-19 Forecasting Hub in terms of averaged Weighted Interval Score. Additionally, we introduce a novel uncertainty quantification methodology designed for the onset of emerging epidemics. Developing versatile approaches that do not rely on historical data and can maintain high accuracy in the face of novel pandemics is critical for enhancing public health decision-making and strengthening preparedness for future outbreaks.

97 MATHEMATICS AND COMPUTING↗

Tensor Decompositions for Count Data that Leverage Stochastic and Deterministic Optimization

There is growing interest to extend low-rank matrix decompositions to multi-way arrays, or tensors. One fundamental low-rank tensor decomposition is the canonical polyadic decomposition (CPD). The challenge of fitting a low-rank, nonnegative CPD model to Poisson-distributed count data is of particular interest. Several popular algorithms use local search methods to approximate the global maximum likelihood estimator from local minima. Simultaneously, a recent trend in theoretical computer science and numerical linear algebra leverages randomization to solve very large, hard problems. The typical approach is to use randomization for a fast approximation and determinism for refinement to yield effective algorithms with theoretical guarantees. Two popular algorithms for Poisson CPD reflect that emergent dichotomy: CP Alternating Poisson Regression is a deterministic algorithm and Generalized Canonical Polyadic decomposition makes use of stochastic algorithms in several variants. This work extends recent work to develop two new methods that leverage randomized and deterministic algorithms for improved accuracy and performance.

97 MATHEMATICS AND COMPUTING↗

Acoustic Impedance Inversions for Offshore CO 2 Storage: South Georgia Embayment, United States

This study builds and expands on previous CO 2 storage resource assessment studies of the southeastern offshore Atlantic margin by providing a detailed evaluation on how rock porosities and permeabilities are distributed across the Upper Cretaceous strata restricted to the South Georgia Embayment (SGE). Using legacy industry two-dimensional seismic reflection and well data, this assessment is the first application of multiple seismic inversion techniques in this area. This workflow provides a reliable and repeatable model-based inversion which gives an improved image to discriminate lithology and predict porosity. The workflow is applicable to future CO 2 storage resource assessment studies elsewhere. The inversion results indicate that distinct porosity and permeability regimes are present and distributed in the Upper Cretaceous strata within the SGE. The impedance and porosity relationships show well-founded and reliable correlation. These relationships reveal low impedance coincident to the high porosity intervals which are proposed as potential reservoir intervals for CO 2 storage. In addition, the result shows that the Upper Cretaceous strata have two main potential reservoirs in the lower part. These are overlain by a thick impermeable interval, mostly shale, which has high impedance, low porosity, and low permeability and extends within the SGE. This result is in agreement with a previous study that also proposed two significant storage reservoirs for CO 2 in the Upper Cretaceous strata. Since porosity distribution is estimated using multiple methods, it follows the trends of seismic signature and structures of the Upper Cretaceous strata. The extracted values of porosity, ranging from 15 to 36%, and permeability, ranging from 1 to 100 mD, are close to the measured values from the well core data at the Upper Cretaceous strata interval.

54 ENVIRONMENTAL SCIENCES↗

Multiple-input experimental modal analysis

The development of experimental modal analysis techniques is reviewed. System and excitation assumptions are discussed. The methods examined include the forced normal mode excitation method, the frequency response function method, the damped complex exponential response method, the Ibrahim time domain approach, the polyreference approach, and mathematical input-output model methods. The current trend toward multiple input utilization in the estimation of system parameters is noted.

Allemang, R. J.↗