Spectral Ratio Research Articles

Spectral ratioing of [formula omitted] to constrain the dust particle size distribution of comets

A numerical model of cometary dust environments is used to gain a deeper understanding of the relationship between the brightness (Afρ) and the dust particle size distribution in the coma. Specifically, the spectral ratio of Afρ(425 nm)/Afρ(900 nm) is modelled for a wide range of parameters and tied to the power-law index. The studied parameters are dust composition, terminal outflow velocity and the dust production rate day–night asymmetry. We find that the spectral ratio of Afρ modelled at 425 nm and 900 nm correlates with the power-law index of the particle size distribution. This method could be used to place constraints on the dust size distributions of comets and as a result improve the use of Afρ as a proxy for cometary activity. Optically red dust indicates that the scattering is dominated by large particles.

Bayesian Network Inference for Low-Magnitude Nonnatural Seismic Event Discrimination

Abstract In response to the gaps in understanding the causal relationship between seismic waveform features and the types of seismic events, this research is focused on seismic events of low magnitude (ML≤3.0) in the North China region. Using the Bayesian network theory, we conduct an analysis to infer event types for natural earthquakes, artificial explosions, and mining collapses, and the outcomes achieved notable efficacy for the discrimination of seismic events. Through the analysis of seismic waveforms from 1818 events, we systematically extracted and quantified 55 features in temporal, spectral, and energy domains, which were then recoded as node variables for subsequent analysis. The new data set was subject to select nodes with strong associations to the node type. Subsequently, Bayesian network topologies were constructed using three different algorithms to reconstruct the custom network, calculating posterior probabilities and marginal probabilities. Simultaneously, an extensive evaluation with precision–recall curves of the network structure was carried out, encompassing accuracy, precision, recall, and F1-score. Ultimately, sensitivity analysis was performed on each node to reveal the extent of the influence of node variations on the inference of the node type. The findings showed that the sensitivity of discrimination of seismic events was notably high for several features, including high-frequency P/S spectral ratio values (11 to ∼20 Hz), central frequency, dominant frequency, average frequency, rise and decay average frequency, the real part of the complex cepstral coefficients, peak ground acceleration, and zero crossing. In the classification of natural earthquakes, artificial explosions, and mining collapses, it was observed that the probability of mining collapses was maximized when peak ground acceleration was less than 1526.08, and concurrently, the P/S spectral ratio (11 to ∼20 Hz) fell within the range of −0.25 to −0.02.

Current spectral norm and phase variation based fault region identification for active distribution network

The paper presents a fault region identification method for the active distribution network (ADN) with limited PMU. First, PMU configuration and region division strategies are proposed based on the network topology. Next, the difference in three-phase current variations between the normal and fault regions of the ADN is analyzed. A multi-dimensional state monitoring matrix is built using the measured current data. The spectral norm ratio coefficient is constructed based on the 2-norm to lower the complexity of the multi-dimensional state monitoring matrix and quantify the difference in state changes before and after the fault occurs in each region. Then, the fault region is identified by combining each region’s spectral norm ratio coefficient and the change of the current phase. Finally, an IEEE 33-node simulation model is created to simulate faults under different working conditions. According to the simulation results, the suggested approach is less impacted by fault type, neutral point grounding mode, and transition resistance. Furthermore, even if the communication does not match the rigorous synchronization requirements, the proposed method can still complete the fault identification of the distribution network correctly and has high robustness.

Seismic vulnerability distribution of the earthquake prone area in Central Aceh

Abstract Central Aceh is known as an earthquake-prone area due to the presence of the Sumatra Fault and secondary fault crossing the area resulting the M 6.1 earthquake in 2013. The objective of the research is to investigate seismic vulnerability of Central Aceh. Microtremor observations were conducted at 25 locations in Central Aceh analysed using the Horizontal to Vertical Spectral Ratio (HVSR). Seismic amplification is represented by the maximum values of H/V. The dominant frequency f0 representing the sediment thickness was obtained when the H/V curve is maximum. The frequency values ranged from 14.9 to 25.9 Hz, and the seismic vulnerability values varied from 0.02 to 3.1. The top soil consists of thin layers of sedimentary material dominated by tertiary rocks. The rocks are classified as Type IV and Class I based on the dominant frequency values. This result is consistent with the geology, consisting of sedimentary and volcanic rocks. Areas with a seismic vulnerability value > 1.0 is considered to have a high risk of earthquake damage because they are likely to experience earthquake amplification. Areas with seismic vulnerability value > 1.0 are found in West Aceh and the Geureudong volcano area.

Revealing crack initiation and extension of brittle rock via time-lapse acoustic attenuation: Insights from active ultrasonic experiments

Predicting brittle rock failure is essential for rock mechanics research and early disaster alerts in deep rock engineering. Successful prediction relies on identifying crack initiation and extension. However, currently, there is no universally accepted method for accurately identifying and tracking crack evolution to predict brittle failure of the surrounding rock. In this study, the propagation parameters of transmission ultrasonic waves were used to characterize crack behaviors at the laboratory scale. The differential acoustic attenuation term was proposed by modifying the spectral ratio method to evaluate the dissipation capability of the corresponding medium to the transmitted ultrasonic waves. In addition to acoustic attenuation, velocity and amplitude attenuation were introduced to characterize the damage evolution during brittle failure. The indicators were analyzed by conducting a uniaxial compression test on a sandstone sample with a single flaw. During compression, active ultrasonic surveys were conducted to attain the propagation parameters. The two-dimensional digital image correlation method was simultaneously applied to capture the full-field strain on the front surface of the sample. The shear and tensile strain profiles indicated that the shear damage initially accumulated at the tips of the flaw, causing a reduction in the velocity and amplitude of the transmitted ultrasonic waves. The acoustic attenuation indicator in the flaw-tip regions increased before the shear strain began. As the damage developed further, the acoustic attenuation increased overall, with numerous local fluctuations in the flaw-tip regions due to stress rearrangement and damage development, while the velocity and amplitude remained consistently low. In contrast, intact regions that were far away from the tips of the flaw exhibited different characteristics. The velocity and amplitude variations were more staged, while the acoustic attenuation changed gradually. The results indicated that the variations in the propagation parameters exhibited specific patterns before the damage was initiated, offering a potential for predicting crack initiation. The variations in the acoustic attenuation, velocity, and amplitude attenuation allow for tracking the crack propagation.

Quantitative electroencephalography in cerebral amyloid angiopathy

ObjectiveWe investigated whether quantitative electroencephalography (qEEG) correlates with cognition and cortical superficial siderosis (cSS) in cerebral amyloid angiopathy. MethodsWe included patients with sporadic (sCAA) and hereditary Dutch-type CAA (D-CAA). Spectral measures and the phase lag index (PLI) were analyzed on qEEG. Cognition was assessed with the MoCA and cSS presence was scored on 3T-MRI. Linear regression analyses were performed to investigate these qEEG measures and cognition. Independent samples T-tests were used to analyze the qEEG measure differences between participants with and without cSS. ResultsWe included 92 participants (44 D-CAA; 48 sCAA). A lower average peak frequency (β[95 %CI] = 0.986[0.252–1.721]; P = 0.009) and a higher spectral ratio (β[95 %CI] = −0.918[−1.761–−0.075]; P = 0.033) on qEEG correlated with a lower MoCA score, irrespective of a history of symptomatic intracerebral hemorrhage (sICH). The PLI showed no correlation to the MoCA. qEEG slowing was not different in those with or without cSS. ConclusionsSpectral qEEG (but not PLI) reflects cognitive performance in patients with CAA with and without a history of sICH. We found no association between qEEG slowing and cSS. SignificanceqEEG could be a valuable biomarker, especially in challenging cognitive testing situations in CAA, and a potential predictive tool in future studies.

Application of Alkaline Deep Eutectic Solvents as a Green Alternative to the Traditional Extractants for the Isolation of Humic Substances.

The presented study focused on the possibility of using alkaline deep eutectic solvents (ADESs) as green extractants for the isolation of humic substances (HSs) from peat and lignite in a process intensified by ultrasound. For this purpose, the extraction procedure was statistically described on the basis of the Box-Behnken design, and the carboxyl group content in the obtained products was optimized due to the ADES composition, ultrasound intensity, and extraction time. For optimal extraction conditions, the experimental carboxyl content in the isolated products was equal to 3.71 and 2.96 mmol g-1 for the HSs extracted from peat and lignite, respectively. These values were similar to the results for the reference samples, which were HSs extracted using 0.1 M NaOH, as well as humic acids and sodium humates purchased from Sigma-Aldrich. The qualitative assessment of the products obtained was based on spectroscopic methods, including FTIR, 1H NMR, and UV-vis. The analyses carried out for the isolated samples revealed the characteristic structures of HSs, including components of aliphatic chains and aromatic core as well as carboxyl, ester, and amino groups. Simultaneously, the results of the spectral ratio of E 280/E 472 showed the significant differences between the relative amount of lignin for the samples tested.

Simulation method for multi-source information fusion space target

Current space target simulation systems suffer from limitations, such as a single simulated spectral band, inability to allocate spectral ratios between bands, and poor imaging quality of multi-spectral simulated images. We propose a multi-source information fusion spatial target simulation method system that features a “series-parallel” hybrid connection. We optimize the design of a multi-spectral common aperture optical system and experimentally verify the proposed method. The experimental results show that the proposed system can simulate a wide range of spectral bands ranging from ultraviolet to long-wave infrared. Furthermore, it achieves precise control over the ratio of simulated radiant energy among spectral bands, with a spectral simulation error of less than 4.8%. Notably, it effectively images stars against the cosmic background and spacecraft across different spectral bands. This system serves as a theoretical cornerstone and provides crucial technological support for performance testing and rapid iterative development of multi-source information fusion space exploration systems.

Quantitative analysis of laser-generated ultrasonic wave characteristics and their correlation with grain size in polycrystalline materials

Abstract The quantitative relationship between nanosecond pulsed laser parameters and the characteristics of laser-generated ultrasonic waves in polycrystalline materials was evaluated. The high energy of the pulsed laser with a large irradiation spot simultaneously generated ultrasonic longitudinal and shear waves at the epicenter under the slight ablation regime. An optimized denoising technique based on wavelet thresholding and variational mode decomposition was applied to reduce noise in shear waves with a low signal-to-noise ratio. An approach for characterizing grain size was proposed using spectral center frequency ratio (SCFR) based on time-frequency analysis. The results have demonstrated that the generation regime of ultrasonic waves is not solely determined by the laser power density; even at high power densities, a high energy with a large spot can generate an ultrasonic waveform dominated by the thermoelastic effect. It is ascribed to the intensification of the thermoelastic effect with the proportional increase in laser irradiation spot area for a given laser power density. Furthermore, both longitudinal and shear wave SCFRs are linearly related to grain size in polycrystalline materials; however, the shear wave SCFR is more sensitive to finer-grained materials. This study holds great significance for evaluating metal material properties using laser ultrasound.

Evaluation of vertical-to-horizontal spectral acceleration ratio for the 2023 Turkey-Syria seismic sequences

This paper aims to evaluate the vertical-to-horizontal, V/H, spectral acceleration ratios from the earthquake ground motions recorded during the February 6, 2023 Turkey-Syria seismic events, which include doublet events with moment magnitudes of Mw7.7 and Mw7.6 that struck the city of Kahramanmaraş and a considerable aftershock with Mw6.7 that mainly struck the city of Nurdağı. Statistics of the observed V/H spectral ratios computed from records that have peak ground acceleration larger than 0.1g showed that most records of the three earthquakes yielded spectral V/H ratios higher than 2/3, which is a typical value recommended by several seismic codes. Additionally, comparison of the observed V/H spectra with respect to the Ground Motion Predictive Models, GMPMs, for V/H ratios introduced by Bozorgnia and Campbell (2016), BC model, and Bommer et al. (2011), BAK model, revealed that both GMPMs mostly underestimate the V/H ratios at long spectral periods. Residual analysis showed that the capability of the BAK model for large source-to-site distances is more than that for small distances for which the ratios may be underestimated specially at long periods, while the estimations of the BC model do not significantly depend on the distance. Additionally, comparison of the observed V/H ratios with those computed based on the current 2018 edition of the Turkish seismic code indicated that the three considered earthquakes mostly yielded V/H ratios considerably higher than the code-based values in the long period range that could have higher damaging potential than those ratios in short period range.

Magnitude estimation and site characterization in southwestern British Columbia: application to earthquake early warning

In this study, we took a close look at the Ocean Networks Canada’s earthquake early warning system in southwestern British Columbia through analysis of the magnitude estimates by this system and characterization of site conditions for both onshore and offshore stations. Using magnitude values estimated at each station, over hundreds of notifications, we provided station terms to correct the magnitudes for stations that systematically generate high or low magnitude values. Moreover, by compiling a rich ground motion amplitude dataset and applying the horizontal-to-vertical spectral ratio method from Fourier amplitude of acceleration and response spectral acceleration, we investigated site characterization through evaluation of non-linear site response behavior and estimation of the site dominant frequency (fpeak) and its peak amplitude (Apeak) for each station. In general, no strong evidence of non-linearity is observed at any stations considering the magnitude-distance distribution of ground motions in this study. Offshore sites show fpeak and Apeak in the range of approximately 1.7–6 Hz and 0.4–1.2 (in base-10 log unit), respectively, whereas onshore sites show approximately 1–6 Hz for fpeak and 0.3–0.7 (in base-10 log unit) for Apeak.

Optimizing alteration mineral detection: A fusion of multispectral and hyperspectral remote sensing techniques in the Sar-e-Chah-e Shur, Iran

The Sar-e-Châh-e-Shur region in the Lut Block of eastern Iran has significant potential for diverse ore mineral deposits, particularly copper and iron mineralization. This study used multispectral remote sensing data, including ASTER, Sentinel 2, and Landsat 8 for detecting alteration zones associated with ore mineralization in the study area. Spectral ratios associated with porphyry copper-iron alterations and advanced argillic-iron alterations were utilized to establish training data, supplemented by hyperspectral data processing derived from Hyperion (spaceborne) and HyMap (airborne) sensors. Several classification algorithms were implemented, and evaluation metrics such as the Kappa coefficient and overall accuracy were calculated using a confusion matrix to assess classification performance. Field validation was conducted through thin mineral section analysis and geochemical stream sediment data, yielding a Normalized Score (NS) for conformity estimation. Notably, the Landsat 8, ASTER, and Sentinel 2 - classes category 2 -artificial neural networks-thresholds 2 (LAS-2-ANN2) class (representing advanced argillic alterations) was achieved a high score of 3.5, corresponding to an 87.5% match to the geological evidence. In addition, hyperspectral data processing using artificial neural network algorithms was yielded an overall score of 2.43, indicating a 61% match. In conclusion, this investigation highlighted the efficacy of fusing multispectral and hyperspectral data for accurate mineral alteration mapping in high potential regions, offering valuable insights for future exploration campaigns.

3D shear wave velocity imaging of the subsurface structure of granite rocks in the arid climate of Pan de Azúcar, Chile, revealed by Bayesian inversion of HVSR curves

Abstract. Seismic methods are emerging as efficient tools for imaging the subsurface to investigate the weathering zone. The structure of the weathering zone can be identified by differing shear wave velocities as various weathering processes will alter the properties of rocks. Currently, 3D subsurface modelling of the weathering zone is gaining increasing importance as results allow the identification of the weathering imprint in the subsurface not only from top to bottom but also in three dimensions. We investigated the 3D weathering structure of monzogranite bedrock near the Pan de Azúcar National Park (Atacama Desert, northern Chile), where the weathering is weak due to the arid climate conditions. We set up an array measurement that records seismic ambient noise, which we used to extract the horizontal-to-vertical spectral ratio (HVSR) curves. The curves were then used to invert for 1D shear wave velocity (Vs) models, which we then used to compile a pseudo-3D model of the subsurface structure in our study area. To invert the 1D Vs model, we applied a transdimensional hierarchical Bayesian inversion scheme, allowing us to invert the HVSR curve with minimal prior information. The resulting 3D model allowed us to image the granite gradient from the surface down to ca. 50 m depth and confirmed the presence of dikes of mafic composition intruding the granite. We identified three main zones of fractured granite, altered granite, and the granite bedrock in addition to the mafic dikes with relatively higher Vs. The fractured granite layer was identified with Vs of 1.4 km s−1 at 30–40 m depth, while the granite bedrock was delineated with Vs of 2.5 km s−1 and a depth range between 10 and 50 m depth. We compared the resulting subsurface structure to other sites in the Chilean coastal cordillera located in various climatic conditions and found that the weathering depth and structure at a given location depend on a complex interaction between surface processes such as precipitation rate, tectonic uplift and fracturing, and erosion. Moreover, these local geological features such as the intrusion of mafic dikes can create significant spatial variations to the weathering structure and therefore emphasize the importance of 3D imaging of the weathering structure. The imaged structure of the subsurface in Pan de Azúcar provides the unique opportunity to image the heterogeneities of a rock preconditioned for weathering but one that has never experienced extensive weathering given the absence of precipitation.

Validity of Acoustic Measures Obtained Using Various Recording Methods Including Smartphones With and Without Headset Microphones.

The goal of this study was to assess various recording methods, including combinations of high- versus low-cost microphones, recording interfaces, and smartphones in terms of their ability to produce commonly used time- and spectral-based voice measurements. Twenty-four vowel samples representing a diversity of voice quality deviations and severities from a wide age range of male and female speakers were played via a head-and-thorax model and recorded using a high-cost, research standard GRAS 40AF (GRAS Sound & Vibration) microphone and amplification system. Additional recordings were made using various combinations of headset microphones (AKG C555 L [AKG Acoustics GmbH], Shure SM35-XLR [Shure Incorporated], AVID AE-36 [AVID Products, Inc.]) and audio interfaces (Focusrite Scarlett 2i2 [Focusrite Audio Engineering Ltd.] and PC, Focusrite and smartphone, smartphone via a TRRS adapter), as well as smartphones direct (Apple iPhone 13 Pro, Google Pixel 6) using their built-in microphones. The effect of background noise from four different room conditions was also evaluated. Vowel samples were analyzed for measures of fundamental frequency, perturbation, cepstral peak prominence, and spectral tilt (low vs. high spectral ratio). Results show that a wide variety of recording methods, including smartphones with and without a low-cost headset microphone, can effectively track the wide range of acoustic characteristics in a diverse set of typical and disordered voice samples. Although significant differences in acoustic measures of voice may be observed, the presence of extremely strong correlations (rs > .90) with the recording standard implies a strong linear relationship between the results of different methods that may be used to predict and adjust any observed differences in measurement results. Because handheld smartphone distance and positioning may be highly variable when used in actual clinical recording situations, smartphone + a low-cost headset microphone is recommended as an affordable recording method that controls mouth-to-microphone distance and positioning and allows both hands to be available for manipulation of the smartphone device.

Site Response in the Walnut Creek–Concord Region of San Francisco Bay, California: Ground-Motion Amplification in a Fault-Bounded Basin

ABSTRACT Thirty-seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault-bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well-recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source-site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source-site spectral inversion method yields the best least-squares fit to site response for a database of ground-motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface-wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close-in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface-wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear-wave velocity using a global search algorithm yields VS30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite-element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin-edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

Maximization of cytochrome C oxidase enzyme activity by optimizing color conversion from red organic light-emitting diodes

Photobiomodulation (PBM) therapy using red and near-infrared (NIR) light can effectively activate cytochrome C oxidase (CCO) enzyme to promote wound healing in living organisms. To maximize this therapeutic effect, the use of a light source with a spectrum absorbed by the CCO enzyme is necessary. A red/NIR light source optimized for CCO enzyme was obtained by partially color-converting a red organic light-emitting diode with copper indium sulfide/zinc sulfide quantum dots (QDs). To increase color conversion efficiency, the dispersion of QDs was improved by using a fibrous color conversion layer (FCCL). By optimizing the QD content in the FCCL, a red/NIR light source similar to the absorption spectrum of CCO enzyme was obtained with a spectral matching ratio of 81.7 %. When this light source was used in PBM therapy, it was experimentally verified that it showed a much faster wound closure effect compared to general red light, and it was proven that this was due to increased activity of the CCO enzyme. This study is a method of obtaining a red/NIR light source for PBM suitable for a specific enzyme performing a specific function, and is particularly useful due to its excellent form factor of being large, thin, and flexible.

Identification of Ground Vulnerability Based on Kg, PGA, and GSS Assessment by Single Station Microtremor in Bendosari Hamlet, Kalasan Village, Yogyakarta

Bendosari Hamlet, Kalasan, Yogyakarta is a research area that is one of the areas affected by the earthquake that occurred on May 27, 2006. The earthquake with a magnitude of 6.4 caused quite severe damage to buildings around the research area. Based on the phenomenon of earthquakes that have occurred, there needs to be further study related to vulnerable areas in the event of another earthquake. Research data is primary data taken using a single statiun microtremor tool with an overall acquisition point of 50 measurement points. This study aims to find out the value of seismic vulnerability index (Kg), Peak Ground Acceleration (PGA), and Ground Shear Strain (GSS). Using the Horizontal to Vertical Spectral Ratio (HVSR) method, the dominant frequency result (f0) and amplification factor value (A0) from data that has been recorded through the microtremor tool. Based on research that has been done, the results are obtained in the f0rm of seismic vulnerability index values that are in the range of 0.49 - 30.18, the value of the acceleration of ground movement is in the range of 0.21 - 0.61 Gal, and the ground shear strain value is in the range of 4 x 10-6 to 38 x 10-5. Point K26 which is the original documentation of buildings affected by the earthquake, analyzed that the value of Vs of 858.112 m / s on the vs ground profile indicates that there is hard soil at a depth after 27.2 m. This shows that the thickness of subsurface sedimentary rocks at point K26 is quite thick. Keywords: Bendosari Hamlet; Mikrotremor; GSS; HVSR; PGA; Seismic Vulnerability

On-Site Earthquake Early Warning Model for Selected Records in the NGA-West2 Dataset Using S- and P-Wave Spectral Ratios

ABSTRACT On-site earthquake early warning (EEW) requires the best estimates of earthquake magnitude and distance parameters within a few seconds after the P-wave arrival for estimating the subsequent S-wave parameters. The errors in the estimates of the earthquake P-wave parameters will propagate into the estimates of the S-wave parameters. To solve this problem, we used the methodology by Zhao and Zhao (2019), which uses the spectral ratio R3P between the response spectral values of the first 3 s of the S wave and that of the first 3 s of the P wave, referred to as the R3P model. The modeling presented here was based on strong-motion records from the Next Generation Attenuation (NGA)-West2 dataset. We also used the spectral ratio RFP between the response spectral values from the full records from the S-wave arrival time to the end of the record and that of the first 3 s P wave to develop a second EEW model (the RFP model). An advantage of these two models is that the magnitude and hypocentral distance are not required in considerable magnitude and distance ranges. This means that the errors in the estimated source and path parameters from the first 3 s of the P wave will not affect the model predictions. A theoretical justification for these results is that the magnitude and the distance scaling rates for the first 3 s of the P wave are similar to those of the first 3 s of the S wave. This may also apply to the full S-wave window within useful EEW magnitude and distance ranges. We also found that when the estimated magnitude and distance for a record are necessary, the effect of the corresponding errors would be smaller than using a ground-motion prediction equation (GMPE), because the magnitude and distance scaling rates from this study are smaller than those of many GMPEs.

Inferring the Focal Depths of Small Earthquakes in Southern California Using Physics-Based Waveform Features

ABSTRACT Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ∼8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2<ML<4) recorded at local distances (<150 km), for which depth estimation is especially challenging. We find that differential magnitudes (Mw/ML–Mc) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡log10(M0)+3log10(fc), and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. We also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

High-resolution spatial elastic wave velocity variation detected using multiple-linear-array microtremor measurements

We propose a method for capturing the near-surface elastic wave velocity variation using microtremors. The proposed method does not directly estimate the velocity of the medium; instead, it measures the spatial velocity change. We apply the stretching method, which is often used to detect the time variation of shear wave velocity based on the monitoring of seismic ambient noise, to microtremor recordings. Because spatially dense ambient noise records are unavailable, we performed microtremor measurements using multiple linear arrays with multiple sensors along the survey line to collect as many recordings as possible. The synchronised microtremor recordings of each array were used to compute the horizontal-to-vertical spectral ratio (H/V), cross-correlation (CC), autocorrelation (AC), and single-station cross-correlation (SC). The H/V ratio and the waveforms of these correlation functions showed significant spatial variation. The velocity variation was measured quantitatively using the stretching method, which correlates a compressed and stretched reference function with each function on the line. The velocity variation of the vertical and radial components of CC compares well with the variation of the H/V peak, indicating that the radial components of CC retrieved the Rayleigh-wave Green’s function. The waveforms of AC and SC also showed spatial variation but the measured value of velocity change was smaller than that for CC. Further investigation is awaited to conclude if this is due to low-quality data or inappropriate measurement conditions. The reliability of the velocity change was confirmed by comparing the results to dispersion curves obtained using the conventional array method.