scholarly journals Application of the Wiener filter to magnetic profiling in the volcanic environment of Mt. Etna (Italy)

1996 ◽  
Vol 39 (1) ◽  
Author(s):  
C. Del Negro
Keyword(s):  
Synthetic Data ◽  
Wiener Filter ◽  
Magnetic Anomalies ◽  
Fault System ◽  
Volcanic Area ◽  
Wiener Filters ◽  
Mt Etna ◽  
Band Pass ◽  
Magnetic Profile ◽  
Mean Square Error Criterion

The frequency-domain Wiener filtering was applied to magnetic anomalies in the volcanic area of Mt. Etna. This filter, under suitable conditions (additive noise, linear processing and mean-square error criterion), can furnish an effective tool for discriminating the geologic feature of interest (the signal) from the noise. The filter was first tested with synthetic data. Afterwards it was applied to a magnetic profile carried out across the principal fault system of the Mt. Etna volcano, that hosted the dykes feeding both the 1989 and the 1991-93 eruptions. The magnetic anomalies linked to the volcanic section and those linked to the contact between the clay basement and the lava coverage show significant spectral overlap. Thus by estimating the power spectrum of the signal, obtained resolving the forward problem, a least-squares Wiener filter has been designed. In such context, it was possible to verify the effectiveness of Wiener filters, whereas traditional band-pass filtering proved inadequate. In fact, analysis of the noise showed that all the meaningful components of the observed magnetic field were resolved. The results put further constraints on location and geometry of the shallow plumbing system of Mt. Etna.

THE DETERMINATION OF DIGITAL WIENER FILTERS BY MEANS OF GRADIENT METHODS

Geophysics ◽  
1973 ◽  
Vol 38 (2) ◽  
pp. 310-326 ◽  
Author(s):  
R. J. Wang ◽  
S. Treitel
Keyword(s):  
Seismic Data ◽  
High Speed ◽  
Gradient Methods ◽  
Wiener Filter ◽  
Gradient Algorithm ◽  
Rapid Convergence ◽  
True Solution ◽  
Peripheral Device ◽  
Wiener Filters

The normal equations for the discrete Wiener filter are conventionally solved with Levinson’s algorithm. The resultant solutions are exact except for numerical roundoff. In many instances, approximate rather than exact solutions satisfy seismologists’ requirements. The so‐called “gradient” or “steepest descent” iteration techniques can be used to produce approximate filters at computing speeds significantly higher than those achievable with Levinson’s method. Moreover, gradient schemes are well suited for implementation on a digital computer provided with a floating‐point array processor (i.e., a high‐speed peripheral device designed to carry out a specific set of multiply‐and‐add operations). Levinson’s method (1947) cannot be programmed efficiently for such special‐purpose hardware, and this consideration renders the use of gradient schemes even more attractive. It is, of course, advisable to utilize a gradient algorithm which generally provides rapid convergence to the true solution. The “conjugate‐gradient” method of Hestenes (1956) is one of a family of algorithms having this property. Experimental calculations performed with real seismic data indicate that adequate filter approximations are obtainable at a fraction of the computer cost required for use of Levinson’s algorithm.

scholarly journals Interferometric SAR Phase Denoising Using Proximity-Based K-SVD Technique

Sensors ◽  
2019 ◽  
Vol 19 (12) ◽  
pp. 2684
Author(s):  
Chandrakanta Ojha ◽  
Adele Fusco ◽  
Innocenzo M. Pinto
Keyword(s):  
Noise Reduction ◽  
Synthetic Data ◽  
Sar Interferometry ◽  
Phase Image ◽  
Redundant Representations ◽  
X Band ◽  
Interferometric Phase ◽  
Phase Images ◽  
Mt Etna ◽  
Interferometric Sar

This paper addresses the problem of interferometric noise reduction in Synthetic Aperture Radar (SAR) interferometry based on sparse and redundant representations over a trained dictionary. The idea is to use a Proximity-based K-SVD (ProK-SVD) algorithm on interferometric data for obtaining a suitable dictionary, in order to extract the phase image content effectively. We implemented this strategy on both simulated as well as real interferometric data for the validation of our approach. For synthetic data, three different training dictionaries have been compared, namely, a dictionary extracted from the data, a dictionary obtained by a uniform random distribution in [ − π , π ] , and a dictionary built from discrete cosine transform. Further, a similar strategy plan has been applied to real interferograms. We used interferometric data of various SAR sensors, including low resolution C-band ERS/ENVISAT, medium L-band ALOS, and high resolution X-band COSMO-SkyMed, all over an area of Mt. Etna, Italy. Both on simulated and real interferometric phase images, the proposed approach shows significant noise reduction within the fringe pattern, without any considerable loss of useful information.

Detection of sharp lateral discontinuities through the analysis of surface-wave propagation

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. EN77-EN90 ◽  
Author(s):  
Paolo Bergamo ◽  
Laura Valentina Socco
Keyword(s):  
Surface Wave ◽  
Fault Location ◽  
Synthetic Data ◽  
Finite Element Method Simulation ◽  
Fault System ◽  
Energy Concentration ◽  
Data Sets ◽  
Data Set ◽  
Shallow Subsurface ◽  
Velocity Models

Surface-wave (SW) techniques are mainly used to retrieve 1D velocity models and are therefore characterized by a 1D approach, which might prove unsatisfactory when relevant 2D effects are present in the investigated subsurface. In the case of sharp and sudden lateral heterogeneities in the subsurface, a strategy to tackle this limitation is to estimate the location of the discontinuities and to separately process seismic traces belonging to quasi-1D subsurface portions. We have addressed our attention to methods aimed at locating discontinuities by identifying anomalies in SW propagation and attenuation. The considered methods are the autospectrum computation and the attenuation analysis of Rayleigh waves (AARW). These methods were developed for purposes and/or scales of analysis that are different from those of this work, which aims at detecting and characterizing sharp subvertical discontinuities in the shallow subsurface. We applied both methods to two data sets, synthetic data from a finite-element method simulation and a field data set acquired over a fault system, both presenting an abrupt lateral variation perpendicularly crossing the acquisition line. We also extended the AARW method to the detection of sharp discontinuities from large and multifold data sets and we tested these novel procedures on the field case. The two methods are proven to be effective for the detection of the discontinuity, by portraying propagation phenomena linked to the presence of the heterogeneity, such as the interference between incident and reflected wavetrains, and energy concentration as well as subsequent decay at the fault location. The procedures we developed for the processing of multifold seismic data set showed to be reliable tools in locating and characterizing subvertical sharp heterogeneities.

scholarly journals Northward migration of the Oregon forearc on the Gales Creek fault

Geosphere ◽  
2020 ◽  
Vol 16 (2) ◽  
pp. 660-684 ◽  
Author(s):  
Ray E. Wells ◽  
Richard J. Blakely ◽  
Sean Bemis
Keyword(s):  
Gravity Data ◽  
Magnetic Anomalies ◽  
Fault System ◽  
Coast Range ◽  
Geologic Mapping ◽  
Willamette Valley ◽  
Potential Field Data ◽  
Columbia River Basalt ◽  
Geophysical Surveys ◽  
Step Over

Abstract The Gales Creek fault (GCF) is a 60-km-long, northwest-striking dextral fault system (west of Portland, Oregon) that accommodates northward motion and uplift of the Oregon Coast Range. New geologic mapping and geophysical models confirm inferred offsets from earlier geophysical surveys and document ∼12 km of right-lateral offset of a basement high in Eocene Siletz River Volcanics since ca. 35 Ma and ∼8.8 km of right-lateral separation of Miocene Columbia River Basalt at Newberg, Oregon, since 15 Ma (∼0.62 ± 0.12 mm/yr, average long-term rate). Relative uplift of Eocene Coast Range basalt basement west of the fault zone is at least 5 km based on depth to basement under the Tualatin Basin from a recent inversion of gravity data. West of the city of Forest Grove, the fault consists of two subparallel strands ∼7 km apart. The westernmost, Parsons Creek strand, forms a linear valley southward to Henry Hagg Lake, where it continues southward to Newberg as a series of en echelon strands forming both extensional and compressive step-overs. Compressive step-overs in the GCF occur at intersections with ESE-striking sinistral faults crossing the Coast Range, suggesting the GCF is the eastern boundary of an R′ Riedel shear domain that could accommodate up to half of the ∼45° of post–40 Ma clockwise rotation of the Coast Range documented by paleomagnetic studies. Gravity and magnetic anomalies suggest the western strands of the GCF extend southward beneath Newberg into the Northern Willamette Valley, where colinear magnetic anomalies have been correlated with the Mount Angel fault, the proposed source of the 1993 M 5.7 Scotts Mills earthquake. The potential-field data and water-well data also indicate the eastern, Gales Creek strand of the fault may link to the NNW-striking Canby fault through the E-W Beaverton fault to form a 30-km-wide compressive step-over along the south side of the Tualatin Basin. LiDAR data reveal right-lateral stream offsets of as much as 1.5 km, shutter ridges, and other youthful geomorphic features for 60 km along the geophysical and geologic trace of the GCF north of Newberg, Oregon. Paleoseismic trenches document Eocene bedrock thrust over 250 ka surficial deposits along a reverse splay of the fault system near Yamhill, Oregon, and Holocene motion has been recently documented on the GCF along Scoggins Creek and Parsons Creek. The GCF could produce earthquakes in excess of Mw 7, if the entire 60 km segment ruptured in one earthquake. The apparent subsurface links of the GCF to other faults in the Northern Willamette Valley suggest that other faults in the system may also be active.

The IDQ curve: A tool for evaluating the direction of remanent magnetization from magnetic anomalies

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. J85-J98
Author(s):  
Shuang Liu ◽  
Xiangyun Hu ◽  
Dalian Zhang ◽  
Bangshun Wei ◽  
Meixia Geng ◽  
...  
Keyword(s):  
Remanent Magnetization ◽  
Field Data ◽  
Eastern China ◽  
Synthetic Data ◽  
Magnetic Anomalies ◽  
Natural Remanent Magnetization ◽  
Magma Intrusion ◽  
Magnetization Direction ◽  
Ore Bodies ◽  
Using Data

Natural remanent magnetization acts as a record of the previous orientations of the earth’s magnetic field, and it is an important feature when studying geologic phenomena. The so-called IDQ curve is used to describe the relationship between the inclination ( I) and declination ( D) of remanent magnetization and the Köenigsberger ratio ( Q). Here, we construct the IDQ curve using data on ground and airborne magnetic anomalies. The curve is devised using modified approaches for estimating the total magnetization direction, e.g., identifying the maximal position of minimal reduced-to-the-pole fields or identifying correlations between total and vertical reduced-to-the-pole field gradients. The method is tested using synthetic data, and the results indicate that the IDQ curve can provide valuable information on the remanent magnetization direction based on available data on the Köenigsberger ratio. Then, the method is used to interpret field data from the Yeshan region in eastern China, where ground anomalies have been produced by igneous rocks, including diorite and basalt, which occur along with magnetite and hematite ore bodies. The IDQ curves for 24 subanomalies are constructed, and these curves indicate two main distribution clusters of remanent magnetization directions corresponding to different structural units of magma intrusion and help identify the lithologies of the magnetic sources in areas covered by Quaternary sediments. The estimated remanent magnetization directions for Cenozoic basalt are consistent with measurements made in paleomagnetism studies. The synthetic and field data indicate that the IDQ curve can be used to efficiently estimate the remanent magnetization direction from a magnetic anomaly, which could help with our understanding of geologic processes in an area.

Interpretation of magnetic anomalies of geological and archaeological origins in a volcanic area (Tusculum site, Lazio, Italy): Methodological proposals

2020 ◽  
Vol 173 ◽  
pp. 103942
Author(s):  
P. Diarte-Blasco ◽  
A.M. Casas ◽  
A. Pocoví ◽  
J.J. Villalaín ◽  
A. Muñoz ◽  
...  
Keyword(s):  
Magnetic Anomalies ◽  
Volcanic Area

scholarly journals Real time retrieval of volcanic cloud particles and SO<sub>2</sub> by satellite using an improved simplified approach

2016 ◽  
Vol 9 (7) ◽  
pp. 3053-3062 ◽  
Author(s):  
Sergio Pugnaghi ◽  
Lorenzo Guerrieri ◽  
Stefano Corradini ◽  
Luca Merucci
Keyword(s):  
Sulfur Dioxide ◽  
Atmospheric Model ◽  
Volcanic Plume ◽  
Volcanic Area ◽  
Simplified Approach ◽  
Mt Etna ◽  
Plume Temperature ◽  
Volcanic Cloud ◽  
Preliminary Preparation ◽  
The Mean

Abstract. Volcanic plume removal (VPR) is a procedure developed to retrieve the ash optical depth, effective radius and mass, and sulfur dioxide mass contained in a volcanic cloud from the thermal radiance at 8.7, 11, and 12 µm. It is based on an estimation of a virtual image representing what the sensor would have seen in a multispectral thermal image if the volcanic cloud were not present. Ash and sulfur dioxide were retrieved by the first version of the VPR using a very simple atmospheric model that ignored the layer above the volcanic cloud. This new version takes into account the layer of atmosphere above the cloud as well as thermal radiance scattering along the line of sight of the sensor. In addition to improved results, the new version also offers an easier and faster preliminary preparation and includes other types of volcanic particles (andesite, obsidian, pumice, ice crystals, and water droplets). As in the previous version, a set of parameters regarding the volcanic area, particle types, and sensor is required to run the procedure. However, in the new version, only the mean plume temperature is required as input data. In this work, a set of parameters to compute the volcanic cloud transmittance in the three quoted bands, for all the aforementioned particles, for both Mt. Etna (Italy) and Eyjafjallajökull (Iceland) volcanoes, and for the Terra and Aqua MODIS instruments is presented. Three types of tests are carried out to verify the results of the improved VPR. The first uses all the radiative transfer simulations performed to estimate the above mentioned parameters. The second one makes use of two synthetic images, one for Mt. Etna and one for Eyjafjallajökull volcanoes. The third one compares VPR and Look-Up Table (LUT) retrievals analyzing the true image of Eyjafjallajökull volcano acquired by MODIS aboard the Aqua satellite on 11 May 2010 at 14:05 GMT.

Sign in / Sign up

Export Citation Format

Share Document