Isostatic Response Research Articles

Constraining the origin of the Norwegian strandflat – The influence of isostatic and dynamic surface changes

The Norwegian strandflat is a prominent low-relief bedrock surface found near sea level along most of the west coast of Norway. Its origin has been discussed throughout the last 130 years but is yet to be resolved. Some studies suggest that the strandflat represent a tropical weathering front of Mesozoic age that has since been buried and re-exhumed, while others relate its origin to Pleistocene periglacial and glacial processes and/or wave-induced weathering and erosion. Previous interpretations of the strandflat have considered postglacial isostatic uplift, but the impacts of isostatic changes due to glacial erosion and deposition, as well as dynamic surface changes driven by mantle convection, have been largely overlooked. Here we examine how geomorphological-driven isostatic changes and dynamic surface changes have influenced the land-surface elevation along the Norwegian coast during late Pliocene-Quaternary (the last ca. 3 million years). We employ quantitative estimates of glacial erosion and deposition to assess the flexural isostatic response from the resulting load changes. Our analyses show that patterns of geomorphic isostatic adjustments and dynamic surface changes are generally not reflected in the present elevation of the strandflat. Only the loading effect from the deposition of the North Sea Fan can clearly be correlated with the submerged strandflat found near Stad (~62 °N). Our results imply that if the strandflat formed synchronously along the Norwegian coast as a flat surface at sea level, the strandflat we observe today must have developed after the majority of late Pliocene-Quaternary glacial erosion took place, but prior to the main deposition of the North Sea Fan. This would place strandflat formation within the last few glacial cycles, but before the Last Glacial Maximum (LGM). This inferred pre-LGM age of the strandflat is generally consistent with cosmogenic nuclide exposure ages and observed striations on the strandflat. Finally, we examine ice cover and land-surface changes relative to sea level during the last 80,000 years and find no extended periods favorable for synchronous strandflat formation across all regions along the Norwegian coast. This implies that either the strandflat is diachronous, or that the processes of formation have either been extremely fast under certain conditions or are independent of sea level, for instance related to glacial erosion.

Flexural isostatic response of continental-scale deltas to climatically driven sea level changes

Abstract. The interplay between climate-forced sea level change, erosional and depositional processes, and flexural isostasy in deep time on passive margin deltas remains poorly understood. We performed a series of conceptual simulations to investigate flexural isostatic responses to high-frequency fluctuations in water and sediment load associated with climatically driven sea level changes. We model a large drainage basin that discharges to a continental margin and produces a large deltaic depocenter, then prescribe synthetic and climatic-driven sea level curves of different frequencies to assess flexural response. Results show that flexural isostatic responses are bidirectional over 100–1000 kyr timescales and are in sync with the magnitude, frequency, and direction of sea level fluctuations and that isostatic adjustments play an important role in driving along-strike and cross-shelf river mouth migration and sediment accumulation. Our findings demonstrate that climate-forced sea level changes produce a feedback mechanism that results in self-sustaining creation of accommodation into which sediment is deposited and plays a major role in delta morphology and stratigraphic architecture.

Tectonics or rebound: Pleistocene fault reactivation in the highest mountains of the Carpathians

The steep morphology, highest elevation within the Carpathians, and distinct fault borders of the Tatra Mts. support the assumption of recent tectonic activity. However, for decades, the Quaternary tectonic activity of the Tatra Mts. has been poorly evidenced, and a late Pleistocene fault rupture was discovered recently. Using the protected environment of the caves, we have utilized 230Th/U dating of damaged speleothems to extend the record of Quaternary deformation of the Tatra Mts. up to 0.5 Ma. The results from ten caves reveal five periods of increased activity at 465–332, 280–260, 204.8–188.8, 127–86.6, and 29.5–10.11 ka ago. Three mechanisms of fault reactivation were identified through paleostress analysis of cave passage offsets: (1) gravitational sliding; (2) NNE–SSW transpression, likely associated with the Alps–Carpathians–Pannonian (ALCAPA) NNE motion; and (3) tectonic-driven SE-trending compression that reactivated the strike–slip and the gently dipping oblique lower-order faults. The geochronological data indicate that these three stress regimes operated during the late Pleistocene; however, our fault-slip data do not allow us to identify the dominant process. Herein, we suggest that the fault slip was caused by the isostatic response of the Tatra block to the unloading of mountain glaciers, coupled with sediment evacuation overlapping with regional-scale tectonic processes. These processes are locally obliterated by debuttressing. We also show that the distances between the caves and possible seismic sources (i.e., the Sub-Tatric and Ružbachy faults) are sufficiently short to break speleothems in response to oscillation. Given that for most caves an earthquake of >Mw6.5 would be destructive, it seems co-seismic deformations in caves are highly likely considering the Mw > 7 potential of the Sub-Tatric Fault.

Gravity evidence for a heterogeneous crust of Mercury

We modeled gravity data to explore Mercury’s internal structure and show the presence of crustal heterogeneities in density. We first evaluated the lithospheric flexure occurring in the spherical harmonic degree range 5–80, according to the flexural isostatic response curve. We thus estimated a mean elastic lithosphere thickness of about 30 pm 10 km and modeled the crust-mantle interface, which varies from 19 to 42 km depth, according to a flexural compensation model. The isostatic gravity anomalies were then obtained as the residual field with respect to the contributions from topography and lithospheric flexure. Isostatic anomalies are mainly related to density variations in the crust: gravity highs mostly correspond to large-impact basins suggesting intra-crustal magmatic intrusions as the main origin of these anomalies. Isostatic gravity lows prevail, instead, above intercrater plains and may represent the signature of a heavily fractured crust.

Tectonic influence of Cenozoic post-rift alkaline magmatic intrusions on the onshore basement of the Campos and Santos basins, SE Brazil

This study focuses on the tectonic control of post-rift Eocene alkaline intrusions that outcrop in the Atalaia Peninsula, SE Brazil, on the onshore limit of both the Campos and Santos basins, the most prolific Brazilian oil basins. To investigate these tectonic-magmatic processes, high-resolution unmanned aerial vehicle (UAV) images, structural field data, and geochronological analyses were performed. This study characterized several alkaline lithotypes (syenites, N–S and NW–SE trachytes, NW–SE lamprophyres, and E–W phonolites dikes), cutting a N–S oriented Paleoproterozoic orthogneiss basement and NE–SW tholeiitic dikes. A new Ar–Ar age of 53.82 ± 0.7 Ma was obtained on phlogopite from one lamprophyre dike. Geometric and dynamic analysis indicate that these alkaline intrusions are synchronous and their emplacements can be explained by a sub-horizontal NNE–SSW intraplate extension linked to a sub-vertical σ1. Magma rising occurred by forceful injection through transcrustal sets of discontinuities that reached the mantle. Strike to oblique-slip and late normal faults support an anticlockwise rotational extension direction after emplacement. These alkaline tectonic-magmatic dynamics may have caused an isostatic response, influencing the evolution of the abutting offshore basins while triggering a large influx of sediments during the Paleocene-Eocene.

Late Oligocene Orogen‐Scale Tilting in Northern Tibet: A Response to Northward Injection of the Tibetan Lower Crust?

AbstractThe East Kunlun Shan (EKLS) in northern Tibet occupies boundaries of the low‐relief topography and lower crustal low‐velocity zone in the interior plateau, making it ideal for exploring the relationship of surface deformation with underlying geodynamic processes. We used previous and new apatite (U‐Th)/He data to analyze the exhumation history and pattern throughout the EKLS and link surface deformation to deep structures. Integrated (U‐Th)/He ages reveal the rapid exhumation at 27–25 Ma, due to the coeval orogen‐scale tilting of the EKLS. Along with the crustal structures beneath the EKLS, it is inferred that orogen‐scale tilting is the isostatic response of the nonuniform crustal thickening related to northward injection of the Tibetan lower crust. This study highlights the role of ductile deformation within the lower crust in mountain building in northern Tibet, which shares a similarity with mountain building pattern in the eastern plateau margin.

Flexural Response to Erosional Unloading of Continental Margins: An Example From the Bering Sea, USA

AbstractThe geodynamic response of continental margins to erosion may reveal current and past changes in the location of the continental slope. We model the three‐dimensional isostatic response to erosion along the retreating Beringian Margin to test the hypothesis that the margin and its canyons were dominantly eroded during the Pleistocene. Our results predict 900 m of flexural uplift, though evidence for this magnitude of uplift is not observed in Pleistocene and younger deposits, suggesting a substantial amount of margin erosion occurred pre‐Pleistocene. Uplift of this magnitude may be represented by a Late Miocene unconformity across the Bering Shelf and upwarped subsurface units at the base of the Beringian Margin. The distribution of uplift from varying canyon incision and margin retreat scenarios has implications for the extent of geodynamic response during erosional periods.

The influence of lithospheric rheology and surface processes on the preservation of escarpments at rifted margins

The formation of rift flank topography along rifted margins is related to the flexural uplift of the footwall on the rift border fault during lithospheric stretching. The preservation of rift flank topography in mature margins can be explained as a feedback mechanism between differential denudation of the margin and regional isostatic response of the lithosphere. However, the contribution of the lithospheric rifting on the escarpment amplitude in both young and mature margins is not easily quantified. In the present work, we used a 2D thermo-mechanical numerical model to simulate lithospheric extension to evaluate the sensitivity of escarpment amplitude and its preservation under different rheological conditions over timescales larger than 100 Myr, comparable with the age of mature rifted margins. We found that the evolution of escarpment uplift and its preservation for tens of millions of years are sensitive to the degree of coupling between the crust and mantle and the regional isostatic response of the lithosphere to surface processes. The flexural uplift as a response to lithospheric stretching are limited to the first 100‐150 km from the rift flank. Without the influence of surface processes, the amplitude of the escarpment monotonically decreases through time due to the lateral flow of the lower crust and thermal cooling of the margin. In the scenarios with surface processes, topographic signature inherited from the rift phase in the post-rift escarpment elevation is possible only in margins where the erosive escarpment retreat is smaller than ∼100 km. In scenarios where the escarpment retreat exceeds this distance, the escarpment elevation is mainly preserved by the combination of the inherited topography (existent before the onset of lithospheric rifting) and flexural rebound of the lithosphere due to the erosion along the margin.

Computational Geodynamics of South American Plate: Contributions and Perspectives

Computational geodynamics is a key tool in Earth Sciences, enabling the simulation of different geodynamic processes at any time scale in Nature or those which cannot be adequately reproduced by analogue modelling in laboratorial conditions. Specifically, the coupling of surface processes of erosion and sedimentation with the internal dynamics of the lithosphere is a complex problem that involves the solution of a set of differential equations that are only adequately solved by numerical models. During the last three decades, different numerical models were developed to explain the importance of the coupling between the surface and internal dynamics of the Earth, both in active margins and in stable tectonic domains, showing how the coupling leads to counter-intuitive results not observed when each process is analyzed separately. One example of this complex coupling is the feedback mechanism between erosion of the landscape and the regional isostatic response of the lithosphere. In this work, we present simple isostatic and flexural elements that highlight the importance of surface processes on the stress and strain pattern in lithospheric plates. Initially, we present a review on the development of computational geodynamics at University of São Paulo. This review is followed by an analysis of the density structure of the Earth and how the high-density contrast at the Earth’s surface creates a major impact on the isostatic equilibrium of the lithosphere when variations on topographic loads are taken into account. Additionally, we show that the wavelength of the denudation of the Earth’s surface due to fluvial dynamics corresponds to the characteristic length scale for flexural bending of the lithosphere, maximizing the flexural stresses in the lithosphere. Finally, we present recent works on the coupling between surface and lithospheric processes and future challenges for the development of computational geodynamics, with possible strategies to solve them.

Total isostatic response to the complete unloading of the Greenland and Antarctic Ice Sheets

The land surface beneath the Greenland and Antarctic Ice Sheets is isostatically suppressed by the mass of the overlying ice. Accurate computation of the land elevation in the absence of ice is important when considering, for example, regional geodynamics, geomorphology, and ice sheet behaviour. Here, we use contemporary compilations of ice thickness and lithospheric effective elastic thickness to calculate the fully re-equilibrated isostatic response of the solid Earth to the complete removal of the Greenland and Antarctic Ice Sheets. We use an elastic plate flexure model to compute the isostatic response to the unloading of the modern ice sheet loads, and a self-gravitating viscoelastic Earth model to make an adjustment for the remaining isostatic disequilibrium driven by ice mass loss since the Last Glacial Maximum. Feedbacks arising from water loading in areas situated below sea level after ice sheet removal are also taken into account. In addition, we quantify the uncertainties in the total isostatic response associated with a range of elastic and viscoelastic Earth properties. We find that the maximum change in bed elevation following full re-equilibration occurs over the centre of the landmasses and is +783 m in Greenland and +936 m in Antarctica. By contrast, areas around the ice margins experience up to 123 m of lowering due to a combination of sea level rise, peripheral bulge collapse, and water loading. The computed isostatic response fields are openly accessible and have a number of applications for studying regional geodynamics, landscape evolution, cryosphere dynamics, and relative sea level change.

Building a Young Mountain Range: Insight Into the Growth of the Greater Caucasus Mountains From Detrital Zircon (U‐Th)/He Thermochronology and 10Be Erosion Rates

AbstractThe Greater Caucasus (GC) Mountains within the central Arabia‐Eurasia collision zone are an archetypal example of a young collisional orogen. However, the mechanisms driving rock uplift and forming the topography of the range are controversial, with recent provocative suggestions that uplift of the western GC is strongly influenced by an isostatic response to slab detachment, whereas the eastern half has grown through shortening and crustal thickening. Testing this hypothesis is challenging because records of exhumation rates mostly come from the western GC, where slab detachment may have occurred. To address this data gap, we report 623 new, paired zircon U‐Pb and (U‐Th)/He ages from seven different modern river sediments, spanning a ∼400 km long gap in bedrock thermochronometer data. We synthesize these with prior bedrock thermochronometer data, recent catchment averaged 10Be cosmogenic exhumation rates, topographic analyses, structural observations, and plate reconstructions to evaluate the mechanisms growing the GC topography. We find no evidence of major differences in rates, timing of onset of cooling, or total amounts of exhumation across the possible slab edge, inconsistent with previous suggestions of heterogeneous drivers for exhumation along‐strike. Comparison of exhumation across timescales highlight a potential acceleration, but one that appears to suggest a consistent northward shift of the locus of more rapid exhumation. Integration of these new datasets with simple models of orogenic growth suggest that the gross topography of the GC is explainable with traditional models of accretion, thickening, and uplift and does not require any additional slab‐related mechanisms.

Oceanic isostasy as a trigger for the rift-to-drift transition

AbstractA long-standing missing link in our understanding of the Wilson cycle is how a continental rift transitions to seafloor spreading. The variety of rift structures and transition timings at rift margins do not easily lend themselves to some specific degree of strain and/or magmatism as the tipping point. Invariably ignored in the process, but a potential key to the conundrum, is the isostatic response that comes with ocean loading during and after inundation. Ocean mass redistribution on variably subsiding crust drives flow in the asthenosphere in much the same way a growing icecap drives a corresponding outward mantle flow. This flow alters mantle tractions of the rift system, with disappearance of basal resistance, and even adds a push to the rifting process. Evidence for ocean inundation facilitating self-sustained seafloor spreading is observed in the Atlantic, around the Afar triple junction, and elsewhere, indicating that the ocean should not be considered simply incidental to the creation of oceanic basins.

Flexural subsidence modelling of post-rift paleobathymetry and sedimentary infill in the northern South China Sea margin

The South China Sea (SCS) is the largest marginal sea on the SE Asian continental margin, located at the junction of the Pacific, Indo-Australia and Eurasian plates, and has undergone the joint effect of the Tethys and Pacific tectonic domains. Currently, the SCS remains debatable for its evolution and origin. However, its northern margin preserves the Cenozoic tectonic and subsidence history, which allow us to restore the past continental margin (0–33.7 Ma), to estimate subsidence, thickness and sediment accumulation correlated with regional tectonics. This paper has employed a 2D flexural backstripping technique and spatio-temporal variations by enhancing the seismic profiles with flexurally interpolated modelled sections. According to our findings, the transition was mainly controlled by three structural stages. (1) Early-Late Oligocene: records the highest subsidence due to the initiation of seafloor spreading, isostatic response of continental breakup, associated uplifts therein, and retreat of the Pacific Plate. (2) Early-Middle Miocene: mature spreading attributed to a slow phase of subsidence during the Early Miocene and a rapid increase of subsidence associated with the closure of the SCS in the Middle Miocene. (3) Thermal subsidence: sediments shifted since the Late Miocene towards mechanically weak and stretched regions (central-southern and eastern) which were triggered by a gravitational pull force. A fault system and differential stretched crust relocated sediment to the regional stress directions of the Pacific Plate.

‘Block and basin’ style rift basins: sedimentological insights from the Mississippian Fell Sandstone Formation

The block and basin tectonostratigraphic framework for the northern Pennine (rift) Basin, within which buoyant granite intrusions core intra-basin fault-bounded blocks, has long held traction. However, many of the elements of this framework are rooted in primitive tectonic models and, perhaps unsurprisingly, corresponding depositional models often reflect this. Using sedimentological and sedimentary provenance approaches, the synrift (Mississippian) fluvio-deltaic Fell Sandstone Formation and age-equivalent strata within the northern Pennine Basin are examined. Highlighted divergences from classically depicted models relate to occurrences of pre-Carboniferous basement domes or monoclines, which are unbounded by major vertically displacing (>100 m) fault systems. Such structures in the northern Pennine Basin are all granite-cored and their origins are associated with their buoyancy and flexural isostatic processes. One such basement dome, the Cheviot Block, confined and deflected the Fell Sandstone fluvio-deltaic system from the west, causing locally elevated net sand content and variations in the dominant palaeodrainage direction. Central parts of the Alston Block, which forms a regional monocline along an east–west axis, were comparatively uplifted because of flexural isostatic responses to granite intrusions. The findings presented are at variance not only with classically depicted depositional models for the region, but also with more general depictions of dominantly normal fault-driven rift basin systems. Supplementary material: Tables of data locations with derivation, trace element data and major element (oxide) data are available at https://doi.org/10.6084/m9.figshare.c.5733257

Modelling sea-level fingerprints of glaciated regions with low mantle viscosity

Abstract. Global patterns of sea-level change – often termed “sea-level fingerprints” – associated with future changes in ice/water mass re-distribution are a key component in generating regional sea-level projections. Calculation of these fingerprints is commonly based on the assumption that the isostatic response of the Earth is dominantly elastic on century timescales. While this assumption is accurate for regions underlain by mantle material with viscosity close to that of global average estimates, recent work focusing on the West Antarctic region has shown that this assumption can lead to significant error where the viscosity is significantly lower than typical global average values. Here, we test this assumption for fingerprints associated with glaciers and ice caps. We compare output from a (1D) elastic Earth model to that of a 3D viscoelastic model that includes low-viscosity mantle in three glaciated regions: Alaska, southwestern Canada, and the southern Andes (Randolph Glacier Inventory (RGI) regions 1, 2, and 17, respectively). This comparison indicates that the error incurred by ignoring the non-elastic response is of the order of 1 mm in most areas (or about 1 % of the barystatic signal) over the 21st century with values reaching the centimetre level in glaciated regions. However, in glaciated regions underlain by low-viscosity mantle, the non-elastic deformation can result in relative sea-level changes with magnitudes of up to several tens of centimetres (or several times the barystatic value). The magnitude and spatial pattern of this non-elastic signal is sensitive to variations in both the projected ice history and regional viscosity structure, indicating the need for loading models with high spatial resolution and improved constraints on regional Earth viscosity structure to accurately simulate sea-level fingerprints in these regions. The anomalously low mantle viscosity in these regions also amplifies the glacial isostatic adjustment signal associated with glacier changes during the 20th century, causing it to be an important (and even dominant) contributor to the modelled relative sea-level changes over the 21st century.

Late Quaternary sea level, isostatic response, and sediment dispersal along the Queen Charlotte fault

Abstract The active Pacific margin of the Haida Gwaii and southeast Alaska has been subject to vigorous storm activity, dramatic sea-level change, and active tectonism since glacial times. Glaciation was minimal along the western shelf margin, except for large ice streams that formed glacial valleys to the shelf break between the major islands of southeast Alaska and Haida Gwaii. Upon deglaciation, sediment discharge was extensive, but it terminated quickly due to rapid glacial retreat and sea-level lowering with the development of a glacio-isostatic forebulge, coupled with eustatic lowering. Glacial sedimentation offshore ended soon after 15.0 ka. The shelf became emergent, with sea level lowering by, and possibly greater than, 175 m. The rapid transgression that followed began sometime before 12.7 ka off Haida Gwaii and 12.0 ka off southeast Alaska, and with the extreme wave-dominated environment, the unconsolidated sediment that was left on the shelf was effectively removed. Temperate carbonate sands make up the few sediment deposits presently found on the shelf. The Queen Charlotte fault, which lies just below the shelf break for most of its length, was extensively gullied during this short period of significant sediment discharge, when sediment was transported though the glacial valleys and across the narrow shelf through fluvial and submarine channels and was deposited offshore as sea level dropped. The Queen Charlotte fault became the western terminus of the glacio-isostatic forebulge, with the fault acting as a hinged flap taking up the uplift and collapse along the fault of 70+ m. This may have resulted in the development of the distinctive fault valley that presently acts as a very linear channel pathway for sediment throughout the fault system.

The Hydro-Isostatic Rebound Related to Megalake Chad (Holocene, Africa): First Numerical Modelling and Significance for Paleo-Shorelines Elevation

Lake Chad, the largest freshwater lake of north-central Africa and one of the largest lakes of Africa, is the relict of a giant Quaternary lake (i.e., Megalake Chad) that developed during the early- to mid-Holocene African Humid Period. Over the drylands of the Sahara Desert and the semi-arid Sahel region, remote sensing (optical satellite imagery and digital elevation models) proved a successful approach to identify the paleo-shorelines of this giant paleo-lake. Here we present the first attempt to estimate the isostatic response of the lithosphere due to Megalake Chad and its impact on the elevation of these paleo-shorelines. For this purpose, we use the open source TABOO software (University of Urbino, Italy) and test four different Earth models, considering different parameters for the lithosphere and the upper mantle, and the spatial distribution of the water mass. We make the simplification of an instantaneous drying-up of Megalake Chad, and compute the readjustment related to this instant unload. Results (i.e., duration, amplitude, and location of the deformation) are then discussed in the light of four key areas of the basin displaying prominent paleo-shoreline morpho-sedimentary features. Whatever the Earth model and simplification involved in the simulations, this work provides a strong first-order evaluation of the impact on hydro-isostasy of Megalake Chad. It demonstrates that a water body similar to this megalake would induce a significant deformation of the lithosphere in the form of a vertical differential uplift at basin-scale reaching up to 16 m in the deepest part of the paleo-lake, and its shorelines would then be deflected from 2 m (southern shorelines) to 12 m (northern shorelines), with a maximum rate of more than 1 cm y−1. As such, any future study related to the paleo-shorelines of Megalake Chad, should integrate such temporal and spatial variation of their elevations.

Modelling Sea-Level Fingerprints of Glaciated Regions with Low Mantle Viscosity

Abstract. Sea-level fingerprints define the spatially varying relative sea-level response to changes in grounded ice distribution. These fingerprints are a key component in generating regional sea-level projections. Calculation of these fingerprints is commonly based on the assumption that the isostatic response of the Earth is dominantly elastic on century time scales. While this assumption is accurate for regions underlain my mantle material with viscosity close to that of global average estimates, recent work focusing on the Antarctic region has shown that this assumption can led to significant error when the viscosity departs significantly from typical average values. Here we test this assumption for fingerprints associated with glaciers and ice caps. We compare output from a (1D) elastic Earth model to that of a 3D viscoelastic model which includes low viscosity mantle in three glaciated regions: Alaska, southwestern Canada and the southern Andes (Randolph Glacier Inventory (RGI) regions 1, 2 & 17, respectively). This comparison indicates that the error incurred by ignoring the non-elastic response is generally less than 1 cm over the 21st century but can reach magnitudes of up to several 10s of centimetres in low viscosity areas. This error can have large spatial gradients where crustal uplift in ice covered (or previously ice covered) areas changes into subsidence when moving away from the loading centres to areas peripheral to the mass loss. The existence of these large gradients indicates the need for loading models with high spatial resolution to accurately simulate sea-level fingerprints in these regions. We conclude that sea-level projections for Alaska, southwestern Canada and the southern Andes should not be based on elastic Earth models.

Quantification of Erosion and Uplift in a Rising Orogen—A Large-Scale Perspective (Late Tortonian to Present): The Case of the Gibraltar Arc, Betic Cordillera, Southern Spain

The present study deals with the morphometric quantification of erosion and illustrates the uplift component triggered by denudation (isostasy) in the growth and evolution of a rising orogeny by the application of Airy isostasy concepts. The Gibraltar Arc, located in the Western–Central sector of the Betic Cordillera, developed an exceptional geological scenario during the Messinian Salinity Crisis since the thin emerged fringe of the uprising Cordillera disconnected the Atlantic and Mediterranean basins, generating a relevant misbalance and asymmetry in the fluvial erosion between the two slopes of the emergent orogeny. Our analysis was applied to 50 individual drainage basins (spatial isostatic units) in the Western–Central Betic Cordillera, allowing us to obtain individual and bulk estimates for these isostatic parameters. GIS-based numerical estimations were obtained using LiDAR Digital Elevation Models (DEMs) provided by the Spanish Geographical Institute and reconstructed pre-incision surface models obtained from proxy paleo-elevation data, estimated from stratigraphic and geomorphological littoral to shallow marine markers. The obtained values for geophysical relief, denudation plates, erosion/uplift rates and computed accumulated uplift (245–407 ±20 m) are higher for the ancient Mediterranean slope of the orogen. On the contrary, the Atlantic slope presents an accumulated uplift of only 138–236 ±20 m, indicating the strong control of the ancient Messinian Atlantic–Mediterranean water divide. The temporal study of erosion indicates that most of the difference in uplift in the Mediterranean slope was achieved during or soon after the Messinian Salinity Crisis, resulting in mean uplift rates of 0.21 mm/y, but practically null (0.01 mm/y) for the Atlantic slope. The comparison of the geophysical relief models with proxy paleo-elevation data allowed us to assess the current state of the denudation process in the range. The results indicate that, towards the west of the range denudation compensated elevation, and is actively back-feeding isostatic rebound. Therefore, the contribution of external processes to mountain range elevation through isostasy is quantitatively estimated using elevation data. In this case, a relevant part of the surface uplift (50-55%) is undertaken by the orogen. Ultimately, the Messinian Salinity Crisis-related isostatic response to differential denudation may be behind the quaternary westward tilting of Iberia, causing more than 70% of the Peninsula to drain towards the Atlantic.

Implications of glacial isostatic adjustment on petroleum reservoirs in the Grand banks of Newfoundland

Glacial Isostatic Adjustment (GIA), or postglacial rebound refers to the crustal response to glacial loading and unloading processes. In this investigation, we run numerical simulations using different rheological parameters and ice loading histories to obtain models of the glacial isostatic response for the Grand Banks of Newfoundland. The objective is to assess in general the potential impact of GIA on hydrocarbon reservoirs and trapping structures in terms of tilting, deformation, fault activity or geochemical implications on the trapped fluids. Comparisons may be drawn to related studies, including analogous impacts on reservoirs in the Barents Sea, as well as validation of GIA model predictions of vertical crustal motion using GNSS stations, levelling networks and relative sea level (RSL) observations at numerous sites across eastern Canada. By running a series of GIA models, a range of vertical motion rates are output through time to determine displacement amounts, and calculations are performed on these outputs to disseminate other valuable data attributes such as the displacement gradient and differences in time and space. In general, the vertical motion rate is seen to change drastically over time for certain points, as the lag in response of regions near the ice periphery leads to a distinct pattern of uplift and subsidence. The Grand Banks generally experienced between 57 m of subsidence to 34 m of uplift since the last glacial maximum, depending on the specific location and the GIA model parameters. At the reservoir scale, model results along a 10 km baseline show up to −1.0–1.5 m of differential vertical displacement, which could have led to reservoir tilting or even hydrocarbon migration. As a result of this study, it is clear that there are implications of GIA on hydrocarbon reservoirs since the Last Glacial Maximum, and these effects are ongoing.