We statistically investigate the Coronavirus Disease 19 (COVID-19) pandemic, which became particularly invasive in Italy in March 2020. We show that the high apparent lethality or case fatality ratio (CFR) observed in Italy, as compared with other countries, is likely biased by a strong underestimation of the number of infection cases. To give a more realistic estimate of the lethality of COVID-19, we use the actual (March 2020) estimates of the infection fatality ratio (IFR) of the pandemic based on the minimum observed CFR and analyze data obtained from the Diamond Princess cruise ship, a good representation of a “laboratory” case-study from an isolated system in which all the people have been tested. From such analyses, we derive more realistic estimates of the real extent of the infection as well as more accurate indicators of how fast the infection propagates. We then isolate the dominant factors causing the abnormal severity of the disease in Italy. Finally, we use the death count—the only data estimated to be reliable enough—to predict the total number of people infected and the interval of time when the infection in Italy could end.

The complexity of coseismic slip distributions influences the tsunami hazard posed by local and, to a certain extent, distant tsunami sources. Large slip concentrated in shallow patches was observed in recent tsunamigenic earthquakes, possibly due to dynamic amplification near the free surface, variable frictional conditions or other factors. We propose a method for incorporating enhanced shallow slip for subduction earthquakes while preventing systematic slip excess at shallow depths over one or more seismic cycles. The method uses the classic k⁻² stochastic slip distributions, augmented by shallow slip amplification. It is necessary for deep events with lower slip to occur more often than shallow ones with amplified slip to balance the long-term cumulative slip. We evaluate the impact of this approach on tsunami hazard in the central and eastern Mediterranean Sea adopting a realistic 3D geometry for three subduction zones, by using it to model ~ 150,000 earthquakes with \(M_{w}\) from 6.0 to 9.0. We combine earthquake rates, depth-dependent slip distributions, tsunami modeling, and epistemic uncertainty through an ensemble modeling technique. We found that the mean hazard curves obtained with our method show enhanced probabilities for larger inundation heights as compared to the curves derived from depth-independent slip distributions. Our approach is completely general and can be applied to any subduction zone in the world.

The 2011 Tohoku earthquake produced an unexpected large amount of shallow slip greatly contributing to the ensuing tsunami. How frequent are such events? How can they be efficiently modelled for tsunami hazard? Stochastic slip models, which can be computed rapidly, are used to explore the natural slip variability; however, they generally do not deal specifically with shallow slip features. We study the systematic depth-dependence of slip along a thrust fault with a number of 2D dynamic simulations using stochastic shear stress distributions and a geometry based on the cross section of the Tohoku fault. We obtain a probability density for the slip distribution, which varies both with depth, earthquake size and whether the rupture breaks the surface. We propose a method to modify stochastic slip distributions according to this dynamically-derived probability distribution. This method may be efficiently applied to produce large numbers of heterogeneous slip distributions for probabilistic tsunami hazard analysis. Using numerous M9 earthquake scenarios, we demonstrate that incorporating the dynamically-derived probability distribution does enhance the conditional probability of exceedance of maximum estimated tsunami wave heights along the Japanese coast. This technique for integrating dynamic features in stochastic models can be extended to any subduction zone and faulting style.

Campi Flegrei is a densely populated volcanic area in Italy. Its inner caldera periodically experiences uplift and subsidence, known as bradyseism, also accompanied by seismic activity. In the last decade, with uplift rates up to 2 cm/month, about nine-thousand earthquakes were recorded. Upon request of the local administration, the most updated data were collected and analyzed to evaluate the risk management strategy consisting of structural retrofitting according to the building code. Here it is shown that the reference moment magnitude is in the range \(\left({{\mathrm{4.4,5.1}}}\right)\), based on fault mapping, geomorphological inference, earthquake relocation, stress-drop analysis, and ground motion modelling. Earthquake forecasting enabled computing the exceedance probabilities of these magnitudes. Earthquake engineering showed that the minimum magnitudes expected to cause exceedance of design ground motion, are larger than the reference magnitudes. Finally, the risk reduction implied by the safety levels of new constructions was assessed for reinforced concrete buildings.

SUMMARY: Understanding mechanical processes occurring on faults requires detailed information on the microseismicity that can be enhanced today by advanced techniques for earthquake detection. This problem is challenging when the seismicity rate is low and most of the earthquakes occur at depth. In this study, we compare three detection techniques, the autocorrelation FAST, the machine learning EQTransformer, and the template matching EQCorrScan, to assess their ability to improve catalogues associated with seismic sequences in the normal fault system of Southern Apennines (Italy) using data from the Irpinia Near Fault Observatory (INFO). We found that the integration of the machine learning and template matching detectors, the former providing templates for the cross-correlation, largely outperforms techniques based on autocorrelation and machine learning alone, featuring an enrichment of the automatic and manual catalogues of factors 21 and 7, respectively. Since output catalogues can be polluted by many false positives, we applied refined event selection based on the cumulative distribution of their similarity level. We can thus clean up the detection lists and analyse final subsets dominated by real events. The magnitude of completeness decreases by more than one unit compared to the reference value for the network. We report b-values associated with sequences smaller than the average, likely corresponding to larger differential stresses than for the background seismicity of the area. For all the analysed sequences, we found that main events are anticipated by foreshocks, indicating a possible preparation process for main shocks at subkilometric scales.