Regional moment tensor uncertainty due to mismodeling of the crust

The retrieval of earthquake moment tensor (MT) requires the response of the medium, in which seismic waves travel from the hypocenter to the stations, to be known. In inverting long-period (LP) seismic data (teleseismic and LP regional records), a gross earth model is sufficient; with decreasing periods, a more detailed model is needed. This is the case when waveforms of weak earthquakes at regional distances are to be inverted. Regional moment tensors (RMTs) of mostly Mediterranean earthquakes are determined on a routine basis by the Swiss Seismological Survey (SED) by using averaged models of the earth's crust. By inverting broad-band records of the M_w=4.8 earthquake near Udine, N Italy, on Feb. 14, 2002, we tested the sensitivity of the MT solution with respect to possible errors in the earth model used and in the location of the hypocenter depth. We perturbed the P and S velocities and the thickness in the 1-D earth model in the range from 3% to 30% of the parameter values and constructed estimates of confidence regions of the MT and error bars of the source time function (STF) and scalar moment in three frequency bands. Similarly, these error characteristics were determined assuming a mislocation in the hypocenter depth. We found that, in the band of periods from 25 to 50 s, the mechanism is resolved well (at the confidence level 95% at least) up to an earth model uncertainty of 30%, in the passband 10–25 s up to about 10%, but it is undetermined completely at periods of 5–10 s. An error in hypocenter depth of as much as double the value reported by the location procedure does not destroy the resolution of the mechanism at periods above 10 s. In the RMT catalog of the SED, earthquakes of M_w greater than about 3.5 are processed at periods above 30 s; thus, the solutions for these events are robust with respect to a possible uncertainty in the earth model used. Mechanisms of weaker earthquakes, retrieved from short periods, should be interpreted with caution.     

A complex, Subplinian-type eruption from low-viscosity, phonolitic to tephri-phonolitic magma: the AD 472 (Pollena) eruption of Somma-Vesuvius, Italy

The combined use of field investigation and laboratory analyses allowed the detailed stratigraphic reconstruction of the Pollena eruption (472 AD) of Somma-Vesuvius. Three main eruptive phases were recognized, related either to changes in the eruptive processes and/or to relative changes of melt composition. The eruption shows a pulsating behavior with deposition of pyroclastic fall beds and generation of dilute and dense pyroclastic density currents (PDC). The eruptive mechanisms and transportation dynamics were reconstructed for the whole eruption. Column heights were between 12 and 20 km, corresponding to mass discharge rates (MDR) of 7×10⁶ kg/s and 3.4×10⁷ kg/s. Eruptive dynamics were driven by magmatic fragmentation of a phono-tephritic to tephri-phonolitic magma during Phases I and II, whereas phreatomagmatic fragmentation dominated Phase III. Magma composition varies between phonolitic and tephritic-phonolitic, with melt viscosity likely not in excess of 10³ Pa s. The volume of the pyroclastic fall deposits, calculated by using of proximal isopachs, is 0.44 km³. This increases to 1.38 km³ if ash volumes are extrapolated on a log thickness vs. square root area diagram using one distal isopach and column height.Editorial responsibility: R Cioni     

Kinetic energy–rainfall intensity relationship for Central Cebu, Philippines for soil erosion studies

In the study of soil erosion, specifically on detachment of soil particles by raindrop impact, kinetic energy is a commonly suggested indicator of the raindrop's ability to detach soil particles from the soil mass. Since direct measurement of kinetic energy requires sophisticated and costly instruments, the alternative approach is to estimate it from rainfall intensity. The present study aims at establishing a relationship between rainfall intensity and kinetic energy for rainfalls in Central Cebu, Philippines as a preface of a wider regional investigation.

Drop size distributions of rainfalls were measured using the disdrometer RD-80. There are two forms of kinetic energy considered here. One is kinetic energy per unit area per unit time (KE_R, J m⁻² h⁻¹) and the other is kinetic energy per unit area per unit depth (KE, J m⁻² mm⁻¹). Relationships between kinetic energy per unit area per unit time (KE_R) and rainfall intensity (I) were obtained using linear and power relations. The exponential model and the logarithmic model were fitted to the KE–I data to obtain corresponding relationships between kinetic energy per unit area per unit depth of rainfall (KE) and rainfall intensity (I). The equation obtained from the exponential model produced smaller standard error of estimates than the logarithmic model.     

Evaluation of current statistical approaches for predictive geomorphological mapping

Predictive models are increasingly used in geomorphology, but systematic evaluations of novel statistical techniques are still limited. The aim of this study was to compare the accuracy of generalized linear models (GLM), generalized additive models (GAM), classification tree analysis (CTA), neural networks (ANN) and multiple adaptive regression splines (MARS) in predictive geomorphological modelling. Five different distribution models both for non-sorted and sorted patterned ground were constructed on the basis of four terrain parameters and four soil variables. To evaluate the models, the original data set of 9997 squares of 1 ha in size was randomly divided into model training (70%, n=6998) and model evaluation sets (30%, n=2999).In general, active sorted patterned ground is clearly defined in upper fell areas with high slope angle and till soils. Active non-sorted patterned ground is more common in valleys with higher soil moisture and fine-scale concave topography. The predictive performance of each model was evaluated using the area under the receiver operating characteristic curve (AUC) and the Kappa value. The relatively high discrimination capacity of all models, AUC=0.85–0.88 and Kappa=0.49–0.56, implies that the model's predictions provide an acceptable index of sorted and non-sorted patterned ground occurrence. The best performance for model calibration data for both data sets was achieved by the CTA. However, when the predictive mapping ability was explored through the evaluation data set, the model accuracies of CTA decreased clearly compared to the other modelling techniques. For model evaluation data MARS performed marginally best.Our results show that the digital elevation model and soil data can be used to predict relatively robustly the activity of patterned ground in fine scale in a subarctic landscape. This indicates that predictive geomorphological modelling has the advantage of providing relevant and useful information on earth surface processes over extensive areas, such data being unavailable through more conventional survey methods.     

Crustal structure of the Lofoten–Vesterålen continental margin, off Norway

By compiling wide-angle seismic velocity profiles along the 400-km-long Lofoten–Vesterålen continental margin off Norway, and integrating them with an extensive seismic reflection data set and crustal-scale two-dimensional gravity modelling, we outline the crustal margin structure. The structure is illustrated by across-margin regional transects and by contour maps of depth to Moho, thickness of the crystalline crust, and thickness of the 7+ km/s lower crustal body. The data reveal a normal thickness oceanic crust seaward of anomaly 23 and an increase in thickness towards the continent–ocean boundary associated with breakup magmatism. The southern boundary of the Lofoten–Vesterålen margin, the Bivrost Fracture Zone and its landward prolongation, appears as a major across-margin magmatic and structural crustal feature that governed the evolution of the margin. In particular, a steeply dipping and relatively narrow, 10–40-km-wide, Moho-gradient zone exists within a continent–ocean transition, which decreases in width northward along the Lofoten–Vesterålen margin. To the south, the zone continues along the Vøring margin, however it is offset 70–80 km to the northwest along the Bivrost Fracture Zone/Lineament. Here, the Moho-gradient zone corresponds to a distinct, 25-km-wide, zone of rapid landward increase in crustal thickness that defines the transition between the Lofoten platform and the Vøring Basin. The continental crust on the Lofoten–Vesterålen margin reaches a thickness of 26 km and appears to have experienced only moderate extension, contrasting with the greatly extended crust in the Vøring Basin farther south. There are also distinct differences between the Lofoten and Vesterålen margin segments as revealed by changes in structural style and crustal thickness as well as in the extent of elongate potential-field anomalies. These changes may be related to transfer zones. Gravity modelling shows that the prominent belt of shelf-edge gravity anomalies results from a shallow basement structural relief, while the elongate Lofoten Islands belt requires increased lower crustal densities along the entire area of crustal thinning beneath the islands. Furthermore, gravity modelling offers a robust diagnostic tool for the existence of the lower crustal body. From modelling results and previous studies on- and off-shore mid-Norway, we postulate that the development of a core complex in the middle to lower crust in the Lofoten Islands region, which has been exhumed along detachments during large-scale extension, brought high-grade, lower crustal rocks, possibly including accreted decompressional melts, to shallower levels.