Figure 1.
Comparison of the manufacturer-specified self-noise power spectral density (PSD) of representative accelerometer technologies used in seismic monitoring and Structural Health Monitoring (SHM) applications. The comparison includes conventional MEMS accelerometers (ADXL355, VS1002, and SI1003), QMEMS accelerometers (Epson M-A352 and M-A370), and force-balance accelerometers (Nanometrics Titan, Güralp Fortimus, Kinemetrics EpiSensor ES-T, and Lunitek Triton). The New High Noise Model (NHNM) and New Low Noise Model (NLNM) [
35], together with reference spectra (thick black curves) for earthquakes of increasing magnitude (M1.5–M5.5) at an epicentral distance of 10 km (dashed grey curves), are included for comparison. The PSD curves were derived from the self-noise specifications reported in the corresponding manufacturers’ datasheets.
Figure 1.
Comparison of the manufacturer-specified self-noise power spectral density (PSD) of representative accelerometer technologies used in seismic monitoring and Structural Health Monitoring (SHM) applications. The comparison includes conventional MEMS accelerometers (ADXL355, VS1002, and SI1003), QMEMS accelerometers (Epson M-A352 and M-A370), and force-balance accelerometers (Nanometrics Titan, Güralp Fortimus, Kinemetrics EpiSensor ES-T, and Lunitek Triton). The New High Noise Model (NHNM) and New Low Noise Model (NLNM) [
35], together with reference spectra (thick black curves) for earthquakes of increasing magnitude (M1.5–M5.5) at an epicentral distance of 10 km (dashed grey curves), are included for comparison. The PSD curves were derived from the self-noise specifications reported in the corresponding manufacturers’ datasheets.
Figure 2.
Comparison of signal processing chains between a conventional capacitive MEMS accelerometer (a) and the QMEMS accelerometer (b). The QMEMS approach utilizes direct acceleration-to-frequency conversion and digital frequency measurement, avoiding the conventional analog amplification and voltage-based A/D conversion chain, thereby reducing the associated sources of noise and non-linearity.
Figure 2.
Comparison of signal processing chains between a conventional capacitive MEMS accelerometer (a) and the QMEMS accelerometer (b). The QMEMS approach utilizes direct acceleration-to-frequency conversion and digital frequency measurement, avoiding the conventional analog amplification and voltage-based A/D conversion chain, thereby reducing the associated sources of noise and non-linearity.
Figure 3.
a) External view of the M-A370 accelerometer module. The compact aluminum housing (48×24×16 mm³) contains the 3-axis sensor assembly and digital readout electronics, offering a significant size advantage for dense array deployments. b) Internal structure of the M-A370 sensor module. The exploded view (left) shows the differential arrangement of the six transducer units (two per axis) mounted within the package. The detailed schematic (top right) illustrates the high-density tungsten proof mass attached to the quartz cantilever and DETF resonator. The SEM image (bottom right) shows the fabricated Double-Ended Tuning Fork (DETF) resonator.
Figure 3.
a) External view of the M-A370 accelerometer module. The compact aluminum housing (48×24×16 mm³) contains the 3-axis sensor assembly and digital readout electronics, offering a significant size advantage for dense array deployments. b) Internal structure of the M-A370 sensor module. The exploded view (left) shows the differential arrangement of the six transducer units (two per axis) mounted within the package. The detailed schematic (top right) illustrates the high-density tungsten proof mass attached to the quartz cantilever and DETF resonator. The SEM image (bottom right) shows the fabricated Double-Ended Tuning Fork (DETF) resonator.
Figure 4.
Functional architecture of the multiparametric Smart Sensor Box used in this study. The system integrates a high-sensitivity Epson M-A352/M-A370 QMEMS accelerometer, a GNSS-disciplined Time Base Unit (TBU), a dual-core STM32H7 microcontroller, and an embedded single-board computer providing edge-computing capabilities. The modular architecture supports deterministic timing, local data processing, multi-parameter environmental sensing, and distributed communication through wired and wireless interfaces, providing a scalable framework for seismic monitoring and Structural Health Monitoring applications.
Figure 4.
Functional architecture of the multiparametric Smart Sensor Box used in this study. The system integrates a high-sensitivity Epson M-A352/M-A370 QMEMS accelerometer, a GNSS-disciplined Time Base Unit (TBU), a dual-core STM32H7 microcontroller, and an embedded single-board computer providing edge-computing capabilities. The modular architecture supports deterministic timing, local data processing, multi-parameter environmental sensing, and distributed communication through wired and wireless interfaces, providing a scalable framework for seismic monitoring and Structural Health Monitoring applications.
Figure 5.
Experimental setup adopted for the laboratory validation of the proposed QMEMS-based acquisition systems at the L.E.D.A. (Laboratory of Earthquake Engineering and Dynamic Analysis) Research Centre. (a) Complete Smart Sensor Box mounted on the SPEKTRA APS129 electrodynamic shaker for testing under the operational configuration. (b) Epson M-A352 and M-A370 QMEMS sensing elements mounted side-by-side on a rigid fixture for direct comparison of the sensor response independently of the enclosure. (c) Güralp Fortimus force-balance accelerometer installed on the shaker table a. (d) Frequency–amplitude excitation scheme adopted to evaluate sensor performance over the range of amplitudes and frequencies shown in panel.
Figure 5.
Experimental setup adopted for the laboratory validation of the proposed QMEMS-based acquisition systems at the L.E.D.A. (Laboratory of Earthquake Engineering and Dynamic Analysis) Research Centre. (a) Complete Smart Sensor Box mounted on the SPEKTRA APS129 electrodynamic shaker for testing under the operational configuration. (b) Epson M-A352 and M-A370 QMEMS sensing elements mounted side-by-side on a rigid fixture for direct comparison of the sensor response independently of the enclosure. (c) Güralp Fortimus force-balance accelerometer installed on the shaker table a. (d) Frequency–amplitude excitation scheme adopted to evaluate sensor performance over the range of amplitudes and frequencies shown in panel.
Figure 6.
Measured relative amplitude error as a function of excitation frequency for the three configurations: (a) Sensor Box integrating the Epson M-A370 QMEMS (b) standalone Epson M-A370 QMEMS accelerometer mounted directly on the shake table and (c) Güralp Fortimus force-balance accelerometer. The curves correspond to sinusoidal excitation tests performed at nominal peak acceleration amplitudes of 0.1, 0.5, 1, 3, and 5 m s⁻².
Figure 6.
Measured relative amplitude error as a function of excitation frequency for the three configurations: (a) Sensor Box integrating the Epson M-A370 QMEMS (b) standalone Epson M-A370 QMEMS accelerometer mounted directly on the shake table and (c) Güralp Fortimus force-balance accelerometer. The curves correspond to sinusoidal excitation tests performed at nominal peak acceleration amplitudes of 0.1, 0.5, 1, 3, and 5 m s⁻².
Figure 7.
Experimental setup and representative stabilization diagrams obtained during the dynamic testing campaign on the 2:3 scale masonry building at the L.E.D.A. laboratory (modified from [
41]. The comparison among the M-A352 (a), PCB 356A17 (b), VS1002 (c), and ADXL355 (d) sensors highlights the influence of their self-noise on the identification of structural modal frequencies under ambient-vibration conditions, with the M-A352 and PCB 356A17 providing the clearest stabilization diagrams.
Figure 7.
Experimental setup and representative stabilization diagrams obtained during the dynamic testing campaign on the 2:3 scale masonry building at the L.E.D.A. laboratory (modified from [
41]. The comparison among the M-A352 (a), PCB 356A17 (b), VS1002 (c), and ADXL355 (d) sensors highlights the influence of their self-noise on the identification of structural modal frequencies under ambient-vibration conditions, with the M-A352 and PCB 356A17 providing the clearest stabilization diagrams.
Figure 8.
Relative timing drift between the GNSS-derived PPS and the synthetic PPS during holdover operation. The reported values represent the measured timing difference between the GNSS-derived PPS and the synthetic PPS generated by the TBU as a function of elapsed time after GNSS disconnection. The timing offset was determined from direct oscilloscope measurements of the temporal separation between the rising edges of the two PPS signals. During the approximately 4.4 h experiment, the relative timing error remained within ± 80 μs.
Figure 8.
Relative timing drift between the GNSS-derived PPS and the synthetic PPS during holdover operation. The reported values represent the measured timing difference between the GNSS-derived PPS and the synthetic PPS generated by the TBU as a function of elapsed time after GNSS disconnection. The timing offset was determined from direct oscilloscope measurements of the temporal separation between the rising edges of the two PPS signals. During the approximately 4.4 h experiment, the relative timing error remained within ± 80 μs.
Figure 9.
Example of synchronization validation between two co-located Smart Sensor Boxes (CT02 and CT04) operating under GNSS-disciplined timing and sampling simultaneously at 1000 Hz. The upper panels compare the recorded N–S acceleration waveforms and their Power Spectral Density (PSD), showing an excellent agreement in both the time and frequency domains. The central panels report the global coherence and phase difference between the two signals, demonstrating coherence values close to unity and negligible phase shifts across the frequency band of interest. The lower panels show the evolution of coherence and cross-correlation within selected frequency bands (0.5–5 Hz, 5–10 Hz, and 10–30 Hz), confirming the stability of the synchronization throughout the entire recording. Sub-sample cross-correlation analysis yields a residual relative timing offset of only 0.009 samples, corresponding to approximately 9.1 μs, with a correlation coefficient equal to 1.0 and an equivalent phase difference of only 0.02° at 5 Hz.
Figure 9.
Example of synchronization validation between two co-located Smart Sensor Boxes (CT02 and CT04) operating under GNSS-disciplined timing and sampling simultaneously at 1000 Hz. The upper panels compare the recorded N–S acceleration waveforms and their Power Spectral Density (PSD), showing an excellent agreement in both the time and frequency domains. The central panels report the global coherence and phase difference between the two signals, demonstrating coherence values close to unity and negligible phase shifts across the frequency band of interest. The lower panels show the evolution of coherence and cross-correlation within selected frequency bands (0.5–5 Hz, 5–10 Hz, and 10–30 Hz), confirming the stability of the synchronization throughout the entire recording. Sub-sample cross-correlation analysis yields a residual relative timing offset of only 0.009 samples, corresponding to approximately 9.1 μs, with a correlation coefficient equal to 1.0 and an equivalent phase difference of only 0.02° at 5 Hz.

Figure 10.
Drift rate of the CT02–CT04 timing offset during GNSS holdover operation and comparison with the external ambient temperature measured by the BME280 sensor. The left axis shows the first derivative of the timing drift (μs h⁻¹), whereas the right axis reports the corresponding temperature evolution. Labels indicate the cumulative timing drift (ms) measured at each epoch.
Figure 10.
Drift rate of the CT02–CT04 timing offset during GNSS holdover operation and comparison with the external ambient temperature measured by the BME280 sensor. The left axis shows the first derivative of the timing drift (μs h⁻¹), whereas the right axis reports the corresponding temperature evolution. Labels indicate the cumulative timing drift (ms) measured at each epoch.
Figure 11.
Southern Italy map with the spatial distribution of the urban seismic observatories of CMSU. Colored squares indicate the epicenters of historical earthquakes with magnitude larger than 5.5. In the insets are shown the investigated areas together with the station distribution (yellow triangles) within the urban environments of Catania, Campi Flegrei, Messina–Reggio Calabria, Potenza, port of Vulcano, and Ragusa. In Catania, Potenza, and Ragusa, QMEMS sensors are integrated with force-balance and conventional MEMS accelerometers, whereas only QMEMS deployments are represented for the Messina–Reggio Calabria and Campi Flegrei volcanic area, where more than 50 QMEMS M-A352 accelerometer are currently in operation.
Figure 11.
Southern Italy map with the spatial distribution of the urban seismic observatories of CMSU. Colored squares indicate the epicenters of historical earthquakes with magnitude larger than 5.5. In the insets are shown the investigated areas together with the station distribution (yellow triangles) within the urban environments of Catania, Campi Flegrei, Messina–Reggio Calabria, Potenza, port of Vulcano, and Ragusa. In Catania, Potenza, and Ragusa, QMEMS sensors are integrated with force-balance and conventional MEMS accelerometers, whereas only QMEMS deployments are represented for the Messina–Reggio Calabria and Campi Flegrei volcanic area, where more than 50 QMEMS M-A352 accelerometer are currently in operation.
Figure 12.
Example of continuous seismic recordings acquired at station OCF04 of the Campi Flegrei Urban Seismic Observatory using an Epson M-A352 QMEMS accelerometer. The left panel shows a continuous record collected on 1 September 2025, including the Md 4.0 earthquake together with several smaller local events. The panel on the right illustrates the three-component recording of a Md 1.0 earthquake, highlighting the capability of the QMEMS sensor to detect and accurately record very weak local seismic events.
Figure 12.
Example of continuous seismic recordings acquired at station OCF04 of the Campi Flegrei Urban Seismic Observatory using an Epson M-A352 QMEMS accelerometer. The left panel shows a continuous record collected on 1 September 2025, including the Md 4.0 earthquake together with several smaller local events. The panel on the right illustrates the three-component recording of a Md 1.0 earthquake, highlighting the capability of the QMEMS sensor to detect and accurately record very weak local seismic events.
Figure 13.
Response of the Messina–Reggio Calabria Urban Seismic Observatory (OSU) stations to the deep Mw 6.1 earthquake of 1 June 2026. (a) Overview of the event, showing the vertical (Z) component waveforms recorded at twenty OSU stations, the corresponding regional INGV ShakeMap (
https://shakemap.ingv.it/data/46107472/current/products/pga.jpg), and the geographical distribution of the stations across the Messina Strait. (b) Spatial distribution of the maximum PGV and PGA values. For each station, the reported values correspond to the maximum among the three components (Z, N, and E). Circle size and colour are proportional to the measured amplitudes, while the numerical labels indicate the corresponding PGV (cm s⁻¹) and PGA (% g) values.
Figure 13.
Response of the Messina–Reggio Calabria Urban Seismic Observatory (OSU) stations to the deep Mw 6.1 earthquake of 1 June 2026. (a) Overview of the event, showing the vertical (Z) component waveforms recorded at twenty OSU stations, the corresponding regional INGV ShakeMap (
https://shakemap.ingv.it/data/46107472/current/products/pga.jpg), and the geographical distribution of the stations across the Messina Strait. (b) Spatial distribution of the maximum PGV and PGA values. For each station, the reported values correspond to the maximum among the three components (Z, N, and E). Circle size and colour are proportional to the measured amplitudes, while the numerical labels indicate the corresponding PGV (cm s⁻¹) and PGA (% g) values.
Figure 14.
Example of the products generated by PSAVE at the ME016 station of Messina network during the the Mw 6.1 earthquake of 1 June 2026 (22:12:35 UTC), which occurred at a hypocentral depth of approximately 250 km from Messina–Reggio Calabria Urban Seismic Observatory networks for. (a) PSAVE impact-based EEW summary showing both predictive and observed parameters; (b) pseudo-spectral acceleration curves and three-component acceleration time histories.
Figure 14.
Example of the products generated by PSAVE at the ME016 station of Messina network during the the Mw 6.1 earthquake of 1 June 2026 (22:12:35 UTC), which occurred at a hypocentral depth of approximately 250 km from Messina–Reggio Calabria Urban Seismic Observatory networks for. (a) PSAVE impact-based EEW summary showing both predictive and observed parameters; (b) pseudo-spectral acceleration curves and three-component acceleration time histories.
Figure 15.
Adopted sensor layout at the clock tower during the ambient vibration campaign. M-A352 accelerometers were installed at the ground floor and at the clock level of the tower.
Figure 15.
Adopted sensor layout at the clock tower during the ambient vibration campaign. M-A352 accelerometers were installed at the ground floor and at the clock level of the tower.
Figure 16.
Experimental setup and comparison of the self-noise characteristics of the M-A352 and M-A370 accelerometers at the clock tower of the Durham Castle (United Kingdom). Power spectral densities (PSDs) calculated from one hour (01:00–02:00 UTC) of ambient-noise recordings acquired simultaneously by co-located sensors installed at the top and base of the Clock Tower reveal several structural resonance frequencies between 2 and 10 Hz, with larger amplification observed at the top of the structure.
Figure 16.
Experimental setup and comparison of the self-noise characteristics of the M-A352 and M-A370 accelerometers at the clock tower of the Durham Castle (United Kingdom). Power spectral densities (PSDs) calculated from one hour (01:00–02:00 UTC) of ambient-noise recordings acquired simultaneously by co-located sensors installed at the top and base of the Clock Tower reveal several structural resonance frequencies between 2 and 10 Hz, with larger amplification observed at the top of the structure.
Figure 17.
Structural Health Monitoring of the “Papa Giovanni XXIII” reinforced-concrete arch bridge in Ragusa (southern Italy). (a) General view of the bridge. (b) Stabilization diagram obtained from Operational Modal Analysis (OMA) during the dynamic characterization together with the corresponding first bending mode shape identified at 1.61 Hz. (c) Comparison of power spectral densities computed from one hour of nighttime ambient-vibration recordings acquired by co-located ADXL355 and M-A352 accelerometers, highlighting the lower self-noise and improved detectability of structural resonances provided by the M-A352 sensor. (d) Automatic tracking of the first five modal frequencies over several months for the three components at a representative monitoring node.
Figure 17.
Structural Health Monitoring of the “Papa Giovanni XXIII” reinforced-concrete arch bridge in Ragusa (southern Italy). (a) General view of the bridge. (b) Stabilization diagram obtained from Operational Modal Analysis (OMA) during the dynamic characterization together with the corresponding first bending mode shape identified at 1.61 Hz. (c) Comparison of power spectral densities computed from one hour of nighttime ambient-vibration recordings acquired by co-located ADXL355 and M-A352 accelerometers, highlighting the lower self-noise and improved detectability of structural resonances provided by the M-A352 sensor. (d) Automatic tracking of the first five modal frequencies over several months for the three components at a representative monitoring node.
Figure 18.
Instrumented buildings ED001–ED003 and corresponding stabilization diagrams derived from ambient vibration measurements through the SSI-Data method. Stable poles indicate the identified modal frequencies, revealing two flexural modes and one torsional mode for each structure.
Figure 18.
Instrumented buildings ED001–ED003 and corresponding stabilization diagrams derived from ambient vibration measurements through the SSI-Data method. Stable poles indicate the identified modal frequencies, revealing two flexural modes and one torsional mode for each structure.
Figure 19.
Pseudo-acceleration response spectra and three-component acceleration time histories of the Md 4.4 earthquake of 20 May 2024 recorded at stations OCF03 (top panels, basement level) and OCF02 (bottom panels, rooftop level) installed in building ED001 (see
Figure 18). The Husid plot (green line), representing the cumulative energy release of the signal, is superimposed on the time histories and was used to estimate the significant durations T75 and T90, whose values are indicated by the red arrows. Comparison between the basement and rooftop recordings highlights the amplification effects associated with the dynamic response of the structure.
Figure 19.
Pseudo-acceleration response spectra and three-component acceleration time histories of the Md 4.4 earthquake of 20 May 2024 recorded at stations OCF03 (top panels, basement level) and OCF02 (bottom panels, rooftop level) installed in building ED001 (see
Figure 18). The Husid plot (green line), representing the cumulative energy release of the signal, is superimposed on the time histories and was used to estimate the significant durations T75 and T90, whose values are indicated by the red arrows. Comparison between the basement and rooftop recordings highlights the amplification effects associated with the dynamic response of the structure.
Figure 20.
Automatic tracking of the modal frequencies of building ED001 derived from power spectral density (PSD) analysis during the Md 3.5 and Md 4.4 earthquakes of 20 May 2024. The upper panel shows the vertical acceleration record (HNZ), while the lower panels display the evolution of the spectral peaks identified on the three components (HNZ, HNN, and HNE). Zoomed views (right panels) highlight the temporary decrease of the first modal frequency during the Md 4.4 earthquake and its subsequent recovery; a similar but less pronounced behavior was previously observed during the preceding Md 3.5 event. The conceptual sketch illustrates the temporary reduction of structural stiffness during shaking, associated with soil–structure interaction and nonlinear structural behavior, and the subsequent recovery of the original dynamic properties after the event, indicating the absence of permanent damage.
Figure 20.
Automatic tracking of the modal frequencies of building ED001 derived from power spectral density (PSD) analysis during the Md 3.5 and Md 4.4 earthquakes of 20 May 2024. The upper panel shows the vertical acceleration record (HNZ), while the lower panels display the evolution of the spectral peaks identified on the three components (HNZ, HNN, and HNE). Zoomed views (right panels) highlight the temporary decrease of the first modal frequency during the Md 4.4 earthquake and its subsequent recovery; a similar but less pronounced behavior was previously observed during the preceding Md 3.5 event. The conceptual sketch illustrates the temporary reduction of structural stiffness during shaking, associated with soil–structure interaction and nonlinear structural behavior, and the subsequent recovery of the original dynamic properties after the event, indicating the absence of permanent damage.
Table 1.
Comparison of representative force-balance (FBA), MEMS and QMEMS accelerometers for seismic monitoring and Structural Health Monitoring (SHM) applications. Most of the listed sensors were employed in the experimental activities presented in this work. Specifications are taken from the corresponding manufacturers’ datasheets.
Table 1.
Comparison of representative force-balance (FBA), MEMS and QMEMS accelerometers for seismic monitoring and Structural Health Monitoring (SHM) applications. Most of the listed sensors were employed in the experimental activities presented in this work. Specifications are taken from the corresponding manufacturers’ datasheets.
Table 2.
Comparison of representative IEPE (ICP®) accelerometers commonly used for structural dynamics and laboratory testing. Most of the listed sensors were employed in the experimental activities presented in this work. Specifications are taken from the corresponding manufacturers’ datasheets; dynamic range values marked with (*) were estimated as described in the text.
Table 2.
Comparison of representative IEPE (ICP®) accelerometers commonly used for structural dynamics and laboratory testing. Most of the listed sensors were employed in the experimental activities presented in this work. Specifications are taken from the corresponding manufacturers’ datasheets; dynamic range values marked with (*) were estimated as described in the text.
Table 3.
Typical synchronization requirements for seismic monitoring, Operational Modal Analysis (OMA), and Structural Health Monitoring (SHM) applications.
Table 3.
Typical synchronization requirements for seismic monitoring, Operational Modal Analysis (OMA), and Structural Health Monitoring (SHM) applications.
Table 4.
Measured synchronization performance of the system investigated in this study.
Table 4.
Measured synchronization performance of the system investigated in this study.
Table 5.
Modal parameters identified by Operational Modal Analysis for the first three monitored buildings (ED001–ED003), including modal frequency, damping ratio, and complexity index.
Table 5.
Modal parameters identified by Operational Modal Analysis for the first three monitored buildings (ED001–ED003), including modal frequency, damping ratio, and complexity index.
Table 6.
Indicative engineering guidelines relating accelerometer self-noise and timing synchronization requirements for seismic monitoring, Operational Modal Analysis (OMA), and Structural Health Monitoring (SHM) applications. The reported values represent indicative performance ranges derived from laboratory experiments, field deployments, and representative applications discussed in this work, complemented by evidence from the literature. They should not be interpreted as strict thresholds, since actual requirements depend on the monitored structure, ambient vibration levels, monitoring objectives, and the target application.
Table 6.
Indicative engineering guidelines relating accelerometer self-noise and timing synchronization requirements for seismic monitoring, Operational Modal Analysis (OMA), and Structural Health Monitoring (SHM) applications. The reported values represent indicative performance ranges derived from laboratory experiments, field deployments, and representative applications discussed in this work, complemented by evidence from the literature. They should not be interpreted as strict thresholds, since actual requirements depend on the monitored structure, ambient vibration levels, monitoring objectives, and the target application.