Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface estimation”

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

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

At least 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Infrastructure-Relevant Storms of the Last Century

Historical extreme storm events are widely used to make Probable Maximum Precipitation (PMP) estimates, which form the cornerstone of large water management infrastructure safety. Past studies suggest that extreme precipitation processes can be sensitive to land surface feedback and the planetary warming trend, which makes the future safety of large infrastructures questionable given the projected changes in land cover and temperature in the coming decades. In this study, a numerical modeling framework was employed to reconstruct 10 extreme storms over CONUS that occurred during the past 100 years, which are used by the engineering profession for PMP estimation for large infrastructures such as dams. Results show that the correlation in daily rainfall for such reconstruction can range between 0.4 and 0.7, while the correlation for maximum 3-day accumulation (a standard period used in infrastructure design) is always above 0.5 for post-1948 storms. This suggests that current numerical modeling and reanalysis data allow us to reconstruct big storms after 1948 with acceptable accuracy. For storms prior to 1948, however, reconstruction of storms shows inconsistency with observations. Our study indicates that numerical modeling and data may not have advanced to a sufficient level to understand how such old storms (pre-1948) may behave in future warming and land cover conditions. However, the infrastructure community can certainly rely on the use of model reconstructed extreme storms of the 1948-present period to reassess safety of our large water infrastructures under assumed changes in temperature and land cover.

Chen, Xiaodong↗
DOE OSTI2020

A highly efficient uranium grabber derived from acrylic fiber for extracting uranium from seawater

Acrylic fiber can be chemically converted to an am- idoxime and carboxylate containing chelating adsorbent for extraction of uranium from seawater. At an optimized ratio of amidoxime:carboxylate, the chelating fiber in real seawater shows a higher uranium, shorter saturation time and lower vanadium adsorption capacity relative to similar high-surface area chelating fibers developed recently using a radiation-induced grafting method. This simple and low-cost synthesis method can be scaled up to mass produce the chelating fiber for recovering metals from aquatic environments including production of uranium from seawater.

Pan, H B.↗
DOE OSTI2023

Quantifying Site Effects and Their Influence on Earthquake Source Parameter Estimations Using a Dense Array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

58 GEOSCIENCES↗
DOE OSTI2021

Deforestation strengthens atmospheric transport of mineral dust and phosphorus from North Africa to the Amazon

Phosphorus contained in atmospheric mineral dust aerosol originating from Africa fertilizes tropical forests in Amazonia. However, the mechanisms influencing this nutrient transport pathway remain poorly understood. Here we use the Community Earth System Model to investigate how large-scale deforestation affects mineral dust aerosol transport and deposition in the tropics. We find that the surface biophysical changes that accompany deforestation produce a warmer, drier and windier surface environment that perturbs atmospheric circulation and enhances long-range dust transport from North Africa to the Amazon. Tropics-wide deforestation weakens the Hadley circulation, which in turn, leads to a northward expansion of the Hadley cell and increases surface air pressure over the Sahara Desert. Additionally, the high pressure anomaly over the Sahara, in turn, increases northeasterly winds across North Africa and the tropical North Atlantic Ocean, which subsequently increases dust transport to the South America continent. We estimate that the annual atmospheric phosphorus deposition from dust significantly increases by 27% (P < 0.01) in the Amazon under a scenario of complete deforestation. These interactions exemplify how land surface changes can modify tropical nutrient cycling, which in turn, may have consequences for long-term changes in tropical ecosystem productivity and biodiversity.

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

Molecular chemisorption of N 2 on IrO 2 (110)

We investigated the adsorption of N 2 on stoichiometric and O-rich IrO 2 (110) surfaces using temperature programmed desorption (TPD) experiments and density functional theory (DFT) calculations. TPD shows that N 2 desorbs predominantly from the stoichiometric-IrO 2 (110) surface in a well-defined peak at 270 K for N 2 coverages below about 0.5 ML, and that a shoulder centered near 235 K develops in the N 2 TPD traces as the coverage approaches saturation, indicating that adsorbed N 2 molecules become destabilized at high N 2 coverage. Experiments of N 2 adsorption onto O-rich IrO 2 (110) surfaces provide evidence that N 2 adsorbs exclusively on the coordinatively-unsaturated Ir atoms (Ircus) of the surface, and that pre-adsorbed O-atoms (“on-top” oxygen) stabilize adsorbed N 2 molecules, causing the main N 2 TPD peak to shift toward higher temperature with increasing oxygen coverage. Consistent with prior results, our DFT calculations predict that an N 2 molecule preferentially adsorbs into an upright configuration on an Ircus atom of the IrO 2 (110) surface and achieves a binding energy of about 100 kJ/mol. The computed binding energy agrees well with our experimental estimate of ~90 kJ/mol for low N 2 coverages on stoichiometric IrO 2 (110). Lastly, the DFT calculations also quantitatively reproduce the observed stabilization of N 2 by co-adsorbed on-top O-atoms and predict the destabilization of N 2 on IrO 2 (110) as the N 2 adlayer becomes crowded at high coverage.

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

Predictive Risk Mapping of Schistosomiasis in Madagascar Using Ecological Niche Modeling and Precision Mapping

Schistosomiasis is a neglected tropical disease (NTD) found throughout tropical and subtropical Africa. In Madagascar, the condition is widespread and endemic in 74% of all administrative districts in the country. Despite the significant burden of the disease, high-resolution risk maps have yet to be produced to guide national control programs. This study used an ecological niche modeling (ENM) and precision mapping approach to estimate environmental suitability and disease transmission risk. The results show that suitability for schistosomiasis is widespread and covers 264,781 km 2 (102,232 sq miles). Covariates of significance to the model were the accessibility to cities, distance to water, enhanced vegetation index (EVI), annual mean temperature, land surface temperature (LST), clay content, and annual precipitation. Disease transmission risk is greatest in the central highlands, tropical east coast, arid-southwest, and northwest. An estimated 14.9 million people could be at risk of schistosomiasis; 11.4 million reside in rural areas, while 3.5 million are in urban areas. This study provides valuable insight into the geography of schistosomiasis in Madagascar and its potential risk to human populations. Because of the focal nature of the disease, these maps can inform national surveillance programs while improving understanding of areas in need of medical interventions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Sea spray aerosol production flux retrieval based on Doppler lidar measurements

Supermicron sea spray aerosols (SSA) play a crucial role in stratocumulus cloud microphysics and drizzle formation. However, our understanding of SSA's spatiotemporal distribution, production mechanisms, and fluxes in the marine boundary layer remains limited. This study introduces a novel approach to determine supermicron SSA production flux and number concentration at various heights above the ocean surface. Our method leverages Doppler lidar data, specifically attenuated backscatter and vertical velocity, collected at a 105 m range gate during the Eastern Pacific Clouds and Aerosols for Precipitation Experiment (EPCAPE) at Scripps Pier, La Jolla, California. To focus on nascent SSA, we analyzed data from periods with wind directions between 225° and 315° and wind speeds exceeding 4 m s −1 . Cloud and precipitation events were excluded using relative humidity, backscatter thresholds, and 2D Video Disdrometer observations. Results show that calculated supermicron SSA number concentrations at 105 m ranged from 0.26 to 2.93 cm −3 for surface wind speeds between 4 and 6 ms −1 . When exceeding the method's Limit of Detection, estimated production fluxes ranged from 0.20 to 1.53 cm −2 s −1 . This methodology has potential for application in various oceanic locations, likely enhancing our understanding of SSA's role in atmospheric processes and climate interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Atom Probe Tomography Analysis of Mica

Laser-assisted atom probe tomography (APT) is a relatively new, powerful technique for sub-nanometric mineral and biomineral analysis. However, the laser-assisted APT analysis of highly anisotropic and chemically diverse minerals, such as phyllosilicates, may prove especially challenging due to the complex interaction between the crystal structure and the laser pulse upon applying a high electric field. Micas are a representative group of nonswelling clay minerals of relevance to a number of scientific and technological fields. In this study, a Mg-rich biotite was analyzed by APT to generate preliminary data on nonisotropic minerals and to investigate the effect of the crystallographic orientation on mica chemical composition and structure estimation. Here, the difference in results obtained for specimens extracted from the (001) and ( hk 0) mica surfaces indicate the importance of both experimental parameters and the crystallography. Anisotropy of mica has a strong influence on the physicochemical properties of the mineral during field evaporation and the interpretation of APT data. The promising results obtained in the present study open the way to future innovative APT applications on mica and clay minerals and contribute to the general discussion on the challenges for the analysis of geomaterials by atom probe tomography.

3D reconstructions↗
DOE OSTI2026

Nanoionics Drastically Accelerating Mass Transfer at Elevated Temperatures over 750 °C

Nanoionics were previously considered thermally unstable and infeasible for devices operating above 500 °C. Here, we elucidate the design principle for establishing stable nanoionics from various oxides. We utilized reversible solid oxide cells (SOCs) as the test bed and implemented nanoionics using atomic layer deposition (ALD). We demonstrate a straightforward, interface-controlled, practical approach to render a conformal, ∼15 nm thick ALD film, which initially thermodynamically favors the formation of a solid solution with the substrate into surface nanoionics with single or double layers of nanograins with random crystal orientations. The nanoionics exhibited conductivity estimated to be 7 orders of magnitude higher than that of their bulk-scale counterpart. They demonstrated conformability with uniform grain sizes of ∼15 nm, even after electrochemical operation for ∼500 h at 750 °C and 1000 h at 850 °C. The thermal stability and conductivity of such nanoionics represent a conceptual and technological framework in nanoionics.

atomic layer deposition↗
DOE OSTI2026

As above, not so below: Ion fractionation in planetary analog ices

The geophysical evolution and astrobiological potential of ocean worlds are indelibly linked to the chemical compositions of their oceans and ice shells. In the absence of direct measurements, empirical estimates of subsurface ocean compositions have relied on the assumption that the ionic compositions of ice shell surfaces are representative of their underlying ocean compositions. Here, we present experimental results demonstrating that ion fractionation—the differential entrainment of ion species in forming ices—is likely a prevalent process on ocean worlds, suggesting that planetary ice shell compositions do not directly reflect their underlying ocean compositions. We measure in-ice depletions and amplifications of relative ion concentrations ranging between −40 and +77%, compared to the parent fluid composition. Although this may complicate the interpretation of spacecraft data, ion fractionation provides a mechanism for generating compositionally diverse ices that could help explain the geological complexity of planetary ice shells.

58 GEOSCIENCES↗
DOE OSTI2026

Rapid Analyses of Sparse Seismoacoustic Data Reveals the Timing and Size of the Accurate Energetic Systems Explosion

On 10 October 2025 an explosion occurred at a facility operated by Accurate Energetic Systems in Humphreys County, Tennessee. The incident resulted in 16 fatalities and created a debris field over several square kilometers. To address remaining questions about explosion timing and size, we collected about 20 seismic and 19 acoustic records of the blast from sensors up to hundreds of kilometers away. We then deployed 10 distinct physics-based, reduced order models (ROMs) that used validated geological structure and atmospheric conditions from the time of the event, along with observations of body- and surface-wave energy, as well as acoustic overpressure and phase duration. Each ROM predicted either timing, yield estimates, or both. We binned these estimates and their uncertainties according to each ROMs’ assumptions about confinement (aboveground, buried fully coupled, and buried partially coupled) and combined these estimates with other forensic data to conclude that the event occurred as a single, aboveground explosion on 10 December 2025 12:47:50.8 ±0.1 s with a yield equivalent to 11.8 [2.3,16.5] tons of Trinitrotoluene. Our estimates align with the Bureau of Alcohol, Tobacco, Firearms and Explosives inventory reports of 11–13 tons. This multimethod approach demonstrates the use of remotely observed geophysical data to rapidly aid conventional forensic investigations of accidental explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Compact, high-power superconducting radio-frequency accelerators for environmental applications

Electron-beam irradiation has been proven to destructively reduce or eliminate a wide variety of organic chemicals, viruses and bacteria from wastewater, as well as reducing sulfur and nitrous oxides emission from coal-fired power plants. It is estimated that there are approximately 30,000 such particle accelerators in use worldwide for industrial processes including surface and bulk processing of material, medical sterilization, and environmental remediation. Maximum beam power from accelerators used in these applications is currently limited to less than 500 kW, and a higher beam power is needed to reduce treatment costs. The market availability of continuous-wave (CW) electron-beam accelerators with power of the order of ~100 kW is also very limited. Superconducting radio-frequency (SRF) linear accelerators (linacs) are commonly used at basic research laboratories throughout the world due to their exceptionally high efficiency as compared to current industrial varieties. Recent advances in cryogenics, SRF thin films and high-power magnetrons allow for the design of increasingly compact and efficient CW SRF electron linacs in the energy range 1-10 MeV, with up to 1 MW of beam power. A new class of compact, high-efficiency electron-beam accelerators may provide cost-effective solutions for a range of industrial and environmental remediation applications, particularly with respect to tackling one such class of contaminants so-called “forever chemicals”, including per- and polyfluoroalkyl substances (PFASs), which are ubiquitous in a wide range of products and for which there is currently no effective destruction technology. We have recently demonstrated the key technologies required for such accelerators by: 1) operating a Nb3Sn SRF accelerating cavity cooled by commercial cryocoolers up to an accelerating gradient of 12.4 MV/m and 2) demonstrating the phase-locking as well as a high-efficiency power combining scheme for industrial magnetron transmitters. In this presentation, we will introduce the principles and benefits of MW-class SRF linacs, based on the conduction-cooled SRF technology that we have demonstrated for environmental remediation. We propose the development of a 4 MeV, 20 kW prototype to be built at Jefferson Lab as a first accelerator demonstrator unit. We will also present the results from samples study on the effect of electron-beam irradiation on so-called “forever chemicals” such as 1,4-dioxane and PFAS, using an existing multi-purpose 10 MeV, low-power, CW SRF linac at Jefferson Lab.

Ciovati, Gianluigi↗
DOE OSTI2025

Diffusion kinetics of 3 He in pyroxene and plagioclase and applications to cosmogenic exposure dating and paleothermometry in mafic rocks

In this study, we investigate the diffusivity of cosmogenic 3 He in a variety of plagioclase and pyroxene compositions, and its application to paleothermometry and exposure dating in these minerals, through stepwise degassing experiments. While cosmogenic 3 He has been utilized for exposure dating in pyroxene for decades due to its retentivity, plagioclase, often found along with pyroxene in mafic rocks, is generally less retentive of cosmogenic noble gas. However, the diffusivity of 3 He in either plagioclase or pyroxene has not yet been measured quantitatively. A challenge in measuring diffusion kinetics by step-degassing experiments in poorly retentive minerals is the fact that significant amounts of He can be lost prior to the experiment. To address this issue, we apply a forward “multiple diffusion domain” (MDD) inversion model that includes model predictions of initial gas loss during irradiation and storage of the samples to account for this observation and add constraints to the diffusion parameters. We find that 3 He diffusivity in plagioclase appears to be highly variable. This variability can be explained by the MDD inversion models' inability to constrain the diffusion parameters when significant gas has been lost during irradiation and/or prolonged storage prior to experiment analysis, resulting in an overestimation of 3 He retentivity. Plagioclase samples that were kept frozen after irradiation, to limit this initial gas loss, yielded the most reliable estimate of diffusion kinetics. We find that 3 He in plagioclase is diffusively lost at Earth's surface temperatures on a timescale of 100 years and is therefore unsuitable for surface temperature paleothermometry. On the contrary, we find cosmogenic 3 He in pyroxene to be retentive at Earth's surface temperatures on a 1-million-year timescale.

Bergelin, Marie [Berkeley Geochronology Center (BG↗
DOE OSTI2021

Cutting out the middleman: calibrating and validating a dynamic vegetation model (ED2-PROSPECT5) using remotely sensed surface reflectance

Canopy radiative transfer is the primary mechanism by which models relate vegetation composition and state to the surface energy balance, which is important to light- and temperature-sensitive plant processes as well as understanding land–atmosphere feedbacks. In addition, certain parameters (e.g., specific leaf area, SLA) that have an outsized influence on vegetation model behavior can be constrained by observations of shortwave reflectance, thus reducing model predictive uncertainty. Importantly, calibrating against radiative transfer outputs allows models to directly use remote sensing reflectance products without relying on highly derived products (such as MODIS leaf area index) whose assumptions may be incompatible with the target vegetation model and whose uncertainties are usually not well quantified. Here, we created the EDR model by coupling the two-stream representation of canopy radiative transfer in the Ecosystem Demography model version 2 (ED2) with a leaf radiative transfer model (PROSPECT-5) and a simple soil reflectance model to predict full-range, high-spectral-resolution surface reflectance that is dependent on the underlying ED2 model state. We then calibrated this model against estimates of hemispherical reflectance (corrected for directional effects) from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States. The calibration significantly reduced uncertainty in model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗
DOE OSTI2024

Continuous snow depth, ground interface temperature and shallow soil temperature measurements from 2021-10-1 to 2022-6-14, Seward Peninsula, Alaska

The dataset contains co-located snow depth, ground interface temperature, and shallow soil temperature measured at 98 discrete locations in a watershed located along the Nome-Teller road at mile marker 27 (referred to as T27), and at 53 discrete locations on a hillslope located along the Kougarok road at mile marker 64 (referred to as K64), in Seward Peninsula, Alaska. The dataset aims to understand the local heterogeneity of snow depth, snow temperature, and soil temperature dynamics and their interactions in a discontinuous permafrost region. The dataset is also valuable to train and evaluate machine learning and physical models to predict snow depth or the impact of snow depth on ground surface temperature. At each location, temperatures above and below the ground surface were measured by a pair of vertically deployed distributed temperature profiling probes designed based on Dafflon et al (2022). The probes have high precision temperature sensors spaced at 5 or 10 cm. The mean daily snow depth was estimated by identifying the pair of consecutive sensors with maximum drop of daily temperature high-frequency fluctuations. The ground interface temperature was measured by the sensor located 1-5 cm above the ground surface at 15-minute intervals. The shallow soil temperature was measured by the sensor located 1-5 cm below the ground surface at 15-minute intervals. The dataset includes a description of the probe locations in the "Probe_locations_*.csv" file and the 3 data files (Snow_depths_*.csv, Ground_interface_temperatures_*.csv, Shallow_soil_temperatures_*.csv). * is either T27 or K64, which are the two study sites. In each data file, each column corresponds to a measurement location. Metadata files include data descriptions (_dd.csv) for tabular data. All included files are listed and described in xxxx_flmd.csv. A more detailed description of data processing, along with an updated dataset incorporating multiple seasons and improved snow depth estimation is available at https://doi.org/10.15485/2480365 (Wang et al., 2025a, Wang et al., 2025b).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A New Martian Crater Chronology: Implications for Jezero Crater

Crater chronologies are a fundamental tool to assess the relative and absolute ages of planetary surfaces when direct radiometric dating is not available. Martian crater chronologies are derived from lunar crater spatial densities on terrains with known radiometric ages, and thus they critically depend on the Moon-to-Mars extrapolation. This extrapolation requires knowledge of the time evolution of the impact flux, including contributions from various impactor populations, factors that are not trivially connected to the dynamical evolution of the early Solar System. In this paper, we will present a new Martian crater chronology based on current dynamical models, and consider the main sources of uncertainties (e.g., impactor size–frequency distribution; dynamical models with late and early instabilities, etc.). The resulting “envelope” of Martian crater chronologies significantly differs from previous chronologies. The new Martian crater chronology is discussed using two interesting applications: Jezero crater’s dark terrain (relevant to the NASA Mars 2020 mission) and the southern heavily cratered highlands. Our results indicate that Jezero’s dark terrain may have formed ∼3.1 Ga, i.e., up to 0.5 Gyr older than previously thought. In addition, available crater chronologies (including our own) overestimate the number of craters larger than 150 km on the southern highlands, suggesting either that large craters have been efficiently erased over Martian history or that dynamical models need further refinement. Further, our chronology constrains the age of Isidis basin to be 4.05–4.2 Ga and that of the Borealis basin to be 4.35–4.40 Ga; these are predictions that can be tested with future sample and return missions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Analysis of Seismic Signals Generated by Vehicle Traffic with Application to Derivation of Subsurface Q-Values

Correct identification and modeling of anthropogenic sources of ground motion are of considerable importance for many studies, including detection of small earthquakes and imaging seismic properties below the surface. Here, to understand signals generated by common vehicle traffic, we use seismic data recorded by closely spaced geophones normal to roads at two sites on San Jacinto fault zone. To quantify the spatiotemporal and frequency variations of the recorded ground motions, we develop a simple analytical solution accounting for propagation and attenuation of surface waves. The model reproduces well-observed bell-shaped spectrograms of car signals recorded by geophones close to roads, and it can be used to estimate frequency-dependent Q-values of the subsurface materials. The data analysis indicates Q-values of 3–40, for frequencies up to 150 Hz for road-receiver paths at the two examined sites. The derived Q-values are consistent with attenuation factors of surface waves previously obtained with other methods. The analytical results and analysis procedure provide a highly efficient method for deriving Q-values of shallow subsurface materials.

58 GEOSCIENCES↗