Robust inversion of vertical electrical sounding data using a multiple reweighted least-squares method

The root cause of the instability problem of the least-squares (LS) solution of the resistivity inverse problem is the ill-conditioning of the sensitivity matrix. To circumvent this problem a new LS approach has been investigated in this paper. At each iteration, the sensitivity matrix is weighted in multiple ways generating a set of systems of linear equations. By solving each system, several candidate models are obtained. As a consequence, the space of models is explored in a more extensive and effective way resulting in a more robust and stable LS approach to solving the resistivity inverse problem. This new approach is called the multiple reweighted LS method (MRLS). The problems encountered when using the L ₁- or L ₂-norm are discussed and the advantages of working with the MRLS method are highlighted. A five-layer earth model which generates an ill-conditioned matrix due to equivalence is used to generate a synthetic data set for the Schlumberger configuration. The data are randomly corrupted by noise and then inverted by using L ₂, L ₁ and the MRLS algorithm. The stabilized solutions, even though blurred, could only be obtained by using a heavy ridge regression parameter in L ₂- and L ₁-norms. On the other hand, the MRLS solution is stable without regression factors and is superior and clearer. For a better appraisal the same initial model was used in all cases. The MRLS algorithm is also demonstrated for a field data set: a stable solution is obtained.     

Double-difference tomography of P- and S-wave velocity structure beneath the western part of Java,Indonesia

West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults. In this study, double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java. To infer the geometry of the structure beneath the study area, precise earthquake hypocenter determination was first performed before tomographic imaging. For this, earthquake waveform data were extracted from the regional Meteorological, Climatological, Geophysical Agency (BMKG) network of Indonesia from South Sumatra to Central Java. The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected, the key features being events of magnitude > 3, azimuthal gap < 210° and number of phases > 8. A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters. The hypocenter locations were then relocated using double-difference tomography (tomoDD). A significant reduction of travel-time (root mean square basis) and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java. Double-difference tomography was found to give a clear velocity structure, especially beneath the volcanic arc area, i.e., under Mt Anak Krakatau, Mt Salak and the mountains complex in the southern part of West Java. Low velocity anomalies for the P and S waves as well as the v_P/v_S ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.     

CFD simulation and experimental analysis of flow dynamics and grinding performance of opposed fluidized bed air jet mill

This paper presents a three dimensional Computational Fluid Dynamics (CFD) model to investigate the flow dynamics of solid–gas phases during fine grinding in an air jet mill. Alpine 100AFG fluidized bed air jet mill is considered for the study and the jet milling model is simulated using FLUENT 6.3.2 using a standard k-ε model. The model is developed in GAMBIT 2.3.16 and meshed by tet/hybrid (T-Grid) and Triangular (Pave) meshes. The effects of operating parameters such as solid feed rate, grinding air pressure and internal classifier speed on the performance of the jet mill are analyzed. The CFD simulation results are presented in the forms of dual phase vector plot, volume fraction of phases and particle trajectories during fine grinding process. The mass of ground feed entering and leaving the cyclone (underflow) is also computed by simulation. The proposed model gives realistic predictions of the flow dynamics within the jet mill. Experiments are conducted on the Alpine 100AFG jet mill to study the particle size, morphology and mass of the ground product. The numerical results are found in good agreement with the experimental results.     

The 2015 eruption of Gamalama volcano (Ternate Island–Indonesia): precursor,crisis management,and community response

GeoJournal - Gamalama is an active stratovolcano on Ternate, a small volcanic island in Maluku Utara, Indonesia. Since 1510, a total of 77 eruptions have been recorded, with various impacts on the...     

Multivariate fire risk models using copula regression in Kalimantan,Indonesia

Natural Hazards - Forest fires have become a national issue yearly and elicited serious attention from the government and researchers in Indonesia. Copula-based joint distribution can construct a...     

Landslide hazard and risk assessment and their application in risk management and landuse planning in eastern flank of Menoreh Mountains,Yogyakarta Province,Indonesia

The Menoreh Mountains in Yogyakarta are severely affected by landslides. Due to the high population densities, mass movements are generally damaging and fatal. More than other Javanese mountains, the Menoreh Mountains cumulate several factors causing landslides. Therefore, it is necessary to evaluate the ways to map landslide risk in order to improve the risk mitigation. The objectives of this paper are to provide landslide hazard and risk assessment that will be useful for risk prevention and landuse planning in the Menoreh Mountains. So far, risk management has been developed by the Research Centre for Disasters Gadjah Mada University in collaboration with the Regional Development Planner (BAPPEDA), which carries out fundamental and applied researches. The results of the studies have been integrated in the risk prevention and landuse planning in order to improve the integrated landslide mitigation programme.     

EM2INV — A finite difference based algorithm for two-dimensional inversion of geoelectromagnetic data

The paper presents an efficient finite difference based 2D-inversion algorithm, EM2INV, for geoelectromagnetic data. The special features of the algorithm are

Aquifer parameter estimation from surface resistivity data

This paper is devoted to the additional use, other than ground water exploration, of surface geoelectrical sounding data for aquifer hydraulic parameter estimation. In a mesoscopic framework, approximated analytical equations are developed separately for saline and for fresh water saturations. A few existing useful aquifer models, both for clean and shaley sandstones, are discussed in terms of their electrical and hydraulic effects, along with the linkage between the two. These equations are derived for insight and physical understanding of the phenomenon. In a macroscopic scale, a general aquifer model is proposed and analytical relations are derived for meaningful estimation, with a higher level of confidence, of hydraulic parameter from electrical parameters. The physical reasons for two different equations at the macroscopic level are explicitly explained to avoid confusion. Numerical examples from existing literature are reproduced to buttress our viewpoint.