Matched bandpass filtering separates potential-field data into anomaly components representing different source depths (Phillips 2001). The implementation of the filtering was applied to the Mauritania magnetic data (Fig. 2) using a GX in Oasis Montaj based on a USGS software package (Phillips 1997). Matched filtering is accomplished in three steps. The initialization program first prepares the input data grid for Fourier transform by extending the rows and columns and then computes the Fourier transform and the natural logarithms of the radially-symmetric part (RSP) of the Fourier power spectrum and the non-radially symmetric part of the power spectrum. The matched bandpass filters are designed interactively by fitting equivalent source layers to the log of RSP power (Fig. 4) and also non-linearly adjusts the equivalent layer parameters to better fit the observed spectrum. The actual bandpass filtering (Fig. 5) calculates the inverse Fourier transforms and removes the row and column extensions (Phillips 2001).
The RSP of the power spectrum of the Mauritania data was matched by a three-layer equivalent model (Fig. 4), consisting of two shallow equivalent dipole layers and a single deep equivalent magnetic half space. The three corresponding bandpass filters (Fig. 5) were applied to the new data (Fig. 1) and data merged with the UN (Fig. 2).
Figure 4. The radially symmetric power spectrum of the data (green) and the power spectrum of a matching three- layer equivalent model (smooth blue line).
Figure 5. The matched bandpass filters corresponding to the three equivalent layers.
The results of applying bandpass filters 1 and 2 to the data are shown in Figures 6 and 7 respectively. Bandpass 3 (not shown) was not used in the analysis because it primarily contained very short-wavelength, low-amplitude noise and small anomalies associated with dikes, sills and BIF’s that are contained in the bandpass 2 filtered data. In Mauritania, matched bandpass filtering was used to separate the anomalies produced by the shallow dikes, sills, BIF’s and ophilolites (Fig. 6) from anomalies produced by deeper crystalline basement (Fig. 7). This low frequency magnetic map looks similar to the pseudo gravity map (Fig. 3).
Figure 6. Bandpass filtered image (filter 2, Fig. 5) of the magnetic data.
Figure 7. Bandpass filtered image (filter 1, Fig. 5) of the magnetic data.
Share with your friends: |