Spectral Density Measurements Research Articles

Estimation of Interferences in Magnetoencephalography (MEG) Brain Data Using Intelligent Methods for BCI-based Neurorehabilitation Applications

<p>Brain-Computer Interface (BCI) neurorehabilitation offers the potential to improve recovery and quality of life for stroke survivors. It aims to restore lost physical and mental abilities through motor and cognitive therapies. Magnetoencephalography (MEG) signals are a major advancement in BCI technology as they provide accurate and consistent assessments of brain activity for control and interaction applications. MEG is indispensable for recording the magnetic fields produced in the brain during motor imagery tasks due to its capability to evaluate cerebral activity with remarkable temporal resolution. However, one of the major challenges associated with MEG recording is the loss of signal quality due to physiological artifacts and ambient noise. Additionally, the head movement of the individual during the recording process can result in the introduction of artifacts into the recorded data, which can distort the spatial mapping of brain activity. This, in turn, can jeopardize the reliability and accuracy of the results obtained. This study aims to identify the most effective technique for removing artifacts from MEG signals by conducting a comparative performance analysis of prominent denoising algorithms, such as Infomax, FastICA, SOBI, and SWT. The findings conclude that Infomax is the most effective algorithm for removing physiological artifacts from a signal while maintaining the integrity and essential features of the original data. FastICA was found to be the second most effective algorithm. Infomax outperformed FastICA in Power Spectral Density (PSD) and Percentage Root mean square error Difference (PRD) measurements.</p>

Temperature dependence of the low-frequency noise in AlGaN/GaN fin field effect transistors

Low-frequency (LF) noise measurements are compared for Schottky-gate AlGaN/GaN heterostructure planar and fin field-effect transistors (FinFETs) as functions of gate voltage and measuring temperature. The noise of each device type is consistent with a carrier number fluctuation model. Similar effective defect-energy Eo distributions are derived for each of the two device architectures from measurements of excess drain-voltage noise-power spectral density vs temperature from 80 to 380 K. Defect- and/or impurity-related peaks are observed in the inferred energy distributions for Eo < 0.2 eV, Eo ≈ 0.45 eV, and Eo > 0.6 eV. Significant contributions to the LF noise are inferred for nitrogen vacancies and ON and FeGa impurity complexes. Ga dangling bonds at fin interfaces with gate metal are likely candidates for enhanced noise observed in FinFETs, relative to planar devices. Reducing the concentrations of these defects and impurity complexes should reduce the LF noise and enhance the performance, reliability, and radiation tolerance of GaN-based high electron mobility transistors.

Climatology of the open–closed boundary using ULF wave observations from South Pole, McMurdo, and distributed Antarctic AGOs

It has been shown that a proxy determination of the magnetospheric open–closed magnetic field line boundary (OCB) location can be made by examining the ultra-low-frequency (ULF) wave power in magnetometer data, with particular interest in the Pc5 ULF waves with periods of 3–10 min. In this study, we present a climatology of such Pc5 ULF waves using ground-based magnetometer data from the South Pole Station (SPA), McMurdo (MCM) station, and the Automatic Geophysical Observatories (AGOs) located across the Antarctic continent, to infer OCB behavior and variability during geomagnetically quiet times (i.e., Ap < 30 nT). For each season [i.e., austral fall (20 February 2017–20 April 2017), austral winter (20 May 2017–20 July 2017), austral spring (20 August 2017–20 October 2017), and austral summer (20 November 2017–20 January 2018)], north–south (i.e., H-component) magnetic field line residual power–spectral density (PSD) measurements taken during geomagnetically quiet periods within a 60-day window centered at the austral solstice/equinox are averaged in 10-min temporal bins to form the climatology at each station. These residual PSDs thus enable the analysis of Pc5 activity (and lower period “long-band” oscillations) and, thus, OCB location/variability as a function of season and magnetic latitude. The dawn and dusk transitions across the OCB are analyzed, with a discussion of dawn and dusk variability during nominally quiet geomagnetic periods. In addition, latitudinal dependencies of the OCB and peak Pc5 periods at each station are discussed, along with the empirical Tsyganenko model comparisons to our site measurements.

Graph-Based Electroencephalography Analysis in Tinnitus Therapy.

Tinnitus is the perception of sounds like ringing or buzzing in the ears without any external source, varying in intensity and potentially becoming chronic. This study aims to enhance the understanding and treatment of tinnitus by analyzing a dataset related to tinnitus therapy, focusing on electroencephalography (EEG) signals from patients undergoing treatment. The objectives of the study include applying various preprocessing techniques to ensure data quality, such as noise elimination and standardization of sampling rates, and extracting essential features from EEG signals, including power spectral density and statistical measures. The novelty of this research lies in its innovative approach to representing different channels of EEG signals as new graph network representations without losing any information. This transformation allows for the use of Graph Neural Networks (GNNs), specifically Graph Convolutional Networks (GCNs) combined with Long Short-Term Memory (LSTM) networks, to model intricate relationships and temporal dependencies within the EEG data. This method enables a comprehensive analysis of the complex interactions between EEG channels. The study reports an impressive accuracy rate of 99.41%, demonstrating the potential of this novel approach. By integrating graph representation and deep learning, this research introduces a new methodology for analyzing tinnitus therapy data, aiming to contribute to more effective treatment strategies for tinnitus sufferers.

Power spectral density analysis to investigate morphology of Ni films

Ni thin films grown by thermal evaporation and sputtering under different deposition conditions are characterised by structural and morphological properties using X-ray diffraction (XRD) and atomic force microscopy (AFM) techniques. XRD results suggested the growth of polycrystalline face-centred cubic Ni phase for all the samples. Morphological characteristics of the films were compared by analysing AFM data for root mean square roughness, height–height correlation function and power spectral density (PSD) measurements. Applying fractal and k-correlation fitting models to the PSD data, different morphological parameters are quantified. The study suggested that Ni films grown at higher substrate temperature (∼150 °C) by thermal evaporation and at low Ar pressure (∼0.2 Pa) by sputtering techniques yielded films of small surface roughness with a Brownian fractal self-affine surface.

Fractal Neurodynamics.

The neuronal ongoing electrical activity in the brain network, the neurodynamics, reflects the structure and functionality of generating neuronal pools. The activity of neurons due to their excitatory and inhibitory projections is associated with specific brain functions. Here, the purpose was to investigate if the local ongoing electrical activity exhibits its characteristic spectral and fractal features in wakefulness and sleep across and within subjects. Moreover, we aimed to show that measures typical of complex systems catch physiological features missed by linear spectral analyses. For this study, we concentrated on the evaluation of the power spectral density (PSD) and Higuchi fractal dimension (HFD) measures. Relevant clinical impact of the specific features of neurodynamics identification stands primarily in the potential of classifying cortical parcels according to their neurodynamics as well as enhancing the effectiveness of neuromodulation interventions to cure symptoms secondary to neuronal activity unbalances.

Laboratory measurements of thermal noise

An experimental study was undertaken of the aspect of underwater ambient noise which normally receives only a cursory reference—namely, thermal noise. A brief survey of the theory and literature, and consideration of how to separate amplifier noise from hydrophone noise, is followed by laboratory measurements on piston and spherical hydrophones. Results are presented of impedance and power spectral density measurements undertaken in water and other fluids. Measurements with the piston hydrophones were also performed as a function of distance from tank and water surface boundaries or from another hydrophone. The cross power spectral densities are presented as a function of separation distance between pairs of hydrophones. All hydrophones had their free field receive sensitivities and beam patterns measured, but the spectral densities have not yet been converted to pressures as the more interesting features are seen near resonances, where phase calibrations will also be required.

A broadband seismological observatory at Larsemann Hills, Antarctica: Noise characteristics and data quality

Antarctica, the southernmost ice-covered continent, plays a pivotal role in unraveling Earth's intricate evolutionary processes. Several studies have been conducted to comprehend the crustal structure of West Antarctica. However, our understanding of the sub-surface structure in the East Antarctica region, particularly along Princess Elizabeth Land (PEL), remains limited due to the scarcity of geophysical data. To bridge this knowledge gap, we (NCESS) established a permanent broadband seismological observatory in the captivating region of Larsemann Hills, East Antarctica. In addition, we investigate the noise characteristics at the station, employing power spectral density (PSD) measurements. Results have been promising, indicating that the noise levels are within the range of New High Noise Model and New Low Noise Model. Monthly variations in PSD reveal different noise levels throughout the year, with winter months exhibiting lower levels of microseismic noise attributed to lower temperatures and frozen sea surfaces. The long-period noise is higher in April & May, possibly due to high-speed winds and lower in September & October. The outcomes serve as a testament to the success of our installation and ensure a valuable data set. The data will contribute to a comprehensive investigation of sub-surface structures in the PEL region of East Antarctica, enhancing our understanding of geological processes and tectonic evolution. Furthermore, the findings will serve as a valuable resource for future research and contribute to our knowledge of Earth's dynamic processes.

Low-Profile Capacitive Load Cells for Underground Mining Material and Wear Classification to Promote Worker Safety

This work proposes a capacitive load cell for conical picks to enable underground continuous mining machine operators to perform their roles away from known hazardous regions near the machine. The load cell is embedded in commercially available flexible printed circuit board, integrates with the target tooling, and demonstrates in situ force sensing of vibration signatures for continuous mining cutting tools. Changes in material constitution, tool mass, and tool geometry cause modal variations in vibrational response measurable with force sensors at the cutting interface. Time-series measurements are captured during rock cutting tests using a linear cutting machine. These measurements are segmented into small windows, less than 0.25 s, and are preprocessed using the fast Fourier transform, which highlights the modal variations. The transformed measurements are then classified into different material and wear categories using support-vector machines with the radial basis function kernel. Different normalization schemes and Fourier transform methods are tested for signal preprocessing. Results show that the power spectral density measurements with normally distributed coefficients give good results for material classification, while the normalized time-domain measurements give better results for wear classification. Under laboratory conditions, this technique is shown to classify material and wear categories with F1 score above 0.85 out of 1.0 for our experiment. This technology could be used to assist operators in assessing material and wear conditions from a safer distance. It has applications in the coal mining industry as well as other applications which use conical picks such as road milling.

Direct Measurement of Acoustic Spectral Density and Fractional Topological Charge

LDOS is a fundamental spectral property that plays a central role in various physical phenomena, such as wave-matter interactions and spontaneous emissions. The role of LDOS in acoustics was uncovered only in recent years and the measurement of acoustic LDOS has not yet been achieved. Here, we report on the direct measurement of the LDOS in acoustic systems. The acoustic LDOS is quantified here through the measurement of the volume flow rate and the acoustic pressure with a local excitation-probe configuration. Based on this method, we detect the LDOS in the one-dimensional acoustic Su-Schrieffer-Heeger model and observe the fractional topological charge of the system. Our work unveils the role of the LDOS in acoustic phenomena and paves the way toward characterizing and tailoring the LDOS in acoustic systems, which holds promise for applications in acoustic emission control, energy localization, and wave-matter interaction.

A Pilot Investigation of Visual Pathways in Patients with Mild Traumatic Brain Injury.

In this study, we examined visual processing within primary visual areas (V1) in normal and visually impaired individuals who exhibit significant visual symptomology due to sports-related mild traumatic brain injury (mTBI). Five spatial frequency stimuli were applied to the right, left and both eyes in order to assess the visual processing of patients with sports-related mild traumatic brain injuries who exhibited visual abnormalities, i.e., photophobia, blurriness, etc., and controls. The measurement of the left/right eye and binocular integration was accomplished via the quantification of the spectral power and visual event-related potentials. The principal results have shown that the power spectral density (PSD) measurements display a distinct loss in the alpha band-width range, which corresponded to more instances of medium-sized receptive field loss. Medium-size receptive field loss may correspond to parvocellular (p-cell) processing deprecation. Our major conclusion provides a new measurement, using PSD analysis to assess mTBI conditions from primary V1 areas. The statistical analysis demonstrated significant differences between the mTBI and control cohort in the Visual Evoked Potentials (VEP) amplitude responses and PSD measurements. Additionally, the PSD measurements were able to assess the improvement in the mTBI primary visual areas over time through rehabilitation.

Analysis of seismic noise of broadband seismological stations installed along the Western Ghats

The Western Ghats (WG) is one of the great escarpments that extend ~1500 km parallel to the west coast of India in the direction of NNW-SSE. We deployed a network of seven broadband seismological stations along the WG to decipher its evolution. In the present study, we investigate the characteristics of different kinds of noises at the stations by utilizing the power spectral density (PSD) measurements. Further, the results are compared with the global standard noise models to assess the data quality. The PSD results reveal that the short period (0.1–1 s) cultural noise is more prominent at stations AGMB, SDPR and MGLI when compared to the other stations, especially during day hours of Indian Standard Time (IST) zone since these sites were in the proximity of roads. The seasonal variations are observed especially in the microseismic period range and noise levels are more prominent in the months of July to August since the western part of India experiences peak monsoon during this period. These variations are observed especially at PCH and KNUR stations as their locations were near the coastline. Furthermore, the results indicate that the noise levels are more prominent in vertical components than that of horizontal components in the microseismic period range whereas it is reversed in the short and long period ranges. Later, the results indicate that the noise levels at all the localities of stations are within the global standard noise models, which suggests that our installation of broadband seismological stations has been successful at all seven locations and has good data quality.

Optical receiver characterizations and corrections for ground-based and airborne measurements of spectral actinic flux densities

Abstract. Solar actinic radiation in the ultraviolet and visible range (UV/VIS) perpetuates atmospheric photochemistry by inducing photolysis processes which form reactive radical species. Photolysis frequencies are rate constants that quantify the rates of photolysis reactions and therefore constitute important parameters for quantitative analyses. Photolysis frequencies are usually calculated from modelled or measured solar spectral actinic flux densities. Suitable measurement techniques are available, but measurement accuracy can suffer from non-ideal 2π or 4π solid-angle reception characteristics of the usually employed 2π optical receivers or receiver combinations. These imperfections, i.e. deviations from an angle-independent response, should be compensated for by corrections of the measured data. In this work, the relative angular sensitivities of four commonly used 2π quartz receivers were determined in the laboratory in a range 280–660 nm. Based on this information, the influence of the non-ideal responses on measured spectral actinic flux densities for ground-based and airborne applications was investigated for a wide range of atmospheric conditions. Spectral radiance distributions and contributions of direct, diffuse downward and diffuse upward spectral actinic flux densities were calculated with a radiative transfer model to derive the corrections. The intention was to determine the ranges of possible corrections under realistic measurement conditions and to derive simple parametrizations with reasonable uncertainties. For ground-based 2π measurements of downward spectral actinic flux densities, corrections typically range <10 % dependent on wavelength and solar zenith angle, with 2 %–8 % uncertainties covering all atmospheric conditions. Corrections for 4π airborne measurements were determined for the platforms Zeppelin NT (New Technology) and HALO (High Altitude and Long Range Research Aircraft) in altitude ranges 0.05–2 and 0.2–15 km, respectively. Total, downward and upward spectral actinic flux densities were treated separately. In addition to various atmospheric conditions, different ground albedos and small (<5∘) aircraft attitude variations were considered in the uncertainties, as well as aircraft headings with respect to the sun in the case of HALO. Corrections for total and downward spectral actinic flux densities again typically range <10 % dependent on wavelength, solar zenith angle and altitude, with 2 %–10 % uncertainties covering all atmospheric conditions for solar zenith angles below 80∘. For upward spectral actinic flux densities, corrections were more variable and significantly greater, up to about −50 % at low altitudes and low ground albedos. A parametrization for corrections and uncertainties was derived using uncorrected ratios of upward / downward spectral actinic flux densities as input, applicable independent of atmospheric conditions for a given wavelength, solar zenith angle and altitude. The use was limited to conditions with solar zenith angles <80∘ when direct sun radiation cannot strike upward- and downward-looking receivers simultaneously. Examples of research flights with the Zeppelin and HALO are discussed, as well as other approaches described in the literature.

Galactic Cosmic-Ray Propagation in the Inner Heliosphere: Improved Force-field Model

A key goal of heliophysics is to understand how cosmic rays propagate in the solar system’s complex, dynamic environment. One observable is solar modulation, i.e., how the flux and spectrum of cosmic rays change as they propagate inward. We construct an improved force-field model, taking advantage of new measurements of magnetic power spectral density by Parker Solar Probe to predict solar modulation within the Earth’s orbit. We find that modulation of cosmic rays between the Earth and Sun is modest, at least at solar minimum and in the ecliptic plane. Our results agree much better with the limited data on cosmic-ray radial gradients within Earth’s orbit than past treatments of the force-field model. Our predictions can be tested with forthcoming direct cosmic-ray measurements in the inner heliosphere by Parker Solar Probe and Solar Orbiter. They are also important for interpreting the gamma-ray emission from the Sun due to scattering of cosmic rays with solar matter and photons.

Magnetic Field Dropouts and Associated Plasma Wave Emission near the Electron Plasma Frequency at Switchback Boundaries as Observed by the Parker Solar Probe

The first solar encounters by the Parker Solar Probe revealed the magnetic field to be dominated by short field reversals in the radial direction, referred to as “switchbacks.” While radial velocity and proton temperature were shown to increase inside the switchbacks, ∣B∣ exhibits very brief dropouts only at the switchback boundaries. Brief intensifications in spectral density measurements near the electron plasma frequency, f pe, were also observed at these boundaries, indicating the presence of plasma waves triggered by current systems in the form of electron beams. We perform a correlative study using observations from the Parker FIELDS Radio Frequency Spectrometer and Fluxgate Magnetometer to compare occurrences of spectral density intensifications at the electron plasma frequency (f pe emissions) and ∣B∣ dropouts at switchback boundaries during Parker’s first and second solar encounters. We find that only a small fraction of minor ∣B∣ dropouts are associated with f pe emissions. This fraction increases with ∣B∣ dropout size until all dropouts are associated with f pe emissions. Brief spikes in the differential electron flux measured by the SWEAP Solar Probe Analyzer for Electron sensors also occur in conjunction with nearly all f pe emissions. This suggests that in the presence of strong ∣B∣ dropouts, electron currents that create the perturbation in ∣B∣ along the boundaries are also stimulating plasma waves such as Langmuir waves.

Phase Noise Spectral Density Measurement of Broadband Frequency-Modulated Radar Signals

In continuous-wave (CW) radar systems, such as frequency-modulated (FMCW), frequency-stepped (FSCW), or orthogonal frequency-division multiplexing (OFDM) radar systems, the range and velocity uncertainty are significantly impaired by phase noise decorrelation. Therefore, radar designers require accurate knowledge of their synthesizers’ phase noise profiles to assess and predict radar performance. However, commercial phase noise analyzers cannot determine phase noise during modulation, and this may differ notably from phase noise in the pure CW mode. Recent methods for FMCW phase noise analysis usually require comprehensive <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">${a}$ </tex-math></inline-formula> <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">priori</i> knowledge of modulation parameter and are prone to systematic deviations. To overcome these issues, we propose a new approach based on differential analysis of subsequent time-domain measurements. This method retains, statistical phase noise information while reducing systematic influences. For the first time, less <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">${a}$ </tex-math></inline-formula> <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">priori</i> signal knowledge is required, and the method works for nearly any kind of broadband signal modulation. The concept requires only a digitizer (e.g., an oscilloscope) and some digital signal processing. The proposed method is first experimentally tested with different phase-locked-loop (PLL)-based synthesizer phase noise profiles. The obtained phase noise profiles agree perfectly with the results of an established measurement system. After this proof of basic functionality, the unique phase noise analysis capability for BB modulated signals is demonstrated with PLL-generated FMCW signals. The results reveal a significant phase noise difference between the different setups and clearly show the capability and benefit of the novel phase noise spectral density measurement concept.

Denoising EEG Signals for Real-World BCI Applications Using GANs.

As a measure of the brain's electrical activity, electroencephalography (EEG) is the primary signal of interest for brain-computer-interfaces (BCI). BCIs offer a communication pathway between a brain and an external device, translating thought into action with suitable processing. EEG data is the most common signal source for such technologies. However, artefacts induced in BCIs in the real-world context can severely degrade their performance relative to their in-laboratory performance. In most cases, the recorded signals are so heavily corrupted by noise that they are unusable and restrict BCI's broader applicability. To realise the use of portable BCIs capable of high-quality performance in a real-world setting, we use Generative Adversarial Networks (GANs) that can adopt both supervised and unsupervised learning approaches. Although our approach is supervised, the same model can be used for unsupervised tasks such as data augmentation/imputation in the low resource setting. Exploiting recent advancements in Generative Adversarial Networks (GAN), we construct a pipeline capable of denoising artefacts from EEG time series data. In the case of denoising data, it maps noisy EEG signals to clean EEG signals, given the nature of the respective artefact. We demonstrate the capability of our network on a toy dataset and a benchmark EEG dataset developed explicitly for deep learning denoising techniques. Our datasets consist of an artificially added mains noise (50/60 Hz) artefact dataset and an open-source EEG benchmark dataset with two artificially added artefacts. Artificially inducing myogenic and ocular artefacts for the benchmark dataset allows us to present qualitative and quantitative evidence of the GANs denoising capabilities and rank it among the current gold standard deep learning EEG denoising techniques. We show the power spectral density (PSD), signal-to-noise ratio (SNR), and other classical time series similarity measures for quantitative metrics and compare our model to those previously used in the literature. To our knowledge, this framework is the first example of a GAN capable of EEG artefact removal and generalisable to more than one artefact type. Our model has provided a competitive performance in advancing the state-of-the-art deep learning EEG denoising techniques. Furthermore, given the integration of AI into wearable technology, our method would allow for portable EEG devices with less noisy and more stable brain signals.

Experimental Assessment of a Magnetic Induction-Based Receiver for Magnetic Communication

This paper presents the topology of a novel approach to magnetic communication using a differential magnetic induction (MI)-based receiver and a differential MI receiving sensor. In this paper, a differential MI sensor based on two ferromagnetic cores is proposed as a receiving sensor, unlike the air coil-type MI sensor in the conventional search coil sensor concept. This differential MI sensor has the advantages of ultra-high sensitivity characteristics of the pT/√Hz level, which can detect weak magnetic fields in magnetic communication; moreover, the sensor is smaller than a conventional air coil MI sensor. The proposed differential MI sensor contributes to improving sensor performance by increasing its signal-to-noise ratio. The design and fabrication of the proposed MI sensor were based on a printed circuit board (PCB). The pickup coil of the PCB-based MI sensor directly wound the pickup coil onto a ferromagnetic core composed of Ni-Zn ferrite material. To analyze the key factors that affected the performance of the receiver, the magnetic field-to-voltage conversion ratio (MVCR) and equivalent magnetic spectral density measurements of the proposed PCB-based MI sensor were performed. Wireless digital communication using quadrature phase shift keying (QPSK), which is less sensitive to noise and has a high data rate, was used to evaluate the proposed MI-based receiver. The transmitted and received waveforms were compared to confirm that the transmitted digital data were accurately received as a result of the final demodulation of the receiver. Additionally, several performance metrics, such as constellation and error vector magnitude, were measured. The results of the comprehensive analysis confirmed the applicability of the proposed differential MI-based receiver to a magnetic field.

The Wavefront Power Spectral Density Measurement of Aspheric Lens with Long Focal-Length Using a Computer-Generated Hologram

To ensure the performance of the optical system, the machining accuracy of lens with long focal lengths is required to ensure the image quality. A new method for lens transmission wavefront power spectral density (PSD) in mid-frequency domain measurement using binary phase computer-generated hologram (CGH) is presented. This technique is widely applicable and is particularly useful for measuring large-size lenses with long focal lengths. A comparison experiment of the CGH measurement with results from a Fizeau sphere interferometry method is carried out to verify the accuracy and convenience of the measurement. Furthermore, measurement uncertainty due to CGH fabrication process is analysed. Analysis of the CGH test showed the overall accuracy of less than 1 nm RMS for a sphere lens with over 30 m focal length and Φ410 mm clear aperture. CGH can provide reference spheres with high precision, in the meantime greatly shorten air space, thus reducing the effect of vibration and air turbulence, therefore is of great importance for lens transmission wavefront PSD measurement. The realization of high precision, high efficiency and nondestructive testing of long focal-lens wavefront PSD ensure the ultra-precision and certainty level of machining, hence improving the comprehensive performance of the optical system.

Gisola: A High-Performance Computing Application for Real-Time Moment Tensor Inversion

Abstract Automatic moment tensor (MT) determination is essential for real-time seismological applications. In this article, Gisola, a highly evolved software for MT determination, oriented toward high-performance computing, is presented. The program employs enhanced algorithms for waveform data selection via quality metrics, such as signal-to-noise ratio, waveform clipping, data and metadata inconsistency, long-period disturbances, and station evaluation based on power spectral density measurements in parallel execution. The inversion code, derived from ISOLated Asperities—an extensively used manual MT retrieval utility—has been improved by exploiting the performance efficiency of multiprocessing on the CPU and GPU. Gisola offers the ability for a 4D spatiotemporal adjustable MT grid search and multiple data resources interconnection to the International Federation of Digital Seismograph Networks Web Services (FDSNWS), the SeedLink protocol, and the SeisComP Data Structure standard. The new software publishes its results in various formats such as QuakeML and SC3ML, includes a website suite for MT solutions review, an e-mail notification system, and an integrated FDSNWS-event for MT solutions distribution. Moreover, it supports the ability to apply user-defined scripts, such as dispatching the MT solution to SeisComP. The operator has full control of all calculation aspects with an extensive and adjustable configuration. MT’s quality performance, for 531 manual MT solutions in Greece between 2012 and 2021, was measured and proved to be highly efficient.