Spatial variation law of vertical derivative zero points for potential field data

In this paper we deduce the analytic solutions of the first- and second-order vertical derivative zero points for gravity anomalies in simple regular models with single, double, and multiple edges and analyze their spatial variation. For another simple regular models where it is difficult to obtain the analytic expression of the zero point, we try to use the profile zero points to analyze the spatial variation. The test results show that the spatial variation laws of both first- and second-order vertical derivative zero points are almost the same but the second-order derivative zero point position is closer to the top surface edge of the geological bodies than the first-order vertical derivative and has a relatively high resolution. Moreover, with an increase in buried depth, for a single boundary model, the vertical derivative zero point location tends to move from the top surface edge to the outside of the buried body but finally converges to a fixed value. For a double boundary model, the vertical derivative zero point location tends to migrate from the top surface edge to the outside of the buried body. For multiple boundary models, the vertical derivative zero point location converges from the top surface edge to the outside of the buried body where some zero points coincide and finally vanish. Finally, the effectiveness and reliability of the proposed method is verified using real field data.

A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential field data

Edge detection and enhancement techniques are commonly used in recognizing the edge of geologic bodies using potential field data. We present a new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative which has the functions of both edge detection and enhancement techniques. First, we calculate the total horizontal derivative （THDR） of the potential-field data and then compute the n-order vertical derivative （VDRn） of the THDR. For the n-order vertical derivative, the peak value of total horizontal derivative （PTHDR） is obtained using a threshold value greater than 0. This PTHDR can be used for edge detection. Second, the PTHDR value is divided by the total horizontal derivative and normalized by the maximum value. Finally, we used different kinds of numerical models to verify the effectiveness and reliability of the new edge recognition technology.

Vertical derivative conversion of irregular-range potential fi eld data based on improved projection onto convex sets method

Gravity and magnetic exploration areas are usually irregular,and there is some data defi ciency.Missing data must be interpolated before the vertical derivative conversion in the wavenumber domain.Meanwhile,for improved processing precision,the data need to be edge-padded to the length required by the fast Fourier transform algorithm.For conventional vertical derivative conversion of potential fi eld data(PFD),only vertical derivative conversion is considered,or interpolation,border padding,and vertical derivative conversion are executed independently.In this paper,these three steps are considered uniformly,and a vertical derivative conversion method for irregular-range PFD based on an improved projection onto convex sets method is proposed.The cutoff wavenumber of the filter used in the proposed method is determined by fractal model fi tting of the radial average power spectrum(RAPS)of the potential fi eld.Theoretical gravity models and real aeromagnetic data show the following:(1)The fitting of the RAPS with a fractal model can separate useful signals and noise reasonably.(2)The proposed iterative method has a clear physical sense,and its interpolation,border padding error,and running time are much smaller than those of the conventional kriging and minimum curvature methods.

Identification of the Caroline Plate boundary:constraints from magnetic anomaly

The Caroline Plate is located among the Pacific Plate,the Philippine Sea Plate,and the India Australia Plate,and plays a key role in controlling the spreading direction of the Philippine Sea Plate.The Caroline Submarine Plateau(or Caroline Ridge)and the Eauripik Rise on the south formed a remarkable T-shaped large igneous rock province,which covered the northern boundary between the Caroline Plate and the Pacific Plate.However,relationship between these tectonic units and magma evolution remains unclear.Based on magnetic data from the Earth Magnetic Anomaly Grid(2-arc-minute resolution)(V2),the normalized vertical derivative of the total horizontal derivative(NVDR-THDR)technique was used to study the boundary of the Caroline Plate.Results show that the northern boundary is a transform fault that runs 1400 km long in approximately 28 km wide along the N8°in E-W direction.The eastern boundary is an NNW-SSE trending fault zone and subduction zone with a width of tens to hundreds of kilometers;and the north of N4°is a fracture zone of dense faults.The southeastern boundary may be the Lyra Trough.The area between the southwestern part of the Caroline Plate and the Ayu Trough is occupied by a wide shear zone up to 100 km wide in nearly S-N trending in general.The Eauripik transform fault(ETF)in the center of the Caroline Plate and the fault zones in the east and west basins are mostly semi-parallel sinistral NNW-SSE–trending faults,which together with the eastern boundary Mussau Trench(MT)sinistral fault,the northern Caroline transform fault(CTF),and the southern shear zone of the western boundary,indicates the sinistral characteristics of the Caroline Plate.The Caroline hotspot erupted in the Pacific Plate near the CTF and formed the west Caroline Ridge,and then joined with the Caroline transform fault at the N8°.A large amount of magma erupted along the CTF,by which the east Caroline Ridge was formed.At the same time,a large amount of magma developed southward via the eastern branch of the ETF,forming the northern segment of the Eauripik Rise.Therefore,the magmatic activity of the T-shaped large igneous province is obviously related to the fault structure of the boundary faults between the Caroline Plate and Pacific Plate,and the active faults within the Caroline Plate.

Ground Magnetic, GPR, and Dipole-Dipole Resistivity for Landfill Investigation

Ground magnetic, ground penetrating radar (GPR), and dipole-dipole resistivity were carried out to environmentally investigate a landfill. In this context, these geophysical techniques were conducted to identify the subsurface contents of the landfill, furthermore, specify any possible leakage and/or contamination in the study area. The ground-magnetic survey carried out in the study area comprised 31 profiles each 120 m in length. Different wavelength filters were applied to the measured data. Vertical derivative, downward continuation, apparent susceptibility, band-pass, and analytical signal filters separated successfully the shallow sources. Whereas, upward continuation and low-pass Gaussian filters isolated significantly the deep magnetic sources. 3D Euler deconvolution (SI = 3) remarkably estimated the depths of the shallow sources (0 - 10 m) of the landfill contents. The conducted GPR and dipole-dipole resistivity allocated tangibly the locations and depths of the near surface anomalies. Both techniques didn’t reveal any possible leakage and/or contamination. Noteworthy, integration among magnetic, GPR, and dipole-dipole resistivity confirmed positively the results of each method. Nevertheless, some anomalies were recognized successfully by one technique and not by the others.