Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface 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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Influence of Time and Frequency Domain Wave Forcing on the Power Estimation of a Wave Energy Converter Array

Industry-specific tools for analyzing and optimizing the design of wave energy converters (WECs) and associated power systems are essential to advancing marine renewable energy. This study aims to quantify the influence of phase information on the device power output of a virtual WEC array. We run the phase-resolving wave model FUNWAVE-TVD (Total Variation Diminishing) to generate directional waves at the PacWave South site offshore from Newport, Oregon, where future WECs are expected to be installed for testing. The two broad cases presented correspond to mean wave climates during warm months (March–August) and cold months (September–February). FUNWAVE-TVD time series of sea-surface elevation are then used in WEC-Sim, a time domain numerical model, to simulate the hydrodynamic response of each device in the array and estimate their power output. For comparison, WEC-Sim is also run with wave energy spectra calculated from the FUNWAVE-TVD simulations, which do not retain phase information, and with wave spectra computed using the phase-averaged model Simulating WAves Nearshore (SWAN). The use of spectral data in WEC-Sim requires a conversion from frequency to time domain by means of random superposition of wave components, which are not necessarily consistent because of the linear assumption implicit in this method. Thus, power response is characterized by multiple realizations of the wave climates.

13 HYDRO ENERGY↗
DOE OSTI2020

Sweeping Jet Film Cooling at High Blowing Ratio on a Turbine Vane

Abstract Experimental and numerical investigations were performed to study the effects of high blowing ratios and high freestream turbulence on sweeping jet film cooling. Experiments were conducted on a nozzle guide vane suction surface in a low-speed linear cascade at a range of blowing ratios of 0.5–3.5 and freestream turbulence of 0.6% and 14.3%. Infrared thermography was used to estimate the adiabatic cooling effectiveness. Thermal field and boundary layer measurements were conducted at a cross-plane at x/D = 12 downstream of the hole exit. Results were compared with a baseline 777-shaped hole and showed that the sweeping jet hole has improved cooling effectiveness at high blowing ratios (M > 1). The thermal field data revealed that the coolant separates from the surface at high blowing ratios for the 777-shaped hole while the coolant remains attached for the sweeping jet hole. Boundary layer measurements further confirmed that due to the sweeping motion of the jet, the effective jet momentum of the sweeping jet hole remains much lower than that of a 777-shaped hole. Thus, the coolant remains closer to the wall even at high blowing ratios. Large Eddy simulations (LES) were performed for both the sweeping jet and 777-shaped hole to evaluate the interaction between the coolant jet and the freestream in the near hole regions. Results showed that 777-shaped hole has a strong jetting at high blowing ratio that originates inside the hole breakout edges thus causing the jet to blow-off from the surface. In contrast, the sweeping jet does not show this behavior due to its internal geometry and flapping motion of the jet.

Engineering↗
DOE OSTI2022

Probing the initial stages of iron surface corrosion: Effect of O 2 and H 2 O on surface carbonation

Iron plays a vital role in natural processes such as water, mineral, iron, and nutrient cycles. Iron undergoes reduction-oxidation and catalytic reactions to produce various corrosion films depending on its chemical environment. Near ambient pressure X-ray photoelectron spectroscopy, polarized modulated infrared reflection absorption spectroscopy, and Auger electron spectroscopy were used to study the key reactants, from O 2 (g), H 2 O vapor, Na + and Cl – on the initial stages of iron surface corrosion. Here, with increasing the ratio of O 2 and H 2 O, surface hydrocarbons were shown to oxidize into carbonates, while the Cl – was found to migrate into the interface. The effect of each individual reactant was measured separately and water was shown to have a first order rate dependence on the carbonate growth at low pressures, with little dependence for O 2 . Near ambient pressures, both H 2 O and O 2 were found to increase the carbonate growth, which was estimated using the Langmuir isotherm model, yielding Gibbs energies between –9.8 to –8.5 kJ/mol. A mechanism is suggested to explain the oxidation is catalyzed by NaCl on iron surfaces and the adventitious hydrocarbons served as the source for surface carbonation. These findings have implications for understanding other surface catalytic and redox interface chemistry in complex environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Discerning structure sensitivity of Ni-core Pd-shell nanoparticles for enhanced nitrite reduction

Denitrification using Pd-based catalysts converts N-oxyanions into harmless products. However, the metal cost and scarcity are barriers to application. In this work, we synthesize water-suspended, structure-controlled Ni-core/Pd-shell nanoparticles (Ni@Pd NPs) that show Pd-catalyzed nitrite (NO 2 − ) reduction catalysis promoted by Ni. The NPs immobilized on Al 2 O 3 (0.47 wt % Pd and 0.74 wt % Ni) have higher NH 4 + selectivity and higher catalytic activity than wet-impregnated Pd/Al 2 O 3 (1.0 wt % Pd), i.e., 315.4 vs. 132.1 Lg surface Pd −1 min −1 . X-ray photoelectron spectroscopy (XPS) and CO-diffuse reflectance infrared Fourier transform spectroscopy (CO-DRIFTS) analyses show that Ni increases Pd electron density, which is proposed to lower the bond dissociation energy barrier for nitric oxide surface intermediate. Low-coordinated Pd atoms correlate to high activity and low NH 4 + selectivity. Using CatCost software, we estimate Ni@Pd/Al 2 O 3 to have a ∼40% lower manufacturing cost compared to wet-impregnated Pd/Al 2 O 3 . This rationally designed catalyst illustrates how nickel promotes Pd catalysis and how Pd usage can be lowered for denitrification and other hydrogenation reactions.

CO-DRIFTS↗
DOE OSTI2022

Convective Momentum Transport and Its Impact on the Madden‐Julian Oscillation in E3SM‐MMF

Abstract Convective momentum transport (CMT) is the process of vertical redistribution of horizontal momentum by small‐scale turbulent flows from moist convection. Traditional general circulation models (GCMs) and their multiscale modeling framework (MMF) versions poorly represent CMT due to insufficient information of subgrid‐scale flows at each GCM grid. Here the explicit scalar momentum transport (ESMT) scheme for representing CMT is implemented in the Energy Exascale Earth System Model‐Multiscale Modeling Framework (E3SM‐MMF) with embedded 2‐D cloud‐resolving models (CRMs), and verified against E3SM‐MMF simulations with 3‐D CRMs and observations. The results show that representing CMT by ESMT helps reduce climatological mean precipitation model bias over the western Pacific and the ITCZ regions, which is attributed to the weakened mean easterlies over the Pacific. Also, CMT from simulations with 2‐D and 3‐D CRMs impose a similar impact on Kelvin waves by reducing their variability and slowing down their phase speed, but opposite impacts on the Madden‐Julian Oscillation (MJO) variability. The ESMT scheme readily captures the climatological mean spatial patterns of the zonal and meridional components of CMT and their variability across multiple time scales, but shows some differences in estimating its magnitude. CMT mainly affects the MJO by decelerating its winds in the free troposphere, but accelerates its near‐surface winds. This study serves as a prototype for implementing CMT scheme in the MMF simulations, highlighting its crucial role in reducing model bias in mean state and spatiotemporal variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The advancing wave front on a sloping channel covered by a rod canopy following an instantaneous dam break

The drag coefficient Cd for a rigid and uniformly distributed rod canopy covering a sloping channel following the instantaneous collapse of a dam was examined using flume experiments. The measurements included space x and time t high resolution images of the water surface h(x, t) for multiple channel bed slopes So and water depths behind the dam Ho along with drag estimates provided by sequential load cells. Using these data, an analysis of the Saint-Venant equation (SVE) for the front speed was conducted using the diffusive wave approximation. An inferred Cd=0.4 from the h(x, t) data near the advancing front region, also confirmed by load cell measurements, is much reduced relative to its independently measured steady-uniform flow case. This finding suggests that drag reduction mechanisms associated with transients and flow disturbances are more likely to play a dominant role when compared to conventional sheltering or blocking effects on Cd examined in uniform flow. The increased air volume entrained into the advancing wave front region as determined from an inflow–outflow volume balance partly explains the Cd reduction from unity.

Mechanics↗
DOE OSTI2022

Ac losses in field-cooled type I superconducting cavities

As superconductors are cooled below their critical temperature, stray magnetic flux can become trapped in regions that remain normal. The presence of trapped flux facilitates dissipation of ac current in a superconductor, leading to losses in superconducting elements of microwave devices. In type II superconductors, dissipation is well-understood in terms of the dynamics of vortices hosting a single flux quantum. In contrast, the ac response of type I superconductors with trapped flux has not received much attention. Building on Andreev's early work (Andreev 1967 Sov. Phys. JETP 24 1019), here we show theoretically that the dominant dissipation mechanism is the absorption of the ac field at the exposed surfaces of the normal regions, while the deformation of the superconducting/normal interfaces is unimportant. We use the developed theory to estimate the degradation of the quality factors in field-cooled cavities, and we satisfactorily compare these theoretical estimates to the measured field dependence of the quality factors of two aluminum cavities. We also identify a regime in which the dissipated power depends weakly on the Ginzburg-Landau parameter; this makes it possible to apply our findings to cavities made of other materials, such as niobium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗
DOE OSTI2019

Microwave Meissner screening properties of proximity-coupled topological-insulator / superconductor bilayers

The proximity-coupled topological-insulator / superconductor (TI/SC) bilayer system is a representative system to realize topological superconductivity. In order to better understand this unique state and enable future applications of the TI/SC bilayer, a comprehensive characterization and understanding of the microscopic properties of the bilayer are required. In this work, a microwave Meissner screening study, which exploits a high-precision microwave resonator technique, is conducted on the SmB 6 /YB 6 thin-film bilayers as an example TI/SC system. The study reveals spatially dependent electrodynamic screening response of the TI/SC system that is not accessible to other techniques, from which the corresponding microscopic properties of a TI/SC bilayer can be obtained. The TI thickness dependence of the effective penetration depth suggests the existence of a bulk insulating region in the TI layer. The spatially dependent electrodynamic screening model analysis provides an estimate for the characteristic lengths of the TI/SC bilayer: normal penetration depth, normal coherence length, and the thickness of the surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A Method for Dimensionally Adaptive Sparse Trigonometric Interpolation of Periodic Functions

We present a method for dimensionally adaptive sparse trigonometric interpolation of multidimensional periodic functions belonging to a smoothness class of finite order. This method targets applications where periodicity must be preserved and the precise anisotropy is not known a priori. To the authors' knowledge, this is the first instance of a dimensionally adaptive sparse interpolation algorithm that uses a trigonometric interpolation basis. The motivating application behind this work is the adaptive approximation of a multi-input model for a molecular potential energy surface (PES) where each input represents an angle of rotation. Our method is based on an anisotropic quasi-optimal estimate for the decay rate of the Fourier coefficients of the model; a least-squares fit to the coefficients of the interpolant is used to estimate the anisotropy. Thus, our adaptive approximation strategy begins with a coarse isotropic interpolant, which is gradually refined using the estimated anisotropic rates. The procedure takes several iterations where ever-more accurate interpolants are used to generate ever-improving anisotropy rates. We present several numerical examples of our algorithm where the adaptive procedure successfully recovers the theoretical “best” convergence rate, including an application to a periodic PES approximation. An open-source implementation of our algorithm resides in the Tasmanian UQ library developed at Oak Ridge National Laboratory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Steady-State Ultrasonic Non-Destructive Evaluation for Aerospace & Additive Manufacturing Applications [Slides]

How does ultrasonic inspection work? Rich phenomenology of interactions between ultrasonic waves and features in plate-like structures. Hidden defects/discontinuities create measurable changes in the surface wavefield. Based on these measurements, we can build tools that detect the presence of different kinds of defects and estimate structural/material properties (e.g. local thickness).

42 ENGINEERING↗
DOE OSTI2025

Swept-Lookback Deflectometry for High Performance Concentrating Solar Power Optical Metrology

This report describes an initial investigation into a proposed solution to the important problem of performing a detailed evaluation of heliostat optical performance, in situ in a heliostat field. Our approach is to place digital cameras in a position near the receiver where they look back toward the heliostat mirrors. The pixels of each camera sensor identify a set of small cells on the mirror surface, each corresponding to a “mixel.” By either passing reflected sunbeam over the camera or passing the camera through the reflected sunbeam, the cameras intercept sunlight reflected from each mixel. We then analyze the recorded video data to determine times when each mixel transitions from dark to light, and then back to dark. We then use these transitions to construct vectors from the camera to the mixel, and then from the mixel to the edge of the Sun at that moment. We then compute the surface normal at the mixel, which bisects the angle between these vectors. Performing this analysis for all mixels in the mirror yields a high-resolution map of slope across the mirror surface. We have implemented most of this process, successfully collecting data for an example heliostat facet and computing a preliminary estimated slope map. However, more work remains to complete this calculation, since certain factors and transformations are not yet included. Our observations so far support our hypothesis that such a system is possible, but we have not yet completed our quantitative evaluation of the concept.

14 SOLAR ENERGY↗
DOE OSTI2020

Structural, Criticality, and Radiation Dose Calculations to support SNF Loading into a DOE Standard Canister

The DOE Standard Canister Demonstration Project includes the development of an internal support structure (ISS) for the 4.6-m long and 45.7-cm diameter canister. This study presents structural evaluations, criticality safety assessments, and dose rate calculations conducted within the scope of the ISS design process to support smooth canister loading operations and to ensure safe storage, transportation, and disposal of Peach Bottom 1 Core II (PB2) and Fort St. Vrain (FSV) SNF, currently stored at the CPP-603 facility of the INL. The ISS includes a 316L stainless steel basket that holds twelve PB2 SNF rods. Further, six individual, 78.7-cm long, 316L stainless steel columns are equally spaced and welded to the inner canister wall at the lower end of the shell. These columns represent the FSV basket and can hold one FSV SNF element. An A92014 T6 aluminum spacer disc is bolted to the bottom plate of the PB2 basket to vertically restrain the FSV SNF element after the PB2 basket is placed inside the canister on top of the FSV basket. The structural evaluations of the ISS followed applicable ASME BPVC.III.3 guidelines and included finite element (FE) analyses of the PB2 basket structure; analyses of welds and bolds; buckling analyses of selected components, and acceptability assessments of the expected basket deformations under loading operations. The criticality safety assessments used the Monte Carlo N-Particle (MCNP) software architecture version 6.2 including ENDF/B-V continuous energy cross-section libraries, considering intact SNF in a single storage overpack or two multi-storage overpack configurations, and intact or failed SNF configured for disposal. The dose rate computations are based on source terms taken from the DOE Spent Fuel Database. The isotopic composition was decay corrected for the year 2022 using the ORIGEN module in the SCALE suite. A 19-group photon spectrum and 27-group neutron-source spectra were generated and used in MCNP to calculate estimated dose-equivalent rates, both on DOE Standard Canister contact and at a radial distance of from the canister surface. The results of this study indicate a structurally sound system that can uphold its criticality safety functions throughout its intended operational phases. Further, they increase confidence that sufficient radiological protection is technically achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Insulation Activation and Contamination in Low-Power Experiments

Among the diverse advanced nuclear reactor concepts being developed are designs that expect to operate at much higher temperatures than conventional light-water reactors, and therefore require alternative materials for components capable of withstanding these extreme temperatures. Insulation is one such component where alternative high-temperature industrial composites are being considered. Low-power experiments can provide an opportunity for advancing our understanding of material behavior under irradiation; however, these experiments are not without occupational hazards. The focus of this paper is the potential for activation of insulating composites when exposed to neutron fluence during reactor experiments, and to discuss the subsequent contamination of the testing area. This paper will also discuss practical measures for preventing the creation and inhalation of activated particulate matter (i.e., dust), such as engineering and administrative controls and personal protective equipment. This paper evaluates the suitability of several insulation materials for use in reactor experiments based on their activation from neutron irradiation. The materials considered are Pyrogel XTE by Aspen Aerogel, Cerablanket and Kaowool by Morgan Advanced Materials, Thermo-12 Gold by Johns Manville, and Maxsil CF6-2000 by McAllister Mills. Although these materials are not yet widely used in the nuclear industry, they are common in other industries due to their thermal insulating properties. These materials produce dust while handling that can pose health and safety hazards, which could be further exacerbated by activation. The analysis presented in this paper utilizes a deterministic method developed to predict the activation source term of the provided insulation materials. The method calculates the atomic density of each element within the composites and simulates their exposure to a specified neutron flux density representative of low power experimental conditions. During the irradiation phase and throughout the extensive cooling period, neutron absorptions and nuclear decays are tracked for relevant isotopes. The residual activity within each composite is then assessed to identify the material with the lowest potential radiological hazard to workers. Dust will be created during installation and removal of insulating materials and may also be generated during the experiment. Contamination of surrounding surfaces due to dispersion is highly likely unless preventative measures are taken. The results of this analysis are used to estimate an internal dose from inhalation of activated dust, which is significantly lower than the allowable whole-body dose. However, the risk of internal exposure may increase if the dust becomes a transport vector for nuclear material. The goals of this paper are to enhance the understanding of the potential activation and contamination hazards from the use of these insulating materials in low-power nuclear experiments and to foster safer work environments through informed decision-making and strategic planning. This work ultimately contributes to a broader understanding of how new non-nuclear materials for advanced reactors impact the radiological safety of co-located workers, helping to ensure the responsible management of radioactive materials during installation and decommissioning activities.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2020

Comparison of total emitted solar-induced chlorophyll fluorescence (SIF) and top-of-canopy (TOC) SIF in estimating photosynthesis

We report many studies have shown that solar-induced fluorescence (SIF) has a good potential to predict gross primary production (GPP) of vegetation. What we measured by remote sensing or near-surface platforms is top-of-canopy (TOC) SIF (SIF toc ), which is not necessarily equal or proportional to total emitted SIF (SIF tot ) from the entire canopy due to the (re)absorption and scattering effects. However, photosynthesis, the process that plants use to fix carbon from the atmosphere, occurs at the entire vertical canopy. Here, by using the recollision theory, we calculated SIF tot at 760 nm from the measured SIF toc , hyperspectral reflectance (R), canopy interception (i) and leaf albedo (ω l ). Among them, both SIF toc and R can be obtained from concurrent TOC spectral measurements; i and ω l in the near-infrared region can be estimated from the open access datasets with a good accuracy. Our result confirms that the measured SIF toc only accounts for a small fraction of SIF tot : the above-the-canopy sensor can only “see” on average 22.9% of SIF tot at Harvard Forest. SIF tot has the following advantages over SIF toc in estimating GPP: (1) SIF tot improves the diurnal estimate of canopy GPP, especially capable to capture the midday depression of photosynthesis which may cause the large discrepancies between SIF toc and GPP on a diurnal basis, (2) SIF tot produces a stronger correlation with GPP from plants with complex canopy structure or under sky conditions with more diffuse irradiance, and (3) the SIF tot -GPP relationship shows a stronger resilience to environmental stresses. The fluorescence escape ratio (f esc ), the ratio between SIF toc multiplied by π and SIF tot , is mostly determined by the sun-canopy-sensor geometry and leaf inclination distribution. The effect of LAI and the leaf chlorophyll concentration on f esc is marginal at the 760 nm wavelength. Our results suggest that converting SIF toc into SIF tot provides a better tool to understand and estimate GPP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Neural Density Estimation and Uncertainty Quantification for ChemCam Spectra [Slides]

The ChemCam instrument of Curiosity uses laser-induced breakdown spectroscopy (LIBS). It fires a laser at target and vaporizes rock surfaces, creating a plasma. Three spectrographs divide the plasma light into wavelengths for chemical analysis: ultraviolet, violet, and visible near-infrared. Regression methods (SVR, PCR, CNN) have been employed for calibration (prediction of the elemental composition of samples); however, labeled ChemCam samples are limited. Here, we focus on unsupervised learning and employ generative models from ChemCam analysis. Further, we use labels (supervised) in combination to the generative model to compute uncertainties related to predictions. We report generative modeling can be successfully applied to model real-world data. Normalizing flow models can be efficiently constructed on latent spaces for fast downstream inference. Unsupervised and supervised learning can be combined to form an uncertainty quantification framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Bayesian Inference for the Seismic Moment Tensor Using Regional Waveforms and Teleseismic- P Polarities with a Data-Derived Distribution of Velocity Models and Source Locations

The largest source of uncertainty in any source inversion is the velocity model used in the transfer function that relates observed ground motion to the seismic moment tensor. However, standard inverse procedure often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. Here, we incorporate this uncertainty into an estimation of the seismic moment tensor using a data-derived distribution of velocity models based on complementary geophysical data sets, including thickness constraints, velocity profiles, gravity data, surface-wave group velocities, and regional body-wave travel times. The data-derived distribution of velocity models is then used as a prior distribution of Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional and teleseismic-P waveforms. The use of multiple data sets is important for gaining resolution to different components of the moment tensor. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and interpreted in terms of the most probable source type.

58 GEOSCIENCES↗