A hybrid optimization scheme for self-potential measurements due to multiple sheet-like bodies in arbitrary 2D resistivity distributions

Self-potential is a passive geophysical method that can be applied in a straightforward manner with minimum requirements in the field. Nonetheless, interpretation of self-potential data is particularly challenging due to the inherited non-uniqueness present in all potential methods. Incorporating information regarding the target of interest can facilitate interpretation and increase the reliability of the final output. In the current paper, a novel method for detecting multiple sheet-like targets is presented. A numerical framework is initially described that simulates sheet-like bodies in an arbitrary 2D resistivity distribution. A scattered field formulation based on finite differences is employed that allows the edges of the sheet to be independent of the grid geometry. A novel analytical solution for two-layered models is derived and subsequently used to validate the accuracy of the proposed numerical scheme. Lastly, a hybrid optimization is proposed that couples linear least-squares with particle-swarm optimization in order to effectively locate the edges of multiple sheet-like bodies. Through numerical and real data, it is proven that the hybrid optimization overcomes local minimal that occurs in complex resistivity distributions and converges substantially faster compared to traditional particle-swarm optimization.     

A discriminatory diagram of massive versus stratified deposits based on the sedimentation and bedload transportation rates. Experimental investigation and application to pyroclastic density currents

Dense gas-particle jets similar to collapsing eruption columns were generated by large-scale experiments. The column collapse resulted in a ground-hugging current forming stratified layers with bedding similar to natural pyroclastic density current deposits. At the impact of the collapsing column on the ground, a thick, massive bed was formed due to a high sedimentation rate that dumped turbulence due to high clast concentration. Down-current, flow expansion favoured turbulence and dilute gas-particle current that formed thin rippled layers deposited under traction. Experiments fed with fine ash (median size 0·066 mm) formed deposits without tractional structures, because fine particles, as other sedimentary fine material, is cohesive and exposes a limited surface to the shear stress. Experimental outcomes show that massive beds are formed where the sedimentation rate per unit width S_rw exceeds the bedload transportation rate Q_b by two orders of magnitude. A lower ratio generates traction at the base of the flow and formation of shear structures that increase in wavelength and height with a decreasing flux. This study presents a diagram that provides a useful addition for facies analysis of pyroclastic density currents, provided that deposits representing sustained sedimentation can be identified in the field. In the diagram a decrease in the S_rw/Q_b ratio corresponds to an increase in bedform size. Application of the diagram for hazard assessment purposes allows the reconstruction of the mass eruption rate of the Agnano–Monte Spina eruption at Campi Flegrei, which is the main variable defining the intensity of past eruptions, and of the Bingham rheology of the massive underflow of the Mercato pyroclastic density current at Vesuvius.