Source Time Function Research Articles

Imaging the Northern Los Angeles Basins with Autocorrelations

Abstract We show reflectivity cross sections for the San Gabriel, Chino, and San Bernardino basins north of Los Angeles (LA), California, determined from autocorrelations of ambient noise and teleseismic earthquake waves. These basins are thought to channel the seismic energy from earthquakes on the San Andreas fault to LA, and a more accurate model of their depth is important for hazard mitigation. We use the causal side of the autocorrelation function (ACF) to determine the zero-offset reflection response. To minimize the smoothing effect of the source time function, we remove the common mode from the autocorrelation to reveal the zero-offset reflection response. We apply this to 10 temporary nodal lines consisting of a total of 758 geophones with an intraline spacing of 250–300 m. We also show that the ACF from teleseismic events can provide illumination on the subsurface that is consistent with ambient noise. Both autocorrelation results compare favorably to receiver functions.

Two unequal shots in the same hole: statistical versus deterministic deconvolution

Determination of the source time function and earth impulse response for explosive source seismic reflection data is traditionally performed using predictive deconvolution, which assumes the reflection seismogram is the convolution of a minimum-phase wavelet with a white, random, stationary sequence of impulses. We compare this statistical method with a new deterministic method that uses two unequal shots close together in the same hole and assumes a seismogram is the convolution of a source time function with a Green's function, plus noise. The seismograms from the two shots share essentially the same Green's function to any receiver, because the shots are close together. The Green's function is estimated by deconvolving the two seismograms at a given receiver for the corresponding source time functions, modelled independently; source parameters are obtained by trial-and-error by minimising the difference between the two estimates of the Green's function. We test these competing theories using seismic reflection data acquired with two charges fired separately in the same hole into a receiver spread of 1178 bunched geophone groups at 10 m intervals: a 1 kg charge at 42 m depth and 2 kg at 40 m depth. The criterion for comparison of the theories is the similarity of the recovered impulse responses at any receiver. Each method recovers a wavelet, or source time function, for each shot and two estimates of the impulse response, or Green's function, at each receiver. The Green's functions recovered from the two shots using the deterministic method are almost identical and are measurably more similar than the impulse responses recovered using predictive deconvolution. That is, the deterministic method is better. The resolution of explosive source seismic reflection data can now be significantly improved with the new method and minimal additional field effort of loading one extra charge per shot hole.

Rupture process of the 2019 Mw 6.8 Davao Del Sur earthquake, Philippines from geodetic and seismic data

Because of its complex tectonic setting, the Philippines experiences vigorous seismic activity. On December 15, 2019, an Mw 6.8 earthquake struck Davao del Sur, the Philippines, severely damaging infrastructure and adversely affecting the population. In this study, we used Interferometry Synthetic Aperture Radar (InSAR) and teleseismic broadband data to explore the earthquake’s source parameters, rupture process, and tectonic characteristics through a Bayesian modeling approach and Coulomb stress analysis. An analysis of InSAR suggested that the mainshock occurred on an NW–SE oriented single plane within the Cotabato fault system, with the Makilala-Malungon fault being the causative fault. The Joint inversion revealed that the rupture initiated in the middle crust (depth: approximately 16.5 km), exhibiting a blind, left-lateral strike-slip fault rupture propagation with minimal reverse slip and no surface rupture. The slip was concentrated in a distinct asperity spanning a depth of 3–24 km, with the peak slip reaching approximately 1.8 m. The source time function of this event indicated a single pulse with a total source duration of approximately 17 s. The static stress drop value was low, which may be related to the brittle condition of the crust. Coulomb stress changes calculations suggested earthquake risk in the future.

Isotropic High‐Frequency Radiation in Near‐Fault Seismic Data

AbstractWe compare Fourier Amplitude Spectra of Fault Normal (FN) and Fault Parallel (FP) seismograms at near‐fault sites for seven strike‐slip earthquakes with moment magnitudes Mw ≥ 6. For all events we find large FN/FP ratios at low frequencies consistent with near‐fault S‐wave radiation patterns for strike‐slip earthquakes. However, the difference diminishes with increasing frequency and FN/FP is about 1 above a transition frequency. The results may reflect small tensile/isotropic components in the earthquake rupture zones that homogenize the high‐frequency radiation in different directions at near‐fault sites. The FN/FP ratios at low frequencies and transition frequencies above which FN ∼ FP vary among the analyzed earthquakes and have no clear correlation with the magnitudes. The lack of correlation may signify a characteristic scale (e.g., process zone size, duration of source time function) controlling the isotropic radiation, and/or wave propagation and other effects that mask the source effects.

The Impact of Source Time Function Complexity on Stress-Drop Estimates

ABSTRACT Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time- and frequency-domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016) source time function (STF) database and analyze how rupture complexity drives stress-drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress-drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved.

Physical Characteristics Controlling Radiation from Heterogeneous Ruptures—Finite Faults

ABSTRACT Kinematic simulations of ground motion require representations of the earthquake source: the distribution of final slip, parameters of the source time function, and the velocity of rupture travel. There is a significant ambiguity in prescribing these physical characteristics, causing uncertainty in the resulting motions that needs to be quantified. The representation integral is an appropriate tool: it allows exact calculation of the source effect in both the near and far fields in the frequency band of practical interest. The commonly used distributions of slip have a k-square shape of their wavenumber spectra. Various k-square slips change the slope of the radiated spectra in the range of ∼−2.5 and −4.0 in both the far and near fields. The spectra generated by randomly disturbed constant slip are indistinguishable from those emitted by k-square faults. In both cases, variations in peak values of ground velocity and acceleration between realizations are relatively insignificant: under ∼15% for the same hypocenter position. The slopes of the Fourier spectra produced exclusively by the form of the slip function and the slip heterogeneity are equivalent to using a formal kappa filter with κ ranging from ∼0.025 to 0.045 s. No ad hoc high-frequency filtering (of kappa or fmax type) is required if fault finiteness is accounted for. Geometric irregularity of rupture fronts, at least for the way the front progression is randomized in our case, does not appreciably affect the slopes of the spectra. Its principal effect is in blurring the directivity, reducing the sharpness of radiated pulses. The most influential parameter affecting the peak ground motions for several commonly used slip functions is the maximum velocity of slip: scaling of vmax causes a proportional scaling in peak ground acceleration. This parameter is the most important to constrain to reduce ambiguities in predicted ground motions.

Constraints on the 2013 Saravan intraslab earthquake using permanent GNSS, InSAR and seismic data

SUMMARY On 16 April 2013, an Mw = 7.7 earthquake struck the border of Iran and Pakistan in the central part of the Makran subduction zone with a reported depth of 80 km by USGS. This rare event in this poorly instrumented region helps to shed light on the kinematics of the subducting slab. We investigate source parameters of the Saravan intraslab normal earthquake using RADARSAT-2 SAR images in three ascending tracks, nine permanent GNSS sites and teleseismic data. The maximum coseismic displacement occurred at the SRVN GNSS station with 54.1 mm southeast horizontal and 42.7 mm upward vertical displacements. The coseismic ascending InSAR displacement maps illustrate a continuous and smooth NE-trending elliptical shape deformation pattern with a maximum of ∼29 cm of displacement away from the satellite. We use 25 broad-band teleseismic P-waveforms to estimate the focal mechanism of the main shock. A joint uniform inversion of InSAR, GNSS and teleseismic data reveals a NW-dipping SW-striking fault and a primarily normal-faulting earthquake with a minor right-lateral strike-slip component. The static slip distribution of the InSAR coseismic maps localizes variable slip at depths between 50 and 81 km with a maximum amplitude >3 m at 60–75.5 km depth, rupturing the oceanic crust of the subducted slab. The kinematic slip distribution exhibits a well-constrained slip pattern with a nucleation depth of 65 km. The source time function indicates that the earthquake reaches its maximum moment tensor release at ∼8 and ∼16 s. The NE-trend of the Saravan earthquake slip pattern, the orientation of the volcanic arc, and the distribution of the intraslab intermediate-depth normal earthquakes provide new insights into slab geometry in the central Makran subduction zone. We suggest that the slab bending at the hinge of subducting Arabian Plate is oblique along a NE–SW direction parallel to the volcanic arc rather than the shoreline or deformation front, and it is likely to be the reason for an oblique volcanic arc in the Makran subduction zone. These new constraints on the Makran slab geometry will help further studies in establishing realistic coupling maps for seismic hazard assessment.

Comparison between rupture parameters of intermediate and deep earthquakes at the Peru–Brazil–Bolivia border and northern Chile

We determined the main parameters of the source rupture process of intermediate- and deep-depth earthquakes occurring in the Peru–Brazil–Bolivia border region and northern Chile. The parameters of depth, fault-plane orientation, scalar seismic moment, slip distribution, and radiated seismic energy are obtained from seismograms. We selected 15 intermediate-depth earthquakes (100 km < h < 300 km) and 10 very deep earthquakes (h > 500 km) with magnitudes MW ≥ 6.0. For most events, the slip distribution over the rupture plane shows a single asperity, and the source time function presents a simple pulse. There are differences between intermediate-depth and deep earthquakes. The rupture areas, maximum slip and source time function (STF) duration are larger for intermediate-depth events than for deep events. Additionally, the STF’s show a sharper increase for deep earthquakes. The scaled radiated seismic energy shows larger values for deep depth events. The stress regime pattern derived from the obtained focal mechanism agrees with the geometry of the subduction of the Nazca plate. At intermediate depths, in the northern area up to 12°S, the stress pattern corresponds to a horizontal extension, while in the southern area, the tension axes dip at an angle of 30°. At deep depths, the stress regime corresponds to vertical compression in the north and dips of approximately 45° in the south.

Rupture Model of the 5 April 2024 Tewksbury, New Jersey, Earthquake Based on Regional Lg-Wave Data

Abstract On 5 April 2024, an earthquake of magnitude 4.8 occurred in Tewksbury, New Jersey. It was the largest instrumentally recorded event since 1900 in New Jersey and southern New York. Millions of people around New York City, ∼65 km east-northeast of Tewksbury, felt the shaking from the mainshock, but the epicentral area experienced no known significant property damages. We determine the focal mechanism, which is oblique faulting, and retrieve the Lg-wave relative source time functions (RSTFs) from the stations at regional distances to understand rupture processes and ground motions. Our fault-slip models well explain azimuthal variations of the RSTFs. The models show the rupture propagating toward the east-northeast (∼50° to 60°), not along the fault strike. The slip distribution on the nodal plane striking north and dipping to the east shows a slip area of 1.1 km radius with the rupture propagating down-dip. The down-dip rupture may account for the observed lack of strong shaking in the epicentral area.

Forward modeling of quake’s infrasound recorded in the stratosphere on board balloon platforms

Acoustic waves generated by seismic waves contain information on the internal structure of planets, and can be sensed by pressure sensors onboard high-altitude balloons. To identify the various contributions (infrasound signal, noise, balloon response, etc.) in such pressure records, a full waveform modeling is implemented and completed by infrasound ray tracing and additional data analysis. Here, we analyze the Stratéole-2 pressure data associated with two earthquakes (Garcia et al. Geophys. Res. Lett. 49(15):e98844, 2022) and compared these to full waveform simulations by SPECFEM2D-DG-LNS software. Even if our simulations do not precisely reproduce the waveform observed in the frequency range [0.05, 0.3] Hz, we show that the waveform presents more sensitivity to quake and internal structure parameters than to atmospheric structure, and that seismic surface wave dispersion is observed in balloon pressure records. The long-duration pressure oscillations observed after the main infrasonic signal cannot be fully reproduced by our one-dimensional input model even when source time function complexity and aftershocks are considered. These features are ascribed mainly to the complex vertical ground movements below the balloon and partly to late secondary infrasound arrivals excited by the interactions of seismic waves with the topography. These results enhance the advantages and limitations of quake-related infrasound observations on board terrestrial and planetary balloon platforms.Graphical

Mechanisms and Seismological Signatures of Rupture Complexity Induced by Fault Damage Zones in Fully‐Dynamic Earthquake Cycle Models

AbstractDamage zones are common around faults, but their effects on earthquake mechanics are still incompletely understood. Here, we investigate how damage affects rupture patterns, source time functions (STF) and ground motions in 2D fully‐dynamic cycle models. We find that back‐propagating rupture fronts emerge in large faults and can be triggered by residual stresses left by previous ruptures or by damage‐induced pulse‐to‐crack transitions. Damage‐induced back‐propagating fronts are modulated by slip rate oscillations, amplify high‐frequency radiation, and sharpen the multiple peaks in STF even in the absence of frictional heterogeneity or fault segmentation. Near‐field ground motion is predominantly controlled by stress heterogeneity left by prior seismicity, and further amplified within the damage zone by trapped waves and outside it by secondary rupture fronts. This study refines our knowledge on damage zone effects on earthquake rupture and identifies their potentially observable signatures in the near and far field.

Three‐Dimensional Teleseismic Elastic Reverse‐Time Migration With Deconvolution Imaging Condition and Its Application to Southwest Japan

AbstractWe have developed a novel deconvolution‐based reverse time migration method to image lithospheric structures using teleseismic data recorded by dense seismic arrays. Unlike traditional approaches that rely on the retrieval of Green's functions (or receiver functions), the new method directly utilizes the recorded three‐component (3‐C) seismic waveforms to reconstruct subsurface wavefields, which has the advantage of avoiding the estimation and removal of source time function. Importantly, it enables the asymptotic estimation of P‐to‐S transmission conversion coefficients at the elastic discontinuities and enhances resolving power for the fine‐scale heterogeneities. Taking these improvements, we have obtained a full three‐dimensional (3‐D) high‐resolution image of the subduction zone beneath southwest Japan. This image provides a comprehensive configuration of the severely deformed subducting plate beneath the Kii Channel and Kinki Peninsula.

Empirical Green’s function analysis of some induced earthquake pairs from the Groningen gas field

We have applied the empirical Green’s function (EGF) method to 53 pairs of earthquakes, with magnitudes ranging from M = 0.4 to M = 3.4, induced by gas production from the Groningen field in the Netherlands. For a subset of the events processed, we find that the relative source time functions obtained by the EGF deconvolution show clear indications of a horizontal component of rupture propagation. The earthquake monitoring network used has dense azimuthal coverage for nearly all events such that wavelet duration times can be picked as a function of source-station azimuth and inverted using the usual Doppler broadening model to estimate rupture propagation strike, distance, and velocity. Average slip velocities have also been estimated and found to be in agreement with typical published values. We have used synthetic data, from both a simple convolutional model of the seismogram and more sophisticated finite difference rupture simulations, to validate our data processing workflow and develop kinematic models which can explain the observed characteristics of the field data. Using a measure based on the L1-norm to discriminate results of differing quality, we find that the highest quality results show very good alignment of the rupture propagation with directions of the detailed fault map, obtained from the full-field 3D seismic data. The dip direction rupture extents were estimated from the horizontal rupture propagation distances and catalogue magnitudes showing that, for all but the largest magnitude event (the M = 3.4 event of 8th January 2018), the dip-direction extent is sufficiently small to be contained wholly within the reservoir.

Bayesian Inversion of MDJ Seismograms for Source Parameters and Yields of DPRK Underground Nuclear Explosions 2 to 6

AbstractThe Democratic People’s Republic of Korea (DPRK) conducted six underground nuclear explosions at the Punggye-ri nuclear test site at Mount Mantap, a granite peak. Test 1 was separate from tests 2 to 6, which were within about 1 km of each other. Using seismograms recorded at Mudanjiang (MDJ) seismic station in China, I propose a new approach to obtain source parameters, source time functions and yields of events 2 to 6, assuming they share the same Green’s function from Punggye-ri to MDJ. Each source is modelled as a spherical cavity in a homogeneous isotropic elastic full space, with four independent parameters constrained by published data on the properties of granite and analysis of the recorded MDJ seismograms. The effect of the ground surface is included as a planar reflection that modifies the pressure at the cavity boundary. The Green’s function for each event is estimated by deconvolving the seismogram for the estimated source time function. Very fast simulated annealing (VFSA) is used to search the parameter space to minimise the root-mean square difference among the estimated Green’s functions and their mean. The estimated Green’s functions are similar and differ in amplitude by less than a factor of 2. Green’s functions from Punggye-ri to other seismic stations may be obtained by deconvolving the seismograms for the corresponding source time functions. Mapping the nonlinear zones surrounding each explosion indicates that this part of the site was destroyed by the underground nuclear tests before its official destruction on 24 May 2018.

Experimental full waveform inversion for elastic material characterization with accurate transducer modeling

For quality assurance and in-service safety, sufficient characterization of the materials is essential during the development, manufacture, and processing of a material in any industrial setting. Evaluation of elastic coefficients, material microstructures, morphological features, and related mechanical properties are all part of the process of characterizing a material. Ultrasonic signals are sensitive to material properties such as wave speeds, attenuation, diffusion backscattering, microstructural variation, and elastic characteristics (e.g., elastic modulus, hardness, etc.). Ultrasonic computed tomography (USCT) is an emerging imaging method that can be implemented to obtain a quantitative estimation of material properties. In this study, a source estimation technique was initially proposed to obtain the source time function for accurate forward modeling by constructing a linear inverse problem for the unknown transducer modeling. Finally, a material characterization approach was proposed with accurate source estimation to extract wave speed distribution from an elastic material by employing a wave-based method, known as full waveform inversion (FWI). Systematic performance analysis of the proposed FWI model with accurate source estimation was assessed using experimental and synthetic data obtained from a 6061 aluminum sample. Overall, the proposed FWI technique has successfully reconstructed the wave speed distribution, exhibiting the potential of the proposed method of material characterization in various engineering applications.

The 2020 Mw 7.0 Samos (Eastern Aegean Sea) Earthquake: joint source inversion of multitype data, and tsunami modelling

SUMMARY We present a kinematic slip model and a simulation of the ensuing tsunami for the 2020 Mw 7.0 Néon Karlovásion (Samos, Eastern Aegean Sea) earthquake, generated from a joint inversion of high-rate GNSS, strong ground motion and InSAR data. From the inversion, we find that the source time function has a total duration of ∼20 s with three peaks at ∼4, 7.5 and 15 s corresponding to the development of three asperities. Most of the slip occurs at the west of the hypocentre and ends at the northwest downdip edge. The peak slip is ∼3.3 m, and the inverted rake angles indicate predominantly normal faulting motion. Compared with previous studies, these slip patterns have essentially similar asperity location, rupture dimension and anticorrelation with aftershocks. Consistent with our study, most published papers show the source duration of ∼20 s with three episodes of increased moment releases. For the ensuing tsunami, the eight available gauge records indicate that the tsunami waves last ∼18–30 hr depending on location, and the response period of tsunami is ∼10–35 min. The initial waves in the observed records and synthetic simulations show good agreement, which indirectly validates the performance of the inverted slip model. However, the synthetic waveforms struggle to generate long-duration tsunami behaviour in simulations. Our tests suggest that the resolution of the bathymetry may be a potential factor affecting the simulated tsunami duration and amplitude. It should be noted that the maximum wave height in the records may occur after the decay of synthetic wave amplitudes. This implies that the inability to model long-duration tsunamis could result in underestimation in future tsunami hazard assessments.

Detectability analysis of very low frequency earthquakes: methods and application in Nankai using F-net and DONET broad-band seismometers

SUMMARY For a more quantitative discussion of slow earthquake activity, we evaluated the detectable limits of very low frequency earthquakes (VLFEs), which are seismic slow earthquakes observed in very low-frequency (&amp;lt; 0.05 Hz) bands in the Nankai subduction zone. We performed numerical simulations using a local 3-D model and used the observed noise level of permanent broad-band seismometers. First, we investigated the effects of the source-time functions on the maximum amplitudes of the VLFE signals at a certain station. The maximum amplitudes of the VLFE signals were controlled by the VLFE moment rate. The detectable limit of VLFEs at each source location can be defined as the lowest moment rate of detectable VLFEs, which radiate signals larger than the noise levels of any component at ≥ 3 stations. For inland seismometers only, the detectable limits of VLFEs at deep (30–40 km) and shallow (≤ 10 km) depths were 1012–1012.3 and 1012.7 N·m s−1, respectively. Due to the geometrical spreading of VLFE signals and large noise levels in horizontal components, offshore seismometers improved the detectability of shallow VLFEs in regions where seismometers were densely deployed. Based on our detectability and published catalogues, shallow slow earthquakes are less active south-southwest off the Kii Peninsula, where geodetic studies expect mechanical coupling.

Rapid Estimation of Single-Station Earthquake Magnitudes with Machine Learning on a Global Scale

ABSTRACT The foundation of earthquake monitoring is the ability to rapidly detect, locate, and estimate the size of seismic sources. Earthquake magnitudes are particularly difficult to rapidly characterize because magnitude types are only applicable to specific magnitude ranges, and location errors propagate to substantial magnitude errors. We developed a method for rapid estimation of single-station earthquake magnitudes using raw three-component P waveforms observed at local to teleseismic distances, independent of prior size or location information. We used the MagNet regression model architecture (Mousavi and Beroza, 2020b), which combines convolutional and recurrent neural networks. We trained our model using ∼2.4 million P-phase arrivals labeled by the authoritative magnitude assigned by the U.S. Geological Survey. We tested input data parameters (e.g., window length) that could affect the performance of our model in near-real-time monitoring applications. At the longest waveform window length of 114 s, our model (Artificial Intelligence Magnitude [AIMag]) is accurate (median estimated magnitude within ±0.5 magnitude units from catalog magnitude) between M 2.3 and 7.6. However, magnitudes above M ∼7 are more underestimated as true magnitude increases. As the windows are shortened down to 1 s, the point at which higher magnitudes begin to be underestimated moves toward lower magnitudes, and the degree of underestimation increases. The over and underestimation of magnitudes for the smallest and largest earthquakes, respectively, are potentially related to the limited number of events in these ranges within the training data, as well as magnitude saturation effects related to not capturing the full source time function of large earthquakes. Importantly, AIMag can determine earthquake magnitudes with individual stations’ waveforms without instrument response correction or knowledge of an earthquake’s source-station distance. This work may enable monitoring agencies to more rapidly recognize large, potentially tsunamigenic global earthquakes from few stations, allowing for faster event processing and reporting. This is critical for timely warnings for seismic-related hazards.

Sensing Optical Fibers for Earthquake Source Characterization Using Raw DAS Records

AbstractDistributed Acoustic Sensing (DAS) is becoming a powerful tool for earthquake monitoring, providing continuous strain‐rate records of seismic events along fiber optic cables. However, the use of standard seismological techniques for earthquake source characterization requires the conversion of data in ground motion quantities. In this study we provide a new formulation for far‐field strain radiation emitted by a seismic rupture, which allows to directly analyze DAS data in their native physical quantity. This formulation naturally accounts for the complex directional sensitivity of the fiber to body waves and to the shallow layering beneath the cable. In this domain, we show that the spectral amplitude of the strain integral is related to the Fourier transform of the source time function, and its modeling allows to determine the source parameters. We demonstrate the validity of the technique on two case‐studies, where source parameters are consistent with estimates from standard seismic instruments in magnitude range 2.0–4.3. When analyzing events from a 1‐month DAS survey in Chile, moment‐corner frequency distribution shows scale invariant stress drop estimates, with an average of Δσ = (0.8 ± 0.6) MPa. Analysis of DAS data acquired in the Southern Apennines shows a dominance of the local attenuation that masks the effective corner frequency of the events. After estimating the local attenuation coefficient, we were able to retrieve the corner frequencies for the largest magnitude events in the catalog. Overall, this approach shows the capability of DAS technology to depict the characteristic scales of seismic sources and the released moment.

Sensitivity of the second seismic moments resolution to determine the fault parameters of moderate earthquakes

Second-degree seismic moments provide a simple description of the spatiotemporal extent of the earthquake source. Finite source attributes such as rupture length, width, duration, velocity, and propagation direction can be estimated by computing second-degree seismic moments without the need for a predefined rupture model. This is achieved by analyzing the properties of apparent source time functions (ASTFs) obtained from seismic signals recorded at different stations after eliminating instrument responses and path effects. In this study, to define the limits of its application in the analysis of small earthquakes and to evaluate the sensitivity and reliability of the results to uncertainties due to observations and prior knowledge, we modeled a synthetic seismic source and examined how potential uncertainties in hypocentral depth, velocity model, focal mechanism, source duration, and number of recording stations can affect the inversion results. An accurate ASTF is essential to obtain robust results and our findings show that the mean values of the key source parameters, i.e., fracture size, source duration, and rupture velocity, are generally well reproduced in all sensitivity tests, with some exceptions, within the standard deviation. We also demonstrate that large uncertainties in the hypocentral depth and inaccurate velocity models introduce a significant bias, especially in rupture size and average centroid velocity, indicating the strong influence of ray path calculation in the inversion process. These resolution limits must therefore be taken into account when interpreting the results obtained with this technique.