Temporal variations in the spatial distribution of shoot density in a Posidonia oceanica meadow and patterns of genetic diversity

The meadows of the Mediterranean endemic seagrass Posidonia oceanica exhibit relatively high variations of structural and biometric features at various spatial scales. An investigation performed in 1992 in the meadow off Lacco Ameno (Island of Ischia, Gulf of Naples, Italy) detected peculiar spatial patterns of plant distribution, characterized by nestlike structures with radially increasing or decreasing shoot densities. Eight years later (2000), geo‐referenced collections at selected points were repeated to trace the temporal variations of shoot density and investigate the recurrence of the density cores previously detected. In addition, shoots for molecular analyses were collected to check the hypothesis that nestlike patterns exhibit highest levels of genetic variability, due to the confluence of several genetically distinct stolons. The 2000 survey confirmed the presence of the main density cores detected in 1992, although their spatial distribution was slightly shifted and a general decrease of spatial anisotropy was observed, probably due to an increased disturbance, mainly due to pleasure boat anchoring. Patterns of genetic diversity showed a more complex picture, well related to the shoot density spatial pattern, especially when compared with the previous 1992 survey. Patterns of genetic diversity confirmed our previous hypotheses on the genesis of shoot density cores, suggesting they are produced over long time, due to a slow stolonization process and a convergence of different genotypes. Regression of the meadow and decrease of density may lead, in short periods, to a homogenization of the density patterns, while genetic diversity cores represent a long‐term ‘memory’ of their previous distribution.     

On using Boussinesq-type equations near the shoreline: a note of caution

We briefly analyze some characteristics of the behavior in very shallow waters i.e. near the shoreline of high-order (dispersive-nonlinear) Boussinesq-type equations. By using the Carrier and Greenspan (1958) solution as test flow conditions we illustrate the behavior of both purely dispersive and dispersive-nonlinear contributions near the shoreline. It is also shown that Boussinesq-type equations can be more usefully handled in the swash zone if written in terms of the total water depth.     

Decorrugation and removal of directional trends of magnetic fields by the wavelet transform: application to archaeological areas

The problem of removing directional trends frequently occurs in the processing of magnetic data and also in the subsequent steps of data interpretation. The so-called corrugations are typical directional trends occurring in levelled data, which may be removed in several ways. Classical techniques are based on high-pass filtering of the data and successively filtering these transformed data with directional cosine filters. Other linear features are due to real sources, such as pipelines in shallow surveys or dike swarms in regional surveys. They should, nevertheless, be considered as noise, due to the fact that their effect is strong and tends to hide the field features related to structures of more interest. We deal with both kinds of problem, presenting the results of a study in an archaeological area of southern Italy. Decorrugation of magnetic field anomalies is performed using a method based on the excellent space–frequency localization properties of wavelet bases, allowing a very sharp filtering of the field along a selected direction. We compare this technique with the classical one in a synthetic case and find that the wavelet decorrugation is simpler and produces low distortion maps. Besides the field decorrugation, the wavelet approach was also shown to be useful in the subsequent enhancement of the measured field. In fact, we show that the wavelet analysis offers a unique framework where various filtering problems (directional, isotropic, global or local as well) may be easily solved. As regards the archaeological case, strong noisy effects from elongated sources (pipelines) were successfully removed in a sharp and local way.     

Modelling gravity variations consistent with ground deformation in the Campi Flegrei caldera (Italy)

In recent years (1970–72 and 1982–84) two inflation episodes took place in the Campi Flegrei caldera (Italy), characterized by significant ground uplift and gravity variations. An elastic half-space model with vertical density stratification is employed to compute the displacement field and the gravity variations produced by the deformation of buried layers, following the inflation of a spherically symmetric deformation source. Contributions to gravity variations are produced by dilation/contraction of the medium, by the displacements of density interfaces (the free surface and subsurface layers) and of source boundaries and, possibly, by new mass input from remote distances into the source volume. Three cases were examined in detail: In case I, the magma chamber is identified as the deformation source and volume and pressure increase in the magma chamber is due to input of new magma from remote distances; in case II deformation is due to magma differentiation within the magma chamber (deformation source with constant mass); in case III the geothermal system is identified as the deformation source and a pressure increase, possibly driven by the exsolution of high temperature and high pressure volatiles in the magma chamber, is assumed to play a dominant role. From the comparison between measured and computed gravity residuals (free-air-corrected gravity variations) we can assess that, in case I, an inflation source with constant density would predict gravity residuals compatible with observations, whereas an expansion at constant mass (case II) would predict gravity residuals much lower than observed. The resolving power of gravity data however prevents accurate assessment of the density of the emplaced material. In case III, the pervasive density increase of the geothermal fluids induced by pressure increase is assumed to be the main source of gravity variations. The average porosity value required for this model to match both the ground deformation and the gravity residuals is found to be ˜10%, a value which is compatible with measured porosity values at Campi Flegrei in deep wells. The subsidence phases following both inflation episodes and the gravity residuals during subsidence lead us to consider case III as more plausible, even if a suitable combination of cases I and III cannot be discarded.     

Thermal, seismic and infrasonic evidences of variable degassing rates at Stromboli volcano

During June 1999, we measured the amplitude and rate (number of events per second) at which gas exited the vent at Stromboli volcano as discrete gas bursts or puffs. This allowed us to identify two styles of gas burst (puffing) activity. The first is characterized by frequent, rapidly rising puffs, the second by less frequent, slowly rising puffs. Each style persisted over 5–40-min-long durations and was associated with a high and low number of strombolian explosions per hour, respectively. Each period was also associated with characteristic delay times between the arrival of the infrasonic and thermal signals during strombolian explosions; the delays were longer during vigorous puffing periods. To explain our observations, we propose a model in which the degassing process cycles between vigorous and weak degassing phases. During vigorous degassing phases, bubble layers ascend the conduit at a frequency of 0.5–1.0 s⁻¹. This high degassing level reflects a gas-rich magma column and leads to an increased rate in the formation of shallow foams and, hence, an increase in puffing and explosive activity, as well as a higher free surface level and/or gas jet velocity. During weak phases, bubble layers ascend the conduit at a reduced frequency of 0.2–0.3 s⁻¹. During such times the magma column is poor in gas. This leads to a decreased rate of foam layer formation and hence a reduction in puffing and explosive activity, as well as a lower free surface level and/or gas jet velocity. Variations in puffing activity can thus be used to track changes in the rate at which the shallow system is supplied by fresh, gas-rich magma. Our observations indicate that the two degassing styles last from 5 to 40 min and that the switch from one to the other occurs over a matter of minutes.     

Optimal energy‐based seismic design of non‐conventional Tuned Mass Damper (TMD) implemented via inter‐story isolation

Inter‐story isolation, an effective strategy for mitigating the seismic risk of both new and existing buildings, has gained more and more interest in recent years as alternative to base isolation, whenever the latter results to be impractical, technically difficult or uneconomic. As suggested by the name, the technique consists in inserting flexible isolators at floor levels other than the base along the height of a multi‐story building, thus realizing a non‐conventional Tuned Mass Damper (TMD). Consistent with this, an optimal design methodology is developed in the present paper with the objective of achieving the global protection of both the structural portions separated by the inter‐story isolation system, that is, the lower portion (below the isolation system) and the isolated upper portion (above the isolation system). The optimization procedure is formulated on the basis of an energy performance criterion that consists in maximizing the ratio between the energy dissipated in the isolation system and the input energy globally transferred to the entire structure. Numerical simulations, performed under natural accelerograms with different frequency content and considering increasing isolation levels along the height of a reference frame structure, are used to investigate the seismic performance of the optimized inter‐story isolation systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Inversion of potential field data using the structural index as weighting function rate decay

Nonparametric inverse methods provide a general framework for solving potential‐field problems. The use of weighted norms leads to a general regularization problem of Tikhonov form. We present an alternative procedure to estimate the source susceptibility distribution from potential field measurements exploiting inversion methods by means of a flexible depth‐weighting function in the Tikhonov formulation. Our approach improves the formulation proposed by Li and Oldenburg (1996, 1998) , differing significantly in the definition of the depth‐weighting function. In our formalism the depth weighting function is associated not to the field decay of a single block (which can be representative of just a part of the source) but to the field decay of the whole source, thus implying that the data inversion is independent on the cell shape. So, in our procedure, the depth‐weighting function is not given with a fixed exponent but with the structural index N of the source as the exponent. Differently than previous methods, our choice gives a substantial objectivity to the form of the depth‐weighting function and to the consequent solutions. The allowed values for the exponent of the depth‐weighting function depend on the range of N for sources: 0 ≤N≤ 3 (magnetic case). The analysis regarding the cases of simple sources such as dipoles, dipole lines, dykes or contacts, validate our hypothesis. The study of a complex synthetic case also proves that the depth‐weighting decay cannot be necessarily assumed as equal to 3. Moreover it should not be kept constant for multi‐source models but should instead depend on the structural indices of the different sources. In this way we are able to successfully invert the magnetic data of the Vulture area, Southern Italy. An original aspect of the proposed inversion scheme is that it brings an explicit link between two widely used types of interpretation methods, namely those assuming homogeneous fields, such as Euler deconvolution or depth from extreme points transformation and the inversion under the Tikhonov‐form including a depth‐weighting function. The availability of further constraints, from drillings or known geology, will definitely improve the quality of the solution.     

Toward a full multiscale approach to interpret potential fields

The way potential fields convey source information depends on the scale at which the field is analysed. In this sense a multiscale analysis is a useful method to study potential fields particularly when the main field contributions are caused by sources with different depths and extents. Our multiscale approach is built with a stable transformation, such as depth from extreme points. Its stability results from mixing, in a single operator, the wavenumber low‐pass behaviour of the upward continuation transformation of the field with the enhancement high‐pass properties of n‐order derivative transformations. So, the complex reciprocal interference of several field components may be efficiently faced at several scales of the analysis and the depth to the sources may be estimated together with the homogeneity degrees of the field. In order to estimate the source boundaries we use another multiscale method, the multiscale derivative analysis, which utilizes a generalized concept of horizontal derivative and produces a set of boundary maps at different scales. We show through synthetic examples and application to the gravity field of Southern Italy that this multiscale behaviour makes this technique quite different from other source boundary estimators. The main result obtained by integrating multiscale derivative analysis with depth from extreme points is the retrieval of rather effective information of the field sources (horizontal boundaries, depth, structural index). This interpretative approach has been used along a specific transect for the analysis of the Bouguer anomaly field of Southern Apennines. It was set at such scales, so to emphasize either regional or local features along the transect. Two different classes of sources were individuated. The first one includes a broad, deep source with lateral size of 45∼50 km, at a depth of 13 km and having a 0.5 structural index. The second class includes several narrower sources located at shallowest depths, ranging from 3–6 km, with lateral size not larger than 5 km and structural indexes ranging from 1–1.5. Within a large‐scale geological framework, these results could help to outline the mean structural features at crustal depths.