Understanding the Geology of the Philippines through Gravity Anomalies

The Philippine Archipelago is a complex island arc system, where many regions still lack geopotential studies. This study aims to present a general discussion of the Philippine gravity anomaly distribution. The high-resolution isostatic anomaly digital grid from the World Gravity Map (WGM) was processed and correlated with the Philippines’ established geology and tectonics. This study also investigated the gravity signatures that correspond to the regional features, e.g., geology, structures, sedimentary basins, and basement rocks of the study area. Upward continuation, high-pass, and gradient filters (i.e., first vertical derivative, horizontal gradient) were applied using the Geosoft Oasis Montaj software. The interpreted gravity maps’ results highlighted the known geologic features (e.g., trench manifestation, ophiolite distribution, basin thickness). They revealed new gravity anomalies with tectonic significance (e.g., basement characterization). The isostatic gravity anomaly map delineates the negative zones. These zones represent the thick sedimentary accumulations along the trenches surrounding the Philippine Mobile Belt (PMB). The Philippine island arc system is characterized by different gravity anomaly signatures, which signify the density contrast of subsurface geology. The negative anomalies (< 0 mGal) represent the thick sedimentary basins, and the moderate signatures (0 to 80 mGal) correspond to the metamorphic belts. The distinct very high gravity anomalies (> 80 mGal) typify the ophiolitic basement rocks. The gravity data’s upward continuation revealed contrasting deep gravity signatures; the central Philippines of continental affinity (20 – 35 mGal) was distinguished from the remaining regions of oceanic affinity (45 – 200 mGal). Local geologic features (e.g., limestone, ophiolitic rocks) and structures (e.g., North Bohol Fault, East Bohol Fault) were also delineated downward continuation and gravity gradient maps of Bohol Island. The WGM dataset’s effectiveness for geologic investigation was achieved by comparing the established geologic features and interpreted gravity anomalies. The processed gravity digital grids provided an efficient and innovative way of investigating the Philippines’ regional geology and tectonics.

The BRAT and GUT Couple: Broadview Radar Altimetry and GOCE User Toolboxes

&lt;p&gt;The scope of this work is to showcase the BRAT (Broadview Radar Altimetry Toolbox) and GUT (GOCE User Toolbox) toolboxes.&lt;/p&gt;&lt;p&gt;The Broadview Radar Altimetry Toolbox (BRAT) is a collection of tools designed to facilitate the processing of radar altimetry data from all previous and current altimetry missions, including Sentinel-3A L1 and L2 products. A tutorial is included providing plenty of use cases on Geodesy &amp; Geophysics, Oceanography, Coastal Zone, Atmosphere, Wind &amp; Waves, Hydrology, Land, Ice and Climate, which can also be consulted in &amp;#160;http://www.altimetry.info/radar-altimetry-tutorial/.&lt;/p&gt;&lt;p&gt;BRAT's last version (4.2.1) was released in June 2018. Based on the community feedback, the front-end has been further improved and simplified whereas the capability to use BRAT in conjunction with MATLAB/IDL or C/C++/Python/Fortran, allowing users to obtain desired data bypassing the data-formatting hassle, remains unchanged. Several kinds of computations can be done within BRAT involving the combination of data fields, that can be saved for future uses, either by using embedded formulas including those from oceanographic altimetry, or by implementing ad-hoc Python modules created by users to meet their needs. BRAT can also be used to quickly visualise data, or to translate data into other formats, e.g. from NetCDF to raster images.&lt;/p&gt;&lt;p&gt;The GOCE User Toolbox (GUT) is a compilation of tools for the use and the analysis of GOCE gravity field models. It facilitates using, viewing and post-processing GOCE L2 data and allows gravity field data, in conjunction and consistently with any other auxiliary data set, to be pre-processed by beginners in gravity field processing, for oceanographic and hydrologic as well as for solid earth applications at both regional and global scales. Hence, GUT facilitates the extensive use of data acquired during GRACE and GOCE missions.&lt;/p&gt;&lt;p&gt;In the current version (3.2), GUT has been outfitted with a graphical user interface allowing users to visually program data processing workflows. Further enhancements aiming at facilitating the use of gradients, the anisotropic diffusive filtering, and the computation of Bouguer and isostatic gravity anomalies have been introduced. Packaged with GUT is also GUT's Variance/Covariance Matrix (VCM) tool, which enables non-experts to compute and study, with relative ease, the formal errors of quantities &amp;#8211; such as geoid height, gravity anomaly/disturbance, radial gravity gradient, vertical deflections &amp;#8211; that may be derived from the GOCE gravity models.&lt;/p&gt;&lt;p&gt;On our continuous endeavour to provide better and more useful tools, we intend to integrate BRAT into SNAP (Sentinel Application Platform). This will allow our users to easily explore the synergies between both toolboxes. During 2020 we will start going from separate toolboxes to a single one.&lt;/p&gt;&lt;p&gt;BRAT and GUT toolboxes can be freely downloaded, along with ancillary material, at https://earth.esa.int/brat and https://earth.esa.int/gut.&lt;/p&gt;

The shape of the Variscan Belt in Central Europe: Strike-slip tectonics versus oroclinal bending

&lt;p&gt;The European Variscan belt sharply changes its trend in easternmost Germany and western Poland, where the ENE- to NE-striking structures are replaced by the ESE- to SE-trending ones. The structures of still another, NNE-SSW strike, take the lead, however, along the SE margin of the Bohemian Massif. The Variscan belt seems, thus, to make nearly a U-turn, encircling the Bohemian Massif from the north. This has been explained for almost a century by assuming a 180&amp;#176; oroclinal loop, in which the Rhenohercynian and Saxothuringian tectonostratigraphic zones inarm the core of the Bohemian Massif. According to this classical view, the outermost tectonostratigraphic zone of the Variscan belt, the Rhenohercynian Zone, continues eastward in the deep substratum of the Permian-Mesozoic basin and reappears at the surface along the eastern rim of the Bohemian Massif.&lt;/p&gt;&lt;p&gt;Since the late 1970s an alternative view has gained an increasing attention that postulates a dextral transpressional regime during the final accretion of the Variscan terranes. This transpressional tectonic context is believed to have resulted from sublatitudinal, right-lateral displacements between Gondwana and Laurussia. Near the Carboniferous-Permian boundary, Gondwana decoupled from the newly formed European Variscan belt and proceeded westward, toward the southern edge of the Laurentian segment of Laurussia, owing to the development of the Appalachian subduction system. Concomitantly with the peak of the Alleghanian orogeny during early Permian, the European Variscan belt experienced a crosscut of its major tectonic zones along a set of dextral strike-slip faults.&lt;/p&gt;&lt;p&gt;In this study, we investigate directions and continuity of structural trends in the external zones of the Variscan orogen in Poland and map a foreland extent of Variscan deformations using seismic, gravimetric-magnetic and borehole data. These permit us testing the orocline- vs strike-slip concepts and develop an overall kinematic model for the NE Variscides.&lt;/p&gt;&lt;p&gt;Matched filtering of isostatic gravity, guided by results of spectral analysis, along with other derivatives of gravity and magnetic fields reveal a dominant WNW-ESE-trending pre-Permian structural grain in the external zones of the Variscan belt in Poland. This trend is confirmed by regional distribution of dips in Carboniferous and Devonian strata that were penetrated by boreholes beneath Permian-Mesozoic sediments. Seismic constraints on the position of the Variscan deformation front come from (1) the GRUNDY 2003 seismic experiment, combining wide-angle reflection-refraction measurements with the near-vertical reflection seismics in central Poland and (2) PolandSPAN and POLCRUST-01 deep reflection profiles in SE Poland. The WNW-ESE structural trend in the Variscan foreland is parallel to a set of major strike-slip fault zones in the area that are considered to convey a significant dextral displacement between Laurussia and Gondwana. The revised position of the Variscan deformation front shows a similar, uninterrupted, generally WNW-ESE trend, up to the SE border of Poland, which indicates an initial continuation of the more internal Variscan zones into the area of the present-day Carpathians. The geometry of the Variscan deformation front along with the pattern of the Variscan structural grain are inconsistent with the idea of an oroclinal loop affecting the external, non-metamorphic Variscan belt.&lt;/p&gt;

Late Palaeozoic strike-slip tectonics versus oroclinal bending at the SW outskirts of Baltica: case of the Variscan belt’s eastern end in Poland

Geophysical and geological data from the eastern sector of the Central European Variscan belt are presented and reviewed in the regional tectonic context. Matched filtering of isostatic gravity, guided by results of spectral analysis, along with other derivatives of gravity and magnetic fields reveal a dominant WNW–ESE-trending pre-Permian structural grain in the external zones of the Variscan belt in Poland. This trend is confirmed by regional distribution of dips in Carboniferous and Devonian strata that were penetrated by boreholes beneath Permian-Mesozoic sediments. Based on these data, two alternative concepts explaining the connection of the Variscan belt and its NE foreland, those of strike-slip tectonics versus oroclinal bending, are discussed. The WNW–ESE structural trend in the Variscan foreland is parallel to a set of major strike-slip fault zones in the area, including those of Upper Elbe, Intra-Sudetic, Odra, Dolsk and Kraków-Lubliniec. These faults are considered to convey a significant dextral displacement between Laurussia and Gondwana. The revised position of the Variscan deformation front shows a similar, uninterrupted, generally WNW–ESE trend, up to the SE border of Poland, which indicates an initial continuation of the Variscan belt into the area of the present-day Western Carpathians. The geometry of the Variscan deformation front along with the pattern of the Variscan structural grain are inconsistent with the idea of an oroclinal loop affecting the external, non-metamorphic Variscan belt. However, the data presented do not entirely rule out an oroclinal loop within the Variscan internides. The still possible options are (1) a semi-oroclinal model postulating ~ 90° bending of the Variscan tectonostratigraphic zones into parallelism with the WNW–ESE strike-slip faults or (2) an orocline limited only to the belt linking the Wolsztyn High and Moravo-Silesian non- to weakly-metamorphic fold-and-thrust belt. Regardless of the kinematic model preferred, our data indicate that structural evolution of the Polish Variscides was concluded with the end-Carboniferous NNE–SSW shortening that resulted in the present-day extent of the Variscan deformation front.

Joint estimation of gravity anomalies using second and third order potential derivatives

Satellite gradiometry data provide the framework for estimating and validating Earth's gravity field from second and third order derivatives of the Earth's gravitational potential. Such procedures are especially useful when applied locally, as they relate to local and regional characteristics of the real gravity field. In the present study a joint inversion procedure is proposed for the estimation of gravity anomalies at sea surface level from second and third order potential derivatives, based on a standard Gauss-Markov estimation model. The estimation procedure is applied for a test area stretching over Iran involving simulated grids from GOCE-only model GGM_TIM_R05 at GOCE altitude and gravity anomalies recovered at sea level. In order to validate the proposed estimation three different reductions have been considered independently, namely the removal of the long-wavelength part of the observed field through a global gravity model, the removal of the high-frequency part of the field through the incorporation of a topographic/isostatic gravity model and the application of variance component estimation. The application of a global gravity model leads to an improvement in the individual component estimation of the order of magnitude 3 per cent to 73 per cent, with a significant reduction in bias to 4 mGal. Smoother gradient components can come out according to removing the topography and taking into account for isostasy that improved up results of recovery to 25 per cent for the radial second order derivative. Finally, the implementation of variance component estimation leads to no significant improvement in results of recovered gravity anomalies.