Data to reproduce the figures in the Optica paper "Earth-Scale Precursors of the Kamchatka Earthquake Revealed via Polarization Monitoring in Undersea Mediterranean Optical Fiber Links" by A. Mecozzi et al.
收藏资源简介：
Data and code to reproduce the figures in the Optica paperA. Mecozzi et al., “Earth-Scale Precursors of the Kamchatka Earthquake Revealed via Polarization Monitoring in Undersea Mediterranean Optical Fiber Links,” Optica, https://doi.org/10.1364/OPTICA.586317 This repository provides the datasets required to reproduce the figures in the above Optica paper, together with MATLAB software to analyze and visualize them. Data files Each .mat file contains the measured transmission matrix in Jones space J (referred to as Jones matrix, from now on) and the corresponding Unix time vector time. The file names encode the start and end timestamps of the dataset. 2025-04-29T000000Z-2026-02-11T000000Z_CHA_ATH.matJones matrix J extracted at the Chania receiver of the Athens–Chania link, with corresponding Unix time time. 2025-03-01T000000Z-2026-02-11T000000Z_GEN_MAR.matJones matrix J extracted at the Genova receiver of the Marseille–Genova link, with corresponding Unix time time. 2025-04-20T000000Z-2026-02-11T000000Z_CAT_HAI.matJones matrix J extracted at the Catania receiver of the Haifa–Catania link, with corresponding Unix time time. 2025-03-03T000000Z-2026-02-11T000000Z_CAT_TAV.matJones matrix J extracted at the Catania receiver of the Tel Aviv–Catania link, with corresponding Unix time time. 2025-05-27T000000Z-2026-02-11T000000Z_ROM_GAR.matJones matrix J extracted at the Roma receiver of the Golfo Aranci–Roma link, with corresponding Unix time time. Processing with Jones.m When any of the above files are provided as input to Jones.m, the script returns the following output variables: phi: the three components of the bare rotation vector. phi1: the three components of the rotation vector compensated with respect to a reference rotation vector, computed by averaging over a time window specified (in hours) by the input parameter hours_averaging. pdl: the polarization-dependent loss (PDL), in dB. time: the Unix time vector (identical to the input time). Twin: the averaging time window corresponding to hours_averaging, expressed in seconds. The implementation follows the algorithm described in:A. Mecozzi et al., “Geophysical sensing using Jones matrices extracted from submarine optical cable transceivers carrying live traffic,” Optica 12, 1712–1719 (2025), https://doi.org/10.1364/OPTICA.57288 Plot reproduction with Plot_and_stats.m The outputs produced by Jones.m can be provided as input to Plot_and_stats.m to reproduce the figures in:A. Mecozzi et al., “Earth-Scale Precursors of the Kamchatka Earthquake Revealed via Polarization Monitoring in Undersea Mediterranean Optical Fiber Links,” Optica, https://doi.org/10.1364/OPTICA.586317
本仓库提供复现Optica论文（A. Mecozzi等人所著《基于地中海海底光纤链路偏振监测揭示堪察加地震的全球尺度前兆》，Optica, https://doi.org/10.1364/OPTICA.586317）中所有图表所需的数据与代码。 本仓库提供复现上述Optica论文中所有图表所需的数据集，以及用于分析和可视化这些数据的MATLAB软件。 数据文件 每个.mat格式文件均包含琼斯空间（Jones space）中测得的传输矩阵J（下文简称琼斯矩阵），以及对应的Unix时间向量time。文件名编码了该数据集的起始与结束时间戳。 2025-04-29T000000Z-2026-02-11T000000Z_CHA_ATH.mat：雅典-干尼亚链路的干尼亚接收端提取的琼斯矩阵J，以及对应的Unix时间向量time。 2025-03-01T000000Z-2026-02-11T000000Z_GEN_MAR.mat：马赛-热那亚链路的热那亚接收端提取的琼斯矩阵J，以及对应的Unix时间向量time。 2025-04-20T000000Z-2026-02-11T000000Z_CAT_HAI.mat：海法-卡塔尼亚链路的卡塔尼亚接收端提取的琼斯矩阵J，以及对应的Unix时间向量time。 2025-03-03T000000Z-2026-02-11T000000Z_CAT_TAV.mat：特拉维夫-卡塔尼亚链路的卡塔尼亚接收端提取的琼斯矩阵J，以及对应的Unix时间向量time。 2025-05-27T000000Z-2026-02-11T000000Z_ROM_GAR.mat：阿尔安奇湾-罗马链路的罗马接收端提取的琼斯矩阵J，以及对应的Unix时间向量time。 使用Jones.m进行处理 将上述任意文件作为输入传入Jones.m脚本时，该脚本将返回以下输出变量： phi：裸旋转矢量的三个分量。 phi1：相对于参考旋转矢量补偿后的旋转矢量的三个分量，参考旋转矢量通过输入参数hours_averaging指定的时间窗口（单位：小时）内取平均计算得到。 pdl：偏振相关损耗（PDL，单位：dB）。 time：Unix时间向量（与输入的time完全一致）。 Twin：与hours_averaging对应的平均时间窗口，单位为秒。 本实现遵循下述文献中描述的算法：A. Mecozzi等人，"基于从承载实时流量的海底光缆收发器提取的琼斯矩阵开展地球物理传感"，Optica 12, 1712–1719 (2025), https://doi.org/10.1364/OPTICA.57288 使用Plot_and_stats.m复现图表 将Jones.m生成的输出作为输入传入Plot_and_stats.m，即可复现下述论文中的图表：A. Mecozzi等人，"基于地中海海底光纤链路偏振监测揭示堪察加地震的全球尺度前兆"，Optica, https://doi.org/10.1364/OPTICA.586317
