Search NASA⌕ Search

SEARCH · Search NASA

Results for “Inverse source location”

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 145 records · Page 8

Characterization of Transport Errors in Chemical Forecasts from a Global Tropospheric Chemical Transport Model

We propose a new methodology to characterize errors in the representation of transport processes in chemical transport models. We constrain the evaluation of a global three-dimensional chemical transport model (GEOS-CHEM) with an extended dataset of carbon monoxide (CO) concentrations obtained during the Transport and Chemical Evolution over the Pacific (TRACE-P) aircraft campaign. The TRACEP mission took place over the western Pacific, a region frequently impacted by continental outflow associated with different synoptic-scale weather systems (such as cold fronts) and deep convection, and thus provides a valuable dataset. for our analysis. Model simulations using both forecast and assimilated meteorology are examined. Background CO concentrations are computed as a function of latitude and altitude and subsequently subtracted from both the observed and the model datasets to focus on the ability of the model to simulate variability on a synoptic scale. Different sampling strategies (i.e., spatial displacement and smoothing) are applied along the flight tracks to search for systematic model biases. Statistical quantities such as correlation coefficient and centered root-mean-square difference are computed between the simulated and the observed fields and are further inter-compared using Taylor diagrams. We find no systematic bias in the model for the TRACE-P region when we consider the entire dataset (i.e., from the surface to 12 km ). This result indicates that the transport error in our model is globally unbiased, which has important implications for using the model to conduct inverse modeling studies. Using the First-Look assimilated meteorology only provides little improvement of the correlation, in comparison with the forecast meteorology. These general statements can be refined when the entire dataset is divided into different vertical domains, i.e., the lower troposphere (less than 2 km), the middle troposphere (2-6 km), and the upper troposphere (greater than 6 km). The best agreement between the observations and the model is found in the lower and middle troposphere. Downward displacements in the lower troposphere provide a better fit with the observed value, which could indicate a problem in the representation of boundary layer height in the model. Significant improvement is also found for downward and southward displacements in the upper troposphere. There are several potential sources of errors in our simulation of the continental outflow in the upper troposphere which could lead to such biases, including the location and/or the strength of deep convective cells as well as that of wildfires in Southeast Asia.

Bey, I.↗

Prominences: The Key to Understanding Solar Activity

Prominences are spectacular manifestations of both quiescent and eruptive solar activity. The largest examples can be seen with the naked eye during eclipses, making prominences among the first solar features to be described and catalogued. Steady improvements in temporal and spatial resolution from both ground- and space-based instruments have led us to recognize how complex and dynamic these majestic structures really are. Their distinguishing characteristics - cool knots and threads suspended in the hot corona, alignment along inversion lines in the photospheric magnetic field within highly sheared filament channels, and a tendency to disappear through eruption - offer vital clues as to their origin and dynamic evolution. Interpreting these clues has proven to be contentious, however, leading to fundamentally different models that address the basic questions: What is the magnetic structure supporting prominences, and how does so much cool, dense plasma appear in the corona? Despite centuries of increasingly detailed observations, the magnetic and plasma structures in prominences are poorly known. Routine measurements of the vector magnetic field in and around prominences have become possible only recently, while long-term monitoring of the underlying filament-channel formation process also remains scarce. The process responsible for prominence mass is equally difficult to establish, although we have long known that the chromosphere is the only plausible source. As I will discuss, however, the motions and locations of prominence material can be used to trace the coronal field, thus defining the magnetic origins of solar eruptions. A combination of observations, theory, and numerical modeling must be used to determine whether any of the competing theories accurately represents the physics of prominences. I will discuss the criteria for a successful prominence model, compare the leading models, and present in detail one promising, comprehensive scenario for prominence formation and evolution that could answer the two questions posed above.

Karpen, Judy T.↗

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗

Comparison of inversion models using AIRSAR data for Death Valley, California

Polarimetric Airborne Synthetic Aperture Radar (AIRSAR) data were collected for the Geologic Remote Sensing Field Experiment (GRSFE) over Death Valley, California, USA, in September 1989. AIRSAR is a four-look, quid-polarizaiton, three frequency instrument. It collects measurements at C-band (5.66 cm), L-band (23.98 cm), and P-band (68.13 cm), and has a GIFOV of 10 meters and a swath width of 12 kilometers. Because the radar measures at three wavelengths, different scales of surface roughness are measured. Also, dielectric constants can be calculated from the data. The scene used in this study is in Death Valley, California and is located over Trail Canyon alluvial fan, the valley floor, and Artists Drive alluvial fan. The fans are very different in mineralogic makeup, size, and surface roughness. Trail Canyon fan is located on the west side of the valley at the base of the Panamint Range and is a large fan with older areas of desert pavement and younger active channels. The source for the material on southern part of the fan is mostly quartzites and there is an area of carbonate source on the northern part of the fan. Artists Drive fan is located at the base of the Black Mountains on the east side of the valley and is a smaller, young fan with its source mostly from volcanic rocks. The valley floor contains playa and salt deposits that range from smooth to Devil's Golf course type salt pinnacles.

Kierein-Young, Kathryn S.↗

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

Abstract We develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Geochemistry & Geophysics↗

A Mesoscale Analysis of Column-Integrated Aerosol Properties in Northern India During the TIGERZ 2008 Pre-Monsoon Period and a Comparison to MODIS Retrievals

The Indo-Gangetic Plain (IGP) of the northern Indian subcontinent produces anthropogenic pollution from urban, industrial and rural combustion sources nearly continuously and is affected by convection-induced winds driving desert and alluvial dust into the atmosphere during the premonsoon period. Within the IGP, the NASA Aerosol Robotic Network (AERONET) project initiated the TIGERZ measurement campaign in May 2008 with an intensive operational period from May 1 to June 23, 2008. Mesoscale spatial variability of aerosol optical depth (AOD, tau) measurements at 500mn was assessed at sites around Kanpur, India, with averages ranging from 0.31 to 0.89 for spatial variability study (SVS) deployments. Sites located downwind from the city of Kanpur indicated slightly higher average aerosol optical depth (delta Tau(sub 500)=0.03-0.09). In addition, SVS AOD area-averages were compared to the long-tenn Kanpur AERONET site data: Four SVS area-averages were within +/- 1 cr of the climatological mean of the Kanpur site, while one SVS was within 2sigma below climatology. For a SVS case using AERONET inversions, the 440-870mn Angstrom exponent of approximately 0.38, the 440-870mn absorption Angstrom exponent (AAE) of 1.15-1.53, and the sphericity parameter near zero suggested the occurrence of large, strongly absorbing, non-spherical aerosols over Kanpur (e.g., mixed black carbon and dust) as well as stronger absorption downwind of Kanpur. Furthermore, the 3km and lOkm Terra and Aqua MODIS C005 aerosol retrieval algorithms at tau(sub 550) were compared to the TIGERZ data set. Although MODIS retrievals at higher quality levels were comparable to the MODIS retrieval uncertainty, the total number of MODIS matchups (N) were reduced with subsequent quality levels (N=25, QA>=0; N=9,QA>=l; N=6, QA>=2; N=1, QA=3) over Kanpur during the premonsoon primarily due to the semi-bright surface, complex aerosol mixture and cloud-contaminated pixels. The TIGERZ 2008 data set provided a unique opportunity to measure the spatial and temporal variations of aerosol loading in the IGP. The strong aerosol absorption derived from ground-based sun/sky radiometer measurements suggested the presence of a predominately black carbon and dust mixture during the pre-monsoon period. Consistent with the elevated heat-pump hypothesis, these absorbing aerosols found across Kanpur and the greater IGP region during the pre-monsoon period likely induced regional atmospheric warming, which lead to a more rapid advance of the southwest Asian monsoon and above normal precipitation over northern India in June 2008.

Giles, D. M.↗

Search for the 1/2 + intruder state in 35 P

The excitation energy of deformed intruder states (specifically the 2p2h bandhead) as a function of proton number Z along N = 20 is of interest both in terms of better understanding the evolution of nuclear structure between spherical Ca 40 and the Island of Inversion nuclei, and for benchmarking theoretical descriptions in this region. At the center of the N = 20 Island of Inversion, the n p n h (where n = 2 , 4 , 6 ) neutron excitations across a diminished N = 20 gap result in deformed and collective ground states, as observed in Mg 32 . In heavier isotones, n p n h excitations do not dominate in the ground states but are present in the relatively low-lying level schemes. With the aim of identifying the expected 2 p 2 h x s 1 / 2 + state in P 35 , the only N = 20 isotone for which the neutron 2p2h excitation bandhead has not yet been identified, the S 36 ( d , He 3 ) P 35 reaction has been revisited in inverse kinematics with the HELical Orbit Spectrometer (HELIOS) at the Argonne Tandem Linac Accelerator System (ATLAS). While a candidate state has not been located, an upper limit for the transfer reaction cross section to populate such a configuration within a 2.5 to 3.6 MeV energy range provides a stringent constraint on the wave function compositions in both S 36 and P 35 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermodynamic measurements in a high pressure hydrogen-oxygen flame using Raman scattering from a broadband excimer laser

Raman scattering is an inelastic molecular scattering process in which incident radiation is reemitted at a fixed change in frequency. Raman spectroscopy can be used to measure the number density and temperature of the irradiated species. The strength of the Raman signal is inversely proportional to the wavelength raised to the fourth power. Consequently, high signal to noise ratios are obtained by using ultraviolet (UV) excitation sources. Using UV sources for Raman Spectroscopy in flames is complicated by the fact that some of the primary constituents in hydrogen-oxygen combustion absorb and reemit light in the UV and these fluorescence processes interfere with the Raman signals. This problem has been handled in atmospheric pressure flames in some instances by using a narrowband tunable excimer laser as a source. This allows for detuning from absorption transitions and the elimination of interfering fluorescence signals at the Raman wavelengths. This approach works well in the atmospheric pressure flame; however, it has two important disadvantages. First, injection-locked narrowband tunable excimer lasers are very expensive. More importantly, however, is the fact that at the high pressures characteristic of rocket engine combustion chambers, the absorption transitions are broadened making it difficult to tune to a spectral location at which substantial absorption would not occur. The approach taken in this work is to separate the Raman signal from the fluorescence background by taking advantage of the fact that Raman signal has nonisotropic polarization characteristics while the fluorescence signals are unpolarized. Specifically, for scattering at right angles to the excitation beam path, the Raman signal is completely polarized. The Raman signal is separated from the fluorescence background by collecting both horizontally and vertically polarized signals separately. One of the polarizations has both the Raman signal and the fluorescence background while the other has only the fluorescence signal. The Raman scatter is the difference between the signals. By choosing an appropriate optical setup, both signals can be obtained simultaneously with the same monochromator; hence, time resolved measurements are possible using this approach.

Hartfield, Roy, Jr.↗

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point-source and extended-source studies rely on the modeling of this diffuse emission for accurate characterization. Here, we describe the development of the Galactic Interstellar Emission Model (GIEM),which is the standard adopted by the LAT Collaboration and is publicly available. This model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse-Compton emission produced in the Galaxy. In the GIEM, we also include large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20deg and we observe an enhanced emission toward their base extending in the north and south Galactic directions and located within approximately 4deg of the Galactic Center.

gamma rays: diffuse background – gamma rays: gen↗

Progress in Modeling Global Atmospheric CO2 Fluxes and Transport: Results from Simulations with Diurnal Fluxes

Progress in better determining CO2 sources and sinks will almost certainly rely on utilization of more extensive and intensive CO2 and related observations including those from satellite remote sensing. Use of advanced data requires improved modeling and analysis capability. Under NASA Carbon Cycle Science support we seek to develop and integrate improved formulations for 1) atmospheric transport, 2) terrestrial uptake and release, 3) biomass and 4) fossil fuel burning, and 5) observational data analysis including inverse calculations. The transport modeling is based on meteorological data assimilation analysis from the Goddard Modeling and Assimilation Office. Use of assimilated met data enables model comparison to CO2 and other observations across a wide range of scales of variability. In this presentation we focus on the short end of the temporal variability spectrum: hourly to synoptic to seasonal. Using CO2 fluxes at varying temporal resolution from the SIB 2 and CASA biosphere models, we examine the model's ability to simulate CO2 variability in comparison to observations at different times, locations, and altitudes. We find that the model can resolve much of the variability in the observations, although there are limits imposed by vertical resolution of boundary layer processes. The influence of key process representations is inferred. The high degree of fidelity in these simulations leads us to anticipate incorporation of realtime, highly resolved observations into a multiscale carbon cycle analysis system that will begin to bridge the gap between top-down and bottom-up flux estimation, which is a primary focus of NACP.

Collatz, G. James↗

GeV–TeV Counterparts of SS 433/W50 from Fermi -LAT and HAWC Observations

The extended jets of the microquasar SS 433 have been observed in optical, radio, X-ray, and recently very-high-energy γ-rays by High Altitude Water Cerenkov (HAWC). The detection of HAWC γ-rays with energies as great as 25 TeV motivates searches for high-energy γ-ray counterparts in the Fermi-Large Area Telescope (LAT) data in the 100 MeV–300 GeV band. In this paper, we report on the first joint analysis of Fermi-LAT and HAWC observations to study the spectrum and location of γ-ray emission from SS 433. Our analysis finds common emission sites of GeV-to-TeV γ-rays inside the eastern and western lobes of SS 433. The total flux above 1 GeV is $\sim 1\times {10}^{-10}\,{\mathrm{cm}}^{-2}\,{{\rm{s}}}^{-1}$ in both lobes. The γ-ray spectrum in the eastern lobe is consistent with inverse Compton emission by an electron population that is accelerated by jets. To explain both the GeV and TeV flux, the electrons need to have a soft intrinsic energy spectrum, or undergo a quick cooling process due to synchrotron radiation in a magnetized environment.

79 ASTRONOMY AND ASTROPHYSICS↗

Using Dark Fiber and Distributed Acoustic Sensing to Characterize a Geothermal System in the Imperial Valley, Southern California

The Imperial Valley, CA, is a tectonically active transtensional basin located south of the Salton Sea; the area hosts numerous geothermal fields, including significant hidden hydrothermal resources without surface manifestations. Development of inexpensive, rugged, and highly sensitive exploration techniques for undiscovered geothermal systems is critical for accelerating geothermal power deployment as well as unlocking a low-carbon energy future. We present a case study utilizing distributed acoustic sensing (DAS) and ambient noise interferometry for geothermal reservoir imaging, utilizing unlit fiber-optic telecommunication infrastructure (dark fiber). The study exploits two days of passive DAS data acquired in early November 2020 over a ~28-km section of fiber from Calipatria, CA to Imperial, CA. We apply ambient noise interferometry to retrieve coherent signals from DAS records and develop a bin stacking technique to attenuate the effects from persistent localized noise sources and to enhance retrieval of coherent surface waves. As a result, we are able to obtain high-resolution two-dimensional (2D) S wave velocity ($V_s$) structure to 3 km depth, based on joint inversion of both the fundamental and higher overtones. We observe a previously unmapped high $V_s$ and low $V_p$/$V_s$ ratio feature beneath the Brawley geothermal system, which we interpret to be a zone of hydrothermal mineralization and lower porosity. This interpretation is consistent with a host of other measurements including surface heat flow, gravity anomalies, and available borehole wireline data. These results demonstrate the potential utility of DAS deployed on dark fiber for geothermal system exploration and characterization in the appropriate geological settings.

Vp/Vs imaging↗

Development of Multi-Sensor Global Cloud and Radiance Composites for DSCOVR EPIC Imager with Subpixel Definition

The Deep Space Climate Observatory (DSCOVR) enables analysis of the daytime Earth radiation budget via the onboard Earth Polychromatic Imaging Camera (EPIC) and National Institute of Standards and Technology Advanced Radiometer (NISTAR). EPIC delivers adequate spatial resolution imagery but only in shortwave bands (317-780 nm), while NISTAR measures the top-of-atmosphere (TOA) whole-disk radiance in shortwave and longwave broadband windows. Accurate calculation of albedo and outgoing longwave flux requires a high-resolution scene identification such as the radiance observations and cloud properties retrievals from low earth orbit (LEO, including NASA Terra and Aqua MODIS, Suomi-NPP VIIRS, and NOAA AVHRR) and geosynchronous (GEO, including GOES east and west, METEOSAT, INSAT-3D, MTSAT-2, and Himawari-8) satellite imagers. The cloud properties are derived using the Clouds and the Earth's Radiant Energy System (CERES) mission Cloud Subsystem group algorithms. These properties have to be co-located with EPIC pixels to provide the scene identification and to select anisotropic directional models (ADMs), which are then used to adjust the NISTAR-measured radiance and subsequently obtain the global daytime shortwave and longwave fluxes. This work presents an algorithm for optimal merging of selected radiance and cloud property parameters derived from multiple satellite imagers to obtain seamless global hourly composites at 5-km resolution. Selection of satellite data for each 5-km pixel is based on an aggregated rating that incorporates five parameters: nominal satellite resolution, pixel time relative to the EPIC time, viewing zenith angle, distance from day/night terminator, and probability of sun glint. To provide a smoother transition in the merged output, in regions where candidate pixel data from two satellite sources have comparable aggregated rating, the selection decision is defined by the cumulative function of the normal distribution so that abrupt changes in the visual appearance of the composite data are avoided. Higher spatial accuracy in the composite product is achieved by using the inverse mapping with gradient search during reprojection and bicubic interpolation for pixel resampling.

Khlopenkov, Konstantin V.↗

Comparing Theoretically Scaled Biomechanical Models

BACKGROUND An investigation into incorporating space suit aspects required a 50th-percentile male model which was not part of the existing dataset of models. Previously, theoretical models for 5th-percentile female and 95th-percentile male were created utilizing a scaled test subject as close as possible to the target height and weight. Scaling factors were created by taking the height ratio and applying it uniformly to the model body segments, then fine tuning to ensure the theoretical model’s height is as expected. This study was initiated to examine the existing method of isometrically scaling in OpenSim [1,2] and to create alternative methods which do not rely on an existing subject being close in height and weight to the theoretical model of interest. Generating a theoretical biomechanical model provides additional abilities without the reliance on available OpenSim models or real subjects. Gained abilities include creating different percentile models and attaining representative anthropometry for any target height/weight, such as targeting specific crew populations or gaps within the current dataset. METHODS AND RESULTS This study includes 3 methods for generating the theoretical model; utilizing the modified unscaled OpenSim Full Body Rajagopal Model (FBRM) [3,4], a dataset of scaled OpenSim models, and existing scaling factors from Dumas et al. [5]. All methods use the Anthropometric Survey of US Army Personnel (ANSUR II) [6] as an input of necessary anthropometric measurements. The following measurements are retrieved from the ANSUR II collection: mass, stature, cervicale height, acromial height, axilla height, waist height, trochanterion height, lateral femoral epicondyle height, lateral malleolus height, acromion-radiale length, radiale-stylion length, palm length, ball of foot length, bicristal breadth, bimalleolar breadth. Some of the measurements are direct segment lengths and others are utilized to derive segment lengths. Height, weight, and age ranges are the required inputs to parse the collection for a mean value of the measurements. A simple iterative process may be necessary for the output of mean mass and stature consistent with the theoretical model of interest. In some cases, one may specify exact values instead of a range but that is dependent upon whether those exact values pertain to a single subject from the ANSUR II collection. Equations from Dumas et al. were used to calculate scaling factors for the two methods that utilize OpenSim scaled and unscaled models. In the third method, Dumas’ scaling factors were used directly instead of generating our own. After scaling factors are achieved, they are applied along with the mass and segment lengths retrieved from ANSUR II in order to calculate the segment mass, Center of Mass (CoM), mass Moment of Inertia (MoI), and the joint location in the parent frame. The theoretical models from the different methods were compared to the previous existing method of isometrically scaling in OpenSim, as well as the standards found in NASA STD-3000 [7] and the NASA Human Integration Design Handbook [8]. The whole-body center of mass is the main comparison performed between the models and the standards. The body segment mass properties were also compared when possible. In some cases, the standards do not provide complete data, or the center of mass location is not provided with respect to the appropriate joint coordinate system. For a 50th-percentile male subject, results from isometric scaling in OpenSim and the method utilizing the unscaled FBRM were compared to the standards in NASA STD-3000. According to the standard, the whole-body vertical CoM location was taken from the head vertex to the CoM. The standards provide 80.2 cm for the CoM vertical component whereas isometric scaling in OpenSim and utilizing the unscaled FBRM provide 82.2 and 82.3 cm, respectively. Through definition, the whole-body CoM in the medial/lateral direction agreed across the models. Investigation is on-going to clearly identify the reference point pertaining to dorsal/ventral whole-body CoM position in standards and relate that to the theoretical models. REFERENCES [1] Delp, S.L., et al., “OpenSim: Open-source Software to Create and Analyze Dynamic Simulations of Movement”, IEEE Transactions on Biomedical Engineering, (2007). [2] Seth, A., Hicks J.L., Uchida, T.K., Habib, A., Dembia, C.L., Dunne, J.J., Ong, C.F., DeMers, M.S., Rajagopal, A., Millard, M., Hamner, S.R., Arnold, E.M., Yong, J.R., Lakshmikanth, S.K., Sherman,n M.A., Delp, S.L. OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. Plos Computational Biology, 14(7). (2018) [3] R. K. Huffman, W. K. Thompson, C. A. Gallo, L. J. Quiocho, “Improvement of Scaling and Inverse Kinematic Results with Additional Upper Body Joints Added to Opensim Rajagopal Model”, NASA Human Research Program Investigator’s Workshop, (2019). [4] Huffman R.K., Thompson W., Gallo C., “Modified OpenSim Rajagopal Full Body Model”, New Technology Report, MSC-26872-1. (August 2020). [5] Dumas R., Cheze L. and Verriest J., 2007. “Adjustments to McConville et al. and Young et al. body segment inertial parameters”. Journal of Biomechanics 40, pp. 543-553 [6] Gordon C.C., 2012 “Anthropometric Survey of U.S. Army Personnel: Methods and Summary Statistics”, US Army Natick Soldier RD&E Center. [7] NASA STD-3000/REV-B, 1995. “The Man-System Integration Standards” [8] NASA/SP-2010-3407/REV-1, 2014. “The Human Integration Design Handbook (HIDH)”

F N Matari↗

Quantifying the impact of temporal analysis of products reactor initial state uncertainties on kinetic parameters

The temporal analysis of products (TAP) reactor, a transient kinetic tool, provides users with information as the catalyst state evolves. However, the state of the art for TAP uncertainty quantification only considers the experimental noise present in the outlet flow signal. Additional sources of uncertainty, including, but not limited to, initial surface coverages, the catalyst zone location, the inert void fraction, and the gas pulse intensity and pulse delay, are not included. For this reason, a framework for quantifying all uncertainty sources present in TAP experiments is presented and applied to a carbon monoxide oxidation case study. Herein, two methods for quantifying these sources of uncertainty are introduced. The first utilizes initial state sensitivities to approximate the parameter variances, as well as to provide insights into the structural certainty of the model. The second generates parameter confidence distributions through an ensemble-based sampling algorithm. This initial state covariance matrix can ultimately be merged with the experimental noise covariance matrix, providing a unified description of the parameter uncertainties for a single TAP experiment.

36 MATERIALS SCIENCE↗

Ammonia downstream from HH 80 North

HH 80-81 are two optically visible Herbig-Haro (HH) objects located about 5 minutes south of their exciting source IRAS 18162-2048. Displaced symmetrically to the north of this luminous IRAS source, a possible HH counterpart was recently detected as a radio continuum source with the very large array (VLA). This radio source, HH 80 North, has been proposed to be a member of the Herbig-Haro class since its centimeter flux density, angular size, spectral index, and morphology are all similar to those of HH 80. However, no object has been detected at optical wavelengths at the position of HH 80 North, possibly because of high extinction, and the confirmation of the radio continuum source as an HH object has not been possible. In the prototypical Herbig-Haro objects HH 1 and 2, ammonia emission has been detected downstream of the flow in both objects. This detection has been intepreted as a result of an enhancement in the ammonia emission produced by the radiation field of the shock associated with the HH object. In this Letter we report the detection of the (1,1) and (2,2) inversion transitions of ammonia downstream HH 80 North. This detection gives strong suppport to the interpretation of HH 80 North as a heavily obscured HH object. In addition, we suggest that ammonia emission may be a tracer of embedded Herbig-Haro objects in other regions of star formation. A 60 micrometer IRAS source could be associated with HH 80 North and with the ammonia condensation. A tentative explanation for the far-infrared emission as arising in dust heated by their optical and UV radiation of the HH object is presented.

Girart, Jose M.↗

Characterizing Baselines and Change in Gas Hydrate Systems using EM Methods

The objective of this project was to advance our understanding of gas hydrate systems in nature by characterizing their electrical properties in the field and in the laboratory. In the laboratory measurements, methane hydrate was synthesized from pure water ice and flash frozen seawater, with varying amounts of sand or silt added. Electrical conductivity was determined by impedance spectroscopy, using equivalent circuit modeling to separate the effects of electrodes and to gain insight into conduction mechanisms. Silt and sand increase the conductivity of pure hydrate, inferred to be contaminant NaCl contributing to conduction in hydrate, to a peak conductivity in agreement with peak resistivities observed in well logs through massive hydrate (3,000--10,000 Ωm). The addition of silt and sand lowers the conductivity of hydrate synthesized from seawater, by an amount consistent with Archie's Law. All samples were characterized using cryogenic scanning electron microscopy and energy dispersive spectroscopy, which shows good connectivity of salt and brine phases. Electrical conductivity measurements of pure hydrate and hydrate mixed with silt during pressure-induced dissociation supports previous conclusions that sediment increases dissociation rate. In order to characterize gas hydrate systems in the field, we collected 360 line kilometers of controlled-source electromagnetic data on Walker Ridge 313, Orca Basin (WR100), Mad Dog (GC781), and Green Canyon 955 in the Gulf of Mexico, all areas with known or seismically inferred gas hydrate deposits and which have be drilled or targeted for future drilling. We deep-towed an EM transmitter that generates an alternating electric field which propagates through the seafloor geology. Data were recorded on 6 receivers towed behind the transmitter at distances between 550 and 1550 m. In the presence of conductive geology, the electric fields will be attenuated, and conversely, in resistive geology the fields will be preserved. Our data were inverted using a 2D inversion method that first optimizes the model-data misfit, then finds the smoothest model fitting the data. This ensures that resistivity structures present in the final model are likely necessary. At each of the proposed drilling sites we found increased resistivity, interpreted as increased hydrate concentrations. However, not only were the primary drilling sites not always more resistive than the alternate sites, at WR313 the strongest resistors were not at the locations targeted for drilling.

03 NATURAL GAS↗

Inferring the Focal Depths of Small Earthquakes in Southern California Using Physics-Based Waveform Features

Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ~8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2 < M L < 4) recorded at local distances (<150 km), for which depth estimation is especially challenging. We find that differential magnitudes (M w /M L –M c ) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡log 10 (M 0 )+3log 10 (f c ), and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. Here, we also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

58 GEOSCIENCES↗