Adjacency selection in Markov Random Fields for high spatial resolution hyperspectral data

 Markov Random Fields, implemented for the analysis of remote sensing images, capture the natural spatial dependence between band wavelengths taken at each pixel, through a suitable adjacency relationship between pixels, to be defined a priori. In most cases several adjacency definitions seem viable and a model selection problem arises. A BIC-penalized Pseudo-Likelihood criterion is suggested which combines good distributional properties and computational feasibility for analysis of high spatial resolution hyperspectral images. Its performance is compared with that of the BIC-penalized Likelihood criterion for detecting spatial structures in a high spatial resolution hyperspectral image for the Lamar area in Yellowstone National Park. Received: 9 March 2001 / Accepted: 2 August 2001     

Fluid circulation at Stromboli volcano (Aeolian Islands, Italy) from self-potential and CO2 surveys

This work addresses the study of fluid circulation of the Stromboli island using a dense coverage of self-potential (SP) and soil CO₂ data. A marked difference exists between the northern flank and the other flanks of the island. The northern flank exhibits (1) a typical negative SP/altitude gradient not observed on the other flanks, and (2) higher levels of CO₂. The general SP pattern suggests that the northern flank is composed of porous layers through which vadose water flows down to a basal water table, in contrast to the other flanks where impermeable layers impede the vertical flow of vadose water. In the Sciara del Fuoco and Rina Grande–Le Schicciole landslide complexes, breccias of shallow gliding planes may constitute such impermeable layers whereas elsewhere, poorly permeable, fine-grained pyroclastites or altered lava flows may be present. This general model of the flanks also explains the main CO₂ patterns: concentration of CO₂ at the surface is high on the porous north flank and lower on the other flanks where impermeable layers can block the upward CO₂ flux. The active upper part of the island is underlain by a well-defined hydrothermal system bounded by short-wavelength negative SP anomalies and high peaks of CO₂. These boundaries coincide with faults limiting ancient collapses of calderas, craters and flank landslides. The hydrothermal system is not homogeneous but composed of three main subsystems and of a fourth minor one and is not centered on the active craters. The latter are located near its border. This divergence between the location of the active craters and the extent of the hydrothermal system suggests that the internal heat sources may not be limited to sources below the active craters. If the heat source strictly corresponds to intrusions at depth around the active conduits, the geometry of the hydrothermal subsystems must be strongly controlled by heterogeneities within the edifice such as craters, caldera walls or gliding planes of flank collapse, as suggested by the correspondence between SP–CO₂ anomalies and structural limits. The inner zone of the hydrothermal subsystems is characterized by positive SP anomalies, indicating upward movements of fluids, and by very low values of CO₂ emanation. This pattern suggests that the hydrothermal zone becomes self-sealed at depth, thus creating a barrier to the CO₂ flux. In this hypothesis, the observed hydrothermal system is a shallow one and it involves mostly convection of infiltrated meteoric water above the sealed zone. Finally, on the base of CO₂ degassing measurements, we present evidence for the presence of two regional faults, oriented N41° and N64°, and decoupled from the volcanic structures.     

Morphological methods and dynamic modelling in landslide hazard assessment of the Campania Apennine carbonate slope

Landslide risk of the Campanian carbonate slopes covered by pyroclastic deposits is mainly connected with the occurrence of high-velocity debris avalanches and debris flows. Analyses show that flows initiate as small translational slides in the pyroclastics. The failure process is controlled by the interaction of both natural and human-induced factors. Geomorphological settings play a decisive role in locating the source failures. Therefore, the crucial aspects in landslide hazard and risk assessment are: (a) recognise the geomorphological control factors, (b) determine parameters defining landslide intensity (velocity, volume, depth of deposit) and (c) predict landslide runout distance. An approach combining geomorphology and numerical analysis has been adopted in the work reported here. Potential future landslide intensity scenarios are simulated predicting the runout behaviour of potential instabilities by using a dynamic model previously calibrated by back-analysing observed events of similar scale and type. The selected area is a sector of the Avella Mountains having the same geomorphological environment as the 1998 Sarno landslides (Campania, Southern Italy).     

Spatial Relations between Classes as Integrity Constraints

The quality and integrity of spatial data is very important to support interoperability among different systems. To reach this aim integrity rules defined by the application play an important role (for example, constraints between object classes). In this article, we propose a methodology to define integrity constraints using user level spatial relations between classes of individuals. We will also provide mapping rules from user level relations to geometric level operators to allow the computation of relations. As a case study, we will define the constraints for the class of rivers and some of its specializations.     

Is the Linné impact crater morphology influenced by the rheological layering on the Moon's surface? Insights from numerical modeling

Linné is a simple crater, with a diameter of 2.23 km and a depth of 0.52 km, located in northwestern Mare Serenitatis. Recent high‐resolution data acquired by the Lunar Reconnaissance Orbiter Camera revealed that the shape of this impact structure is best described by an inverted truncated‐cone. We perform morphometric measurements, including slope and profile curvature, on the Digital Terrain Model of Linné, finding the possible presence of three subtle topographic steps, at the elevation of +20, ?100, and ?200 m relative to the target surface. The kink at ?100 m might be related to the interface between two different rheological layers. Using the iSALE shock physics code, we numerically model the formation of Linné crater to derive hints on the possible impact conditions and target physical properties. In the initial setup, we adopt a basaltic projectile impacting the Moon with a speed of 18 km s^?1. For the local surface, we consider either one or two layers, in order to test the influence of material properties or composite rheologies on the final crater morphology. The one‐layer model shows that the largest variations in the crater shape take place when either the cohesion or the friction coefficient is varied. In particular, a cohesion of 10 kPa marks the threshold between conical‐ and parabolic‐shaped craters. The two‐layer model shows that the interface between the two layers would be exposed at the observed depth of 100 m when an intermediate value (~200 m) for the upper fractured layer is set. We have also found that the truncated‐cone morphology of Linné might originate from an incomplete collapse of the crater wall, as the breccia lens remains clustered along the crater walls, while the high‐albedo deposit on the crater floor can be interpreted as a very shallow lens of fallout breccia. The modeling analysis allows us to derive important clues on the impactor size (under the assumption of a vertical impact and collision velocity equal to the mean value), and on the approximate, large‐scale preimpact target properties. Observations suggest that these large‐scale material properties likely include some important smaller scale variations, disclosed as subtle morphological steps in the crater walls. Furthermore, the modeling results allow advancing some hypotheses on the geological evolution of the Mare Serenitatis region where Linné crater is located (unit S14). We suggest that unit S14 has a thickness of at least a few hundreds of meters up to about 400 m.     

Assessment of the dynamic response of unreinforced masonry structures using a macroelement modeling approach

The seismic performance of unreinforced masonry structures is strongly associated with the interaction between in‐plane and out‐of‐plane mechanisms. The seismic response of these structures has been thoroughly investigated by means of experimental testing, analytical procedures, and computational approaches. Within the framework of the numerical simulations, models based on the finite element method provide a good prediction of the seismic performance of unreinforced masonry structures. However, they usually require a high computational cost and advanced user expertise to define appropriate mechanical properties and to interpret the numerical results. Because of these limitations, simplified models for practical applications have been developed during the last decades. Despite this, a great number of these models focus mostly on the evaluation of the in‐plane response, assuming box (or integral) behavior of the structure. In this paper, a simplified macroelement modeling approach is used to simulate the seismic response of 2 masonry prototypes taking into consideration the combined in‐plane and out‐of‐plane action. The numerical investigations were performed in the static and dynamic fields by using pushover analyses and nonlinear dynamic analyses respectively. The latter is a novel implementation of a model previously developed for static analysis. The results obtained from this study are in good agreement with those provided by a detailed nonlinear continuum FE approach, demonstrating the applicability of this macroelement model with a significant reduction of the computational cost.     

The use of the slope–area function to analyse process domains in complex badland landscapes

This paper explores the effectiveness of the widely-used functional relationship between drainage area (A in m²) and slope (S in m/m) to identify local process domains and aid interpretation of process interactions in a complex badland landscape. In order to perform this investigation, a series of sub-basins tributary to the Formone River in the Orcia catchment (central Italy) were selected as a suitable study area within which to explore our questions, given these basins' general representativeness of local terrain, the availability of a high resolution digital terrain model and previous extensive geomorphological research. Eroding basins containing both calanchi and landslides are common in the sub-humid badland landscape of central Italy, where field observation identifies a complex pattern of erosive processes associated with a history of uplift, despite which parts of the local landscape appear disconnected. Results reveal that the shape of all S–A curves (plotted using S data binned on log A) is comparable with that described in the literature, although sub-basins containing calanchi generally plot with higher S values than non-calanchi ones, except in the ‘fluvial’ section of the plots. Second, when viewed on total data (non-binned) S–A plots, landslide source area domains and calanchi domains are entirely coincident in all basins, supporting a cause–effect relationship. Additional plotting of the frequency characteristics of the raw data in a new way supports the interpretation that calanchi frequently initiate in landslide scars. In general though, although the S–A plots can contribute to the disentanglement of geomorphological behaviour in some complex erosional landscapes, it became apparent that in this landscape, process domains do not separate out with clarity along the A axis as suggested by theory. Despite this, an alternative, broader-scale morphoevolutive model can be proposed for the development of within-landslide calanchi, driven by changes to basin connectivity to the base channel. © 2018 John Wiley & Sons, Ltd.     

Neo-deterministic seismic hazard scenarios for India—a preventive tool for disaster mitigation

Current computational resources and physical knowledge of the seismic wave generation and propagation processes allow for reliable numerical and analytical models of waveform generation and propagation. From the simulation of ground motion, it is easy to extract the desired earthquake hazard parameters. Accordingly, a scenario-based approach to seismic hazard assessment has been developed, namely the neo-deterministic seismic hazard assessment (NDSHA), which allows for a wide range of possible seismic sources to be used in the definition of reliable scenarios by means of realistic waveforms modelling. Such reliable and comprehensive characterization of expected earthquake ground motion is essential to improve building codes, particularly for the protection of critical infrastructures and for land use planning. Parvez et al. (Geophys J Int 155:489–508, 2003) published the first ever neo-deterministic seismic hazard map of India by computing synthetic seismograms with input data set consisting of structural models, seismogenic zones, focal mechanisms and earthquake catalogues. As described in Panza et al. (Adv Geophys 53:93–165, 2012), the NDSHA methodology evolved with respect to the original formulation used by Parvez et al. (Geophys J Int 155:489–508, 2003): the computer codes were improved to better fit the need of producing realistic ground shaking maps and ground shaking scenarios, at different scale levels, exploiting the most significant pertinent progresses in data acquisition and modelling. Accordingly, the present study supplies a revised NDSHA map for India. The seismic hazard, expressed in terms of maximum displacement (Dmax), maximum velocity (Vmax) and design ground acceleration (DGA), has been extracted from the synthetic signals and mapped on a regular grid over the studied territory.     

Improved evaluation of inelastic displacement demands for short‐period masonry structures

This work discusses the simplified estimation of earthquake‐induced nonlinear displacement demands as required by nonlinear static procedures, with particular attention on short‐period masonry structures. The study focuses on systems with fundamental periods between 0.1 and 0.5 s, for which inelastic amplification of the elastic displacement demand is more pronounced; hysteretic force‐displacement relationships characteristic of masonry structures are adopted, because these structures are more commonly found within the considered period range. Referring to the results of nonlinear dynamic analyses of single‐degree‐of‐freedom oscillators, some limitations of the Eurocode 8 and Italian Building Code formulations are first discussed, then an improved equation is calibrated that relates inelastic and elastic displacement demands. Numerical values of the equation parameters are obtained, considering the amount of hysteretic energy dissipation associated with various damage mechanisms observed in masonry structures. Safety factors are also calculated to determine several percentiles of the displacement demand. It is shown that the proposed equation can be extended to more dissipative systems. Finally, the same formulation is adapted to the estimation of seismic displacements when elastic analysis procedures are employed. Copyright © 2017 John Wiley & Sons, Ltd.