Synthesis of a seismic virtual reflector*

We describe a method to process the seismic data generated by a plurality of sources and registered by an appropriate distribution of receivers, which provides new seismic signals as if in the position of the receivers (or sources) there was an ideal reflector, even if this reflector is not present there. The data provided by this method represent the signals of a virtual reflector. The proposed algorithm performs the convolution and the subsequent sum of the real traces without needing subsurface model information. The approach can be used in combination with seismic interferometry to separate wavefields and process the reflection events. The application is described with synthetic examples, including stationary phase analysis and with real data in which the virtual reflector signal can be appreciated.     

Forest Fire Alert System: a Geo Web GIS prioritization model considering land susceptibility and hotspots – a case study in the Carajás National Forest,Brazilian Amazon

To increase the monitoring potential of forest fires, an alert classification methodology using satellite-mapped hotspots has been established to help forest managers in the prioritization of which hotspot to be verified in the field, thus potentially improving the distribution of fire-fighting resources. A computer application was developed based on web-distributed geographical information technology whose main function is to interact automatically generated satellite hotspots and risk areas indicated in fire-susceptibility maps and classify them into five alert levels. The location of the hotspots is available continuously every 4 h, and a susceptibility map is produced daily through map algebra algorithm, which uses static (topography, vegetation and land use) and dynamic (weather) variables. Every process runs through automated geoprocessing routines. The methodology was tested during the dry period of 2007 in the Carajás National Forest, in the Brazilian Amazon, within an area of 400,000 ha. It is a critical area constantly threatened by fires caused by invasions and deforestation owing to intense agribusiness advances and mining activities in its surroundings. This situation results in observations of many hotspots inside the study area for the same day and almost the same time period, in places of extreme opposites, demanding complex rapid analysis and hindering the decision of the displacement of fire-fighting teams. Further, a major mining company operates within the National Forest area, maintaining actions of protection as part of its environmental mining license. Results are presented under three aspects: (i) the credibility of the daily susceptibility map (algorithm), which showed strong correlation between areas of greatest risks and the confirmed forest fires; (ii) the reliability of hotspots (alert levels), confirming 71% of fires; (iii) accuracy in the decision of which hotspot to be checked, which revealed the same number of verifications at different alert levels, 82% confirmed alert 5 hotspots (maximum) and only 50% from alert 1 (minimum), resulting in faster fire-fighting actions, minimizing burned areas and, in some cases, allowing fire control before its spreading. Therefore, the methodology demonstrated that GIS routines are able to determine the relationship between a reality-based, interpreted susceptibility map of the area and satellite-generated hotspots, highlighting the ones of highest hazard level through the alert classification, becoming an important tool to help decisions from the fire-control center, especially for high-risk regions. The methodology may be extrapolated to other forested areas.     

The SHARE European Earthquake Catalogue (SHEEC) 1000–1899

In the frame of the European Commission project “Seismic Hazard Harmonization in Europe” (SHARE), aiming at harmonizing seismic hazard at a European scale, the compilation of a homogeneous, European parametric earthquake catalogue was planned. The goal was to be achieved by considering the most updated historical dataset and assessing homogenous magnitudes, with support from several institutions. This paper describes the SHARE European Earthquake Catalogue (SHEEC), which covers the time window 1000–1899. It strongly relies on the experience of the European Commission project “Network of Research Infrastructures for European Seismology” (NERIES), a module of which was dedicated to create the European “Archive of Historical Earthquake Data” (AHEAD) and to establish methodologies to homogenously derive earthquake parameters from macroseismic data. AHEAD has supplied the final earthquake list, obtained after sorting duplications out and eliminating many fake events; in addition, it supplied the most updated historical dataset. Macroseismic data points (MDPs) provided by AHEAD have been processed with updated, repeatable procedures, regionally calibrated against a set of recent, instrumental earthquakes, to obtain earthquake parameters. From the same data, a set of epicentral intensity-to-magnitude relations has been derived, with the aim of providing another set of homogeneous Mw estimates. Then, a strategy focussed on maximizing the homogeneity of the final epicentral location and Mw, has been adopted. Special care has been devoted also to supply location and Mw uncertainty. The paper focuses on the procedure adopted for the compilation of SHEEC and briefly comments on the achieved results.     

Modeling the Permeability Loss of Metallic Iron Water Filtration Systems

Over the past 30 years the literature has burgeoned with in situ approaches for groundwater remediation. Of the methods currently available, the use of metallic iron (Fe⁰) in permeable reactive barrier (PRB) systems is one of the most commonly applied. Despite such interest, an increasing amount of experimental and field observations have reported inconsistent Fe⁰ barrier operation compared to contemporary theory. In the current work, a critical review of the physical chemistry of aqueous Fe⁰ corrosion in porous media is presented. Subsequent implications for the design of Fe⁰ filtration systems are modeled. The results suggest that: (i) for the pH range of natural waters (>4.5), the high volumetric expansion of Fe⁰ during oxidation and precipitation dictates that Fe⁰ should be mixed with a non‐expansive material; (ii) naturally occurring solute precipitates have a negligible impact on permeability loss compared to Fe⁰ expansive corrosion; and (iii) the proliferation of H₂ metabolizing bacteria may contribute to alleviate permeability loss. As a consequence, it is suggested that more emphasis must be placed on future work with regard to considering the Fe⁰ PRB system as a physical (size‐exclusion) water filter device.     

Volcanic ash detection by GPS signal

We investigate the ability of GPS to detect volcanic plumes at Mt. Etna, Italy. We use a robust statistical approach to highlight whether the presence of a volcanic plume in the atmosphere may really affect the GPS undifferenced post-fit phase residuals. The proposed method has been tested for the September 4–5, 2007 activity of Mt. Etna. This eruption produced powerful lava fountains forming a weak, a few kilometers high plume for several hours, representing typical activity at Etna over the last 5 years. We analyzed data from nineteen Etna permanent GPS stations located on the volcano flanks at different heights and applied a statistical test based on four main steps: (a) realization of a simplified model representing the volcanic plume in atmosphere; (b) evaluation of the GPS satellite and station couples intersecting the plume; (c) calculation of the volcanic plume region crossed by the GPS signal; (d) application of a robust statistical test in order to see whether the volcanic plume affected the GPS signals. Results show that during the September 4–5, 2007 explosive activity, the GPS residuals definitely include the contribution of the volcanic plume. Our analysis shows that values of the GPS residuals are ten times smaller than those found for the Miyakejima eruption (Japan), highlighting a likely relationship between residuals and eruption intensity. In the future, data derived from the GPS stations located on Etna’s flanks could be used to improve the alerting system of volcanic ash, already operating at the Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Etneo.     

Seismic interferometry experiment in a shallow cased borehole using a seismic vibrator source‡

We present the results of a seismic interferometry experiment in a shallow cased borehole. The experiment is an initial study for subsequent borehole seismic surveys in an instrumented well site, where we plan to test other surface/borehole seismic techniques. The purpose of this application is to improve the knowledge of the reflectivity sequence and to verify the potential of the seismic interferometry approach to retrieve high‐frequency signals in the single well geometry, overcoming the loss and attenuation effects introduced by the overburden. We used a walkaway vertical seismic profile (VSP) geometry with a seismic vibrator to generate polarized vertical and horizontal components along a surface seismic line and an array of 3C geophones cemented outside the casing. The recorded traces are processed to obtain virtual sources in the borehole and to simulate single‐well gathers with a variable source‐receiver offset in the vertical array. We compare the results obtained by processing the field data with synthetic signals calculated by numerical simulation and analyse the signal bandwidth and amplitude versus offset to evaluate near‐field effects in the virtual signals. The application provides direct and reflected signals with improved bandwidth after vibrator signal deconvolution. Clear reflections are detected in the virtual seismic sections in agreement with the geology and other surface and borehole seismic data recorded with conventional seismic exploration techniques.     

Holocene mass‐wasting events in Lago Fagnano,Tierra del Fuego (54°S): implications for paleoseismicity of the Magallanes‐Fagnano transform fault

High‐resolution seismic imaging and coring in Lago Fagnano, located along a plate boundary in Tierra del Fuego, have revealed a dated sequence of Holocene mass‐wasting events. These structures are interpreted as sediment mobilizations resulting from loading of the slope‐adjacent lake floor during mass‐flow deposition. More than 19 mass‐flow deposits have been identified, combining results from 800 km of gridded seismic profiles used to site sediment cores. Successions of up to 6‐m thick mass‐flow deposits, pond atop the basin floor and spread eastward and westward following the main axis of the eastern sub‐basin of Lago Fagnano. We developed an age model, on the basis of information from previous studies and from new AMS‐¹⁴C ages on cored sediments, which allows us to establish a well‐constrained chronologic mass‐wasting event‐catalogue covering the last ～12 000 years. Simultaneously triggered, basin‐wide lateral slope failure and the formation of multiple debris flow and postulated megaturbidite deposits are interpreted as the fingerprint of paleo‐seismic activity along the Magallanes‐Fagnano transform fault that runs along the entire lake basin. The slope failures and megaturbidites are interpreted as recording large earthquakes occurring along the transform fault since the early Holocene. The results from this study provide new data about the frequency and possible magnitude of Holocene earthquakes in Tierra del Fuego, which can be applied in the context of seismic hazard assessment in southernmost Patagonia.     

On the explicit determination of stability constants for linearized geodetic boundary value problems

The theory of GBVPs provide the basis to the approximate methods used to compute global gravity models. A standard approximation procedure is least squares, which implicitly assumes that data, e.g. gravity disturbance and gravity anomaly, are given functions in L ²(S). We know that solutions in these cases exist, but uniqueness (and coerciveness which implies stability of the numerical solutions) is the real problem. Conditions of uniqueness for the linearized fixed boundary and Molodensky problems are studied in detail. They depend on the geometry of the boundary; however, the case of linearized fixed boundary BVP puts practically no constraint on the surface S, while the linearized Molodensky BVP requires the previous knowledge of very low harmonics, for instance up to degree 25, if we want the telluroid to be free to have inclinations up to 60^°.     

Space-Wise approach for airborne gravity data modelling

Regional gravity field modelling by means of remove-compute-restore procedure is nowadays widely applied in different contexts: it is the most used technique for regional gravimetric geoid determination, and it is also used in exploration geophysics to predict grids of gravity anomalies (Bouguer, free-air, isostatic, etc.), which are useful to understand and map geological structures in a specific region. Considering this last application, due to the required accuracy and resolution, airborne gravity observations are usually adopted. However, due to the relatively high acquisition velocity, presence of atmospheric turbulence, aircraft vibration, instrumental drift, etc., airborne data are usually contaminated by a very high observation error. For this reason, a proper procedure to filter the raw observations in both the low and high frequencies should be applied to recover valuable information. In this work, a software to filter and grid raw airborne observations is presented: the proposed solution consists in a combination of an along-track Wiener filter and a classical Least Squares Collocation technique. Basically, the proposed procedure is an adaptation to airborne gravimetry of the Space-Wise approach, developed by Politecnico di Milano to process data coming from the ESA satellite mission GOCE. Among the main differences with respect to the satellite application of this approach, there is the fact that, while in processing GOCE data the stochastic characteristics of the observation error can be considered a-priori well known, in airborne gravimetry, due to the complex environment in which the observations are acquired, these characteristics are unknown and should be retrieved from the dataset itself. The presented solution is suited for airborne data analysis in order to be able to quickly filter and grid gravity observations in an easy way. Some innovative theoretical aspects focusing in particular on the theoretical covariance modelling are presented too. In the end, the goodness of the procedure is evaluated by means of a test on real data retrieving the gravitational signal with a predicted accuracy of about 0.4 mGal.     

Laboratory determination of mechanical and hydraulic properties of chemically grouted sands

The purpose of this paper is to investigate mechanical and hydraulic properties of sands treated with mineral-based grouts through the results of a laboratory test programme consisting of unconfined compression tests (UCS), triaxial bender element tests (BeT) and constant flow permeability tests in triaxial apparatus. An improved apparatus was set up for obtaining high quality, multiple grouted specimens from a single column. Two selected natural sands having different grain sizes were grouted with two mineral-based silica grouts, resulting in different levels of improvement. The behaviour of the sands treated by mineral grouts, in terms of strength, initial stiffness and permeability, was compared with that exhibited by more traditional silicate grouts. The results of this study indicate that sands treated with mineral grouts result in higher strengths, higher initial shear modulus and lower permeability values than the sands treated with the silicate solution. The effect of grout type, effective confining pressure, and sand particle-size on small-strain shear modulus of grouted sand specimens was evaluated. Based on test results, the small strain shear modulus increment from treated to untreated specimens has been correlated with the unconfined compressive strength, obtaining a unique relationship regardless of grout type and grain-size of tested sands.