## 4 Results  
The likelihood of the refitted models declined with each step change in the controlled parameter away from its optimal value, as shown for New Zealand in Figure 6. The results for California were similar. The log-likelihood of the refitted model is plotted against the controlled spatial scaling factor in Figure 6a and against the temporal scaling factor in Figure 6b. An order of magnitude change in each scaling factor induced a modest reduction in the log-likelihood. The maximum reduction of about 34 units corresponded to an information loss per earthquake of about 0.2 relative to the overall optimal fit.  
The refitted mixing parameter \(\mu\) tended to increase as the controlled parameter shifted further away from its optimal value, as shown for New Zealand in Figure 7. Again, the results were similar for California. The variation of \(\mu\) with the spatial scaling factor is shown in Figure 7a and against the temporal scaling factor in Figure 7b. The values of \(\mu\) increased from about 0.15 at the optimal fit to greater than 0.5 when the temporal or spatial scaling factors were changed by an order of magnitude. The \(\mu\) value represents the  
Figure 6: Log-likelihood of EEPAS model fitted with controlled values of (**a**) \(\sigma_{A}\) (Table 2) and (**b**) \(a_{T}\) (Table 3) to the New Zealand earthquake catalogue.  
\begin{table}
\begin{tabular}{c c} \hline \hline
**NZ EEPAS-1F** & **California EEPAS-1F** \\ \hline
2.49 & 3.16 \\
2.34 & 3.01 \\
2.19 & 2.86 \\
2.04 & 2.71 \\
1.89 & 2.56 \\
1.74 & 2.41 \\
1.59 & 2.26 \\
1.44 \({}^{\dagger}\) & 2.11 \({}^{\ddagger}\) \\
1.29 & 1.96 \\
1.14 & 1.81 \\
0.99 & 1.66 \\
0.84 & 1.51 \\
0.69 & 1.36 \\
0.54 & 1.21 \\
0.39 & 1.06 \\ \hline \hline \end{tabular} \({}^{\dagger}\) Fitted.  
\end{table}
Table 3: Controlled values of \(a_{T}\) in EEPAS_1F models for NZ and California.  
proportional contribution of the background model to the total EEPAS model rate density. Higher \(\mu\) values thus indicate a greater contribution of the background component and a smaller contribution of the time-varying component. In other words, higher \(\mu\) values indicate that there were fewer target earthquakes with precursors matching the changed spatial and temporal distributions.  
As the controlled parameter was changed, the refitted values of the other parameters changed in a way that was consistent with the notion of a space-time trade-off. The results are shown for New Zealand in Figure 8a and for California in Figure 8b.  
Figure 8: Trade-off of spatial and temporal scaling factors \({\sigma_{A}}^{2}\) and \(10^{\mu\tau}\), respectively, revealed by the fit of the EEPAS model with controlled values of \(\sigma_{A}\) (blue triangles) and \(a_{T}\) (black squares). The straight line with a slope of \(-1\) represents an even trade-off between space and time. (**a**) New Zealand. (**b**) California.  
Figure 7: Fitted values of mixing parameter \(\mu\) (\(0\ \leq\ \mu\ \leq\ 1\)) of the EEPAS model fitted with controlled values of (**a**) \(\sigma_{A}\) and (**b**) \(a_{T}\) to the New Zealand earthquake catalogue.  
In each plot, the pairs of scaling factors resulting from controlling \(\sigma_{A}\) are shown as blue triangles, and those resulting from controlling \(a_{T}\) are shown as black squares. The temporal scaling factor decreased as the controlled spatial scaling factor increased, and the spatial scaling factor decreased as the controlled temporal scaling factor increased. However, the curves had different slopes depending on whether \(\sigma_{A}\) or \(a_{T}\) was the controlled variable. An even trade-off line with a slope of \(-1\) is drawn through the intersection of the two curves (straight blue line in Figure 8a,b). Its slope lies between the average slopes of the two controlled fitting curves.