Search NASA⌕ Search

SEARCH · Search NASA

Results for “atmospheric acoustics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

Advances in geophysical forensic event monitoring

Forensic analysis of man-made, non-nuclear events (such as industrial accidents, explosion experiments and mine collapses) has become more frequent and detailed owing to advancements in geophysical monitoring. Here, in this Technical Review, we demonstrate how geophysical forensic monitoring using seismic, infrasound and hydroacoustic recordings provides insights on events in the solid earth, atmosphere and underwater. Advanced techniques, including machine-learning-based models, have been developed to detect, identify and investigate these events, providing information on location, subevents, sources and explosive yield. The increase in data availability, application of advanced methods and computation and the growth of multitechnology approaches have increased the accuracy of forensic event analysis and enabled more realistic characterization of uncertainties. For example, the 2020 Beirut explosion in Lebanon demonstrated that various seismic, acoustic and other methods could be used to estimate explosive yield (and yield uncertainties) of about 1 ktonne, providing confidence in the application of these methods to smaller events where data are available. However, forensic investigations remain largely limited to known events with identified sources. Increased access to data, sophisticated analysis methods and high-resolution earth models will improve forensic event analysis further, enabling civil and scientific applications, such as localization in the search for the lost ARA San Juan submarine.

geophysics↗

Topographically Scattered Infrasound Waves Observed on Microbarometer Arrays in the Lower Stratosphere

When an acoustic wave strikes a topographic feature, some of its energy is scattered. Sensors on the ground cannot capture these scattered signals when they propagate at high angles. We report observations of upwardly-scattered acoustic waves prior to refraction back to the ground, intercepting them with a set of balloon-borne infrasound microbarometers in the lower stratosphere over northern Sweden. We show that these scattered waves generate a coda whose presence can be related to topography beneath balloons and low-altitude acoustic ducts. The inclination of the coda signals changes systematically with time, as expected from waves arriving from scatterers successively closer to receivers. The codas are present when a temperature inversion channels infrasound from a set of ground chemical explosions along the ground, but are absent following the inversion's dissipation. Since scattering partitions energy away from the main arrival, these observations imply a mechanism of amplitude loss that had previously been inaccessible to measurement. As such, these data and results allow for a better comprehension of interactions between atmospheric infrasound propagation and the solid earth

58 GEOSCIENCES↗

Prediction of regional infrasound produced by supersonic sources using a ray-based Mach cone source

The geometry of the Mach cone produced by a supersonic source is analyzed and mapped into initial conditions used in acoustic ray tracing. The resulting source model is combined with spherical geometry ray tracing methods to enable propagation simulations for infrasonic signals produced by bolides, space debris, rockets, aircraft, and other fast-than-sound sources out to typical infrasonic observation distances of hundreds or thousands of kilometers. Idealized linear and parabolic trajectories typical of bolides and rockets, respectively, are used to demonstrate the calculation of regional infrasonic signals produced by such sources and characteristics of the radiated infrasonic waves are found to vary strongly with the geometry of the trajectory and atmospheric structure. Predicted regional infrasonic signals are compared with those observed from a November 2020 bolide that passed over Scandinavia using a combination of institutionally maintained infrasound stations and “citizen scientist” data from the Raspberry Shake data repository.

47 OTHER INSTRUMENTATION↗

Prediction of regional infrasound produced by supersonic sources using a ray-based Mach cone source

The geometry of the Mach cone produced by a supersonic source is analyzed and mapped into initial conditions used in acoustic ray tracing. The resulting source model is combined with spherical geometry ray tracing methods to enable propagation simulations for infrasonic signals produced by bolides, space debris, rockets, aircraft, and other fast-than-sound sources out to typical infrasonic observation distances of hundreds or thousands of kilometers. Idealized linear and parabolic trajectories typical of bolides and rockets, respectively, are used to demonstrate the calculation of regional infrasonic signals produced by such sources and characteristics of the radiated infrasonic waves are found to vary strongly with the geometry of the trajectory and atmospheric structure. Predicted regional infrasonic signals are compared with those observed from a November 2020 bolide that passed over Scandinavia using a combination of institutionally maintained infrasound stations and “citizen scientist” data from the Raspberry Shake data repository.

42 ENGINEERING↗

Cosmological constraints on dark matter interactions with ordinary matter

Dark matter interactions with electrons or protons during the early Universe leave imprints on the cosmic microwave background and the matter power spectrum, and can be probed through cosmological and astrophysical observations. These interactions lead to momentum and heat exchange between the ordinary and dark matter components, which in turn results in a transfer of pressure from the ordinary to the dark matter. Here, we explore these interactions using a diverse suite of data: cosmic microwave background anisotropies, baryon acoustic oscillations, the Lyman-α forest, and the abundance of Milky-Way subhalos. We derive constraints using model-independent parameterizations of the dark matter–electron and dark matter–proton interaction cross sections and map these constraints onto concrete dark matter models. Our constraints are complementary to other probes of dark matter interactions with ordinary matter, such as direct detection, big bang nucleosynthesis, various astrophysical systems, and accelerator-based experiments. They exclude sufficiently large cross sections for a large range of dark matter masses, which cannot be accessed by direct-detection experiments due to the overburden from the Earth’s atmosphere or crust.

79 ASTRONOMY AND ASTROPHYSICS↗

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗

Electro-optic transduction in silicon via gigahertz-frequency nanomechanics

Interfacing electronics with optical fiber networks is key to the long-distance transfer of classical and quantum information. Piezo-optomechanical transducers enable such interfaces by using gigahertz-frequency acoustic vibrations as mediators for converting microwave photons to optical photons via the combination of optomechanical and piezoelectric interactions. However, despite successful demonstrations, efficient quantum transduction remains out of reach due to the challenges associated with hybrid material integration and increased loss from piezoelectric materials when operating in the quantum regime. Here, we demonstrate an alternative approach in which we actuate 5-GHz phonons in a conventional silicon-on-insulator platform. In our experiment, microwave photons resonantly drive a phononic crystal oscillator via the electrostatic force realized in a charge-biased narrow-gap capacitor. The mechanical vibrations are subsequently transferred via a phonon waveguide to an optomechanical cavity, where they transform into optical photons in the sideband of a pump laser field. Operating at room temperature and atmospheric pressure, we measure a microwave-to-optical photon conversion efficiency of 1.72±0.14×10 −7 in a 3.3 MHz bandwidth. Our results mark a stepping stone towards quantum transduction with integrated devices made from crystalline silicon, which promise efficient high-bandwidth operation and integration with superconducting qubits. Additionally, the lack of need for piezoelectricity or other intrinsic nonlinearities makes our approach applicable to a wide range of materials for potential applications beyond quantum technologies.

42 ENGINEERING↗

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗

Multi‐Sensor Trajectory Reconstruction of the 24 April 2025 Alaska Fireball and Implications for Planetary Defense

On 24 April 2025 at 18:30:57 UTC, a bright daytime fireball over Southcentral Alaska was detected by 37 seismic stations, 16 single infrasound sensors, and four infrasound arrays, yielding 30 ballistic and multiple fragmentation arrivals. Here, the unprecedented density of seismoacoustic coverage enabled detailed reconstruction of the event using acoustic signals, with fragmentation source locations further guiding the identification of Doppler weather radar signatures of a meteorite fall. Incorporation of a radar-derived terminal point yielded a final trajectory solution, which agreed closely with an independent optical trajectory solution from video analysis. The reconstructed entry parameters from seismoacoustic analysis indicate a velocity of 25.3 km/s, an entry angle of 19°, and an energy release of ∼38 t TNT equivalent. Assuming a chondritic composition, the pre-entry object diameter was ∼0.7 m. Using orbital parameters from the optical solution, we estimate meteoroid composition as most likely an L-type ordinary chondrite. The event occurred in the sub-Arctic, where space-based optical systems face challenges in detection, demonstrating the critical role of dense ground-based seismoacoustic networks in characterizing high-latitude atmospheric entries. This uniquely well-recorded event demonstrates the capability of dense seismoacoustic networks to constrain bolide trajectories, energetics, and fragmentation, with radar and optical data providing critical confirmation and complementary perspectives. These results bridge the methodological gap between planetary-defense monitoring of natural impactors and space-traffic analyses of artificial reentries, illustrating how multi-sensor integration can deliver calibration-grade trajectories even for unpredicted events.

Fireball↗

Search for Infrasound Signals in InSight Data Using Coupled Pressure/Ground Deformation Methods

The unprecedented quality and sampling rate of seismometer and pressure sensors of the InSight Mars mission allow us to investigate infrasound through its pressure and ground deformation signals. This study focuses on compliance effects induced by acoustic waves propagating almost horizontally close to the surface. The compliance of acoustic waves is first estimated using the compliance estimates from pressure perturbations moving at wind speed. Then, a marker of compliance events is used to detect events of ground deformation induced by pressure variations, in three frequency bands from 0.4 to 3.2 Hz, from InSight sol 180 to 690. Additional selection criteria are imposed on the detected events to focus on acoustic waves and to remove various noise sources (e.g., wind effects or seismometer artifacts). After an automated selection, the visual inspection of the records allows us to validate two infrasound candidates that cannot be related to pressure perturbations moving at wind speed nor to known noise sources. For our highest quality infrasound candidate, the relation between this event and a convective vortex occuring 10 s later is tested. The azimuth of the vortex position at the time of infrasound detection is not consistent with the arrival azimuth of the suspected infrasound inferred from the polarization of seismometer records, thus the link between these two phenomena cannot be demonstrated. Further investigations would require a better understanding of wind-related noise impacting InSight sensors and of the effects of lateral variations of subsurface mechanical properties on the ground deformations induced by atmospheric pressure variations.

58 GEOSCIENCES↗

Regional infrasonic observations from surface explosions – influence of atmospheric variations and realistic terrain

SUMMARY A pair of 1-ton, conventional surface explosions were conducted at the Nevada National Security Site in the fall of 2020 producing seismoacoustic signatures observable hundreds of kilometres from the source location. Regional infrasonic observations include tropospheric ducting at large distances to the south, a wide stratospheric waveguide with signals observed more than 700 km to the east, and anomalous arrivals in the stratospheric shadow zone. Notable differences in propagation between the events are identified despite the explosions being conducted just two days apart due to a sharp temporal shift in the tropospheric winds as well as structural changes in the stratospheric winds. Propagation simulations of the two events have been completed using a combination of ray tracing and parabolic equation (PE) methods. Simulations have been conducted to quantify the impact of the temporal variations in the atmosphere as well as the influence of terrain on propagation. Temporal variations in reflection locations are found to produce notable changes in downrange propagation due to spatially varying terrain features. Finite frequency effects modelled by the PE are found to predict ensonification not included in corresponding 2D ray tracing simulations. Notable variations in predicted signal amplitude are found due to focusing by along-path and cross-path terrain gradients; though, the later of these is only modelled using fully 3D ray tracing analysis.

58 GEOSCIENCES↗

Equation of State of TiN at High Pressures and Temperatures: A Possible Host for Nitrogen in Planetary Mantles

Nitrogen, the most abundant element in Earth's atmosphere, is also a primary component of solid nitride minerals found in meteorites and on Earth's surface. If they remain stable to high pressures and temperatures, these nitrides may also be important reservoirs of nitrogen in planetary interiors. We used synchrotron X-ray diffraction to measure the thermal equation of state and phase stability of titanium nitride (TiN) in a laser-heated diamond anvil cell at pressures up to ~70 GPa and temperatures up to ~2500 K. TiN maintains the cubic B1 (NaCl-type) crystal structure over the entire pressure and temperature range explored. It has K 0 = 274(4) GPa, K 0 ' = 3.9(2), and γ 0 = 1.39(4) for a fixed V 0 = 76.516(30) Å 3 (based on experimental measurements), q = 1, and θ 0 = 579 K. Additionally, we collected Raman spectra of TiN up to ~60 GPa, where we find that the transverse acoustic (TA), longitudinal acoustic (LA), and transverse optical (TO) phonon modes exhibit mode Grüneisen parameters of 1.66(17), 0.54(15), and 0.93(4), respectively. Based on our equation of state, TiN has a density of ~5.6–6.4 g/cm 3 at Earth's lower mantle conditions, significantly more dense than both the mantle of the Earth and the estimated densities of the mantles of other terrestrial planets, but less dense than planetary cores. Furthermore, we find that TiN remains stable against physical decomposition at the pressures and temperatures found within Earth's mantle, making it a plausible reservoir for deep planetary nitrogen if chemical conditions allow its formation.

58 GEOSCIENCES↗

Ionospheric Disturbances in GNSS TEC Data: SpaceX Falcon 9 Deorbit Maneuvers Over CONUS in April–May 2024

Traveling ionospheric disturbances (TIDs) driven by a large number of internal and external sources are detectable with dense networks of ground‐based Global Navigation Satellite System (GNSS) receivers' measurements of total electron content (TEC). We present the newly developed System for Rapid Analysis of Ionospheric Dynamics (S‐RAID), providing data and visual GNSS TEC products of TIDs of periods ~2–120 min and horizontal resolution up to tens of kilometers for years 2017–2024. The S‐RAID data reveal myriad natural and anthropogenic TIDs from meteorology and space weather, and from human spaceflight activities. In this report, we focus on new prominent disturbances found during SpaceX second stage deorbit maneuvers in April–May 2024. Signatures include waves emanating from the Falcon's trajectory above California and depletions following its deorbit and passage over Arizona. These findings suggest further opportunities to detect and quantify small‐scale events in the ionosphere as well as to understand the responses of the atmosphere‐ionosphere system to known inputs.

58 GEOSCIENCES↗

The Utility of Infrasound in Global Monitoring of Extraterrestrial Impacts: A Case Study of the 2008 July 23 Tajikistan Bolide

Abstract Among the various observational techniques used for the detection of large bolides on a global scale is a low-frequency sound known as infrasound. Infrasound, which is also one of the four sensing modalities used by the International Monitoring System, offers continuous global monitoring and can be leveraged for planetary defense. Infrasonic records can provide an additional dimension for event characterization and a distinct perspective that might not be available through any other observational method. This paper describes the infrasonic detection and characterization of the bolide that disintegrated over Tajikistan on 2008 July 23. This event was detected by two infrasound stations at distances of 1530 and 2130 km. Propagation paths to one of the stations were not predicted by the model, despite being clearly detected. The presence of the signal is attributed to the acoustic energy being trapped in a weak but leaky stratospheric AtmoSOFAR channel. The infrasound signal analysis indicates that the shock originated at the point of the main breakup at an altitude of 35 km. The primary mode of the shock production of the signal detected at the two stations was a spherical blast resulting from the main gross fragmentation episode. The energy estimate, based on the signal period, is 0.17–0.51 kt of TNT equivalent, suggesting a mass of 6.6–23.5 tons. The corresponding object radius, assuming the chondritic origin, was 0.78–1.18 m.

79 ASTRONOMY AND ASTROPHYSICS↗

Real-time gross actinide monitoring in extreme (10 6 :1) beta-gamma background fluid streams using mini-ATMFDs

Real-time, trace alpha activity monitoring in extreme mixed radiation fields is a well-known challenge in fission product-rich streams containing mixtures of beta-gamma-alpha radiation-emitting isotopes, arising due to the extreme 10 6 :1 ratio of background beta-gamma vs. alpha radiation activity levels. This paper presents results from a novel approach to addressing this challenge, based on Tensioned Metastable Fluid Detector (TMFD) sensing technology. An acoustically driven technique was developed and tested at Purdue University (using dissolved Rn-Po and U-chain alpha emitters), and validated (with dissolved Pu) at Savannah River National Laboratory for real-time monitoring of alpha activity in the 100–10,000 dpm/mL range.

Environmental Monitoring↗

Seismic Observations of the OSIRIS-REx Sample Return Capsule Reentry: Deployment, Signal Characteristics, and Wavefield Phenomenology

Controlled spacecraft reentries from interplanetary trajectories provide rare, well-characterized hypersonic sources for advancing seismoacoustic observation techniques. Here we present seismic observations of the OSIRIS-REx sample return capsule (SRC) reentry on 24 September 2023, recorded by 16 three-component nodal seismometers deployed near Eureka, Nevada, at ground distances of 7–20 km from the capsule trajectory. Air-to-ground coupled signals are detected at all stations, exhibiting impulsive onsets consistent with ballistic shock arrivals from the descending Mach cone. We characterize the seismic wavefield through signal amplitude, period, waveform cross-correlation, and array processing. Signal periods decrease systematically with increasing distance from the trajectory within the airport array, indicating that higher-frequency content becomes more prominent at greater offsets, opposite to expectations from geometric spreading and atmospheric absorption. Seismic array processing identifies frequency-dependent back-azimuth variations whose origin remains unresolved; possible contributing factors include source geometry, scattering by fine-scale layered structure in the stratosphere, and near-surface effects. These observations document a spatially complex seismic wavefield from a well-characterized hypersonic line source and provide constraints for future modeling of atmospheric propagation and air-to-ground coupling.

OSIRIS-REx↗

Earthquake detection in a simulated lunar regolith using distributed acoustic sensing

Current models of inner lunar geology have largely been inferred from the seismic experiments and observations performed during the Apollo missions that comprised a relatively small number of seismic instruments. Refining constraints on fundamental lunar relationships such as crust-mantle and mantle-core boundaries in the future will require seismic arrays spanning larger epicentral distances. A promising technology for installing dense seismic arrays with minimal human effort is distributed acoustic sensing (DAS), an approach that allows a single length of fiber optic cable to act as hundreds or thousands of sensors when coupled with a DAS interrogator. While terrestrial uses of DAS technology for seismic monitoring rely on burying the cable to maximize fidelity of seismic signal transmission to the fiber, digging meters of trench to bury optical fiber on lunar or planetary surfaces is logistically infeasible. To evaluate DAS signal attenuation due to surface deployment of cable in lunar regolith, we completed earthquake detection analyses that evaluated the sensitivity of an optic-fiber DAS system to seismic signals at different burial depths. We deployed a single-mode fiber in a 10-m open-bottom wooden box filled with a lunar regolith simulant (LRS) with fiber buried at different depths within the LRS and recorded signals for four regional and local earthquakes. The results were used to identify and evaluate signal attenuation in surface-deployed fiber compared to buried fiber in the LRS. Burial depth responses to active-source signals were also evaluated similar to previous studies characterizing DAS sensitivity of surface-deployed fiber. Atmospheric noise was minimal as the cable was deployed in an indoor environment; however, where observed, atmospheric and anthropogenic noise was filtered out using the same bandpass filtering used to identify earthquake events. We found that signal attenuation of the surface-deployed fiber compared to buried fiber was relatively high in active-source experiments but was not consistently observed in earthquake signals. That burial depth is not highly correlated to attenuation of the observed earthquake signals indicates that in a noise-limited environment, placing DAS-interrogated fiber directly at the regolith surface may be a promising deployment strategy to consider for sensing remote seismic signals during lunar exploration.

58 GEOSCIENCES↗

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗