Basement Geometry Research Articles

3D deformation velocity field analysis and TEM method to detect the tectonic influence on the land subsidence in Zamora, Mexico

The city of Zamora is situated in a tectonic basin and, according to previous interferometric studies, has been experiencing sinking up to 13 cm/yr (2007-2011). Although the reported subsidence pattern (WNW-ESE) is similar to the orientation of the regional fault system, there is a lack of detailed studies to establish the connection between this phenomenon, the thickness of sediments, and the basement geometry. Therefore, this paper presents a 3D deformation velocity field analysis for the period of 2014–2021, using 130 images from the Sentinel-2 satellite to calculate the vertical and east-west components of the subsidence. Likewise, a campaign of 28 Transient Electromagnetic soundings was carried out around the city to determine the thickness of sediments and bedrock geometry. The results of the interferometric analysis reveal a WNW-ESE sinking pattern with the maximum vertical velocities observed near the Zamora Museum, amounting to 108 mm/yr (2014–2017) and 102 mm/yr (2018–2021). Meanwhile, the horizontal component has a heterogeneous behavior in both periods with rates reaching up to 1 cm/year. The Transient Electromagnetic survey results show the presence of two grabens, a horst, and a series of half-grabens. The thickness of sediments varies between 94 and 314 m, with the latter corresponding to the main graben located beneath the city. The axis of this structure and the largest sediment deposits are linked to the area of maximum sinking and the spatial pattern of subsidence.

A speculative ridge within the Kidson Sub-basin – integrated interpretation from geophysical data

This paper presents an analysis of the basement geometry and its implications on the Kidson Sub-basin in the southern Canning Basin. Integration of seismic data primarily near the rim of the sub-basin with airborne electromagnetic data across its central portion reveals an east-northeast oriented ridge at Permian level, suggesting a possible elevated structure in the Ordovician and basement. The proposed ridge effectively separates the Kidson Sub-basin into two distinct parts, providing an explanation for facies and thickness variations observed in the Nambeet and Goldwyer formations, as well as the absence of the Willara Formation in the southeastern portion of the sub-basin. The presence of shallowing trends observed in seismic profiles along the southern flank of the ridge suggests the proximity of the basement high. The ridge is speculative and requires further work to verify its existence, potentially including passive and reflection seismic surveys across the structure. If confirmed, it might be a significant feature with far-reaching implications for the prospectivity of resources in the Kidson Sub-basin. The hydrocarbon sourced from the Nambeet, Goldwyer and Bongabinni formations in the NNW depocentre may have migrated via the extension of the Parallel Range Fault and be trapped in the footwall block over the ridge. Compared to the structures in the Fitzroy Trough, the traps within the Kidson Sub-basin are expected to maintain reasonable integrity and have good potential for petroleum accumulation and carbon sequestration in the thick Paleozoic succession.

The Impact of the Bohemian Spur on the Cooling and Exhumation Pattern of the Eastern Alpine Wedge of the European Alps

AbstractFold and thrust belt architecture may be influenced by basement geometry of the downgoing plate. This influence is notoriously difficult to assess due to a common lack of subsurface constraints and low resolution of exhumation estimates in space and time. The Bohemian Spur is a basement high at the transition from the Alps to the Carpathians. It coincides with narrowing of the foreland basin and an orogen‐scale change of strike. Its location in one of the best‐studied orogens in the world makes it an ideal case for understanding how basement topography influences fold and thrust belt tectonics. However, since thermochronological studies were mainly focused on the core of the Alps, timing and amount of exhumation remain poorly constrained in these peripheral parts of the orogen. We present new apatite (U‐Th)/He and fission track data from the wedge above the Bohemian Spur. Thermally reset ages monitor a so far un(der)appreciated phase of prominent Late Oligocene to Miocene cooling, associated with crustal thickening, uplift and erosion during wedge propagation. Pronounced exhumation on the order of 3–4.5 km can be related to basement steps beneath the advancing wedge. The spur acted as a buttress for foreland‐propagating thrusting, pinning deformation and nucleating antiformal stacking and duplexing and thus exhumation above it. We illustrate how along‐ and across‐strike changes of sub‐detachment topography impact wedge propagation and control fold and thrust belt geometries. The buttressing effect accounts for most of the exhumation, while deep‐seated slab dynamics are of subordinate importance for wedge uplift.

Assessing drivers of high exhumation magnitudes and young cooling ages in the eastern central Andes, southern Peru (13–18°S)

The potential structural controls on exhumation across the southern Peruvian Andes are not well understood, in part due to limited structural studies that co-locate with thermochronometric datasets. We integrate these two datasets and evaluate the relative contribution that fault geometry, magnitude, and shortening rate have on predicted cooling ages. Here we present a balanced cross-section constructed using new structural observations. This section, combined with existing thermochronometer data and a thermokinematic model, investigates the drivers of high exhumation and young canyon thermochronometric ages along the deeply incised Marcapata canyon in southern Peru. Together, these approaches constrain the timing and magnitude of exhumation in this portion of the southern Peruvian Andes and provide a mechanism for documenting how the internal architecture changes along strike.The balanced cross-section (oriented N30E) covers the Subandean Zone to the northeast, the Marcapata canyon on the eastern flank of the southern Peruvian Andes, and the Altiplano-Eastern Cordillera boundary to the southwest (13–18° S). Exhumation is constrained by four low-temperature thermochronometer systems, including apatite and zircon (U-Th)/He (AHe and ZHe, respectively) and fission-track (AFT and ZFT, respectively). The youngest AHe (∼1–3 Ma), AFT (∼3–7 Ma), ZHe (∼4–7 Ma), and ZFT (∼14–17 Ma) ages are located in the center and valley bottom of the Marcapata canyon. The thermokinematically modeled cross-section produces cooling ages determined by fault geometry and kinematics. Reset ZFT ages require burial of Ordovician rocks in excess of 5.5 km above the original 6.5 km depositional depth. We find that the ZFT and ZHe ages in the Eastern Cordillera are sensitive to the history and magnitude of burial, age and location of uplift, and canyon incision. Canyon incision is required to reproduce the youngest canyon thermochronometric ages while slow shortening rates from ∼10 Ma to Present are required to reproduce interfluve thermochronometric ages. Shortening is accommodated by basement faults that feed slip up through three different décollement levels before reaching the surface. The proposed stacked basement geometry sets the first-order cooling signal seen in modeled ages. We determined that the total shortening in this section from the Subandean Zone to the Altiplano is 147.5 km, similar to shortening estimates in an adjacent thermo-kinematically modeled section in the San Gabán canyon 50 km to the southeast. Both the ZHe and ZFT ages in the Marcapata section (4–5 and 14 Ma) are noticeably younger than cooling ages from the San Gabán section (16 and 29 Ma). The Marcapata section's higher magnitude of exhumation is due to a repetition of basement thrusts that continues to elevate the Eastern Cordillera while active deformation occurs in the Subandean Zone. The youngest thermochronometric ages in all four systems are co-located with the overlapping basement thrust geometry. This basement geometry, kinematic sequence of deformation, and canyon incision co-conspire to produce the young cooling ages observed in the Eastern Cordillera.

Numerical Investigation of Existing Tunnel Deformation Induced by Basement Excavation Considering the Unloading Ratio

Basement excavation may induce deformations of the adjacent tunnels. The response of existing tunnels to basement excavation considering the critical unloading ratio is rarely studied. In this study, a three-dimensional numerical model is established to investigate basement–tunnel interaction. Then, the numerical model is validated by simulating the centrifuge model test. Thereafter, the influences of basement geometry and tunnel location relative to the basement on the vertical deformation of the tunnel are studied. The results show that the vertical deformation of the tunnel increases linearly with the unloading ratio, which describes the degree of excavation depth above the tunnel. But there exists a critical unloading ratio of 0.6, beyond which the vertical deformation of the tunnel increases significantly. On this basis, an empirical model is proposed to predict the vertical deformation of the tunnel considering the unloading ratio.

Imaging an Enigmatic Paleovalley with Passive Seismic Methods (Unaweep Canyon, Colorado, United States)

Abstract High-resolution passive seismic imaging of shallow subsurface structures is often challenged by the scarcity of coherent body-wave energy in ambient noise recorded at surface stations. We show that the autocorrelation (AC) of teleseismic P-wave coda extracted from just one month of continuous recording at 5 Hz geophones can overcome this limitation. We apply this method to investigate the longitudinal subsurface bedrock structure of Unaweep Canyon—a paleovalley in western Colorado (United States) with complex evolution. Both fluvial and glacial processes have been proposed to explain the canyon’s genesis and morphology. The teleseismic P-wave coda AC retrieves zero-offset reflections from the shallow (200–500 m depth) basement interface at 120 stations along a 5 km long profile. In addition, we invert interferometrically retrieved surface-wave dispersion for the shear-wave structure of the sedimentary fill. Combined interpretation of these results and other geophysical and well data suggests an overdeepened basement geometry most consistent with glacial processes.

Research on basin effect based on 2-D geophysical prospecting profile

The basin effects alter the frequency content, amplitude and duration of seismic waves with different depth and shape basin. In this paper, we selected 2-D geophysical profiles spanning the Weihe Basin with large area and the Luoquanli Basin with small area. Based on the velocity structure data of the two geophysical profiles, 2-D finite element models were established to analyze the influences of the sedimentary layer in a broad basin and the basement geometry of a small basin on the basin effects using the spectral element method. The results showed that the sedimentary deposits in the basin could significantly amplify the ground motion, and the diffraction waves generated by the lateral inhomogeneity of the basin could prolong the duration of the ground motion. The highlight of this paper, our results showed that the amplification characteristics were controlled by the velocity mean in the vertical direction of sedimentary deposits (vp < 4.5 km/s) for a broad basin with small depth-to-width ratio as Weihe Basin, while the amplification characteristics were dominated by basement geometry of the small basin with large depth-to-width ratio as Luoquanli Basin. In this study, a high-resolution 2-D geophysical refraction profile was directly used to construct a finite element model to study the basin effects, which provided a new method of studying the basin effects.

Metamorphic bedrock geometry of Santorini using HVSR information and geophysical modeling of ambient noise and active-source surface-wave data

Santorini is the most active volcano of the Southern Aegean Volcanic Arc, with volcanic material that has accumulated during the last 600+ Kyrs on top of the pre-volcanic Santorini island. The geometry of the pre-volcanic basement not only provides constraints on the volcanic history, but is also fundamental for seismic hazard assessment, as it partly controls site-effects on strong ground motions. We investigate the geometry of the metamorphic (Cycladic) basement of Santorini using information from both passive (noise) and active source data. We performed a large number (∼280) of single-station ambient noise measurements and collected additional data from 42 previous studies to compute the Horizontal-to-Vertical Spectral Ratio (HVSR) curves. The majority of the recovered HVSR curves show prominent spectral ratio peaks (locally larger than 7–8) indicating a large impedance contrast between pyroclastic volcanic strata and bedrock metamorphic formations. We also performed active (MASW) and passive (Noise Array) surface wave measurements at 4 selected sites to determine local dispersion curves. To find the 1D Vs structure from the surface down to the metamorphic bedrock we performed a joint Monte Carlo inversion of the dispersion curves and the fundamental frequency (f0HVSR) values from the HVSR curves. These models, as well as deep borehole information and HVSR data along the caldera cliff, where the thickness and stratigraphy of the volcanic formations above the bedrock can be directly observed, allowed us to determine a linear relationship between the thickness of the pyroclastic formations and their fundamental resonance period for central-southern Santorini. The resulting geometry of pre-volcanic island shows deep basins around the pre-Alpine bedrock outcrop (Profitis Ilias) and generally agrees with recent results from larger-scale tomographic and gravity data but with some significant local differences. The basins are 200+ meters thick in the Kamari and Perissa area (southeastern Santorini) and >400 m thick in the central (Fira-Imerovigli) area, while there is a thick (∼70-130 m) pyroclastic layer on top of a metamorphic basement “plateau” in the Megalochori area.

4D reconstruction of the Doldenhorn nappe-basement system in the Aar massif: Insights into late-stage continent-continent collision in the Swiss Alps

The inversion of passive margins and their transportation into fold-and-thrust belts is a critical stage of mountain-building processes. In this study, we selected the Doldenhorn Nappe and Aar/Gastern Massifs (Central Swiss Alps) system as an ideal laboratory to document the impact of inherited structures, their along-strike variations, and basement tectonics on the evolution of a fold-and-thrust belt. Three-dimensional geological modelling and cross-section restoration allowed us to reconstruct the 4D evolution of the investigated area during the late-stage Alpine orogeny (30 to 0 Ma). Our results demonstrate that: (i) the Doldenhorn Nappe is the product of the inversion of an asymmetric half-graben basin; (ii) variations in incipient basin sediment thicknesses correlate directly with the along-strike variation of the deformation of the Doldenhorn Nappe; and (iii) the multiphase thick-skinned deformation that overprinted the Doldenhorn Nappe from 22 Ma until today changed the shape of the Doldenhorn Nappe and Aar/Gastern Massifs. This reconstruction shows how thin-skinned nappe formation mechanisms and the nappe geometries are controlled by the initial basement geometry and by the rheological strength contrasts between basement and cover sediments. Basement-involved uplift and shortening controls then the late-stage collisional 3D overprint and mechanics of the fold-and-thrust belts.

Analysis of high-resolution airborne-magnetic data for hydrocarbon generation and preservation potential evaluation of Yola sub-basins, northern Benue Trough, northeastern Nigeria

High-resolution airborne-magnetic data acquired over the Yola extension of the northern Benue Trough, northeastern Nigeria, were analyzed to delineate rock distribution, define basement geometry, determine sedimentary thicknesses, and highlight structural framework among others for the hydrocarbon prospect evaluation of sub-basins in the extension. Filtered and refined magnetic intensity data were subjected to Analytic Signal, Horizontal and Vertical Tilt Derivatives and Euler deconvolution depth weighting analyses to enhance various geological features requisite for petroleum prospect evaluation. Sedimentary rock units flanked the northeast and south widely occur in the central part with volcanic rock intrusions delineated within the sedimentary terrain. Extracted lineaments (less than 5 km and up to 45 km) were observed to dominantly follow the established Benue Trough and Yola sub-basin evolutionary trend along the NE–SW and ESE–WNW directions, respectively. Three prominent sub-basins with sedimentary thicknesses in excess of 3 km and width of more than 60 km, separated by basement highs of mostly volcanic intrusions origin, were identified. The occurrence of relatively thick and extensive sub-basins which are underlain by structurally active basement rocks suggests relatively high geothermal gradient for generation of hydrocarbon from organic matter rich source rocks while the generated oil and gas are preserved within intra-sedimentary structural traps created when the basement fractures were reactivated.

Crustal structure on west-central Baja California peninsula, México, inferred by 3D linear inverse modeling of gravity data

Gravity data from the west-central portion of the Baja California Peninsula, in NW Mexico, were processed on the ellipsoidal reference frame to obtain isostatic residual gravity anomalies. This is a nearly flat region close to sea level and blanked by sediments, where a gravity high bounded by two prominent minima occurs. The gravity high is oriented in a NE-SW direction, and it is possibly the most prominent gravity anomaly transverse to the Baja California Peninsula. The gravity minima, located towards the NW and SE from the gravity high, are respectively related with the Sebastián Vizcaíno and Purísima sedimentary basins. These basins were developed since Early Cretaceous times during the subduction of the Farallon Plate under the North American plate. To infer the crustal structure and the basement geometry in this region, inverse linear modeling is carried on a multilevel data set of upward continued isostatic residual gravity data. Assuming that compaction is the primary diagenetic process on the sediments filling the basins, the initial model is constrained by the surface geology expressed as a density contrast that decreases with depth. The outcome of the inversion is density contrast in a cuboid model, constituted by 67,379 prisms, and its three-dimensional representation suggests that these two basins accommodate as much as five km-thick sedimentary sequences –formed by shallow siliciclastic and deep volcaniclastic rocks penetrated by wells drilled in the region– lying above the inferred basement, which shallows beneath the Lagunitas gravity high, consistently with depth to crystalline rocks penetrated by wells.

Iterative 2D Gravity Modelling for Basement Depth Estimation

In gravity modelling, iterative refinement of 2D basement geometry is a classical problem where the density contrast has to be determined in advance. We present a relatively simple method to estimate 2D basement depth along with its density contrast. Known basement depths at several points and their associated gravity anomalies determine the gravity - depth relationship that is used to obtain the first estimate of the overall basement depth. The initial model response is calculated with a unitary density. The linear regression between calculated model response and actual gravity data results in the density contrast estimation. The misfit or residuals and the gravity - depth relationship are used to update the model iteratively. The proposed technique is relatively stable and results in a recovered model quite similar to the synthetic model. We also explore the possibility to include estimation of the constant the regional anomaly.

Two-Dimensional Gravity Inversion of Basement Relief for Geothermal Energy Potentials at the Harrat Rahat Volcanic Field, Saudi Arabia, Using Particle Swarm Optimization

We invert gravity and magnetic anomalies for basement relief at the Harrat Rahat Volcanic Field (HRVF) for the purpose of evaluating its geothermal energy prospects. HRVF is dominated by basaltic scoria cones and other volcanic rocks overlying the Proterozoic basement. The area considered for this study is located within the northern HRVF and consists mainly of alkali basalts with lesser amounts of benmoreite, mugearite, hawaiite, and trachyte. Our approach adopts a global optimization technique using Particle Swarm Optimization with automated parameter selection, and a two-dimensional gravity-magnetic (GM) forward modeling procedure. The results of the PSO-based approach indicate a depth to the basement at 0.10–624 m, with greater depths within the central region of a solitary anomalous density body in the HRVF. The obtained basement geometry is corroborated by the depth estimates obtained from other potential field inversion methods. The regions with higher prospects are mapped for a targeted future geothermal energy exploration at the HRVF, based on our inversion results.

Imaging the Kanto Basin seismic basement with earthquake and noise autocorrelation functions

SUMMARYSedimentary basins can strongly amplify seismic waves from earthquakes. To better predict future ground motions, detailed knowledge of the sediment thickness and internal structure of basins is required. We image the sediment-to-bedrock interface of the Kanto Basin in Japan using the P-wave reflectivity response from earthquake and ambient seismic noise autocorrelation functions (ACFs) at 286 shallow borehole stations. Earthquake ACFs are computed using P-wave records from 50 Mw 6+ teleseismic events. Noise ACFs are obtained using 1 month of continuous data. Both methods are used to retrieve P-wave traveltimes between the surface and the bedrock interface and map the basin basement geometry. Our prediction of the basement depth agrees generally well with that from a reference velocity model, except for smoother variations in the central part of the basin. Using full-wavefield simulations, we show that the nature of the autocorrelated wavefield has a significant impact on the shape of the ACF waveforms and that earthquake ACFs yield more accurate results in the Kanto Basin.

The Karla impact structure (Russia) explored by potential‐field investigations

AbstractWith no significant topographic expression and limited bedrock exposure, the ~10 km diameter Karla impact structure (Tatarstan, Russia) is poorly known. The age of the impact is also poorly constrained stratigraphically to between 4 and 60 Ma, even if an upper Miocene age is more likely. Targeted gravity and magnetic field surveys were conducted over Karla to explore its size and structure in 2019. Bouguer gravity anomaly data reveal a central positive (+2 mGal) peak ~2 km in diameter surrounded by a concentric negative (−1 mGal) anomaly extending to ~3 km radius; a more irregular, outward‐decreasing (+1 to −1 mGal) positive anomaly extends to 6–8 km radius. A complex impact structure with diameter of 8–10 km is consistent with the Bouguer anomalies. Magnetic field data show 1 to several km‐wavelength anomalies with amplitude variation from +150 to −150 nT and little concentric structure, although the impact feature broadly corresponds to a magnetic low with a weak central high. A 2‐D numerical model of the structure was built using these potential‐field data and petrophysical properties measured on collected samples. It confirms a central uplift composed of Paleozoic sediments and Archean crystalline basement up to 1 km of depth. A 500 m deep collapsed disruption cavity filled by breccia and lacustrine deposits accounts for the Bouguer negative ring. The reversely polarized and weak central magnetic anomalies are controlled by the geometry of the crystalline basement associated with the deformation during the central uplift.

Basement structure investigation using 3-D forward modeling and inversion of geomagnetic data of the Zeit basin area, Gulf of Suez, Egypt

The Zeit basin is located on the western coast of the Gulf of Suez. It is a major, elongated half-graben structure with faults dipping steeply to the southwest. The Zeit basin contains the thickest Miocene sediments and evaporites in the Gulf of Suez. Many exploration activities have been conducted in and around the study area due to the known presence of hydrocarbon resources. In petroleum exploration, the estimation of basement depth (or sedimentary thickness) is a primary exploration risk parameter, in particular, for estimation of heat flow. The magnetic method is one of the few methods that can be used to study the basement structure and estimate the top basement depth. In this paper, the top basement geometry, distribution of the magnetic susceptibility and the basement depth of the Zeit basin were investigated using 3-D forward magnetic modeling and the 3-D inversion of geo-magnetic data. The 3-D forward modeling of the aeromagnetic field data from the Zeit basin area revealed a NW-SE trending half-graben structure, consistent with the known geology of the study area. Magnetic susceptibility cutoff 0.005 SI was applied to show the subsurface relationships among the basement sedimentary units. The calculated depth of the Zeit basin is around 4–5 km. The magnetic field data, calculated using the Gauss Legendre quadrature method, generally agree with the measured data as well as data calculated using the inversion process. The depth and magnetization intensity values used in the 3-D forward modeling, as well as petroleum borehole data, provide a structural basis for the 3-D inversion process.

Unwrapping reworked crust at the Columbia supercontinent margin within central southern Amazon Craton using multi-source geophysics and geochronology data synergy

The geometry and evolution of pre-existing basement in accretionary belts bordering supercontinents are often unclear. Integrative interpretation of long-wavelength potential field satellite data can image deep crust structure, improving our understanding of lithospheric processes that formed these margins bottom-up. Here, we present a multidisciplinary interpretation of the lithospheric architecture of the central southern Amazon Craton, a fragment of an accretionary belt at the southwestern Columbia supercontinent margin. Satellite-borne gravity and magnetic data, airborne magnetic data, passive seismic (Vp/Vs ratio, crustal thickness) and seismic tomography data reveals that basement terranes from the interior of the craton extend into the accretionary margin of Columbia. We demonstrate a vertically heterogeneous structure with an underlying strongly reworked pre-Columbia tectonic wedge that sustained prolonged modification during the supercontinent assembly as corroborated by Nd isotope and geochronology data. Nd isotope data suggest that the protracted orogenic wedge was influenced by subduction angle shifts over time, including addition of substantial juvenile material during slab retreat events. This interplay promoted Craton growth at the supercontinent margin while keeping a subtle record of the pre-existing framework. Our findings point to the possible misrepresentation of basement extension and geometry of supercontinent margins elsewhere.

Basement-decoupled hyperextension rifting: The tectono-stratigraphic record of the salt-rich Pyrenean necking zone (Arzacq Basin, SW France)

AbstractHalf grabens and supra-detachment basins correspond to end-member basin types of magma-poor rift settings, each of them showing a characteristic stratigraphic architecture. The occurrence of a basement-cover décollement has been shown to drastically change the stratigraphic architecture of half graben basins, however, the effect of such basement-cover décollement remains to be documented in supra-detachment basins formed during hyper-extension. We investigate the tectono-stratigraphic record of the Arzacq Basin (SW France) recording the formation of a salt-rich Cretaceous hyperextended rift system. Combining 2-D and 3-D seismic reflection calibrated from well data, we show that this basin is an asymmetric syn-rift extensional syncline growing above a pre-kinematic salt layer. By mapping the sub-salt basement, we show that the formation of this syncline is controlled by the South-Arzacq Fault (SAF), soling in the sub-salt basement. Based on crosscutting relationships and the observed southward migration of syn-rift depocenters, this N110°-striking, 20°-dipping structure accommodates >10 km of thick-skinned extension. The overlying supra-salt cover coherently glided, following the basement geometry. The 3-D segmentation of the SAF and the sub-salt stratigraphic architecture of the Arzacq Basin suggest a roughly dip-slip kinematic. A post-kinematic kilometer-scale uplift is documented on the southern side of the Arzacq Basin. It may result from the increasing lithospheric thinning and thermal support at the end of asymmetric hyperextension. As salt commonly occurs in extensional settings, we believe that our description of the tectono-stratigraphic record of a basement-decoupled supra-detachment basin has global applicability to unleash the tectono-stratigraphic evolution of worldwide hyper-extended rifted margins.

2-dimensional modeling of high-resolution aeromagnetic data over Yola arm of the upper Benue trough, Northeastern Nigeria

2-D modeling of high-resolution aeromagnetic data over the area investigated was carried out with the aid of Oasis MontajTM software in order to get the total magnetic intensity map as well as the residual map from where the modeling was carried out. The data was acquired for Nigeria Geological Survey Agency (NGSA) in 2010 by Fugro international of Netherlands. The specifications of the data were terrain clearance of 80 m, flight line spacing of 500 m and a tie line spacing of 5000 m. The area investigated in this work is located between longitude 120 00' - 130 30' E and Latitude 80 00' – 100 00'N. GM-SYS module of the Oasis montaj was used for this research, with depocentres, structural patterns, basement geometry and morphology as well as estimates of the physical character of the cause to an investigated anomaly determined. The models reveal the horst and graben structures of the basement and the various faults that segmented the area into block patterns. The results revealed sediment thickness in the range of 2-5 km for the models which are in agreement with the results obtained by various researchers in the area. Horst and graben were identified in the models which confirmed the rifting nature of the Benue Trough. The authors of this research are of the view that these depocentres, structural patterns and basement geometry may be a promising site for hydrocarbon traps and reservoirs.
 Keywords: Aeromagnetic Data, GM-SYS, Modeling, Horst, Graben.

Impact of simplifications on numerical modelling of the shallow subsurface at city-scale and implications for shallow geothermal potential

Anthropogenic infrastructures in the shallow subsurface, such as heated basements, tunnels or shallow geothermal systems, are known to increase ground temperatures, particularly in urban areas. Numerical modelling helps inform on the extent of thermal influence of such structures, and its potential uses. Realistic modelling of the subsurface is often computationally costly and requires large amounts of data which is often not readily available, necessitating the use of modelling simplifications. This work presents a case-study on the city centre of Cardiff, UK, for which high resolution data is available, and compares modelling results when three key modelling components (namely ground elevation, hydraulic gradient distribution and basement geometry) are implemented either ‘realistically’, i.e. with high resolution data, or ‘simplified’, utilising commonly accepted modelling assumptions. Results are presented at a point (local) scale and at a domain (aggregate) scale to investigate the impacts such simplifications have on model outputs for different purposes. Comparison to measured data at individual locations shows that the accuracy of temperature outputs from numerical models is largely insensitive to simplification of the hydraulic gradient distribution implemented, while changes in basement geometry affect accuracy of the mean temperature predicted at a point by as much as 3.5 °C. At the domain scale, ground temperatures within the first 20 m show a notable increase (approximately 1 °C volume-averaged and 0.5 °C surface-averaged), while the average heat flux over the domain is about 0.06 W/m2 at 20 m depth. These increased temperatures result in beneficial conditions for shallow geothermal utilisation, producing drilling cost savings of around £1700 per typical household system or about 9% increase in thermal energy potential. Simplifications of basement geometry and (to a lesser degree) the hydraulics can result in an overestimation of these temperatures and therefore over-predict geothermal potential, while the elevation simplification showed little impact.