Earthquake forecasting probability by statistical correlations between low to moderate seismic events and variations in geochemical parameters&#160;
Authors: Pierotti Lisa , Fidani Cristiano , Facca Gianluca , Gherardi Fabrizio 
Published: March 2022

&lt;p&gt;Since late 2002, a network of six automatic monitoring stations is operating in Tuscany, Central Italy, to investigate possible geochemical precursors of earthquakes. The network is operated by the Institute of Geosciences and Earth Resources&amp;#160;(IGG), of&amp;#160;the National Research Council of Italy&amp;#160;(CNR), in collaboration and with the financial support of the Government of the Tuscany Region. The areas of highest seismic risk of the region, Garfagnana, Lunigiana, Mugello, Upper Tiber Valley and Mt. Amiata, are currently investigated. The monitoring stations are equipped with multi-parametric sensors to measure temperature, &lt;em&gt;pH&lt;/em&gt;, electric conductivity, redox potential, dissolved CO&lt;sub&gt;2&lt;/sub&gt; and CH&lt;sub&gt;4&lt;/sub&gt; concentrations in spring waters. The elaboration of long-term time series allowed for an accurate definition of the geochemical background, and for the recognition of a number of geochemical anomalies in concomitance with the most energetic seismic events occurred during the monitoring period (Pierotti et al., 2017).&lt;/p&gt;&lt;p&gt;In an attempt to further exploit data from the geochemical network of Tuscany in a seismic risk reduction perspective, here we present a new statistical analysis that focuses on the possible correlation between low to moderate seismic events and variations in the chemical-physical parameters detected by the monitoring network. This approach relies on the estimate of a conditional probability for the forecast of earthquakes from the correlation coefficient between seismic events and signals variations (Fidani, 2021).&lt;/p&gt;&lt;p&gt;Seismic events (&lt;em&gt;EQ&lt;/em&gt;) are classified according to a magnitude threshold, &lt;em&gt;M&lt;sub&gt;o&lt;/sub&gt;&lt;/em&gt;. We set &lt;em&gt;EQ&lt;/em&gt; = 0, if no seismic events were observed with &lt;em&gt;M &lt; M&lt;sub&gt;o&lt;/sub&gt;&lt;/em&gt;, and &lt;em&gt;EQ&lt;/em&gt; = 1, if at least a seismic event was observed with &lt;em&gt;M &gt; M&lt;sub&gt;o&lt;/sub&gt;&lt;/em&gt;. Chemical-physical (&lt;em&gt;CP&lt;/em&gt;) events were defined based on their appropriate amplitudes threshold &lt;em&gt;A&lt;sub&gt;o&lt;/sub&gt;&lt;/em&gt;, being &lt;em&gt;CP&lt;/em&gt; = 0 if the amplitude &lt;em&gt;A &lt; A&lt;sub&gt;o&lt;/sub&gt;&lt;/em&gt;, and &lt;em&gt;CP&lt;/em&gt; = 1 if &lt;em&gt;A &gt; A&lt;sub&gt;o&lt;/sub&gt;&lt;/em&gt;. Digital time series were elaborated from data collected over the last 10 years, where &lt;em&gt;EQs&lt;/em&gt; were declustered and &lt;em&gt;CPs&lt;/em&gt; detrended for external influences. The couples of events with the same time differences &lt;em&gt;T&lt;sub&gt;EQ&lt;/sub&gt; &amp;#8211; T&lt;sub&gt;CP&lt;/sub&gt;&lt;/em&gt;, between &lt;em&gt;EQs&lt;/em&gt; and &lt;em&gt;CPs&lt;/em&gt;, were summed in a histogram. Then, a Pearson statistical correlation coefficient &lt;em&gt;corr(EQ,CP)&lt;/em&gt; was obtained starting from the covariance definition.&lt;/p&gt;&lt;p&gt;A conditional probability for &lt;em&gt;EQ&lt;/em&gt; forecasting is estimated starting from the correlation coefficient in an attempt to use data from &lt;em&gt;CP&lt;/em&gt; network of Tuscany in a seismic risk reduction framework. The approach consists in an evaluation of &lt;em&gt;EQ&lt;/em&gt; probability in a defined area, given a &lt;em&gt;CP&lt;/em&gt; detection by the station in the same area. The conditional probability &lt;em&gt;P(EQCP)&lt;/em&gt;, when a correlation between &lt;em&gt;EQs&lt;/em&gt; and &lt;em&gt;CPs&lt;/em&gt; exists and time difference is that evidenced by the correlation, is increased by a term proportional to the correlation coefficient as&lt;/p&gt;&lt;p&gt;&lt;img src=&quot;https://contentmanager.copernicus.org/fileStorageProxy.php?f=gnp.9e880d2041e161891512461/sdaolpUECMynit/22UGE&amp;app=m&amp;a=0&amp;c=d3b4ea5ab3bdd75edb9db103b582241c&amp;ct=x&amp;pn=gnp.elif&amp;d=1&quot; alt=&quot;&quot;&gt;&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;with respect to the unconditioned probability &lt;em&gt;P(EQ)&lt;/em&gt; when a &lt;em&gt;CP&lt;/em&gt; event is detected, where &lt;em&gt;P(CP)&lt;/em&gt; is the unconditioned probability of &lt;em&gt;CP&lt;/em&gt;.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;Fidani, C. (2021). Front. Earth Sci. 9:673105.&lt;/p&gt;&lt;p&gt;Pierotti, L. et al. (2017).&amp;#160;Physics and Chemistry of the Earth, Parts A/B/C,&amp;#160;98, 161-172.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;