Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Base stations”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Los Angeles Air Force Base Vehicle-to-Grid Demonstration (Final Project Report)

Electrification of non-tactical vehicle fleets represents a key efficiency and energy security objective for the United States Department of Defense. To achieve electrification, the department targeted vehicle-to-grid services as a way to decrease the overall cost of operating the vehicle fleet and achieve rough parity with traditional internal combustion engine vehicle fleets. This report describes efforts to aggregate a fleet of bi-directional electric vehicles and charging stations to provide regulation up and regulation down in the California Independent System Operator ancillary services market. A 29-vehicle electric vehicle demonstration fleet, consisting of mixed purpose and duty vehicles such as sedans, pickups, vans, and medium-duty trucks, was deployed at the Los Angeles Air Force base. The fleet provided frequency regulation to the California Independent System Operator’s wholesale electricity market to determine the capability of recouping some of the additional costs of procuring electric vehicles and their supporting infrastructure. Lawrence Berkeley National Laboratory, with its partner Kisensum, LLC, developed the fleet scheduling, optimization, and control software to allow the vehicle fleet at the air force base to participate in the ancillary services markets. This report focuses on the control software and market interactions, the significant challenges faced and solutions devised to address them, and examines the potential of using the electric vehicle fleet as an energy storage resource for the base buildings, an application known as vehicle-tobuilding, in providing demand response and emergency backup power. The report discusses key findings related to providing frequency regulation to the California Independent System Operator market, electric vehicle fleet performance, compatibility of varying resource parameters of vehicle fleet aggregation, the need for automated methods for communicating hour-ahead energy bidding, challenges related to battery capacity and charge/discharge rates, and monthly settlement revenue.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Dynamic PRA Methods to Evaluate the Impact on Accident Progression of Accident Tolerant Fuels

Accident tolerant fuels (ATFs) are new nuclear fuels developed in response to the accident at the Fukushima power station in March 2011. The goal of ATFs is to withstand accident scenarios through better performance compared to currently employed fuels (e.g., small-scale hydrogen generation). This paper targets a method for evaluating and comparing ATF performance from a probabilistic risk assessment (PRA) perspective by employing a newly developed combination of event trees and dynamic PRA methods. Compared to classical PRA methods based on event trees and fault trees, dynamic PRA can evaluate with higher resolution the safety impacts of physics dynamics and the timing/sequencing of events on the accident progression without the need to introduce overly conservative modeling assumptions and success criteria. In this paper, we analyze the impact on the accident progression of three different cladding configurations for two initiating events [a large break loss-of-coolant accident (LB-LOCA) and a station blackout (SBO)] by employing dynamic PRA methods. The goal is to compare the safety performance of ATFs (FeCrAl and Cr-coated cladding) and the currently employed Zr-based clad fuel. We employ two different strategies. The first focuses on the identification of success criteria discrepancies between the accident sequences generated by the classical PRA model and the set of simulation runs generated by dynamic PRA using ATF. The second one, on the other hand, directly uses dynamic PRA to evaluate the impact of timing of events (e.g., recovery actions) on accident progression. In conclusion, by applying these methods to the LB-LOCA and SBO initiating events, we show how dynamic PRA methods can provide analysts with detailed and quantitative information on the safety impact of ATFs.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Synergies between repositioning and charging strategies for shared autonomous electric vehicle fleets

The emergence of on-demand shared autonomous electric vehicle (SAEV) service requires careful charging station planning and a joint charging and repositioning strategy to mitigate empty travel. This study couples charging and repositioning events as a means of improving service quality (rider wait times), reducing empty travel due to repositioning or charging, and improving fleet utilization (average daily trips per vehicle and charging queues). This synergy is explored for the Austin, Texas region using POLARIS, an agent-based model. On average, wait times were 39% lower, and average daily trips served per SAEV increased up to 6.4 (or 28%) compared to SAEV repositioning with heuristic charging. Coupling repositioning with charging decreased the fleet's percent empty travel on average by 1.6%, relative to the scenario treating them as independent events (which varies by charging station design). Sparser charging stations reduce investment costs, and operators can leverage this framework to keep traveler wait times low.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Joint optimization of land carbon uptake and albedo can help achieve moderate instantaneous and long-term cooling effects

Abstract Both carbon dioxide uptake and albedo of the land surface affect global climate. However, climate change mitigation by increasing carbon uptake can cause a warming trade-off by decreasing albedo, with most research focusing on afforestation and its interaction with snow. Here, we present carbon uptake and albedo observations from 176 globally distributed flux stations. We demonstrate a gradual decline in maximum achievable annual albedo as carbon uptake increases, even within subgroups of non-forest and snow-free ecosystems. Based on a paired-site permutation approach, we quantify the likely impact of land use on carbon uptake and albedo. Shifting to the maximum attainable carbon uptake at each site would likely cause moderate net global warming for the first approximately 20 years, followed by a strong cooling effect. A balanced policy co-optimizing carbon uptake and albedo is possible that avoids warming on any timescale, but results in a weaker long-term cooling effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗
DOE OSTI2023

SMALL UNMANNED MARINE HYDROKINETIC PLATFORMS FOR POWER GENERATION IN COASTAL AND TIDAL WATERS

The feasibility and optimization of small unmanned mobile marine hydrokinetic (MHK) energy platforms for harvesting marine current energy in coastal and tidal waters are examined. A case study of a platform based on the use of a free-surface waterwheel (FSWW) mounted on an autonomous unmanned surface vehicle (USV) was conducted. Such platforms can serve as recharging stations for aerial drones (UAVs), enabling extension of the UAVs’ autonomous operating time. An unmanned MHK platform potentially meets this need with sustainable power harvested from water currents. For the case study, six different waterwheel configurations were field-tested in the Intracoastal Waterway of South Florida in support of determining the configuration that produced the most power. Required technologies for unmanned operations of the MHK platform were developed and tested. The data from the field-testing were analyzed to develop an empirical relation between the wheel’s theoretical hydrokinetic power produced and the mechanical power harnessed by the MHK platform with various waterwheel configurations during field-testing. The field data was also used to determine the electrical power generated by the FSWW configurations during field-testing. The study has led to the development of standardized testing procedures. The empirical relation is used to examine predicted power production through scaling up different physical aspects of the waterwheel.

McKinney, Adriana↗
DOE OSTI2025

Source Characterization at Low Yields from Physics Experiment One A

The detonation of a well-recorded coupled explosion at the well-characterized Nevada National Security Site, as part of Physics Experiment One, has provided us with an opportunity to compare and contrast methods and results of yield estimation and other event characterization methods. Yield estimates are made using several methods, including local/regional magnitudes, waveform envelopes, acoustic amplitudes, joint seismoacoustic yield, and moment tensors. We include in our analysis a discussion and comparison of applicability, calibrations, assumptions, and uncertainties involved in each. With the exception of magnitude-based methods, all the techniques were more or less successful in recovering the known yield. It is, however, useful to be cognizant of requirements and calibrations, such as station coverage, emplacement conditions, propagation characterization, and uncertainties, to apply the methods to other areas with less complete information. We demonstrate that we can be successful at monitoring low-yield explosions using a variety of methods assuming that the signal-to-noise ratio is adequate at recording stations, and that the material properties are taken into account. Depth estimates are also important, particularly if one allows the possibility of a surface explosion.

Pasyanos, Michael E. [Lawrence Livermore National ↗
DOE OSTI2021

Complex Fault Geometry of the 2020 M ww 6.5 Monte Cristo Range, Nevada, Earthquake Sequence

On 15 May 2020 an M ww 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enabled detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8), which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km-long east–northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure that extends north from the eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest toward the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. We infer that shallow aftershocks throughout the system highlight fault-fracture meshes that connect mapped fault systems at depth. Here, comparing earthquake data with surface ruptures and a simple geodetic fault model sheds light on the complexity of this recent M 6.5 Walker Lane earthquake.

58 GEOSCIENCES↗
DOE OSTI2024

Oak Ridge National Laboratory Neutron Spin Echo Beamline on NB-2

Neutron Spin Echo (NSE) spectroscopy uniquely measures the Q-dependence of slow relaxation dynamics, and having such a machine at HFIR will advance polarized neutron spectroscopy and promote the study of biophysical and energy materials using neutrons. The 2018 Instrument Advisory Board (IAB) advises Neutron Sciences Directorate, ORNL to upgrade the cold neutron delivery system at HFIR, addressing geometrical challenges in neutron transfer from the bright cold source. This is being optimized now: the entrance to the guide system is moving closer to the source, and more guides are being added. However, there are still significant challenges to making the cold flux at the sample world-class. With the proposed guide system, NB-2 can deliver more than 10 6 polarized neutrons per cm 2 /sec below 8Å, and probably 10Å with further optimization. With no potential for running user experiments with neutrons longer than 15Å, the utility of such a machine for biological studies is limited. However, there is a large demand for this class of spectrometer in materials science, chemical engineering, and nanotechnology. Workshops held over the past decade and the three-source vision have highlighted the community's need for a low-angle, high-energy resolution spectrometer at HFIR. A detailed design study of this spectrometer is needed, focusing on optimizing the neutron delivery system for SANS-type studies. This neutron spectrometer would complement studies already performed on both SNSNSE and BASIS by overlapping and extend the dynamic range accessible, while acknowledging the proposed machine, EXPANSE, at the Second Target Station. Two technologies currently exist for such a spectrometer: traditional DC solenoids like those on BL-15 at SNS and IN15 at the ILL, and RF-flipper-based machines similar to RESEDA at FRMII. Either instrument would satisfy the scientific justification. However, there are strong business and scientific cases to utilize the thermal flux, the resonant expertise developed at ORNL, and the potential to utilize an entangled beam of neutrons to probe quantum matter by adopting the neutron resonant spin-echo configuration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

An absolute recoil-carbon polarimeter for polarized light-ions at the BNL Booster

We describe a recoil-carbon polarimeter for the BNL Booster synchrotron capable of measuring the transverse polarization of both the polarized proton beam and the polarized 3 He ion beam throughout the acceleration cycle. The instrument addresses a critical gap in the BNL polarized beam program: no independent polarization measurement currently exists in the injector chain between the 200 MeV Linac polarimeter and the AGS. The proposed new polarimeter will provide continuous, absolute polarization measurements at the Booster stage, supplying additional anchor points for the polarization transmission through the accelerator chain. The polarimeter 9 is based on elastic pC and 3 He C scattering from a thin internal carbon target, with detection of the recoil 12 C nucleus in a fixed six-station silicon detector ring. A single detector geometry provides continuous kinematic coverage from injection to extraction for both beam species without mechanical adjustment. For pC scattering, the available data and the Bonin parametrization provide a well-established basis for estimating the analyzing power and figure of merit over the Booster energy range; for 3 He C the sole experimental anchor is a measurement at 443 MeV, and the Booster polarimeter itself is identified as the instrument to map the analyzing power across the remainder of the acceleration ramp by a ramp-and-return calibration that anchors the absolute 3 He polarization scale to a few percent. Provision for future deuteron polarimetry is incorporated in the chamber design without modification to the existing geometry. With appropriate detector upgrades, the same recoil-carbon technique can be extended to deuteron beams, for which analyzing power data are already available, and, with future analyzing-power measurements, also to 6 Li and 7 Li beams, making the polarimeter a versatile instrument for the full range of light polarized ion species anticipated at BNL.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

A field programmable gate array-based timing and control system for the dynamic compression sector

A field programmable gate array (FPGA)-based timing and trigger control system has been developed for the Dynamic Compression Sector (DCS) user facility located at the Advanced Photon Source (APS) at Argonne National Laboratory. The DCS is a first-of-its-kind capability dedicated to dynamic compression science. All components of the DCS laser shock station—x-ray choppers, single-shot shutter, internal laser triggers, and shot diagnostics—must be synchronized with respect to the arrival of x rays in the hutch. An FPGA synchronized to the APS storage ring radio frequency clock (352 MHz) generates trigger signals for each stage of the laser and x-ray shutter system with low jitter. The developed FPGA-based control system was the first system used to control the laser and the shutter system since its commissioning, and it has been developing since then to improve the timing jitter. The system is composed of a Zynq FPGA, a debug card, line drivers, and a power supply. The delay and offsets of the trigger signals can be adjusted by using a user-friendly graphical user interface with high precision. Here, the details of the system architecture, timing requirements, firmware, and software implementation along with the performance evaluation are presented in this paper. The system offers low timing jitter (15.5 ps rms) with respect to the APS 352 MHz clock, suitable for the 100 ps (FWHM) x-ray bunch duration at the APS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Differences in urban plant community compositions across an urban-rural gradient in Knoxville, TN

Urban forests, or vegetation in areas under heavy human influence, provide many ecosystem services to urban residents such as localized cooling via evapotranspiration, shade, filtering of air pollution, and the associated health benefits of natural spaces. In order to quantify the magnitude of localized cooling by trees growing in varying levels of urbanization (based on % impervious surfaces, e.g., buildings, pavement), urban forest species composition, tree size, and tree density must be characterized. As a part of Oak Ridge National Laboratory’s (ORNL) urban forest temperature study, we conducted tree censuses in five Knoxville city parks where ORNL meteorological stations are deployed. Moreover, we measured every woody plant ≥ 5 cm diameter at breast height (DBH) within a 50 m radius of each site’s meteorological station for its DBH and species identification. When possible, individuals were identified down to species. Certain genera (Quercus spp., Carya spp., Pinus spp.) were identified down to genera in interest of time. Individual and total site basal area were calculated from measured DBH data. Results show notable differences in urban plant community compositions and total woody plant basal area across sites, with more urban sites closer to downtown (West View and SEEED) having lower tree basal area than the more suburban sites (West Hills, Cumberland Estates, and Victor Ashe). We identified 54 species across all sites, with West Hills and Victor Ashe having the highest species diversity. Our results show differences in forest compositions and sizes across Knoxville, which are currently informing ORNL’s evapotranspiration estimates for each site. Data Summary: Census data for West Hills (WH), Cumberland Estates (CE), Victor Ashe (VA), West View (WV), and Socially Equal Energy Efficient Development or SEEED (SD) urban forests in Knoxville, TN, USA, including tree size based on diameter at breast height (DBH; 1.3 m), species identification (Latin and common names), and basal area per stem (BA=π×[.5*DBH]^2). Field data are summarized in this file: “Community_Composition_Data.CSV”. Site-specific data detailing each site’s coordinates, number of stems measured at DBH, average tree DBH, α-diversity (number of species present), and total site basal area (sum of individual basal areas per site) are in this file: “Site_Comparisons.CSV”.

Warren, Jeffrey [ORNL] (ORCID:0000000206804697)↗
DOE OSTI2022

Network Constraints Consideration for Grid-Edge Energy Management System

Increased deployment of distributed energy resources (DER) in distribution system is bringing need for enhanced grid intelligence, control, and flexibility. This is significant at the edge of the grid where DERs, loads or microgrids are located. Integrated DERs in a distribution system need to follow grid codes to avoid violations that results in DER disconnection. To comply with grid code requirements (e.g. IEEE 1547-2018) at the grid edge level, network constrained grid edge energy management system (GEEMS) is proposed in this paper. The objective of GEEMS is to provide economic solution for active and reactive power dispatch set-points at each interval and ensure voltage regulation to support secure interconnection of the grid edge segment to the distribution system. To evaluate the proposed GEEMS framework, four DERs are included into IEEE 13 bus system. GEEMS outperforms the existing economic dispatch-based energy management system by reducing the voltage violation.

Electric vehicle depot charging station↗
DOE OSTI2020

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are the most important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a direct comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. Here, we used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a millimetre VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects, compared to the addition of thermal noise only. Using synthetic data based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M 87, we performed case studies to assess the image quality that can be obtained with the current and future EHT array for different weather conditions. Our synthetic observations show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of our GRMHD source models can be recovered robustly with the EHT2017 array after performing calibration steps, which include fringe fitting, a priori amplitude and network calibration, and self-calibration. With the planned addition of new stations to the EHT array in the coming years, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Integrative investigation of dust emissions by dust storms and dust devils in North Africa

Dust aerosols in North Africa account for >50% of the global total; however dust emission areas are still unclear. Based on the analysis of dust storms simulated with the numerical Weather Research and Forecast (WRF) model, satellite aerosol index (AI), and the dust data observed at 300 meteorological stations over 20 years, the spatio-temporal distribution characteristics of dust storm, dust devil and AI are compared and analyzed. From this study, the results show that: 1) There are two dust emission mechanisms: the dynamically-dominated dust storm and thermally-dominated dust devil; 2) Dust storms occur most frequently in Spring and are concentrated in the areas of Grand Erg Occidental Desert to the Erg Chech-Adrar Desert, the northern part of Grand Erg Oriental, the Atouila Desert to the Ouarane Desert, the Mediterranean coast, the eastern side of Nubian Desert and Bodélé Depression; 3) Dust devils occur most frequently from April to August and are mainly concentrated in the central part of North Africa, especially in the southwest of Hoggar Mountains to the west of Air Mountains, the border area of Egypt - Sudan - Libya and the vicinity of Tibesti Plateau; 4) The spatio-temporal distribution of AI is correlated more with the dust devils emission whereas the annual average contributions by dust storms and dust devils are 61.3% and 38.7%, respectively. This study discovers a new area of dust emissions by dust devils, and provides a better explanation for the spatio-temporal distribution of AI in North Africa.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A DNN surrogate unsteady aerodynamic model for wind turbine loads calculations

A recurrent deep-neural network (DNN) surrogate model capable of modeling the unsteady aerodynamic response and dynamic stall behavior of wind turbine blades has been developed and validated for use in engineering design codes. The model is trained using a subset of the oscillating airfoil experiments conducted at the Ohio State University wind tunnel. The predictions from our DNN model show excellent agreement with the measured data and, in all cases, a marked improvement over the state-of-the-art unsteady aerodynamic models. The DNN-based unsteady aerodynamics model was integrated with OpenFAST to perform full-turbine load computations for the NREL-5MW rotor. The largest differences are observed for the inboard stations, particularly in the pitching moment response, when using the new surrogate model compared to the other models available in OpenFAST.

17 WIND ENERGY↗
DOE OSTI2023

Emulation of seismic-phase traveltimes with machine learning

SUMMARY We present a machine learning (ML) method for emulating seismic-phase traveltimes that are computed using a global-scale 3-D earth model and physics-based ray tracing. Accurate traveltime predictions based on 3-D earth models are known to reduce the bias of event location estimates, increase our ability to assign phase labels to seismic detections and associate detections to events. However, practical use of 3-D models is challenged by slow computational speed and the unwieldiness of pre-computed lookup tables that are often large and have prescribed computational grids. In this work, we train a ML emulator using pre-computed traveltimes, resulting in a compact and computationally fast way to approximate traveltimes that are based on a 3-D earth model. Our model is trained using approximately 850 million P-wave traveltimes that are based on the global LLNL-G3D-JPS model, which was developed for more accurate event location. The training-set consists of traveltimes between 10 393 global seismic stations and randomly sampled event locations that provide a prescribed, distance-dependent geographic sample density for each station. Prediction accuracy is dependent on event-station distance and whether the station was included in the training set. For stations included in the training set the mean absolute deviation (MAD) of the difference between traveltimes computed using ray tracing through the 3-D model and the ML emulator for local, regional, and teleseismic distances are 0.090, 0.125 and 0.121 s, respectively. For tested station locations not included in the training set, MAD values for the three distance ranges increase to 0.173, 0.219 and 0.210 s, respectively. Empirical traveltime residuals for a global reference data are indistinguishable when ML emulation or the 3-D model is used to compute traveltimes. This result holds regardless of whether the recording station is used in ML training or not.

58 GEOSCIENCES↗