The use of «native» wavelet transform for determining lateral density variation of the Volgo-Uralian subcraton

In the last two decades in conjunction with the development of satellite gravimetry, the techniques of regional-scale inverse and forward gravity modeling started to be more actively incorporated in the construction of crustal and lithospheric scale models. Such regional models are usually built as a set of layers and bodies with constant densities. This approach often leads to a certain difference between the initially used measured gravity field and a gravity field that is produced by the model. One of the examples of this kind of models is a recent lithospheric model of the Volgo-Uralian subcraton. In the current study, we are applying the method of «native» wavelet transform to the residual gravity anomaly for defining the possible lateral density variations within the lithospheric layers of Volgo-Uralia. Keywords: wavelet transform; gravity field inversion; forward gravity modeling; Volgo-Uralian subcraton; satellite gravimetry.

New large subglacial lake in Princess Elizabeth Land, East Antarctica, detected by airborne geophysical observations

Knowledge of subglacial lakes is important for understanding the stability of the Antarctica Ice Sheet (AIS) and its contribution to the global sea-level change. We designed an intensified airborne campaign to collect geophysical data in Princess Elizabeth Land (PEL), East Antarctica, during the 2015–2019 CHINARE expeditions. We developed an innovative method to build a set of evidence of a newly detected subglacial lake, Lake Zhongshan. Adaptive RES data analysis allowed us to detect the lake surface and extent. We quantified the lake depth and volume via gravity modeling. Another dataset collected at Lake Vostok provided the ground truth. The results revealed that Lake Zhongshan, located at 73°26'53"S, 80°30'39"E and ~3,603 m below surface, has an area of 328 ± 1 km2, making it the only one in PEL and the fifth largest in Antarctica. These findings are important for understanding subglacial hydrodynamics in PEL, as well as the stability of the AIS.

Integration of magnetic, gravity, and well-logging data interpretation to delineate the structural framework and formation evaluation of Bahariya Formation, North Diyur area, northern Western Desert, Egypt

The Bahariya Formation in Egypt's Western Desert is a major source for minerals and hydrocarbon accumulation. It is also characterized by a relatively high radiation content because it contains iron oxide deposits that attract radioactive elements. The main objectives of our study are to establish depth to basement, basement configuration and related structural elements, and thickness and configuration of the overlain sedimentary section. In addition to the analysis of well-logging data, many advanced techniques have been applied to analyze magnetic and gravity data, including depth estimation, 2D magnetic and gravity modeling, and 3D inversion of potential field data. By integrating all available data, we can determine the structural control of the study area and evaluate the subsurface parameters. Well logging has been used for interpretation of porous and permeable zones, water saturation calculation, and basic lithology identification. The depth to basement in our study ranges from −1700 to −4500 m. The basement is shallow in the northern parts of the study area and deeper in the southern parts. The main clay minerals of the formation are montmorillonite, chlorite, and a mixed clay layer. The Bahariya Formation is composed mainly of sandy clay and sandstone, and therefore it is considered an excellent reservoir.

Basement Mapping of the Fucino Basin in Central Italy by ITRESC Modeling of Gravity Data

Sediments infilling in intermontane basins in areas with high seismic activity can strongly affect ground-shaking phenomena at the surface. Estimates of thickness and density distribution within these basin infills are crucial for ground motion amplification analysis, especially where demographic growth in human settlements has implied increasing seismic risk. We employed a 3D gravity modeling technique (ITerative RESCaling—ITRESC) to investigate the Fucino Basin (Apennines, central Italy), a half-graben basin in which intense seismic activity has recently occurred. For the first time in this region, a 3D model of the Meso-Cenozoic carbonate basement morphology was retrieved through the inversion of gravity data. Taking advantage of the ITRESC technique, (1) we were able to (1) perform an integration of geophysical and geological data constraints and (2) determine a density contrast function through a data-driven process. Thus, we avoided assuming a priori information. Finally, we provided a model that honored the gravity anomalies field by integrating many different kinds of depth constraints. Our results confirmed evidence from previous studies concerning the overall shape of the basin; however, we also highlighted several local discrepancies, such as: (a) the position of several fault lines, (b) the position of the main depocenter, and (c) the isopach map. We also pointed out the existence of a new, unknown fault, and of new features concerning known faults. All of these elements provided useful contributions to the study of the tectono-sedimentary evolution of the basin, as well as key information for assessing the local site-response effects, in terms of seismic hazards.

Crustal structure of the Volgo-Uralian subcraton revealed by inverse and forward gravity modeling

Volgo-Uralia is a Neoarchean easternmost part of the East European craton. Recent seismic studies of the Volgo-Uralian region provided new insights into the crustal structure of this area. In this study, we combine satellite gravity and seismic data in a common workflow to perform a complex study of Volgo-Uralian crustal structure which is useful for further basin analysis of the area. In this light, a new crustal model of the Volgo-Uralian subcraton is presented from a step-wise approach: (1) inverse gravity modeling followed by (2) 3D forward gravity modeling. First, inversion of satellite gravity gradient data was applied to determine the Moho depth for the area. Density contrasts between crust and mantle were varied laterally according to the tectonic units present in the region, and the model is constrained by the available active seismic data. The Moho discontinuity obtained from the gravity inversion was consequently modified and complemented in order to define a complete 3D crustal model by adding information on the sedimentary cover, upper crust, lower crust, and lithospheric mantle layers in the process of forward gravity modeling where both seismic and gravity constraints were respected. The obtained model shows crustal thickness variations from 32 to more than 55 km in certain areas. The thinnest crust with a thickness below 40 km is found beneath the Pericaspian basin, which is covered by a thick sedimentary layer. The thickest crust is located underneath the Ural Mountains as well as in the center of the Volga-Uralian subcraton. In both areas the crustal thickness exceeds 50 km. At the same time, initial forward gravity modeling has shown a gravity misfit of ca. 95 mGal between the measured Bouguer gravity anomaly and the forward calculated gravity field in the central area of the Volga-Uralian subcraton. This misfit was interpreted and modeled as a high-density lower crust which possibly represents underplated material. Our preferred crustal model of the Volga-Uralian subcraton respects the gravity and seismic constraints and reflects the main geological features of the region with Moho thickening in the cratons and under the Ural Mountains and thinning along the Paleoproterozoic rifts, Pericaspian sedimentary basin, and Pre-Urals foredeep.

Joint Seismic and Gravity Data Inversion to Image Intra-Crustal Structures: The Ivrea Geophysical Body Along the Val Sesia Profile (Piedmont, Italy)

We present results from a joint inversion of new seismic and recently compiled gravity data to constrain the structure of a prominent geophysical anomaly in the European Alps: the Ivrea Geophysical Body (IGB). We investigate the IGB structure along the West-East oriented Val Sesia profile at higher resolution than previous studies. We deployed 10 broadband seismic stations at 5 km spacing for 27 months, producing a new database of ∼1000 high-quality seismic receiver functions (RFs). The compiled gravity data yields 1 gravity point every 1–2 km along the profile. We set up an inversion scheme, in which RFs and gravity anomalies jointly constrain the shape and the physical properties of the IGB. We model the IGB’s top surface as a single density and shear-wave velocity discontinuity, whose geometry is defined by four, spatially variable nodes between far-field constraints. An iterative algorithm was implemented to efficiently explore the model space, directing the search toward better fitting areas. For each new candidate model, we use the velocity-model structures for both ray-tracing and observed-RFs migration, and for computation and migration of synthetic RFs: the two migrated images are then compared via cross-correlation. Similarly, forward gravity modeling for a 2D density distribution is implemented. The joint inversion performance is the product of the seismic and gravity misfits. The inversion results show the IGB protruding at shallow depths with a horizontal width of ∼30 km in the western part of the profile. Its shallowest segment reaches either 3–7 or 1–3 km depth below sea-level. The latter location fits better the outcropping lower crustal rocks at the western edge of the Ivrea-Verbano Zone. A prominent, steep eastward-deepening feature near the middle of the profile, coincident with the Pogallo Fault Zone, is interpreted as inherited crustal thickness variation. The found density and velocity contrasts of the IGB agree with physical properties of the main rock units observed in the field. Finally, by frequency-dependent analysis of RFs, we constrain the sharpness of the shallowest portion of the IGB velocity discontinuity as a vertical gradient of thickness between 0.8 km and 0.4 km.

Landscape Connectivity Analysis and Optimization of Qianjiangyuan National Park, Zhejiang Province, China

As natural ecosystems in most parts of the world come under increasing human influence, fragmentation is becoming the major driving factor of the global biodiversity crisis. Therefore, connectivity between habitat patches is becoming even more important. China began building national parks with the primary purpose of protecting nationally representative natural ecosystems and maintaining the integrity of their structure, processes and functions. Research is necessary to improve the internal connectivity of national parks and to propose suggestions for existing functional zoning and biological corridors. In this study, Qianjiangyuan National Park was selected as an example park, and landscape fragmentation was evaluated exponentially and simulated visually. The habitat characteristics of protected species in the region, morphological spatial pattern analysis and the delta of the probability of connectivity were used together to identify key habitat patches and their importance levels in the study area. Potential habitat corridors in the region were then obtained using least-cost path analysis and gravity modeling methods based on the distribution of key habitat and the migration costs of target species. The results of this study show that the disturbed landscape of the study area is dominated by tea plantations and drylands, with central roads being an important factor affecting the overall landscape connectivity. In terms of the distribution of key habitat patches, the mountains have a high value. In terms of area, their size is not directly proportional to their importance for maintaining landscape connectivity in the region, but large area patches are generally of higher importance. In terms of distance, key habitats that are closer to each other have a stronger correlation and a greater possibility for species migration. Combined with the functional zoning of Qianjiangyuan National Park, the setting of strictly protected areas and recreational areas is reasonable, and traditional use areas and ecological conservation areas could be appropriately adjusted according to the distribution of key habitats. The important corridor in the middle of the ecological conservation area is crucial for the overall connectivity of the national park, and the connectivity between strict protected areas will depend on successful protection of the ecological conservation area.