Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local regression”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Observations of Coastal Wind Momentum Flux: Dependence on Fetch and Waves with Comparisons to COARE

Here, using observations from the Martha’s Vineyard Coastal Observatory, this paper investigates how momentum flux in the marine atmospheric surface layer over the coastal ocean varies with fetch, wave age, and wave slope and assesses the performance of the COARE 3.5 bulk parameterization. Long-fetch (at least 300 km) and short-fetch (3–6 km away from land) conditions have very similar momentum flux, with the latter being just 15% higher. The COARE 3.5 wind speed–dependent formulation closely matches the observations. The sea state dependence of wave age and wave slope is analyzed by considering both peak frequency and mean frequency in a wave spectrum. The slope of regression lines between normalized roughness and wave age is sensitive to wind speed ranges and the scatter of momentum flux, potentially explaining why earlier studies did not find a universal formula to characterize the wave dependence. Although the observed momentum flux exhibits an obvious dependence on wave age, a robust quadratic fit between momentum flux and the 10-m neutral wind speed exists only for young waves. For the short-fetch conditions, the momentum flux does not increase as wave slope increases. This may be because waves generated by local wind are still weak and swell that propagates from other areas dominates the wave spectra. In other words, the wave slope computed by considering the spectra does not properly reflect the local wind–wave interaction.

Air-sea interaction↗
NASA NTRS2009

Simulation of RCC Crack Growth Due to Carbon Oxidation in High-Temperature Gas Environments

The carbon wall oxidation technique coupled with a CFD technique was employed to study the flow in the expanding crack channel caused by the oxidation of the channel carbon walls. The recessing 3D surface morphing procedure was developed and tested in comparison with the arcjet experimental results. The multi-block structured adaptive meshing was used to model the computational domain changes due to the wall recession. Wall regression rates for a reinforced carbon-carbon (RCC) samples, that were tested in a high enthalpy arcjet environment, were computationally obtained and used to assess the channel expansion. The test geometry and flow conditions render the flow regime as the transitional to continuum, therefore Navier-Stokes gas dynamic approach with the temperature jump and velocity slip correction to the boundary conditions was used. The modeled mechanism for wall material loss was atomic oxygen reaction with bare carbon. The predicted channel growth was found to agree with arcjet observations. Local gas flow field results were found to affect the oxidation rate in a manner that cannot be predicted by previous mass loss correlations. The method holds promise for future modeling of materials gas-dynamic interactions for hypersonic flight.

Titov, E. V.↗
DOE OSTI2023

Validating the Use of Gaussian Process Regression for Adaptive Mapping of Residual Stress Fields

Probing the stress state using a high density of measurement points is time intensive and presents a limitation for what is experimentally feasible. Alternatively, individual strain fields used for determining stresses can be reconstructed from a subset of points using a Gaussian process regression (GPR). Results presented in this paper evidence that determining stresses from reconstructed strain fields is a viable approach for reducing the number of measurements needed to fully sample a component’s stress state. The approach was demonstrated by reconstructing the stress fields in wire-arc additively manufactured walls fabricated using either a mild steel or low-temperature transition feedstock. Effects of errors in individual GP reconstructed strain maps and how these errors propagate to the final stress maps were assessed. Implications of the initial sampling approach and how localized strains affect convergence are explored to give guidance on how best to implement a dynamic sampling experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamical Downscaling of Earth System Model Data for Energy System Analysis

Assessing energy resources (e.g., solar, wind, and hydro) under future scenarios requires datasets with sufficient spatial and temporal detail to capture variability and extreme events. While global-scale Earth System Model (ESM) projections are widely used, their coarse resolution limits direct application to regional energy system analyses. Dynamical downscaling offers a robust approach to generate physically consistent, fine-scale datasets that better represent local atmospheric processes impacting energy resources. In this work, we present a two-stage approach for producing high-resolution historical and future projections over the contiguous United States (CONUS). First, we optimize the Weather Research and Forecasting (WRF) model configuration for energy-relevant variables - solar irradiance, wind speed, and precipitation - by conducting ERA5-driven simulations at 8-km and 28-km resolution. Multiple physics schemes and model configurations within the WRF are evaluated against observational datasets including the National Solar Radiation Database (NSRDB), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the Stage IV multi-radar/multi-sensor precipitation product for the CONUS domain. Using the best-performing configuration, we dynamically downscale MPI-ESM1-2-HR simulations for 2000-2060 under SSP2-4.5 and SSP5-8.5 scenarios at 4-km spatial and hourly temporal resolution. This presentation will provide a comprehensive analysis of the results from multiple numerical experiments and high-resolution ESM projections. In addition, we will discuss potential applications of our high-resolution datasets within the energy sector and outline future research avenues dedicated to evaluating how extreme weather events influence system performance and resilience.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Machine learning bridges local static structure with multiple properties in metallic glasses

A long-standing challenge in the metallic glass (MG) community has been how to quantitatively gauge the influence of the intricate local packing environment on the response (such as the propensity for atomic rearrangement) of the atomic configuration to external stimuli. Here we establish this structure–property relation by representing the complex amorphous structure using a single, flexibility-orientated structural quantity. This structural flexibility (SF) couples to a bona fide structural representation, the pair distribution function (PDF) of individual atoms, through a weighting function that reflects what matters in the static atomic configuration to dynamic responses. Machine learning is used, employing microscopic flexibility volume as the supervisory signal, to establish via direct regression an optimized weighting vector, which is proven robust for all quenching rates, deformation conditions, and different compositions in a given (e.g., Cu x Zr 100-x ) alloy system. Additionally, the SF is evaluated solely from the particle positions (PDF), for any structure variation, from the atomic scale up to sample average. Strong correlations are demonstrated between SF and a broad range of properties, including vibrational, diffusional, as well as elastic and plastic relaxation responses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of quantitative synchrotron radiation micro-X-ray fluorescence in rice grain

Concentrations of nutrients and contaminants in rice grain affect human health, specifically through the localization and chemical form of elements. Methods to spatially quantify the concentration and speciation of elements are needed to protect human health and characterize elemental homeostasis in plants. Here, an evaluation was carried out using quantitative synchrotron radiation microprobe X-ray fluorescence (SR-µXRF) imaging by comparing average rice grain concentrations of As, Cu, K, Mn, P, S and Zn measured with rice grain concentrations from acid digestion and ICP-MS analysis for 50 grain samples. Better agreement was found between the two methods for high-Z elements. Regression fits between the two methods allowed quantitative concentration maps of the measured elements. These maps revealed that most elements were concentrated in the bran, although S and Zn permeated into the endosperm. Arsenic was highest in the ovular vascular trace (OVT), with concentrations approaching 100 mg kg –1 in the OVT of a grain from a rice plant grown in As-contaminated soil. Quantitative SR-µXRF is a useful approach for comparison across multiple studies but requires careful consideration of sample preparation and beamline characteristics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

Effects of 40Ar and 56Fe ions on retinal photoreceptor cells of the rabbit: implications for manned missions to Mars

Losses of photoreceptor cells (rods) from the retinas of New Zealand white (NZW) rabbits were detectable within 2 years after localized acute irradiation of optic and proximal tissues with > or = 7 Gy of 530 MeV u-1 40Ar ions or > or = 2 Gy of 465 MeV u-1 56Fe ions in the Bragg plateau region of energy deposition. Those limits were determined only from an analysis of variance of dose groups because the shapes of the dose response curves at early post-irradiation times are not known, a concern being addressed by experiments in progress. Losses of photoreceptor cells for the period 0.5-2.5 years post-irradiation, determined by provisional linear regression analysis, were approximately 1.7% Gy-1 and 2.5% Gy-1 for 40Ar and 56Fe ions, respectively.

NASA Discipline Radiation Health↗
NASA NTRSDate not supplied

The INSTEP Monitoring Network: Merging High-and-Low Cost Measurements to Characterize California Wildfires

Despite challenges with data quality and scope, low-cost sensor networks have skyrocketed in popularity over the last 15 years, making air quality data available on refined spatial scales. More recently, studies have leveraged both high and low-quality instruments to create stronger “hybrid” models, with most studies focusing on particulate matter. Low-cost measurements typically represent ground-level emissions only, providing context for human health issues from climate change-driven events such as wildfires. Since low-cost sensors’ capabilities are localized, daily events and microclimates tend to dominate the data rather than larger regional or atmospheric trends. Likewise, their low cost explains their high uncertainty. In contrast, some regulatory-grade instruments produce column measurements as well, providing reliable information on a broader scope. To bridge this gap while expanding into gas-phase measurements, we deployed 12 air quality sensor packages in California, USA during the 2022 wildfire season. These INSTEP (Inexpensive Network Sensor Technology Exploring Pollution) monitors measure carbon monoxide (CO), carbon dioxide (CO2), ozone (O3), nitrogen dioxide (NO2), and several hydrocarbons including methane (CH4) and formaldehyde (HCHO). Half of the monitors were co-located with remote sensing spectrometers: NASA Pandora and Total Column Carbon Observing Network (TCCON). The overlap in pollutants includes NO2, O3, and HCHO between the INSTEP monitors and the Pandora column measurements. TCCON covers column CO, CO2, and CH4, rounding out our comparison. Most of the monitors were distributed throughout the San Francisco Bay area, and an additional three were located within 100 km of Los Angeles. The sites ranged in geographic and population characteristics, including desert, mountainous, coastal, and urban locations. Since varying environmental conditions such as temperature and pressure are known to challenge sensor performance, we will apply newer sensor “calibration” techniques meant to combat this. We will normalize our sensor signals by z-scoring them prior to applying a single calibration model in the form of multivariate linear regression or an artificial neural network. While this technique has been validated for the hydrocarbon and ozone sensor types (metal oxide), it has not yet been tested on electrochemical and non-dispersive infrared sensors, which are also used in the INSTEP monitors. This will serve as a test to see if this normalization technique – or another – is most effective in accounting for environmental differences among sensors. Related data analysis efforts have found success with a variety of geospatial analysis techniques, including weighted network models in which high-quality instruments are given higher weights than their low-cost counterparts. Our preliminary analysis will focus on kriging, which uses a Gaussian algorithm to assign weights, providing estimated pollution levels at locations between monitors. Smoke trajectory and evolution will also be considered using both measurement types. We also aim to baseline subtract our emission estimates from each region to determine which portion of emissions are regional and local, further characterizing burn differences in northern and southern California fires. Future directions include using INSTEP jointly with TEMPO satellite data, and mobile deployments on aircraft and uncrewed aerial vehicles (UAV).

Low-cost sensors↗
DOE OSTI2025

What regulates decomposition in agroecosystems? Insights from reading the tea leaves

Litter decomposition is a critical Earth process, recycling nutrients and setting a portion of plant tissue on a path toward soil organic matter. Despite this importance, we still lack a good understanding of local factors that regulate decomposition, especially in agroecosystems where management plays an outsized role. Using a narrow range of climate and soils, we buried 1,308 pre-manufactured “litter bags” of differing residue quality (i.e., green and rooibos tea leaves) in 109 plots across several management practices to (1) explore the local controls on decomposition in agroecosystems and (2) test the robustness of the Tea Bag Index (TBI). We found that management practices intended to increase soil ecosystem services, that is, soil health, altered the decomposition of both teas. For example, adding nitrogen fertilizer and implementing perennial cropping decreased the extent of green tea decomposition (carbon-to-nitrogen ratio, or C:N = 12.8). No-tillage increased, but perennial cropping decreased, the rate of rooibos tea decomposition (C:N = 50.1). Cropped prairie accelerated green tea decomposition and increased the extent of red tea decomposition. A random forest regression model showed that soil temperature was the strongest predictor of green tea decomposition, but a soil health score also played a significant role in predicting the mass remaining. Soil texture and nutrient availability best predicted rooibos tea decomposition. Finer textured soils seemed to decelerate rooibos decomposition but increased the extent of decomposition. Furthermore, we demonstrated that the TBI metrics correlated somewhat well with empirically derived decomposition constants and were similarly sensitive to the effects of management. Still, the green tea stabilization factor had a substantial prediction bias. Our study increased our basic understanding of what regulates decomposition in agroecosystems. It also showed that the TBI can be a scientifically rigorous citizen science approach to monitoring changes in soil health.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Data-Driven Surrogate Modeling with Microstructure-Sensitivity of Viscoplastic Creep in Grade 91 Steel

Abstract To support the development of advanced steel alloys tailored to withstand extreme conditions, it is imperative to account for the mechanical performance of components, while considering the influence of local microstructure on the macroscopic response. To this end, this study focuses on the development of microstructure-sensitive constitutive models for the mechanical response of Grade 91 steel exposed to extreme thermo-mechanical environments. Polynomial chaos expansion (PCE) surrogates are used to emulate high-fidelity polycrystal simulations of the viscoplastic response of Grade 91 steel as a function of the microstructure fingerprint (e.g., dislocations and precipitates). To cover a wide temperature–stress domain, two separate PCE surrogates—one that captures softening and the other that captures hardening behavior—are combined using another (sparse) Gaussian process regression model. The resulting constitutive creep surrogate model is integrated within the MOOSE finite element framework to simulate the intricate effects of microstructure, in particular MX-phase precipitates, on a component with a graded microstructure. Surrogate sensitivity analysis is applied to quantify the relevant impact of spatially varying microstructure on the creep response in a test-case involving a Grade 91 alloy with a prototypical weld.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Improved Bayesian regularization of inverse problems in vibrations and acoustics using noise-only measurements

Here, this paper studies Tikhonov regularization (ridge regression) parameter selection for problems in vibrations and acoustics. The selection method is based on a popular Bayesian method, but it incorporates measurements of sensor noise. The regularization parameter is closely related to the ratio of system input energy to noise energy, so noise measurements inform the inference procedure and improve parameter identification. In cases where standard Bayesian regularization identifies zero as the optimal regularization parameter, noise measurements guarantee a unique nonzero optimum. Sufficient theoretical criteria are developed for this guarantee. The method is verified in even-determined and under-determined configurations in an acoustic source localization simulation and a vibration load identification experiment. It is shown to yield significant improvements over existing empirical Bayesian regularization. Improvements are larger in the even-determined case and smaller in the under-determined case, wherein the inverse solution is less sensitive to the regularization parameter.

42 ENGINEERING↗
DOE OSTI2022

Understanding Predictability of Daily Southeast U.S. Precipitation Using Explainable Machine Learning

Abstract We investigate the predictability of the sign of daily southeastern U.S. (SEUS) precipitation anomalies associated with simultaneous predictors of large-scale climate variability using machine learning models. Models using index-based climate predictors and gridded fields of large-scale circulation as predictors are utilized. Logistic regression (LR) and fully connected neural networks using indices of climate phenomena as predictors produce neither accurate nor reliable predictions, indicating that the indices themselves are not good predictors. Using gridded fields as predictors, an LR and convolutional neural network (CNN) are more accurate than the index-based models. However, only the CNN can produce reliable predictions that can be used to identify forecasts of opportunity. Using explainable machine learning we identify which variables and grid points of the input fields are most relevant for confident and correct predictions in the CNN. Our results show that the local circulation is most important as represented by maximum relevance of 850-hPa geopotential heights and zonal winds to making skillful, high-probability predictions. Corresponding composite anomalies identify connections with El Niño–Southern Oscillation during winter and the Atlantic multidecadal oscillation and North Atlantic subtropical high during summer.

Pegion, Kathy↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2023

In-Situ Detection and Prediction of WAAM Cross Feature Geometry

Abstract Wire arc additive manufacturing (WAAM) is increasingly used by manufacturers due to its relatively low cost and high deposition rate compared to other metal AM methods, but the parts produced by WAAM can be subject to localized variations in part quality. One such variation is the cross-feature defect, whereby a localized part height increase occurs due to the crossing of deposition toolpaths. Mitigation of this defect is typically achieved using manual path planning strategies, but closed-loop control is underutilized. Since the nature of this defect and of the WAAM process is such that the previous layer’s geometry influences that of the subsequent layer’s, the cross-feature defect geometry changes throughout the deposition. Therefore, any closed-loop control strategy will need to incorporate the dynamic trait of this defect. The present work seeks to implement an in-situ process modeling approach where a regression model can be continuously updated to predict the defect geometry of the subsequent deposition layer based on the historical process data. Several multi-layer cross-feature geometries are deposited and current, voltage, and optical camera data is taken for each layer. The resulting cross-feature geometries are characterized using 3D scanning and the performance and accuracy of the in-situ modeling approach is evaluated.

Thien, Austen↗
DOE OSTI2025

Evaluating the Trustworthiness of Explainable Artificial Intelligence (XAI) Methods Applied to Regression Predictions of Arctic Sea Ice Motion

Abstract Recent advances in explainable artificial intelligence (XAI) methods show promise for understanding predictions made by machine learning (ML) models. XAI explains how the input features are relevant or important for the model predictions. We train linear regression (LR) and convolutional neural network (CNN) models to make 1-day predictions of sea ice velocity in the Arctic from inputs of present-day wind velocity and previous-day ice velocity and concentration. We apply XAI methods to the CNN and compare explanations to variance explained by LR. We confirm the feasibility of using a novel XAI method [i.e., global layerwise relevance propagation (LRP)] to understand ML model predictions of sea ice motion by comparing it to established techniques. We investigate a suite of linear, perturbation-based, and propagation-based XAI methods in both local and global forms. Outputs from different explainability methods are generally consistent in showing that wind speed is the input feature with the highest contribution to ML predictions of ice motion, and we discuss inconsistencies in the spatial variability of the explanations. Additionally, we show that the CNN relies on both linear and nonlinear relationships between the inputs and uses nonlocal information to make predictions. LRP shows that wind speed over land is highly relevant for predicting ice motion offshore. This provides a framework to show how knowledge of environmental variables (i.e., wind) on land could be useful for predicting other properties (i.e., sea ice velocity) elsewhere. Significance Statement Explainable artificial intelligence (XAI) is useful for understanding predictions made by machine learning models. Our research establishes trustability in a novel implementation of an explainable AI method known as layerwise relevance propagation for Earth science applications. To do this, we provide a comparative evaluation of a suite of explainable AI methods applied to machine learning models that make 1-day predictions of Arctic sea ice velocity. We use explainable AI outputs to understand how the input features are used by the machine learning to predict ice motion. Additionally, we show that a convolutional neural network uses nonlinear and nonlocal information in making its predictions. We take advantage of the nonlocality to investigate the extent to which knowledge of wind on land is useful for predicting sea ice velocity elsewhere.

Hoffman, Lauren [Scripps Institution of Oceanograp↗
NASA NTRS1990

Study of the crust and mantle using magnetic surveys by Magsat and other satellites

This article summarizes some of the methodology and results for studies of the earth's mantle and crust. Mantle conductivity studies can be made either by studying signals impressed on the earth from outside, e.g., the ionosphere or magnetosphere, or by studing signals originating in the core and transmitted through the mantle. Crustal field studies begin with a careful selection of the data and subsequent removal of core and external fields by some sort of filtering. Average maps from different local times sometimes differ, presumably due to the remaining presence of fields of external origin. Several techniques for further filtering are discussed. Where large-area aeromagnetic maps are available, crustal maps derived from satellite data can be compared with upward continued data. In general, the comparisons show agreement, with some differences, particularly in and near the auroral belts. The satellite data are further reduced by various methods of inverse and forward modeling, sometimes including reduction to the pole. These techniques are generally unstable at the equator. Common methods of stabilizing the inversions include principle components analysis and ridge regression. Because of the presence of the core field, the entire crustal contribution from the field is not known. Also, there is a basic nonuniqueness to the inverse solutions. Nevertheless, magnetizations that are interpretable can be derived.

Langel, Robert A.↗
DOE OSTI2025

Investigating the effects of local environment on nitrogen vacancies in high-entropy metal nitrides

High-entropy metal nitrides are an important material class in a variety of applications, and the role of nitrogen vacancies is of great importance for understanding their stability and mechanical properties. Here, we study six different high-entropy nitrides with eight different metal species to build a predictive model of the nitrogen-vacancy formation energy. We construct sets of supercells that maximize the number of unique nitrogen environments for a given chemistry, and then use density-functional theory to calculate the energy density for all nitrogen sites, and the vacancy formation energies for the highest, lowest, and a median subset based on the energy densities. The energy density of nitrogen sites correlates with the vacancy formation energies, for binary, ternary, and high-entropy nitrides. A linear regression model predicts the vacancy formation energies using only the nearest-neighbor composition; across our eight metals, we find the largest vacancy formation energies next to Hf, then Zr, Ti, V, Cr, Ta, Nb, and the lowest near Mo. Additionally, we see that binary nitride data show qualitatively similar vacancy formation energy trends for high-entropy nitrides; however, the binary data alone are insufficient to predict the complex nitride behavior. Our model is both predictive and easily interpretable, and correlates with experimental data.

DeSilva, Charith R. [Univ. of Illinois at Urbana-C↗
DOE OSTI2024

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↗