Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RAYLEIGH SCATTERING”

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.

719 records · Page 40

Onset conditions for equatorial spread F

The problem of day-to-day variability in the occurrence of equatorial spread F (ESF) is addressed using multidiagnostic observations and semiempirical modeling. The observational results are derived from a two-night case study of ESF onset conditions observed at Kwajalein Atoll (Marshall Islands) using the ALTAIR incoherent scatter radar and all-sky optical imaging techniques. The major difference between nights when ESF instabilities did not occur (August 14, 1988) and did occur (August 15, 1988) in the Kwajalein sector was that the northern meridional gradient of 6300-A airglow was reduced on the night of limited ESF activity. Modeling results suggest that this unusual airglow pattern is due to equatorward neutral winds. Previous researchers have shown that transequatorial thermospheric winds can exert a control over ESF seasonal and longitudinal occurrence patterns by inhibiting Rayleigh-Taylor instability growth rates. Evidence is presented to suggest that this picture can be extended to far shorter time scales, namely, that 'surges' in transequatorial winds acting over characteristic times of a few hours to a day can result in a stabilizing influence upon irregularity growth rates. The seemingly capricious nature of ESF onset may thus be controlled, in part, by the inherent variability of low-latitude thermospheric winds.

Mendillo, Michael↗

Correlation function metrology for warm dense matter: Recent developments and practical guidelines

X-ray Thomson scattering (XRTS) has emerged as a valuable diagnostic for matter under extreme conditions, as it captures the intricate many-body physics of the probed sample. Recent advances, such as the model-free temperature diagnostic of Dornheim et al. [Nat. Commun. 13 , 7911 (2022)], have demonstrated how much information can be extracted directly within the imaginary-time formalism. However, since the imaginary-time formalism is a concept often difficult to grasp, we provide here a systematic overview of its theoretical foundations and explicitly demonstrate its practical applications to temperature inference, including relevant subtleties. Furthermore, we present recent developments that enable the determination of the absolute normalization, Rayleigh weight, and density from XRTS measurements without reliance on uncontrolled model assumptions. Finally, we outline a unified workflow that guides the extraction of these key observables, offering a practical framework for applying the method to interpret experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear ultrasonic technique for the characterization of microstructure in additive materials

This study employs nonlinear ultrasonic techniques to track microstructural changes in additively manufactured metals. The second harmonic generation technique based on the transmission of Rayleigh surface waves is used to measure the acoustic nonlinearity parameter, β. Stainless steel specimens are made through three procedures: traditional wrought manufacturing, laser-powder bed fusion, and laser engineered net shaping. The β parameter is measured through successive steps of an annealing heat treatment intended to decrease dislocation density. Dislocation density is known to be sensitive to manufacturing variables. In agreement with fundamental material models for the dislocation-acoustic nonlinearity relationship in the second harmonic generation, β drops in each specimen throughout the heat treatment before recrystallization. Geometrically necessary dislocations (GNDs) are measured from electron back-scatter diffraction as a quantitative indicator of dislocations; average GND density and β are found to have a statistical correlation coefficient of 0.852 showing the sensitivity of β to dislocations in additively manufactured metals. Moreover, β shows an excellent correlation with hardness, which is a measure of the macroscopic effect of dislocations.

42 ENGINEERING↗

Atmospheric Correction for Satellite Ocean Color Radiometry

This tutorial is an introduction to atmospheric correction in general and also documentation of the atmospheric correction algorithms currently implemented by the NASA Ocean Biology Processing Group (OBPG) for processing ocean color data from satellite-borne sensors such as MODIS and VIIRS. The intended audience is graduate students or others who are encountering this topic for the first time. The tutorial is in two parts. Part I discusses the generic atmospheric correction problem. The magnitude and nature of the problem are first illustrated with numerical results generated by a coupled ocean-atmosphere radiative transfer model. That code allow the various contributions (Rayleigh and aerosol path radiance, surface reflectance, water-leaving radiance, etc.) to the topof- the-atmosphere (TOA) radiance to be separated out. Particular attention is then paid to the definition, calculation, and interpretation of the so-called "exact normalized water-leaving radiance" and its equivalent reflectance. Part I ends with chapters on the calculation of direct and diffuse atmospheric transmittances, and on how vicarious calibration is performed. Part II then describes one by one the particular algorithms currently used by the OBPG to effect the various steps of the atmospheric correction process, viz. the corrections for absorption and scattering by gases and aerosols, Sun and sky reflectance by the sea surface and whitecaps, and finally corrections for sensor out-of-band response and polarization effects. One goal of the tutorial-guided by teaching needs- is to distill the results of dozens of papers published over several decades of research in atmospheric correction for ocean color remote sensing.

MODIS↗

Resonance scattering from interstellar and interplanetary helium.

The intensity and spatial distribution of the diffuse 584 A glow due to resonance scattering of solar He I, 584 A line radiation by neutral interstellar and interplanetary helium is derived. Two sources of neutral helium are considered: the local interstellar medium interacting with the solar system and dust deionization of the He(++) component of the solar wind. Which of these two sources will dominate the radiation field at 1 AU is dependent upon a number of parameters which have not yet been determined. Results are presented as a function of the width of the solar line, and the temperature and speed of approach of the interstellar medium. The calculated intensities for the most extreme cases vary from approximately .000001 Rayleighs to almost 100 Rayleighs. The relevance of these results to the study of the interplanetary and interstellar medium is discussed.

Paresce, F.↗

Hubble Space Telescope UV spectral observations of Io passing into eclipse

Time-resolved spectra of Io have been obtained with the Faint Object Spectrograph on the Hubble Space Telescope in January 1992 at times centered on the passage of Io into Jupiter's shadow. Two different eclipse observations covered 1100-1600 A and 2250-3300 A. In the far ultraviolet(far-UV) range, emission lines of atomic sulfur and oxygen from Io's atmosphere (similar to those previously detected with the International Ultraviolet Explorer (IUE)) have been observed from Io in sunlight, and the spatial extent of the emitting region has been resolved for the first time: this is 0.5-1 Io radii (R(sub Io)) above the surface. The emission lines are typically 1kR in brightness while Io is in sunlight, and decrease to a few hundred Rayleighs within 20 min or less of Io's passing into shadow. If the emissions are produced in Io's ionosphere, the decrease in shadow appears consistent with the collisional slowing and recombination of photoelectrons in 100-1000 s, with recombination an important quenching process if the dominant ion is molecular (i.e., SO2(+)) condensation, with the residual emission in shadow due either to plasma impact of gas above the hot volcanic calderas or electron impact on S and O. In the near-UV range, we have not detected any airglow emissions from Io's atmosphere in shadow, with the main limitation being a high level of scattered light from Jupiter. We derive a 3 sigma upper limit to the 2560 A SO emission feature of 1 KR, which is close to what is expected from electron impact on SO2 based on the obs erved brightness of the FUV S and O lines in shadow. A high signal-to-noise spectrum of Io's albedo in sunlight reveals a spectral shape similar to laboratory spectra of SO2 frost reflectivity, and the relative albedo spectrum changed as Io passed into eclipse and part of the disk was in shadow. No specific SO2 gas absorption features appear in the albedo spectrum, although there could be substantial gas absorption near 2800 A if the individual lines are narrow and saturated. Finally, we present preliminary models for the near-UV spectrum of Io as functions of SO2 frost areal coverage and SO2 gas density.

Clarke, John T.↗

ARCLITE: the Arctic Lidar Technology Facility at Sondre Stromfjord, Greenland

SRI International is presently developing an Arctic lidar (ARCLITE) facility at the incoherent-scatter radar site in Sondre Stromfjord, Greenland (67.0 degrees N, 209.2 degrees W). The project, funded by NSF through the CEDAR (Coupling Energetics and Dynamics of Atmospheric Regions) initiative, is to perform high-altitude molecular/aerosol backscatter measurements in the Arctic covering the stratosphere/mesosphere region during nighttime and daytime conditions. The lidar system employs an injection seeded, high-energy GCR-5 Spectra Physics Nd:YAG laser for primary operation at 532 nm (550 mJ per pulse at 30 Hz). The receiver system is designed around a 36-inch Cassegrainian telescope of astronomical quality with receiver components and photon-counting electronics for day and night operations. The Rayleigh backscatter measurements will conduct investigations concerning long-term studies of the basic density and temperature structure of the Arctic stratosphere/mesosphere region; effects of stratospheric warming on the structure of the stratosphere/mesosphere region; occurrence and properties of polar stratospheric clouds; occurrence of noctilucent clouds and their properties at high latitudes; and dynamic and thermodynamic coupling between the stratosphere/mesosphere region and the lower thermosphere.

Thayer, J. P.↗

Soft x ray properties of the Geminga pulsar

The ROSAT soft x ray spectrum and pulse profile of the Geminga pulsar are analyzed and interpreted in terms of thermal emission from the surface of the neutron star. The x ray spectrum appears to consist of two blackbody components with T(sub 1) = (5.2 +/- 1.0) x 10 (exp 5) K and T(sub 2) approximately 3 x 10(exp 6) K, respectively. The inferred ratio of surface areas, A(sub 2)/A(sub 1), is approximately 3 x 10(exp -5). Both components are highly modulated at the pulsar rotation period, but the harder x ray pulse is narrower, and leads the main (soft) x ray pulse by about 105 deg of phase. The soft x ray component is interpreted as photospheric cooling of much of the neutron star's surface area, while the small, hot region could be part of the much smaller polar cap heated by energetic particles flowing inward from the magnetospheric accelerator which is responsible for the production of Geminga's gamma rays. Geminga's gamma ray emission is consistent with outer-magnetosphere accelerator models for highly inclined dipoles. These predict the beaming of energetic gamma rays close enough to the star to give copious e(+/-) production in the stellar magnetic field and a large circumstellar pair density from pair inflow toward the surface. These pairs may quench radio emission, and also reflect most of the hard polar cap x rays back to the stellar surface by cyclotron resonance scattering. They are then reemitted from that much larger area at the lower temperature T(sub 1). The single-peaked nature of the x ray pulse and its energy-dependent phase suggest an off-center dipole geometry for the surface magnetic field. Under the assumption that the soft x ray emission comes from the full surface of a neutron star of radius R = 10 km, a distance estimate of (150-400) pc is derived. This range is consistent with the fit interstellar column density of (1.5 +/- 0.5) x 10(exp 20) cm(exp -2). Distances less than 150 pc are probably ruled out both by the lower limit on the column density, and also by the requirement that the Rayleigh-Jeans extrapolation of the soft x ray spectrum not exceed the observed blue flux of the faint optical counterpart. This distance estimate implies that Geminga's efficiency for converting spindown power into gamma-rays is near unity, and that there may be significant beaming of the gamma rays as well. These results tend to bolster the prospect that most of the unidentified high-energy gamma ray sources in the Galactic plane are pulsars, some of which may be radio quiet.

Halpern, J. P.↗

Modeling composite and fluffy grains: The effects of porosity

Recent studies of interplanetary and interstellar dust provide evidence that cosmic dust grains are fluffy, composite objects, highlighting the need for models of the electromagnetic scattering by these grains. Effective medium theory (EMT) with Mie-type series solutions has been used to explore the effects of porosity which would be important in composite dust particles. While this indirect approach is both flexible and computationally efficient, it is not necessarily a good approximation. The need for EMT and its rather restrictive assumptions may be circumvented through a direct computation of the scattering properties via finite element methods, such as the discrete dipole approximation (DDA). Recently, the utility of the DDA method has been advanced significantly through improvements in theory, in numerical algorithms, and in computer hardware. Extensive calculations with the DDA method are used here to examine more directly the effects of porosity. A particular emphasis is placed upon developing a valid methodology. For both solid and porous targets we establish both numerical and physical convergence properties over the range of size parameter that is required for our study. DDA cross sections for grains with a range of porosity are compared to those computed by the EMT/series expansion technique to examine the applicability of several mixing rules, including two extensions of the Bruggeman rule. We show that for particles with Rayleigh vacuum inclusions, the extension proposed by Rouleau & Martin is quite successful. We also investigate the effects of larger, non-Rayleigh vacuum inclusions for various levels of porosity and find that they can be significant.

Wolff, Michael J.↗

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Boundary layer roll circulations during FIRE

The probable mechanism underlying the development of boundary layer roll circulations are studied using wind and temperature profiles measured by the National Center for Atmospheric Research (NCAR) Electra during the stratocumulus phase of the First ISCCP Regional Experiment (FIRE). The expected, or preferred, roll orientations, horizontal wavelengths, and propagation periods are determined by finding the minimum values of the dynamic and thermodynamic forcing parameters, which here are the eddy Reynolds number (Re) and moist Rayleigh number (Ra sub m). These minimum values depend on the height z sub T of the capping temperature inversion and on the values of the Fourier coefficients of the background height-dependent vector wind profile. As input to our nonlinear spectral model, descent and ascent runs by the Electra provide for initial estimates of the inversion height and the wind profiles. In the first phase of the investigation presented here, a mechanism is said to be a probable contributor to the development of roll circulations within the stratocumulus-topped boundary layer if the modeled roll orientation and wavelengths agree with their observed values. Preliminary results using the 14-coefficient model of Haack-Hirschberg (1988) are discussed for the 7 July 1987 Electra Mission 188-A (Flight 5). This mission was flown across a sharp cloud boundary that was within a LANDSAT/SPOT scene. The stratocumulus deck was relatively solid in the eastern part of the scene, while there was a rapid decrease in cloud cover to scattered cumulus clouds aligned in streets to the west. These cloud streets were oriented nearly parallel to the mean wind direction in the layer, which was approximately 340 degrees. The hypothesis that roll circulations occurred in both the relatively clear and the cloudy regions is investigated using as model input a descent profile obtained in the relatively clear air and an ascent profile obtained in the cloudy air. Initial results for the clear air case are that the pure inflection point mode is not possible and the pure thermal mode was oriented 35 degrees to the right of the mean wind direction. The origin of this unacceptably large discrepancy between the observed and modeled results will be investigated further and the conclusions reported at the next FIRE workshop.

Shirer, Hampton N.↗

Novel High Resolution High Temperature In-Line Sensors for Steel Manufacturing

To support efficiency, productivity, yield improvements, and future industry 4.0 and SMART manufacturing objectives for a competitive and prosperous steel industry, this program has embarked on an effort to develop, demonstrate and deploy several novel sensing technologies based on fiber optics for use in production steel facilities. Over the duration of this program, our team has successfully demonstrated three of these technologies: (1) near continuous Rayliegh Scattering optical frequency domain reflectometry (OFDR) sensing with single mode silica fibers for temperature sensing to 700 C, (2) semi-distributed fiber Bragg grating (FBG) sensing with multimode sapphire fibers for temperature sensing to >1600 C, and (3) remote in-situ Raman analysis for slag and flux chemistry analysis at steelmaking temperatures >1550 C. Each of these sensor and interrogation systems was developed, refined, and tested in our labs at Missouri S&T and then successfully deployed at SSAB’s production facilities at two sites, one in Montpelier, IA and one in Mobile, AL.

36 MATERIALS SCIENCE↗

Multiplexable intrinsic Fabry-Pérot interferometers inscribed by femtosecond laser for vibration measurement in high temperature environments

This article reports fabrication technique and demodulation algorithm developments to use multiplexable intrinsic Fabry-Pérot interferometer (IFPI) fiber sensor array for distributed vibration measurements at high temperatures. Using femtosecond laser direct writing scheme, IFPI array were fabricated through laser-induced Type II scattering points in single mode fiber cores. Reflection spectra of IFPI array were demodulated in real-time using modified Bunemen frequency analysis. The demodulation algorithm, which was implemented using a photodetector-array spectrometer, achieves 64-nε dynamic strain resolution at 1-kHz spectral acquisition rate and 2-kHz maximum vibration bandwidth. Performance and stabilities of IFPI sensor array were characterized from room temperatures to 800 °C. The static strain resolution was 0.6 με and the minimum detectable dynamic strain amplitude was 23 nε/√Hz at 800 °C. The multiplex performance was also tested by measuring dynamic strain in time domain through six IFPI sensor array fabricated in one fiber. This paper presents an integrated and low-cost sensing solution to perform distributed vibration measurements in harsh environments.

42 ENGINEERING↗

The Time Series Technique for Aerosol Retrievals over Land from MODIS: Algorithm MAIAC

Atmospheric aerosols interact with sun light by scattering and absorbing radiation. By changing irradiance of the Earth surface, modifying cloud fractional cover and microphysical properties and a number of other mechanisms, they affect the energy balance, hydrological cycle, and planetary climate [IPCC, 2007]. In many world regions there is a growing impact of aerosols on air quality and human health. The Earth Observing System [NASA, 1999] initiated high quality global Earth observations and operational aerosol retrievals over land. With the wide swath (2300 km) of MODIS instrument, the MODIS Dark Target algorithm [Kaufman et al., 1997; Remer et al., 2005; Levy et al., 2007] currently complemented with the Deep Blue method [Hsu et al., 2004] provides daily global view of planetary atmospheric aerosol. The MISR algorithm [Martonchik et al., 1998; Diner et al., 2005] makes high quality aerosol retrievals in 300 km swaths covering the globe in 8 days. With MODIS aerosol program being very successful, there are still several unresolved issues in the retrieval algorithms. The current processing is pixel-based and relies on a single-orbit data. Such an approach produces a single measurement for every pixel characterized by two main unknowns, aerosol optical thickness (AOT) and surface reflectance (SR). This lack of information constitutes a fundamental problem of the remote sensing which cannot be resolved without a priori information. For example, MODIS Dark Target algorithm makes spectral assumptions about surface reflectance, whereas the Deep Blue method uses ancillary global database of surface reflectance composed from minimal monthly measurements with Rayleigh correction. Both algorithms use Lambertian surface model. The surface-related assumptions in the aerosol retrievals may affect subsequent atmospheric correction in unintended way. For example, the Dark Target algorithm uses an empirical relationship to predict SR in the Blue (B3) and Red (B1) bands from the 2.1 m channel (B7) for the purpose of aerosol retrieval. Obviously, the subsequent atmospheric correction will produce the same SR in the red and blue bands as predicted, i.e. an empirical function of 2.1. In other words, the spectral, spatial and temporal variability of surface reflectance in the Blue and Red bands appears borrowed from band B7. This may have certain implications for the vegetation and global carbon analysis because the chlorophyll-sensing bands B1, B3 are effectively substituted in terms of variability by band B7, which is sensitive to the plant liquid water. This chapter describes a new recently developed generic aerosol-surface retrieval algorithm for MODIS. The Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm simultaneously retrieves AOT and surface bi-directional reflection factor (BRF) using the time series of MODIS measurements.

Lyapustin, Alexei↗

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

A Unified Interpretation of Variability in Precipitation Isotope Ratios

Abstract Several mechanisms have been proposed to explain why the isotope ratios of precipitation vary in space and time and why they correlate with other climate variables like temperature and precipitation. Here, we argue that this behavior is best understood through the lens of radiative transfer, which treats the depletion of atmospheric vapor transport by precipitation as analogous to the attenuation of light by absorption or scattering. Building on earlier work by Siler et al., we introduce a simple model that uses the equations of radiative transfer to approximate the two-dimensional pattern of the oxygen isotope composition of precipitation ( δ p ) from monthly mean hydrologic variables. The model accurately simulates the spatial and seasonal variability in δ p within a state-of-the-art climate model and permits a simple decomposition of δ p variability into contributions from gradients in evaporation and the length scale of vapor transport. Outside the tropics, δ p is mostly controlled by gradients in evaporation, whose dependence on temperature explains the positive correlation between δ p and temperature (i.e., the temperature effect). At low latitudes, δ p is mostly controlled by gradients in the transport length scale, whose inverse relationship with precipitation explains the negative correlation between δ p and precipitation (i.e., the amount effect). This suggests that the temperature and amount effects are both mostly explained by the variability in upstream rainout, but they reflect distinct mechanisms governing rainout at different latitudes. Significance Statement The isotopic composition of precipitation has long been used to make inferences about past climates based on its observed relationship with precipitation in the tropics and with temperature at higher latitudes. These relationships—known as the “amount effect” and “temperature effect,” respectively—have been attributed to many different mechanisms, most of which are thought to operate at either high or low latitudes but not both. Here, we present a unified framework for interpreting the isotope variability that can explain the latitude dependence of the temperature and amount effects despite making no distinction between high and low latitudes. Although our results are generally consistent with certain interpretations of the amount effect, they suggest that the temperature effect is widely misunderstood.

54 ENVIRONMENTAL SCIENCES↗