Recently, ring laser technology has provided the first consistent observations of rotational ground motions around a vertical axis induced by earthquakes. "Consistent," in this context, implies that the observed waveforms and amplitudes are compatible with collocated recordings of translational ground motions. In particular, transverse accelerations should be in phase with rotation rate and their ratio proportional to local horizontal phase velocity assuming plane-wave propagation. The ring laser installed at the Fundamentalstation Wettzell in the Bavarian Forest, Southeast Germany, is recording the rotation rate around a vertical axis, theoretically a linear combination of the space derivatives of the horizontal components of motion. This suggests that, in principle, rotation can be derived from seismic-array experiments by "finite differencing." This has been attempted previously in several studies; however, the accuracy of these observations could never be tested in the absence of direct measurements. We installed a double cross-shaped array of nine stations from December 2003 to March 2004 around the ring laser instrument and observed several large earthquakes on both the ring laser and the seismic array. Here we present for the first time a comparison of array-derived rotations with direct measurements of rotations for ground motions induced by the M 6.3 Al Hoceima, Morocco, earthquake of 24 February 2004. With complete 3D synthetic seismograms calculated for this event we show that even low levels of noise may considerably influence the accuracy of the array-derived rotations when the minimum number of required stations (three) is used. Nevertheless, when using all nine stations, the overall fit between direct and array-derived measurements is surprisingly good (maximum correlation coefficient of 0.94).

No abstract

No abstract

Accurate precipitation observations are crucial for water resources management and as inputs for a gamut of hydrometeorological applications. Precipitation data from Integrated Multi-Satellite Retrievals for Global Precipitation Measurement (GPM) (IMERG) have recently been widely used to complement traditional rain gauge systems. However, the satellite precipitation data needs to be validated before being widely used in the applications and this is still missing over the Indonesian maritime continent (IMC). We conducted a validation of the IMERG product version 6 for this region. The evaluation was carried out using gauge data in the period from 2016 to 2020 for three types of IMERG: Early (E), Late (L), and Final (F) from annual, monthly, daily and hourly data. In general, the annual and monthly data from IMERG showed a good correlation with the rain gauge, with the mean correlation coefficient (CC) approximately 0.54–0.78 and 0.62–0.79, respectively. About 80% of stations in the IMC area showed a very good correlation between gauge data and IMERG-F estimates (CC = 0.7–0.9). For the daily assessment, the CC value was in the range of 0.39 to 0.44 and about 40% of stations had a correlation of 0.5–0.7. IMERG had a fairly good ability to detect daily rain in which the average probability of detection (POD) for all stations was above 0.8. However, the false alarm ratio (FAR) value is quite high (<0.5). For hourly data, IMERG’s performance was still poor with CC around 0.03–0.28. For all assessments, IMERG generally overestimated rainfall in comparison with rain gauge. The accuracy of the three types of IMERG in IMC was also influenced by season and topography. The highest and lowest CC values were observed for June–July–August and December–January–February, respectively. However, categorical statistics (POD, FAR and critical success index) did not show any clear seasonal variation. The CC value decreased with higher altitude, but with slight difference for each IMERG type. For all assessments conducted, IMERG-F generally showed the best rainfall observations in IMC, but with slightly difference from IMERG-E and IMERG-L. Thus, IMERG-E and IMERG-L data that had a faster latency than IMERG-F show potential to be used in rainfall observations in IMC.

No abstract

No abstract

Identification of internal structures in an active volcano is mandatory to quantify the physical processes preceding eruptions. We propose a fully unsupervised Bayesian inversion method that uses the point compound dislocation model as a complex source of deformation, to dynamically identify the substructures activated during magma migration. We applied this method at Piton de la Fournaise. Using 7‐day moving trends of Global Navigation Satellite System (GNSS) data preceding the June 2014 eruption, we compute a total of 15 inversion models of 2.5 million forward problems each, without a priori information. Obtained source shapes (dikes, prolate ellipsoids, or pipes) show magma migration from 7–8 km depth to the surface, drawing a mechanical “tomography” of the magma pathway. Our results also suggest source geometries compatible with observed eruptive fissures and seismicity distribution. In case of finite magma volume involved in final dike injection, source volume estimates from this method allow forecasting volumes of erupted lava.

On 2006 May 26 at 23:54 UTC, a moderate shallow crustal earthquake with a moment magnitude of 6.3 occurred in the southern part of Yogyakarta in Java, Indonesia. The earthquake caused severe damages in the area in addition to over 5700 casualties. The cause of this earthquake was initially believed to have been a rupture on the northeast-southwest trending Opak Fault; however, the role of this fault in the earthquake continues to be debated. Therefore, this study presents a subsurface model constructed to characterize the fault geometry associated with the earthquake. We used previously reported aftershock data to image subsurface velocity variations through seismic tomographic inversion of primary waves, shear waves and their velocity ratio (V p /V s). Using data from 10 stations around the hypothetical fault, 588 aftershock events were mostly located 10-15 km east of the Opak River Fault with a maximum depth of approximately 20 km. The seismic tomographic inversion results indicated that severe damage during the earthquake occurred in areas with larger V p /V s ratios associated with unconsolidated sediments, in accordance with previous findings. Furthermore, the configuration of an unnamed fault that was activated during the earthquake is delineated by a velocity anomaly with a depth of up to 5-7 km. This structure is interpreted as a strike-slip fault with a reverse component dipping to the east, striking northeast-southwest.