This work shows the Bouguer and aeromagnetic anomaly maps of the Iberian Chain (NE Spain) and surroundings in relation to its main geological and structural units. In general, both gravity and ...aeromagnetic anomalies are aligned following the same trends in both maps highlighting the major role played by the structural configuration of the Iberian Chain on their origin. In detail most anomalies do not coincide, pointing to different potential field sources to account for them. Further studies would be needed to decipher the origin of most magnetic anomalies observed in the Iberian Chain.
SUMMARY
By utilizing the addition theorems of the arctangent function and the logarithm, we developed a new expression of Bessel’s exact formula to compute the prismatic gravitational field using the ...triple difference of certain analytic functions. The use of the new expression is fast since the number of transcendental functions required is significantly reduced. The numerical experiments show that, in computing the gravitational potential, the gravity vector, and the gravity gradient tensor of a uniform rectangular parallelepiped, the new method runs 2.3, 2.3 and 3.7 times faster than Bessel’s method, respectively. Also, the new method achieves a slight increase in the computing precision. Therefore, the new method can be used in place of Bessel’s method in any situation. The same approach is applicable to the geomagnetic field computation.
SUMMARY
Within the iron metallogenic province of southeast Missouri, USA, there are several mines that contain not only economic iron resources, magnetite and/or hematite, but also contain rare earth ...elements, copper and gold. An area including three major deposits, Pea Ridge, Bourbon and Kratz Spring, was selected for detailed modelling for the upper crustal magnetic susceptibility and density structures. For the study area, ground gravity and high-resolution airborne magnetic and gravity gradiometry data sets are available. An efficient and novel joint inversion algorithm for the simultaneous inversion of these multiple data sets is presented. The Gramian coupling constraint is used to correlate the reconstructed density and magnetic susceptibility models. The implementation relies on the structures of the sensitivity matrices and an efficient minimization algorithm to achieve significant reductions in the memory requirements and computational costs. Consequently, it is feasible to use a laptop computer for the inversion of multiple data sets, each containing thousands of data points, for the recovery of models on the study area, each including approximately one million model parameters. This is the first time that these multiple data sets have been simultaneously inverted for this area. The L1-norm stabilizer is used to provide compact and focused images of the ore deposits. For contrast, independent inversions of each data set are also discussed. In general, our results provide new insights about the concealed ore deposits in the Mesoproterozoic basement rocks of southeast Missouri. Both short- and long-wavelength anomalies exist in the recovered models; these provide a high-resolution image of the subsurface. The geometry and physical properties of the known deposits are determined very well. Additionally, some unknown concealed deposits are revealed; these could be economically valuable and should be considered in future geophysical and geological investigations.
SUMMARY
We have developed a linear 3-D gravity inversion method capable of modelling complex geological regions such as subduction margins. Our procedure inverts satellite gravity to determine the ...best-fitting differential densities of spatially discretized subsurface prisms in a least-squares sense. We use a Bayesian approach to incorporate both data error and prior constraints based on seismic reflection and refraction data. Based on these data, Gaussian priors are applied to the appropriate model parameters as absolute equality constraints. To stabilize the inversion and provide relative equality constraints on the parameters, we utilize a combination of first and second order Tikhonov regularization, which enforces smoothness in the horizontal direction between seismically constrained regions, while allowing for sharper contacts in the vertical. We apply this method to the nascent Puysegur Trench, south of New Zealand, where oceanic lithosphere of the Australian Plate has underthrust Puysegur Ridge and Solander Basin on the Pacific Plate since the Miocene. These models provide insight into the density contrasts, Moho depth, and crustal thickness in the region. The final model has a mean standard deviation on the model parameters of about 17 kg m–3, and a mean absolute error on the predicted gravity of about 3.9 mGal, demonstrating the success of this method for even complex density distributions like those present at subduction zones. The posterior density distribution versus seismic velocity is diagnostic of compositional and structural changes and shows a thin sliver of oceanic crust emplaced between the nascent thrust and the strike slip Puysegur Fault. However, the northern end of the Puysegur Ridge, at the Snares Zone, is predominantly buoyant continental crust, despite its subsidence with respect to the rest of the ridge. These features highlight the mechanical changes unfolding during subduction initiation.
SUMMARY
Computation of gravimetric terrain corrections (TCs) is a numerical challenge, especially when using very high-resolution (say, ∼30 m or less) digital elevation models (DEMs). TC computations ...can use spatial or/and spectral techniques: Spatial domain methods are more exact but can be very time-consuming; the discrete/fast Fourier transform (D/FFT) implementation of a binomial expansion is efficient, but fails to achieve a convergent solution for terrain slopes >45°. We show that this condition must be satisfied for each and every computation-roving point pair in the whole integration domain, not just at or near the computation points. A combination of spatial and spectral methods has been advocated by some through dividing the integration domain into inner and outer zones, where the TC is computed from the superposition of analytical mass-prism integration and the D/FFT. However, there remain two unresolved issues with this combined approach: (1) deciding upon a radius that best separates the inner and outer zones and (2) analytical mass-prism integration in the inner zone remains time-consuming, particularly for high-resolution DEMs. This paper provides a solution by proposing: (1) three methods to define the radius separating the inner and outer zones and (2) a numerical solution for near-zone TC computations based on the trapezoidal and Simpson's rules that is sufficiently accurate w.r.t. the exact analytical solution, but which can reduce the computation time by almost 50 per cent.
Accurate bathymetry estimation is made possible by combining depth data with free-air gravity anomalies on the sea surface recovered from the geoidal heights that are equivalent to the mean sea ...surface derived from satellite radar altimetry. The residual gravity anomalies that represent the short-wavelength effect are required to accurately estimate bathymetry by combining satellite altimetry-derived free-air gravity anomalies and shipborne data including depth and gravity anomalies. In this study, the optimized ensemble model of machine learning techniques was applied to the residual gravity anomalies to estimate bathymetry by the gravity–geologic method (GGM) from various geospatial information including shipborne depth, shipborne gravity anomalies, and satellite altimetry-derived free-air gravity anomalies, in the Ulleung Basin in the East Sea. From the results, the GGM bathymetry predicted using the optimized ensemble model of machine learning was improved by 32.3 m over the GGM bathymetry estimated using the original depth and gravity anomalies. The method presented in this study is for estimating deep-water bathymetry using machine learning, and it has been proven to have superior performance compared with conventional methods.
We present a revised model for the opening of the South Atlantic Ocean founded on a remapping of the continent—ocean boundaries and Aptian salt basins, the chronology of magmatic activity in and ...around the ocean basin and on the timing and character of associated intraplate deformation in Africa and South America. The new plate tectonic model is internally consistent and consistent with globally balanced plate motion solutions. The model includes realistic scenarios for intraplate deformation, pre-drift extension and seafloor spreading. Within the model, Aptian salt basins preserved in the South American (Brazilian) and African (Angola, Congo, Gabon) continental shelves are reunited in their original positions as parts of a single syn-rift basin in near subtropical latitudes (10°S–27°S). The basin was dissected at around 112 Ma (Aptian—Albian boundary) when the model suggests that seafloor spreading commenced north of the Walvis Ridge—Rio Grande Rise.
Salt domes are favorable geological structures associated with hydrocarbon reservoirs and critical mineral resources. Currently, imaging salt diapir subsurface structures have received widespread ...attention. The gravity method has been widely used to interpret salt structures due to the significant density contrast between salts and surrounding sediments. However, complex mass sources at different depths challenge the interpretation of the gravity data over salt domes. In this study, we put forward the mixed L1- and L2-norm regularized inversion of the gravity data observed on the surface to image the 3-D subsurface structure of salt domes. The synthetic tests demonstrate that the mixed-norm regularization inversion of gravity data can effectively recover subsurface complex density anomalies similar to salt dome structures. The practical application in the Vinton Salt Dome reveals a shallow high-density caprock and a deep low-density salt core beneath Vinton Dome, greatly enriching the underground structure information of the study area. It demonstrates that the mixed-norm regularized inversion is a valuable interpretation technique for detecting and delineating salt structures.
Gravity anomaly over shallow waters is one of the fundamental data sources for studying sea level change, ocean currents, and water exchanges between coastal areas and open seas. However, the ...acquirement of gravity data over shallow waters faces multiple challenges due to the degraded quality of satellite altimetry data and the scarcity of surveyed gravimetric observations. To alleviate this problem, we establish a framework for marine gravity anomaly refinement by using satellite-derived bathymetry (SDB). We use a cosine-tapered bandpass filter to extract high-frequency gravity signals from the SDB data, which compensate for the unresolved signals in satellite altimetric gravity data. Numerical experiments over the Discovery Reef and an offshore region near the Port Hedland, Western Australia, demonstrate that the utilization of SDB effectively strengths marine gravity anomaly. By combining the SDB data, the fits between the enhanced gravity anomaly models and surveyed airborne gravity data are improved, by 5.3%-15.7% in comparison with an altimetric gravity model DTU21GRA. The SDB calculated from the linear band model has slightly better performances in gravity anomaly modeling than that computed from the band ratio model and physical-based approach, agreeing well with the SDB validation results. Our results verify the feasibility of using the SDB computed from the physical-based approach for gravity anomaly augmentation, which is of great value in areas devoid of ground-truth depths. This study cements a way for the augmentation of marine gravity anomaly worldwide, especially in remote regions characterized by the scarcity of ground-based gravity data.
SUMMARY
The traditional gravity forward modelling methods for solving partial differential equations (PDEs) only can yield second-order accuracy. When computing the gravity field vector and gradient ...tensor from the obtained potential, those numerical differentiation approaches will inevitably lose accuracy. To mitigate this issue, we propose an efficient and accurate 3-D forward modelling algorithm based on a fourth-order compact difference scheme. First, a 19-point fourth-order compact difference scheme with general meshsizes in x-, y- and z-directions is adopted to discretize the governing 3-D Poisson’s equation. The resulting symmetric positive-definite linear systems are solved by the pre-conditioned conjugate gradient algorithm. To obtain the first-order (i.e. the gravity field vector) and second-order derivatives (i.e. the gravity gradient tensor) with fourth-order accuracy, we seek to solve a sequence of tridiagonal linear systems resulting from the above mentioned finite difference discretization by using fast Thomas algorithm. Finally, two synthetic models and a real topography relief are used to verify the accuracy of our method. Numerical results show that our method can yield a nearly fourth-order accurate approximation not only to the gravitational potential, but also to the gravity field vector and its gradient tensor, which clearly demonstrates its superiority over the traditional PDE-based methods.