Fact-checked by Grok 2 weeks ago

Post-glacial rebound

Post-glacial rebound, also known as glacial isostatic adjustment, is the ongoing vertical motion of the Earth's crust resulting from the viscoelastic relaxation of the mantle in response to the melting of massive ice sheets that depressed the lithosphere during the Pleistocene epoch. This process involves the uplift of formerly glaciated interiors, where the crust rebounds after the removal of ice loads equivalent to several kilometers thick, while peripheral forebulges—regions that bulged upward due to displaced mantle material—now subside as the adjustment propagates outward. Initiated primarily after the Last Glacial Maximum around 21,000 years ago, the phenomenon persists with measurable rates, reaching up to 9–12 mm per year near Hudson Bay in Canada and 5–10 mm per year in Fennoscandia. These motions significantly influence regional relative sea levels, with uplift outpacing eustatic rise in core areas to produce apparent sea level fall, while subsidence in forebulge zones like the U.S. East Coast exacerbates local inundation risks. Observations from GPS, tide gauges, and satellite gravimetry confirm the adjustment's role in geodynamic modeling, mantle viscosity inference, and interpreting Holocene paleoenvironments, underscoring its importance for understanding solid Earth-ice sheet interactions.

Fundamental Principles

Isostatic Equilibrium and Viscoelastic Response

Isostasy describes the gravitational equilibrium achieved when the Earth's lithosphere "floats" atop the denser mantle, following Archimedes' principle of buoyancy in which the thickness and density of crustal material determine its equilibrium depth. This balance assumes a layered Earth structure, with the rigid lithosphere overlying the asthenosphere, which exhibits fluid-like behavior over long timescales due to its viscous properties, enabling slow adjustment to perturbations in surface mass distribution. In the context of post-glacial rebound, isostatic disequilibrium arises from the imposition and subsequent removal of massive ice loads, prompting the mantle to redistribute material viscously to restore equilibrium. The initial response to glacial loading involves predominantly elastic deformation of the lithosphere, resulting in near-instantaneous flexural subsidence under the weight of the ice sheet, with the underlying mantle providing buoyant support. Upon deglaciation, the unloading triggers viscoelastic relaxation, where the combined elastic and viscous properties of the mantle lead to time-dependent crustal uplift as displaced mantle material flows back beneath the depressed region. This viscoelastic behavior contrasts with purely elastic rebound, which would be rapid and complete; instead, viscous flow introduces irreversibility and rate-dependence, with deformation persisting until isostatic equilibrium is reattained. The timescale of this viscoelastic response spans millennia, characterized by an initial rapid phase of uplift followed by exponential decay toward equilibrium, governed by the mantle's rheological parameters such as viscosity, typically on the order of $10^{20} to $10^{21} Pa·s in the upper mantle. Relaxation times range from several thousand to over 10,000 years, reflecting the exponential nature of viscous dissipation, where uplift rates diminish as the driving stress from disequilibrium decreases. This process underscores the mantle's capacity for causal adjustment to surface forcing, with the asthenosphere's low effective viscosity facilitating flow while deeper layers contribute to longer-term stability.

Glacial Loading, Unloading, and Mantle Viscosity

During the Last Glacial Maximum, spanning approximately 26,000 to 19,000 years before present, continental ice sheets imposed immense loads on the Earth's lithosphere, with the Laurentide Ice Sheet attaining maximum thicknesses of 3.4–3.6 km over central regions near Hudson Bay, resulting in viscoelastic crustal depressions exceeding 1 km beneath the thickest ice accumulations. These depressions arose from the displacement of mantle material beneath the load, governed by the density contrast between ice (approximately 917 kg/m³) and the underlying mantle (around 3300 kg/m³), which induced flexural subsidence and outward flow in a viscoelastic Earth. Similar loading occurred under the Fennoscandian Ice Sheet, depressing the crust by hundreds of meters to over a kilometer in Scandinavia. Deglaciation commencing around 19,000 years ago progressively unloaded these regions through ice melt, redistributing mass equivalent to vast ocean volumes and initiating the collapse of peripheral forebulges—uplifted zones 500–1000 km beyond the former ice margins, where mantle material had accumulated during loading, reaching heights up to 60–100 m. This unloading reversed the flow of mantle material, allowing the depressed centers to rebound while peripheral areas subsided as forebulges relaxed, a process persisting today due to the time-dependent nature of viscous relaxation. The spatial pattern of mass redistribution, including eustatic sea-level rise from meltwater, further modulated the isostatic response, with ongoing adjustments reflecting incomplete equilibration. The rate and extent of rebound are profoundly influenced by mantle viscosity, empirically constrained by inverting glacial isostatic adjustment histories to yield upper mantle values typically around 0.5–2 × 10^{21} Pa·s, increasing to 10^{22} Pa·s or higher in the lower mantle to match observed relaxation timescales. These radially varying profiles assume a Maxwell viscoelastic rheology, where shorter-term responses probe shallower, lower-viscosity layers, while long-term signals inform deeper structure; however, debates persist over laterally heterogeneous viscosity, with seismic tomography indicating variations tied to subducting slabs or plumes that could locally accelerate or retard flow by factors of 2–10. Such heterogeneities challenge purely one-dimensional models, as three-dimensional simulations reveal amplified effects on rebound patterns in regions like Fennoscandia.

Historical Discovery and Early Observations

Pre-20th Century Evidence from Raised Shorelines

In 19th-century , geologists documented extensive raised shorelines, known as strandlines, preserved at elevations up to 285 above present in , serving as proxies for post-glacial land . These features, observed along former coastlines, exhibited consistent patterns correlating with the extent of past ice coverage, with higher elevations toward the centers of former ice domes. geologist De Geer mapped these uplifts in 1888, producing the first isobase depicting total post-glacial rebound in across , thereby establishing a glacial isostatic origin over alternative explanations like widespread . De Geer's work resolved debates by linking shoreline displacements directly to ice sheet recession, rejecting uniform global sea-level fall in favor of localized crustal response to deglaciation. Similar evidence emerged in Scotland, where raised beaches indicated differential uplift tied to glacial loading. In 1865, Thomas Francis Jamieson analyzed these shorelines, proposing that the land's depression under ice weight followed by rebound explained their elevation, rather than eustatic sea-level changes alone. Jamieson's observations in eastern Scotland highlighted a "Main Postglacial Shoreline" varying in height, with patterns diminishing away from glaciated highlands, supporting isostatic adjustment over tectonic or erosional causes. This interpretation gained traction as field data showed no corresponding global marine regression sufficient to account for the anomalies. Around Hudson Bay, Canadian geologist Robert Bell reported raised beaches during explorations in the 1880s, noting their prevalence along straight coastal stretches and funnel-shaped bays, with elevations evidencing ongoing land rise post-ice retreat. Bell's findings linked these features to the former Laurentide Ice Sheet's influence, observing wave-like uplift progress from peripheral areas inward, consistent with isostatic recovery rather than localized erosion. Across these regions, early hypotheses favoring erosion-driven buoyancy were supplanted by isostatic models, as shoreline gradients aligned precisely with glacial maxima, providing empirical validation without reliance on contemporary tectonics.

Key Milestones and Researchers in the 20th Century

In 1935, Norman A. Haskell introduced an early mathematical framework for modeling post-glacial rebound as a viscoelastic process, applying the correspondence principle between elastic and viscous responses to isostatic adjustment and deriving a mantle viscosity estimate of approximately $10^{21} Pa·s from observed uplift rates in formerly glaciated regions. This approach marked a shift from purely elastic isostasy toward accounting for time-dependent mantle flow beneath the lithosphere. Building on this, Felix Andries Vening Meinesz in 1937 analyzed post-glacial uplift data from Scandinavia to quantify the Earth's plasticity, modeling the adjustment as viscous relaxation of the mantle under a receding ice load and estimating effective viscosities consistent with slow, ongoing deformation rates of several meters per century in the region. These models highlighted the role of a thin elastic lithosphere overlying a viscous asthenosphere in controlling rebound patterns. In 1941, Beno Gutenberg synthesized observations of post-glacial uplift and eustatic sea-level changes, using Fennoscandian and North American to infer mantle viscosities ranging from $10^{21} to $10^{22} ·s, thereby providing empirical constraints on internal and linking rebound rates—up to 1 cm/year in peripheral zones—to glacial unloading timescales of thousands of years. The 1960s and 1970s saw integration of post-glacial rebound models with emerging , enabling researchers to differentiate isostatic viscoelastic responses from rigid plate motions and active rifting; for instance, uplift in stable cratons like the Canadian was attributed primarily to glacial unloading rather than tectonic . A pivotal advancement occurred in the mid-1970s with W. Richard Peltier's development of a global glacial isostatic adjustment (GIA) framework, which solved the forward problem of surface deformation on a radially stratified, viscoelastic spheroid using Green's functions and normal mode analysis, while addressing the inverse problem to infer mantle viscosity structure from sea-level and uplift records. This approach incorporated realistic ice-load histories and the sea-level equation, formalizing rebound as a coupled system of crustal motion, geoid changes, and relative sea-level variations observable worldwide.

Modern Observational Evidence

Vertical Crustal Motion from Tide Gauges and GPS

Tide gauges provide long-term records of relative sea-level change, which combine absolute sea-level variations with vertical crustal motion; positive relative sea-level fall indicates land uplift exceeding sea-level rise. In regions affected by post-glacial rebound, such as Hudson Bay, tide gauge data from Churchill record an uplift component of approximately 9.0 ± 0.8 mm/yr when corrected using satellite altimetry. Uplift rates derived from tide gauges in the Hudson Bay region range from 1 to 10 mm/yr, decreasing radially from the former center of the Laurentide Ice Sheet. Global Positioning System (GPS) measurements offer direct quantification of absolute vertical crustal motion, isolating land uplift from oceanographic effects and correcting for minor tectonic influences. The BIFROST GPS network, operational in Fennoscandia since the 1990s, records maximum uplift rates of about 10 mm/yr near the Gulf of Bothnia, confirming and refining tide gauge inferences while revealing spatial gradients that diminish to near zero at peripheral distances. In Scandinavia, BIFROST data show uplift rates of 5-8 mm/yr in mid-regions, aligning with tide gauge trends but providing higher precision over baselines of 10-20 years. Empirical validation from dated raised shorelines demonstrates temporal decay in rebound rates, with early Holocene uplift exceeding 10 mm/yr—often several cm/yr immediately post-deglaciation—contrasting present-day mm/yr scales due to viscoelastic relaxation. GPS and tide gauge consistency across these regions underscores the ongoing isostatic response, with rates slowing predictably from ice sheet centers outward.

Horizontal Deformation and Strain Patterns

Horizontal deformation during post-glacial rebound manifests as lateral crustal motions driven by the viscoelastic relaxation of the mantle following glacial unloading, inducing intra-plate strain fields observable through geodetic techniques. In central regions overlying former ice sheet domes, such as around Hudson Bay, the vertical uplift causes radial outward divergence as the lithosphere expands laterally. Conversely, in peripheral zones associated with the collapsing forebulge, horizontal velocities converge inward toward the load centers due to subsidence and mantle flow adjustments. Global Positioning System (GPS) observations quantify these motions, revealing horizontal velocities typically ranging from 0.5 to 2 mm/yr, with far-field rates of 1–2 mm/yr directed toward the former Laurentide Ice Sheet in North America. Interferometric Synthetic Aperture Radar (InSAR) complements GPS by mapping strain gradients, confirming divergent patterns near uplift maxima and convergent flows in annular peripheral belts. Resulting strain patterns exhibit extension beneath central domes, with rates up to 4 × 10^{-9} yr^{-1}, and contractional strain in surrounding semiannular zones at similar magnitudes, reflecting the differential response to isostatic disequilibrium. Data from extensive GPS networks, including those archived by UNAVCO, highlight non-uniformity in these patterns, influenced by lateral variations in mantle viscosity that amplify or attenuate local flows beyond predictions from laterally homogeneous models.

Relative Sea Level Changes and Fingerprinting

Relative sea level (RSL) changes due to glacial isostatic adjustment (GIA) encompass not just vertical crustal deformation but also perturbations to the geoid and solid Earth loading from the viscoelastic response to late Pleistocene deglaciation, creating non-uniform spatial patterns distinct from eustatic sea level variations. The characteristic "fingerprint" of GIA arises from self-gravitation, deformation, and rotational effects, whereby the removal of ice-sheet mass redistributes ocean water away from former load centers, lowering local sea surface heights relative to the geoid, while distant regions experience compensatory rises. Near centers of past glaciation, such as Hudson Bay or Fennoscandia, the combined uplift (often 1–10 mm/yr) and gravitational unloading dominate, yielding ongoing RSL falls exceeding global eustatic rise by up to several mm/yr in paleorecords and modern proxies. This fingerprint manifests as amplified RSL decline proximate to ice margins, with rates influenced by mantle viscosity and load history; for instance, modeling indicates gravitational effects alone contribute negative RSL trends of approximately -0.2 to -0.5 mm/yr in peripheral zones around from historical unloading, compounded by current uplift. In contrast, the migration and collapse of the peripheral forebulge—mantle upwelling displaced during loading—induces subsidence in forebulge regions, elevating RSL. Along the U.S. East Coast, this subsidence adds 0.3–0.8 mm/yr to contemporary RSL rise, with higher values (up to ~1 mm/yr) in the mid-Atlantic where forebulge relaxation is most active, as evidenced by tide gauge records corrected for eustatic and steric components. Empirical isolation of the GIA fingerprint in RSL datasets relies on multi-technique integration, including satellite gravimetry from GRACE (operational since 2002) and GOCE (2009–2013), which map time-variable gravity anomalies to constrain viscoelastic models and subtract contemporary mass signals (e.g., ice melt or hydrology). When fused with radar altimetry for absolute sea surface heights and GPS for crustal motion, these data reveal GIA contributions by forward-modeling predicted fingerprints and inverting residuals, achieving uncertainties below 0.5 mm/yr for regional trends while highlighting model-dependent viscosity variations. Such approaches confirm the dominance of GIA in modulating RSL variability, with gravimetric observations validating modeled bulge subsidence patterns against uncorrelated modern forcings.

Ancillary Effects: Gravity Anomalies and Regional Tilt

Gravity anomalies associated with post-glacial rebound arise from the viscoelastic redistribution of mass following glacial unloading, manifesting as measurable perturbations in Earth's gravity field. Satellite gravimetry from the GRACE mission has identified gravity rate changes of -1 to -5 μGal/yr in actively rebounding regions such as Fennoscandia, Svalbard, and southeastern Alaska, where these negative signals primarily reflect the free-air effect of ongoing crustal uplift outweighing subsidiary mass flow contributions. Ground-based absolute gravimetry corroborates these rates, with measurements in Alaska yielding -3.5 to -5.6 μGal/yr over multi-year campaigns, consistent with GIA predictions that account for mantle viscosity and load history. Regional tilt variations, detectable through high-precision leveling surveys, represent another ancillary signal of rebound, as differential uplift gradients alter the orientation of the local plumb line and influence lake or groundwater level slopes. In the Great Lakes basin, for instance, GIA-driven tilting manifests as a hinge line separating uplift from subsidence, with observed level differences implying tilt rates on the order of 10^{-6} rad/yr, derived from water gauge data and geodetic networks spanning the region. These tilts, measured via repeated leveling loops, perturb hydraulic gradients and have been quantified in North American stable interiors, where GIA contributions exceed tectonic noise by factors of several, as evidenced by GPS-augmented surveys since the late 20th century. Post-glacial rebound also perturbs Earth's rotation, inducing secular polar motion and modulating the Chandler wobble through mass redistribution effects on the planet's moments of inertia. The observed 20th-century polar wander trend, at approximately 3.5 milliarcseconds per year toward the Hudson Bay region, is predominantly attributed to GIA from Laurentide Ice Sheet unloading, with astronomical records dating back to 1900 confirming this drift's consistency with viscoelastic models. GIA further influences Chandler wobble parameters, including its period (around 433 days) and amplitude variations, by altering rotational feedback on a non-hydrostatic Earth, as quantified in analyses separating GIA from atmospheric and oceanic excitations over multi-decadal spans.

Regional Variations and Case Studies

Fennoscandia and Baltic Sea Uplift

The retreat of the Weichselian Ice Sheet from , completing around years before present in the region, initiated pronounced post-glacial rebound centered over the former ice load maximum. This process has resulted in differential vertical crustal motion, with the highest rates persisting in the northern at approximately 9-11 mm per year, tapering to lower values southward and becoming negative south of the zero isobase in the southern . These rates, derived from GPS and leveling observations, reflect ongoing viscoelastic relaxation of the mantle beneath the lithosphere. Total cumulative uplift since deglaciation in the central zone reaches approximately 500 meters, as inferred from raised shoreline elevations and modeling constrained by geological records. Isobase mapping, which delineates contours of equal uplift, originated in the 19th century through observations of emerging coastlines and ancient shorelines, evolving into systematic geodetic frameworks by the early 20th century. These maps reveal an elliptical pattern elongated north-south, with the maximum uplift axis aligned along the Gulf of Bothnia. Rebound has profoundly influenced Baltic Sea evolution, notably contributing to the Ancylus Lake phase around 10,700 to 9,800 years ago, when rapid northern uplift elevated drainage sills, transforming the basin into a large freshwater lake isolated from Atlantic inflows. This stage exemplifies how isostatic adjustment interacted with eustatic changes to dictate paleogeography, with subsequent outlet breaching leading to marine transgressions. Ongoing uplift continues to expose new land, particularly in the Kvarken archipelago, at rates exceeding local sea-level rise. Rebound stresses interact with intraplate tectonics, inducing neotectonic activity including postglacial faults primarily in northern Fennoscandia, where compressive stresses from asymmetric unloading promote thrusting. In the Gulf of Finland, differential uplift contributes to subtle extensional features and crustal strain, modulating minor rifting and seismicity amid the broader isostatic signal. These effects highlight how glacial unloading perturbs lithospheric stress fields, though magnitudes remain subordinate to the primary rebound.

Laurentide Ice Sheet Rebound in North America

The , which covered approximately 13 million square kilometers of during the around 21,000 years ago, exerted maximum load over the region, resulting in the highest present-day post-glacial rates there. GPS measurements indicate vertical uplift rates of approximately 10 per year near , decreasing radially outward to near-zero or negative values farther south. This pattern contrasts with the more and extensive uplift in , where rates exceed 10 per year over a broader peripheral due to differences in ice sheet configuration and mantle viscosity estimates. In the U.S. Midwest, peripheral forebulge collapse manifests as subsidence at rates of 1-2 mm per year, primarily driven by viscous relaxation of the mantle following unloading of the former Laurentide load. This subsidence contributes to relative water level rises in the Great Lakes, altering regional hydrology and necessitating adjustments in benchmarks like the International Great Lakes Datum of 2020, which incorporates glacial isostatic adjustment models to account for ongoing crustal motion. Along the eastern margin, such as in Newfoundland, forebulge migration and collapse have produced complex relative sea-level histories, with initial post-glacial emergence followed by transgression due to subsidence. Radiocarbon dating of marine clays and sediments indicates that sea levels fell below present datum around 10,000-12,000 calibrated years before present before rising, reflecting the northerly collapse of the forebulge at rates tied to mantle flow. Recent GPS arrays across North America, including networks in Canada and the U.S., reveal a deceleration in rebound rates compared to early Holocene predictions, with vertical velocities in central regions stabilizing around 8-10 mm per year in the 2010s-2020s, influenced by viscoelastic Earth models. These observations inform hydrological management in the Great Lakes basin, where isostatic tilt exacerbates water level variability beyond eustatic or climatic factors.

Antarctic and Peripheral Bulge Dynamics

In Antarctica, glacial isostatic adjustment (GIA) manifests as differential crustal uplift, with rates in West Antarctica exceeding those in the East due to greater historical and ongoing ice unloading in the former region. GPS observations indicate present-day uplift rates of up to 5 mm yr⁻¹ across much of West Antarctica, particularly in the Amundsen Sea Embayment where marine-based ice sectors have thinned significantly since the Last Glacial Maximum. In contrast, East Antarctic uplift rates are lower, typically 1–2 mm yr⁻¹, reflecting a thicker and more stable ice cover with minimal net mass loss over the Holocene. This asymmetry arises from viscoelastic mantle flow beneath thinner West Antarctic lithosphere, amplifying rebound responses to deglaciation. Peripheral bulge dynamics extend GIA effects beyond Antarctica, causing subsidence in far-field regions such as southern South America and Australia, where forebulge collapse continues from the Pleistocene unloading. In these areas, subsidence rates reach 0.5–1 mm yr⁻¹, driven by the inward relaxation of the erstwhile peripheral bulge that formed under glacial loading. Three-dimensional Earth models predict enhanced subsidence in these peripheral zones compared to one-dimensional approximations, as lateral viscosity variations influence bulge migration and decay. Contemporary ice mass fluctuations superimpose on long-term GIA signals, notably in East Antarctica where satellite gravimetry detected a mass gain of approximately 130 Gt yr⁻¹ from 2021 to 2022, primarily from anomalous precipitation. This gain induces localized subsidence of 1–2 mm yr⁻¹, partially offsetting the ongoing GIA uplift and complicating isolation of deglacial signals in geodetic data. In West Antarctica, persistent mass loss accelerates uplift beyond Holocene baselines, stabilizing marine ice margins against ocean-driven retreat. Holocene paleodata from raised marine sediments and emergence curves in East Antarctica, such as Prydz Bay, reveal non-linear rebound patterns, with initial rapid uplift rates of 5–10 mm yr⁻¹ in the early Holocene decelerating due to viscoelastic relaxation over regionally thinner lithosphere. These curves indicate peak emergence around 8–6 ka BP, followed by stabilization, consistent with a thinner mantle facilitating faster initial forebulge collapse and load redistribution. Such non-linearity underscores the influence of lateral Earth structure variations on Antarctic GIA, distinct from thicker-lithosphere responses elsewhere. ![Paulson (2007) model of lithospheric uplift rates due to post-glacial rebound, highlighting Antarctic differentials][float-right]

Theoretical Frameworks and Modeling

Formulation of the Sea-Level Equation

The sea-level equation (SLE) constitutes the foundational mathematical relation in glacial isostatic adjustment (GIA) theory, linking past surface mass redistributions—primarily from ice-sheet melting—to observable relative sea-level variations. It derives from conservation of mass, viscoelastic constitutive relations for Earth's mantle, and Poisson's equation for gravitational self-consistency, while incorporating centrifugal effects from rotational perturbations. The equation captures how deglaciation induces crustal uplift beneath former ice loads, peripheral subsidence from forebulge collapse, geoid perturbations (falling locally under load removal, rising distally), and feedback from redistributed ocean water altering the load itself. ![{\displaystyle S=N-U,}}(./assets/16087122c0305721d5080fb484a6581f23f89ef5.svg)[float-right] Fundamentally, relative sea level S (measured with respect to the deforming solid Earth surface) at colatitude θ, longitude λ, and time t equals the geoid height perturbation N minus the radial crustal displacement U (positive upward): S = N – U. Both N and U arise from the linear superposition of responses to surface loads, quantified via load Love numbers hl(s,t) for deformation and kl(s,t) for potential, where l denotes spherical harmonic degree, s is load epicentral distance, and time dependence reflects viscoelastic relaxation under Maxwell rheology (elastic lithosphere over viscous mantle). Direct eustasy from ice-volume equivalent water contributes positively to S, but GIA modifies it regionally: load removal elevates U (reducing local S), while geoid fall (N decrease) further lowers S near former ice centers, yielding net emergence rates up to 1 cm/year in Hudson Bay today. The complete SLE is a Fredholm integral equation of the second kind, self-consistent for ocean-covered sites where water depth equals S: ![{\displaystyle S(\theta ,\lambda ,t)={\frac {\rho {i}}{\gamma }}G{s}\otimes _{i}I+{\frac {\rho {w}}{\gamma }}G{s}\otimes _{o}S+S^{E}-{\frac {\rho {i}}{\gamma }}{\overline {G{s}\otimes _{i}I}}-{\frac {\rho {w}}{\gamma }}{\overline {G{o}\otimes _{o}S}},}(./assets/b75169c5733d7d8430215f0e02f8f14766be42cf.svg)[center] Here, convolution ⊗i integrates ice thickness history I(θ,λ,t') (t' ≤ t) over glaciated regions, while ⊗o does so over oceans; overbars denote global spherical averages ensuring mass conservation; Gs = Gs(h,k) is the sea-level Green's function; ρi ≈ 917 kg/m³ (ice density); ρw ≈ 1028 kg/m³ (seawater); γ incorporates mantle density contrasts; and SE(t) adds non-GIA eustasy (e.g., from minor melt sources). The right-hand feedback term requires iterative solution, as ocean loading S reinforces deformation. For land sites, S = 0 by definition, confining water to oceans. Empirical grounding of the SLE occurs through forward modeling of deglaciation scenarios, with ice histories and Earth parameters inverted against proxy records. Models such as ICE-6GC (VM5a) specify I(θ,λ,t) from 26 ka to present, yielding ~120–130 m eustatic rise calibrated to ~1000 relative sea-level curves from dated isolation basins, beach ridges, and peat (e.g., ~104 14C dates), plus GPS uplift (0.5–12 mm/year) and tide-gauge data. This fits imply lithospheric thickness ~90 km, upper-mantle viscosity ~0.4–1.0 × 1021 Pa·s, and lower-mantle ~1–2 × 1021 Pa·s, with residuals <1–2 m for Holocene S. Discrepancies arise from lateral viscosity variations or ice-model errors, but the framework robustly predicts "fingerprints" of ongoing adjustment.

Numerical Models, Viscosity Estimates, and Uncertainties

Numerical models for glacial isostatic adjustment (GIA), such as (Sea Level Equation solver), numerically solve the viscoelastic response of a layered to deglacial ice loads, incorporating gravitational self-consistency and topographic effects. Developed in since the mid-2000s, version 4.0, released in 2019, enables simulations of surface deformation, gravity anomalies, and relative sea-level changes on a spherical, non-rotating with rheology. These global forward models, often paired with ice histories like ICE-6G_C, predict present-day vertical uplift rates that typically align with GPS observations to within 20-30% in formerly glaciated regions, though discrepancies arise in peripheral zones due to unmodeled lateral heterogeneities. Lateral variations in lithospheric thickness and mantle are critical parameters, as uniform radial structures fail to capture observed strain patterns, necessitating 3D models derived from seismic tomography for improved fits. Mantle viscosity estimates from GIA inversions constrain upper mantle values around $10^{21} Pa·s and lower mantle values 10-30 times higher, inferred by minimizing misfits between modeled and observed uplift or geoid anomalies. These derive from Bayesian or least-squares optimizations against relative sea-level curves and GPS data, but trade-offs between viscosity profiles and ice load amplitude limit resolution, particularly for the transition zone where data sensitivity diminishes. Recent studies incorporating GRACE gravity data refine deep mantle estimates, suggesting lower-than-traditional viscosities in the top third of the lower mantle (around $10^{21.5} Pa·s) to match observed uplift rates near ice centers. Uncertainties in GIA predictions stem primarily from ice-sheet reconstruction errors, estimated at ±15-25% in deglacial volume due to sparse proxy data and chronological ambiguities in Holocene records, propagating to ±1-2 mm/yr errors in uplift rates. Additional variance arises from mantle anisotropy, which can alter rebound rates by up to 10-20% in shear-dominated flow, and from unresolvable radial discontinuities, with Monte Carlo ensembles revealing spatially patterned uncertainties peaking in data-sparse regions like Antarctica. Lateral viscosity heterogeneity, informed by tomography, introduces further ambiguity, as 1D models overestimate peripheral bulge collapse by ignoring plume-related low-viscosity zones. Debates persist over model overpredictions in select regions, such as early 2000s discrepancies in North American forebulge areas where predicted subsidence exceeded GPS by 1-2 mm/yr, attributed initially to overestimated ice volumes but later reconciled via refined Holocene sea-level fingerprints and viscoelastic corrections for ongoing glacier melt. Empirical tests against independent datasets, like GRACE-derived mass trends, highlight persistent tensions in deep viscosity inference, where J₂ degree-2 gravity signals demand higher lower-mantle values than local rebound data, underscoring the need for hybrid seismic-GIA constraints to resolve non-uniqueness. Overall, while benchmarks validate core model physics, quantification of parametric uncertainties via ensemble methods remains essential for applications like sea-level budgeting, where GIA errors contribute up to 20% to global mean estimates.

Geological and Tectonic Consequences

Induced Intraplate Earthquakes

Post-glacial rebound, or glacial isostatic adjustment (GIA), induces intraplate earthquakes primarily through perturbations to the crustal stress regime resulting from the removal of Pleistocene ice loads. The unloading reduces overburden pressure, leading to vertical uplift and lateral extension in formerly glaciated regions, which decreases horizontal compressive stresses and favors normal faulting on pre-existing weaknesses within stable continental interiors. This mechanism is distinct from interplate tectonics, as evidenced by focal mechanisms of intraplate events showing extensional regimes uncorrelated with distant plate boundary forces. Estimated stress changes from GIA range from 5 to 20 MPa in peak rebound zones, sufficient to bring critically stressed faults closer to failure when superimposed on ambient tectonic stresses of similar magnitude. In eastern Canada, encompassing the Ottawa Valley and surrounding Canadian Shield, GIA has been linked to elevated intraplate seismicity rates since deglaciation around 10,000–12,000 years ago. Moderate earthquakes (magnitudes typically below 5) cluster along reactivated ancient faults, such as those in the , where ongoing rebound at rates up to 1–2 mm/year correlates with normal fault mechanisms observed in events like the 2010 M4.0 Val-des-Bois earthquake. Seismicity here remains low overall but exceeds background levels in rebound-affected areas, with no causal ties to major plate boundaries; instead, GIA-driven extension modulates slip on inherited structures under a compressional far-field regime. Empirical data from moment tensor inversions confirm that post-glacial stress relaxation contributes to the observed stress orientations, though debates persist on the relative roles of GIA versus distant ridge-push forces. Similar patterns occur in Fennoscandia and northern Europe, where deglaciation-triggered seismicity manifests as clustered normal-fault events, with GIA stress drops promoting instability without invoking plate boundary dynamics. Quantified via viscoelastic models, these perturbations yield strain rates on the order of 10^{-17} to 10^{-16} s^{-1}, comparable to tectonic rates in intraplate settings and verifiable through paleoseismic records of fault scarps active post-10 ka BP. Overall, while GIA does not generate earthquakes de novo, it amplifies failure potential on faults oriented favorably to the induced extensional field, as substantiated by the absence of comparable seismicity in non-rebounding intraplate regions.

Links to Volcanism and Magma Mobilization

Post-glacial rebound induces mantle decompression through the removal of overlying ice loads, reducing lithostatic pressure and enabling partial melting where the solidus temperature falls below ambient mantle temperatures. Petrologic models demonstrate that this unloading effect promotes magma generation by facilitating decompression melting, particularly in regions like Iceland where rapid deglaciation post-Last Glacial Maximum (LGM) exposed the mantle to pressure drops of several megapascals. Quantitative simulations for Iceland estimate that ongoing glacial isostatic adjustment (GIA) uplift rates of 25-29 mm/yr correlate with enhanced melt production rates of 100-135%, yielding an additional 0.21-0.23 km³/yr of mantle-derived melt, with roughly half focused beneath central volcanic zones. Empirical evidence from Iceland highlights pulses of Holocene basaltic volcanism following LGM deglaciation around 12-10 , with eruption rates increasing up to 50-fold compared to glacial periods, as recorded in lava volumes and geochronologic from the Volcanic . In southern , deglaciation triggered enhanced tectonism and subglacial-to-subaerial eruptions, with dyke injections and low-intensity seismicity accompanying early unloading phases at ice margins, linking isostatic perturbations directly to magma mobilization. Similar patterns appear in glaciated , where edifice rates at volcanoes like Puyehue-Cordón Caulle accelerated post-deglaciation (e.g., 19 ), aligning with GIA-driven decompression rather than independent mantle fluxes. GRACE satellite gravimetry detects ongoing gravity lows in rebounding regions like Iceland, attributable to crustal uplift and mass redistribution from GIA, which spatially overlap with zones of post-glacial magmatic activity and support models of melt-induced density reductions. Although correlations between rebound and volcanism are robust in Iceland, causal interpretations face scrutiny for potential biases in eruption preservation and geochronologic uncertainties (±1-10 ka); alternative explanations, such as CO₂ degassing from deeper sources, are undermined by the precise timing of eruptive clusters—often within centuries of unloading—lacking alignment with non-glacial volatile release chronologies. Peer-reviewed assessments conclude that decompression facilitates magmatism in select cases but requires further multiproxy validation to distinguish from baseline tectonic drivers.

Evolution of Lithospheric Stress Fields

Post-glacial rebound induces changes in lithospheric stress fields through the viscoelastic relaxation of the mantle and flexure of the lithosphere following ice unloading, altering both vertical and horizontal stress components. The removal of glacial loads, which imposed vertical stresses exceeding 30 MPa at the Last Glacial Maximum, leads to a reduction in overburden pressure in formerly glaciated regions, while horizontal stresses decay more slowly due to lateral confinement and mantle flow. This differential relaxation increases deviatoric stresses regionally by 1-5 MPa, as modeled in eastern North America, where present-day GIA contributions can modify intraplate stresses by up to 5 MPa. Such perturbations integrate with plate-scale tectonic forces, potentially influencing the orientation and magnitude of principal stresses across cratonic interiors. In rebound centers like the Canadian Shield, the evolution of stress fields reflects a transition from ice-induced compression to post-deglacial extension in the upper crust, superimposed on pre-existing tectonic patterns. Borehole measurements, including hydraulic fracturing and borehole breakouts, reveal stress orientations in the Shield that deviate from ridge-push or slab-pull predictions, with deviations attributable to ongoing GIA effects over wavelengths exceeding 1000 km. Paleostress indicators, such as fault gouge fabrics and calcite twinning, corroborate these changes, showing radial compression near rebound maxima during early Holocene relaxation. These observations suggest that PGR modulates long-term lithospheric strength, contributing to debates on craton stability by either reinforcing compressional regimes that inhibit rifting or amplifying differential stresses that could precondition intracratonic weaknesses. The integration of vertical uplift and into plate-scale es varies with and lithospheric thickness; higher viscosities delay , sustaining perturbations for . Numerical models indicate that deviatoric maxima of several persist in the beneath rebound peripheries, influencing far-field . through inversion of focal mechanisms in low-seismicity cratons like the confirms GIA's in evolving stress regimes, distinct from purely tectonic drivers.

Implications for Contemporary Sea Level and Climate Attribution

Separating Glacial Isostatic Adjustment from Other Signals

To isolate glacial isostatic adjustment (GIA) from contemporary sea level signals, researchers pair satellite altimetry data, which measure absolute sea surface height relative to the geoid, with modeled GIA corrections to account for ongoing viscoelastic deformation of the solid Earth and associated geoid changes. For instance, the TOPEX/Poseidon mission, operational since August 1992, provides along-track altimetry records that, after applying GIA models such as ICE-6G(VM5a), yield global mean sea level trends adjusted upward by approximately 0.3 mm/year to reflect true ocean volume changes excluding post-glacial effects. Tide gauge records, which capture relative sea level changes, are corrected using co-located Global Positioning System (GPS) measurements of vertical land motion, where uplift rates attributable to GIA—often 1–10 mm/year in formerly glaciated regions like Hudson Bay—are subtracted to derive absolute sea level equivalents. Gravity missions such as GRACE (2002–2017) and GRACE-FO (2018–present) further enable separation by isolating mass redistribution signals from GIA-induced gravity anomalies, distinguishing ice-mass loss from rebound effects through forward modeling and inversion techniques. GIA explains a substantial fraction of regional sea level variance, with modeled contributions reaching 2–5 mm/year in high-latitude areas, countering assumptions of spatially uniform rise by highlighting natural, load-history-dependent patterns that amplify or offset other drivers like steric expansion. In regions such as the U.S. Northeast, GIA subsidence signals of 0.5–1 mm/year contribute up to 30% of observed relative trends, necessitating explicit partitioning to avoid conflating rebound dynamics with anthropogenic mass additions. Empirical analyses from 2020–2025 tide gauge networks, post-GIA and steric corrections, indicate relative sea level stability or rates below 2 mm/year in select rebounding locales like parts of Scandinavia and Canada, consistent with dominant natural adjustment over residual ocean dynamic signals. These corrections underscore GIA's role in reconciling apparent discrepancies between local records and global averages, prioritizing geophysical realism over aggregated attributions.

Regional Discrepancies in Observed Sea Level Rise

Observed sea level changes exhibit significant regional variations, with post-glacial rebound, or glacial isostatic adjustment (GIA), contributing to non-uniform patterns superimposed on global eustatic rise. In formerly glaciated regions such as Scandinavia, ongoing land uplift exceeds contemporary eustatic sea level rise, resulting in relative sea level fall. For instance, uplift rates in the Gulf of Bothnia reach up to 9-10 mm/yr, leading to net relative sea level decreases of approximately 1 mm/yr or more when accounting for global mean rise around 3 mm/yr. Conversely, in peripheral regions distant from former ice sheets, subsidence of the collapsed peripheral bulge amplifies relative sea level rise. Along the U.S. East Coast, this subsidence contributes 1-2 mm/yr to observed rates, with mid-Atlantic areas experiencing subsidence of at least 2 mm/yr due to forebulge dynamics. Satellite gravimetry and altimetry data since 2020 have confirmed these GIA-induced fingerprints, revealing spatial patterns in mass redistribution that align with Holocene sea level proxies. GRACE Follow-On observations and models derived from them show uplift in northern high latitudes and subsidence in mid-latitudes, matching empirical reconstructions of past ice loading. These patterns explain why relative sea level rise accelerates in subsiding zones like the U.S. East Coast while decelerating or reversing in rebound areas, independent of recent climatic forcing. Empirically, post-Last Glacial Maximum () sea level approximately 120-125 globally, with deglacial phases featuring pulses at rates up to 4 per century, demonstrating that heterogeneous and pulsed rises are inherent to isostatic responses rather than anomalies. Such historical variability underscores how GIA modulates contemporary regional discrepancies, with uplift zones counteracting eustatic trends and peripheral enhancing them.

Debates on Acceleration: Empirical Data vs. Model Projections

A 2025 peer-reviewed analysis of relative sea-level records from 243 tide gauges worldwide, spanning periods often exceeding 100 years, concluded that approximately 95% of suitable locations exhibited no statistically significant acceleration in the rate of rise. This empirical assessment, focusing on homogeneity-corrected data to account for datum shifts and instrumental changes, found average rates of 1.2 to 1.8 mm/year across sites, consistent with linear trends rather than quadratic acceleration post-1900. Critics of this view, including analyses of combined tide gauge and satellite altimetry data, assert widespread acceleration, citing rates increasing from about 1.4 mm/year in the early 20th century to 3.7 mm/year since 1993, attributed partly to enhanced thermal expansion and ice-mass loss. However, such claims often rely on shorter satellite records prone to calibration drifts and regional biases, whereas long-term tide gauge ensembles, after rigorous homogeneity testing, support subdued, non-accelerating trends when vertical land motion from glacial isostatic adjustment (GIA) is considered. Projections from IPCC assessments, such as those in the 2021 Sixth Assessment Report, forecast accelerating mean sea-level rise driven by dynamic ice-sheet instabilities, projecting 0.28 to 0.55 meters by 2100 under various emissions scenarios beyond contributions from and GIA. These models incorporate GIA corrections derived from viscoelastic Earth models, yet uncertainties in viscosity—estimated to vary by factors of 2–10 in upper and lower layers—can propagate errors of up to 1 mm/year in regional sea-level budgets, potentially inflating attributed acceleration. Independent evaluations highlight that unverified assumptions in Antarctic ice-shelf collapse and Greenland surface-melt parameterization exceed empirical constraints from mass-balance observations, leading to overprojections when benchmarked against tide gauge-derived eustatic rates. Alternative interpretations emphasize natural multidecadal cycles, such as Atlantic Multidecadal Variability, which modulate observed rates without requiring additional anthropogenic forcing beyond established (approximately 1.5 mm/year globally since 1900). In rebound-dominated regions like and , unadjusted tide gauge records reveal relative sea-level falls of 5–10 /year due to ongoing uplift, masking any eustatic signal and underscoring the necessity of precise modeling for attribution. Globally, after deconvolution, empirical reconstructions from pre-satellite era data indicate no robust post-1900 acceleration exceeding linear trends compatible with steric and isostatic components, challenging model-driven narratives of climate-induced speedup. This discrepancy persists despite advances in parameter estimation, where ensemble modeling reveals prediction spreads of 0.5–2 /year in corrected sea-level trends owing to incomplete deglaciation histories and lateral viscosity heterogeneities. Skeptical analyses, prioritizing tide gauge homogeneity over model ensembles, thus favor interpretations where observed rates align with 20th-century baselines, attributing apparent discrepancies to overreliance on uncertain ice-melt projections rather than empirical validation.

Applications and Broader Impacts

Corrections in Geodetic Datums and Vertical Reference Systems

In regions undergoing post-glacial rebound, geodetic datums such as NAD83 and vertical reference systems like NAVD88 experience systematic vertical drifts due to ongoing glacial isostatic adjustment (GIA), requiring explicit corrections to maintain epoch-independent accuracy. Static realizations of these datums, fixed at reference epochs like 1997.0 for NAD83(CSRS), accumulate errors from unmodeled crustal motion, with GIA contributing dominant vertical velocities in glaciated interiors. Dynamic reference frames incorporate GIA-derived velocity grids to propagate positions forward, reconciling GPS observations with modeled rebound. Canada's NAD83(CSRS) exemplifies these adjustments through the NAD83v70VG crustal velocity model, which fuses GPS with GIA predictions to vertical rates of 0-12 /year across the Canadian , peaking near . These , typically 1-10 per in peripheral zones but higher centrally, enable time-dependent transformations for . In , ETRF realizations aligned to ITRF analogous GIA offsets for Fennoscandian uplift rates of 5-10 /year, refined via GPS fields to mitigate datum inconsistencies. Repeated leveling surveys provide empirical of datum drift, documenting vertical discrepancies of several centimeters over decades in rebound areas, as seen in North American networks tied to tide gauges. Such surveys, combined with gravity measurements, benchmark effects against static orthometric heights, proving indispensable for precise altimetry in and sea-level . Failure to apply these propagates errors into infrastructure projects, including differential straining of pipelines in , where geodetic networks have tracked motions since the 1990s to quantify risks from spatially varying uplift.

Exploration for Hydrocarbons and Groundwater

Post-glacial rebound generates differential stresses within the lithosphere that can reactivate faults and induce fracturing, thereby enhancing permeability in hydrocarbon reservoirs. In the northern North Sea, this process contributes to fault reactivation alongside elevated pore pressures and optimal fault orientations, influencing hydrocarbon column heights and leakage potential. Such uplift-related fracturing facilitates secondary migration pathways and improves reservoir connectivity, as documented in late Cenozoic uplift scenarios applicable to post-glacial contexts. Exploration strategies incorporate glacial isostatic adjustment (GIA) models to compensate for ongoing vertical crustal motions in seismic imaging, ensuring accurate subsurface mapping in rebounding basins. For instance, GIA simulations integrated with geophysical data predict differential uplift rates, aiding the delineation of reservoir traps in regions like the Grand Banks where post-glacial deformation alters structural integrity. This correction is critical for North Sea fields, where unaccounted rebound can distort seismic reflections and lead to misestimation of reservoir volumes. In groundwater contexts, post-depletion aquifer rebound in glaciated regions superimposes elastic recovery signals on isostatic uplift, complicating hydrogeological assessments. Techniques such as InSAR, combined with GNSS and GRACE data, enable joint inversion to separate anthropogenic groundwater storage changes from GIA, as demonstrated in constrained Bayesian models that resolve vertical displacements at millimeter precision. This differentiation prevents conflation of extraction-induced subsidence with rebound, informing sustainable aquifer management in areas like Canada's Hudson Bay lowlands where rates exceed 10 mm/year.

Legal and Territorial Disputes in Rebounding Regions

In regions undergoing post-glacial rebound, such as eastern James Bay in subarctic Canada, ongoing crustal uplift complicates indigenous land claims over coastal islands and shorelines. Rates of uplift in the James Bay area measure approximately 1 to 1.5 meters per century, outpacing contemporary global sea-level rise of about 0.3 meters per century, resulting in net emergence of land and alteration of water boundaries. This dynamic has bearing on unresolved negotiations stemming from the 1975 James Bay and Northern Quebec Agreement (JBNQA), where Canada committed to discussions with Quebec Cree and Nunavik Inuit regarding ownership of islands in eastern James Bay. Glacial isostatic adjustment (GIA) models predict that continued rebound, combined with potential future sea-level acceleration, could shift the submerged or emergent status of certain low-lying islands, thereby affecting resource rights and territorial delineations under the JBNQA framework. These changes intersect with indigenous treaty obligations, as rebound-induced shoreline progradation expands habitable or usable land in areas traditionally claimed by Cree First Nations along Hudson and James Bays. Historical treaties, including those incorporated into the JBNQA, often reference fixed water bodies or coastal features for defining hunting, fishing, and trapping territories, but GIA-driven alterations challenge such static interpretations by progressively exposing former marine areas. For instance, in the "Cree Zone" and "Joint Inuit/Cree Zone" offshore eastern James Bay, uplift forecasts suggest long-term stabilization of island configurations, yet short-term variability heightens contention between Cree communities and Inuit groups over emerging land suitable for settlement or resource extraction. Canadian courts have not yet ruled definitively on GIA's role in treaty reinterpretation, but scholarly analyses emphasize the need for dynamic boundary adjustments to reflect geophysical realities rather than assuming permanence. Under the United Nations Convention on the Law of the Sea (UNCLOS), GIA influences continental shelf delineations in rebounding Arctic regions, where states like Canada, Russia, and Denmark submit data for extended shelf claims beyond 200 nautical miles. Article 76 requires foot-of-slope or sediment-thickness criteria for shelf extension, and GIA corrections to bathymetric and geodetic data are incorporated in submissions to account for historical crustal motion affecting baseline positions. Canada's Arctic submissions in the 2010s, for example, integrated GIA modeling to validate shelf morphology, preventing overestimation of sediment deposition due to isostatic effects. While no arbitrations have hinged solely on GIA disputes as of 2023, the process has informed delimitations, such as Russia's 2015 revised claim on the Lomonosov Ridge, where rebound data refines geoid undulations and relative sea-level histories. In non-Arctic rebound zones like the Baltic Sea, where uplift rates reach 8-9 mm per year, UNCLOS baseline ambulatory nature—tied to the low-water line—implies gradual EEZ expansion for states like Sweden and Finland, though bilateral agreements mitigate potential friction over shared gulfs.

References

  1. [1]
    What is glacial isostatic adjustment? - NOAA's National Ocean Service
    Jun 16, 2024 · Glacial isostatic adjustment, the ongoing movement of land once burdened by ice-age glacier, is also called the Mattress Effect.
  2. [2]
    Glacial Isostatic Adjustment | U.S. Geological Survey - USGS.gov
    Sep 22, 2022 · When the ice sheets retreat, the forebulges slowly collapse while the indentation rebounds. This is called post-glacial isostatic adjustment.
  3. [3]
    [PDF] Atmospherically forced sea-level variability in western Hudson Bay ...
    Oct 6, 2021 · Combining satellite altimeter data with the Churchill tide- gauge data gives an uplift rate of about 9.0 ± 0.8 mm yr−1 (Ray, 2015).Missing: Scandinavia | Show results with:Scandinavia
  4. [4]
    [PDF] Observation of glacial isostatic adjustment in ''stable'' North America ...
    Jan 26, 2007 · [2] Postglacial rebound or glacial isostatic adjustment. (GIA) is the response of the solid Earth to the changing surface load brought about by ...
  5. [5]
    The fingerprints of sea level rise
    Sep 29, 2015 · The East Coast is also on the losing end of another important solid-Earth process that affects regional sea levels: post-glacial rebound. After ...
  6. [6]
    [PDF] Glacial Isostatic Adjustment: Global (post-LGM) & SE Alaska (post-LIA)
    Glacial Isostatic Adjustment. (aka “post-glacial rebound”). • Adding weight of ice depresses crust, forces viscous mantle to flow out of the way.
  7. [7]
    Glacial isostatic adjustment reduces past and future Arctic subsea ...
    Apr 15, 2024 · Including GIA significantly reduces present-day subsea permafrost thickness, chiefly because of hydro-isostatic effects as well as deformation related to ...
  8. [8]
    Isostasy and Flexure of the Lithosphere - Physics Today
    Oct 1, 2002 · In its simplest form, lithospheric isostasy is a restatement of Archimedes' principle: The upper parts of the Earth float on its interior. Such ...<|separator|>
  9. [9]
    Isostasy - an overview | ScienceDirect Topics
    Isostasy is defined as the condition of gravitational equilibrium between the Earth's crust and the underlying viscous mantle, explaining the properties ...
  10. [10]
    Modeling Viscoelastic Solid Earth Deformation Due To Ice Age and ...
    Mar 1, 2023 · The solid Earth is deforming in response to past and present ice loading changes at rates determined by elastic and viscous parameters We ...
  11. [11]
    Surface loading of a viscoelastic earth—I. General theory - Free
    These studies of glacial iso-static adjustment have, in turn, found application in diverse disciplines, including mantle rheology, ice sheet dynamics and ...
  12. [12]
    Regional viscosity variations in Earth's mantle - EGU Blogs
    Feb 28, 2024 · The mantle behaves viscoelastically, where elastic deformation is reversible and occurs over very short timescales, while viscous deformation is ...
  13. [13]
    Modelling post-glacial rebound with lateral viscosity variations
    Important constraints on mantle viscosity can be made by comparing observations of post-glacial rebound (PGR) with model predictions of the Earth's response to ...
  14. [14]
    Inferences of mantle viscosity based on ice age data sets: Radial ...
    Sep 19, 2016 · Our forward prediction of the postglacial decay time for a specific site r is based on finding the best fit function of the form 3 through the ...
  15. [15]
    The impact of water loading on postglacial decay times in Hudson Bay
    May 1, 2018 · Postglacial decay time estimates thus provide a powerful datum for constraining the Earth's viscous structure and improving GIA predictions. We ...<|control11|><|separator|>
  16. [16]
    A glacial isostatic adjustment model for the central and northern ...
    In the final model, the last glacial maximum (LGM) thickness of the Laurentide Ice Sheet west of Hudson Bay was ∼3.4–3.6 km. Conversely, east of Hudson Bay, ...
  17. [17]
    Glacial isostatic adjustment and the radial viscosity profile from ...
    Nov 8, 2002 · Viscosities in the upper mantle are increasing from around 4 × 1020 to 2 × 1021 Pa s. The large jump in viscosity between the upper and lower ...
  18. [18]
    The persisting conundrum of mantle viscosity inferred from mantle ...
    Dec 15, 2024 · In this study, we investigate the effects of mantle viscosity structures on observations of the geoid, mantle structures, and present-day crustal motions and ...
  19. [19]
    Effects of the Last Quaternary Glacial Forebulge on Vertical Land ...
    Jul 1, 2025 · Consequently, the location and height of the forebulge may change after hundreds of thousands of years of glacial loading or unloading.
  20. [20]
    GIA & Trends | Data Portal – GRACE Tellus - NASA
    However, in regions located in the periphery of the last ice age's ice masses, the solid earth experiences what is known as 'forebulge collapse'. For example, ...
  21. [21]
    Glacial isostatic adjustment, relative sea level history and mantle ...
    Abstract. Models of the glacial isostatic adjustment process, which is dominated by the influence of the Late Pleistocene cycle of glaciation and deglaciat.Missing: peer- | Show results with:peer-
  22. [22]
    The Viscosity of the Top Third of the Lower Mantle Estimated Using ...
    Mar 14, 2021 · Laurentia is today rising at just ≈12 mm/yr, showing that solid Earth has to a large extent viscously relaxed in response to ice unloading ...
  23. [23]
    The impact of regional-scale upper-mantle heterogeneity on glacial ...
    Jul 3, 2025 · Here, we investigate the effects of incorporating smaller-scale lateral variability in upper-mantle viscosity into 3-D GIA simulations.
  24. [24]
    Sensitivity of glacial isostatic adjustment observations on 3D Earths ...
    Sep 20, 2025 · One aim of this paper is to review and clarify the physics of the Perturbation method and bring out some important aspects of this method that ...
  25. [25]
    [PDF] Past shore-level and sea-level displacements
    More than. 1 000 highest coastline sites, between 20 and 285 m above present sea-level, are reported within Sweden. There are practically no tides along the ...
  26. [26]
    An introduction to the ups and downs of eustasy - GeoScienceWorld
    The first published maps of post-glacial uplift or rebound by Swedish geologist,. Gerard DeGeer (1888). The contours (called isobases) show uplift in meters, ...<|separator|>
  27. [27]
    A concise history of postglacial land uplift research (from its ...
    Aug 7, 2025 · Thanks to a paper of De Geer in 1888, the glacial isostatic origin was established. Fennoscandia became the classic area of glacial isostasy ...
  28. [28]
    Quaternary sea level change in Scotland | Earth and Environmental ...
    Jan 23, 2018 · This paper summarises developments in understanding sea level change during the Quaternary in Scotland since the publication of the Quaternary of Scotland ...
  29. [29]
    The altitude and age of the main postglacial shoreline in Eastern ...
    This paper discusses the patterns of altitude and age variation of the Main Postglacial Shoreline in eastern Scotland from the Firth of Forth to the Dornoch ...
  30. [30]
    [PDF] R. Bell-Rising of the Land around Hudson Bay. 219
    ... raised beaches on the long straight shore out on the open sea. Hud- son Bay is about 1000 miles long and its outline is funnel- shaped, with James Bay ...
  31. [31]
    Wave-Like Progress of an Epeirogenic Uplift
    departure of the ice, the uplifting wave of the earth's crust has raised the basin of Hudson bay 300 to 500 feet since the sea ... ROBERT BELL, Geol. and ...<|separator|>
  32. [32]
    Haskell [1935] revisited - Mitrovica - 1996 - AGU Journals - Wiley
    Jan 10, 1996 · The Haskell [1935] value of 1021 Pa s for mantle viscosity is a classic and enduring constraint on the rheology of the Earth's interior.
  33. [33]
    A computer simulation of post-glacial rebound with partial water load
    Vening-Meinesz, 1937. F.A. Vening-Meinesz. The determination of the Earth's plasticity from the postglacial uplift of Scandinavia: isostatic adjustment. Verh ...
  34. [34]
    Postglacial sea level: energy method - ScienceDirect.com
    Our method is forward modeling on a laterally homogeneous Maxwell viscoelastic earth using a number of admissible `realistic' viscosity models. We demonstrate ...
  35. [35]
    Sea-level change, glacial rebound and mantle viscosity ... - NASA ADS
    Gutenberg, B., 1941. Changes in sea level, postglacial uplift, and mobility of the Earth's interior, Bull. geol. Soc. Am., 52, 721—772. Hafsten, U., 1956 ...
  36. [36]
    Glacial isostatic adjustment modelling: historical perspectives ...
    Glacial isostatic adjustment (GIA) describes the response of the solid Earth, the gravitational field, and the oceans to the growth and decay of the global ice ...
  37. [37]
    Glacial-Isostatic Adjustment—I. The Forward Problem
    Summary. The isostatic adjustment of a radially stratified visco-elastic spheroid is treated using space-time Green functions for the associated surface ma.
  38. [38]
    Continuous GPS measurements of postglacial adjustment in ...
    Aug 10, 2002 · We find that the maximum observed uplift rate (∼10 mm yr−1) and the maximum predicted uplift rate agree to better than 1 mm yr−1. The ...
  39. [39]
    Vertical crustal motion observed in the BIFROST project
    The cut-off test at the end shows that there is still some drift of the rate during the last twelve-month period, and the present estimate of 10±0.2 mm/year ...
  40. [40]
    [PDF] University of Groningen Rates of Holocene isostatic uplift and ...
    The present mean overall apparent uplift rate is of the order of 4–. 5 mm/yr, but immediately after deglaciation the rate of crustal rebound was several times ...
  41. [41]
    Holocene relative sea-level changes and glacial isostatic ...
    Aug 1, 2011 · Rates of RSL change were highest during the early Holocene and have been decreasing over time, due to the exponential form of the GIA process ...Missing: shorelines | Show results with:shorelines
  42. [42]
    A Robust Estimation of the 3‐D Intraplate Deformation of the North ...
    May 5, 2018 · Glacial isostatic adjustment (GIA) is the main cause of deformation in intraplate North America. Here we use up to 3,271 Global Positioning ...
  43. [43]
    Reversal of the Direction of Horizontal Velocities Induced by GIA as ...
    Aug 27, 2018 · In regions undergoing glacial isostatic adjustment present-day horizontal surface motion is observed to point mostly, but not always, away from ...
  44. [44]
    Glacial Isostatic Adjustment, Intraplate Strain, and Relative Sea ...
    May 28, 2019 · Determining the motion of the North American plate is difficult due to horizontal intraplate deformation, caused mainly by glacial isostatic ...<|separator|>
  45. [45]
    Glacial isostatic adjustment on a rotating earth - Oxford Academic
    The peak radial and horizontal velocities are of the order of 0.5 mm yr−1 ... Indeed, the maximum peak-to-peak radial and geoid displacement rates are ∼1 mm yr−1.
  46. [46]
    Sea level fingerprints and regional sea level change - ScienceDirect
    Aug 1, 2021 · Sea level fingerprints explain regional total sea level changes. Effect of sea level fingerprints is significant in sea level trends near coastlines.
  47. [47]
    [PDF] The Moving boundaries of sea level change
    Numerical prediction of the present-day impact of glacial isostatic adjustment on (a) relative sea level measured by tide gauges, (b) change in sea surface ...
  48. [48]
    On Some Properties of the Glacial Isostatic Adjustment Fingerprints
    Since the Glacial Isostatic Adjustment (GIA) caused by the melting of past ice sheets is still contributing to present-day regional and global sea-level ...
  49. [49]
    The Global Fingerprint of Modern Ice‐Mass Loss on 3‐D Crustal ...
    Aug 16, 2021 · We demonstrate that mass changes in the Greenland Ice Sheet and high latitude glacier systems each generated average crustal motion of 0.1–0.4 ...1 Introduction · 2.1 Greenland · 2.2 Glaciers And Ice Caps
  50. [50]
    The contribution of glacial isostatic adjustment to relative sea-level ...
    Thus, peripheral bulge subsidence is a significant contributor to RSL rise on the Atlantic and Pacific coasts of North America. On the eastern side of the North ...
  51. [51]
    Uncertainty Estimation in Regional Models of Long‐Term GIA Uplift ...
    Jun 30, 2020 · Likewise, away from the former ice sheet centers, GIA uncertainty for relative sea level change is inferred to be ~0.3–0.5 mm/yr along the U.S. ...1 Introduction · 3.2 North America · 5.1 U.S. East Coast
  52. [52]
    The contribution of glacial isostatic adjustment to projections of sea ...
    Jul 30, 2016 · We determine the contribution of glacial isostatic adjustment (GIA) to future relative sea-level change for the North American coastline ...
  53. [53]
    TELLUS GRACE Level-3 1.0-degree Glacial Isostatic Adjustment v1 ...
    To isolate signals of contemporary surface mass loss in the cumulative satellite gravimetry measurements, contemporary GIA rates are computed and subtracted ...
  54. [54]
    Sea-level fingerprints emergent from GRACE mission data
    May 9, 2019 · Here we perform a systematic calculation of sea-level fingerprints of on-land water mass changes using monthly Release-06 GRACE Level-2 Stokes coefficients.
  55. [55]
    Gravity measurements in southeastern Alaska reveal negative ...
    Dec 3, 2010 · The gravity rate of change was found to be −3.5 to −5.6 μGal/yr in the gravity network. Furthermore, gravity results obtained during the 3 years ...
  56. [56]
    Past and present‐day ice mass variation on Svalbard revealed by ...
    Nov 17, 2011 · The gravity rate varies through 10 years of observation; −0.23 μGal/yr in 2000–2002, −3.22 μGal/yr in 2002–2005 and −1.10 μGal/yr in 2005–2010.
  57. [57]
    Observation of glacial isostatic adjustment in “stable” North America ...
    Jan 26, 2007 · The vertical velocities show fast rebound (∼10 mm/yr) near Hudson Bay, the site of thickest ice at the last glacial maximum, which changes to ...
  58. [58]
    Climate-driven polar motion: 2003–2015 - PMC
    Apr 8, 2016 · The 20th century linear drift is generally explained by the glacial isostatic adjustment (GIA) ... Chandler wobble frequency) that we are ...
  59. [59]
    Climate‐Induced Polar Motion: 1900–2100 - AGU Journals - Wiley
    Mar 5, 2025 · ... Earth's rotation axis and induces polar motion. Here ... Glacial Isostatic Adjustment (GIA)—could also affect polar motion significantly.
  60. [60]
    Energy balance of glacial isostatic adjustment - Oxford Academic
    Understanding the feedback between the glacial isostatic adjustment (GIA) and the Earth's rotation ... Precise computation of GIA-induced polar motion (PM) ...
  61. [61]
    A land uplift model in Fennoscandia combining GRACE and ...
    They derived a peak uplift of about 9.5 mm/year placed in the middle of the northern Gulf of Bothnia.
  62. [62]
    [PDF] Computation of Historical Shore Levels in Fennoscandia due to ...
    A consistent and accurate set of postglacial uplift rates along the coasts of Fennoscandia has been published by Ekman (1996). These uplift rates are primarily ...
  63. [63]
    The Postglacial Rebound Signal of Fennoscandia Observed by ...
    Nov 29, 2010 · We have compiled GPS uplift rates for a dozen sites in Fennoscandia for which time series of absolute gravimetry also exists, spanning about a ...
  64. [64]
    Late Weichselian and Holocene shoreline displacement in the ...
    shoreline to the Ångermanland and Baltic area shows a relative uplift at 11,800 years B.P. of 400–450 m in the central area of glaciation. The island of Hitra ...Missing: total | Show results with:total
  65. [65]
    Fennoscandian Land Uplift: Past, Present and Future - ResearchGate
    Aug 10, 2025 · In the 19th century the cause of the land uplift and its environmental consequences were unknown and the rise of arable land from the sea ...
  66. [66]
    The origin of the Baltic Sea
    About 10,700 years ago, the connection to the Skagerrak almost completely closed and the Ancylus Lake was formed, a freshwater lake that at times covered an ...
  67. [67]
    (PDF) The late Quaternary development of the Baltic Sea basin
    13). Freshwater conditions were established in the Baltic basin (the Ancylus Lake stage) due to the glacial isostatic uplift of the central Sweden, the ...
  68. [68]
    Postglacial uplift, neotectonics and seismicity in Fennoscandia
    Oct 1, 2000 · In historic time, the rate of uplift along the coasts of Fennoscandia has been so high that its effects were easily observed within one ...Missing: raised | Show results with:raised
  69. [69]
    Development of the coastal systems of the easternmost Gulf of ...
    Jan 1, 2016 · We examine three questions concerning the post-glacial geological history of the eastern Gulf of Finland: (1) the amplitude of the Holocene ...<|control11|><|separator|>
  70. [70]
    [PDF] Glacial Rebound and Crustal Stress in Finland
    Nov 10, 2003 · Erosion and sedimentation redistribute surface loads and stress the crust. Plate tectonic forces add to the regional background stress.
  71. [71]
    A glacial isostatic adjustment model for the central and northern ...
    An improved model of the central and northern Laurentide Ice Sheet is constrained by 24 relative sea level histories and 18 present-day GPS-measured vertical ...
  72. [72]
    Rise of Great Lakes Surface Water, Sinking of the Upper Midwest of ...
    Jul 31, 2020 · We find this sinking to be produced primarily by viscous collapse of the former Laurentide ice sheet forebulge and secondarily by elastic Great Lakes loading.Missing: peripheral | Show results with:peripheral
  73. [73]
    Dynamic Heights for the International Great Lakes Datum of 2020
    Due to changes from Global Isostatic Adjustment and improvements in positioning technology, IGLD 85 must be updated. The International Great Lakes Datum of 2020 ...Missing: GPS arrays rebound
  74. [74]
    Relative sea‐level history and isostatic rebound in Newfoundland ...
    Mapping the time at which sea-level fell below the present level allows estimation of the rate of northerly migration and collapse of the forebulge at between ...Missing: rise | Show results with:rise
  75. [75]
    Relative sea‐level history and isostatic rebound in Newfoundland ...
    Apr 10, 2025 · Following the last glacial maximum, relative sea level in eastern Newfoundland initially fell and subsequently rose to modern level (Liverman, ...
  76. [76]
    Vertical Displacements and Sea‐Level Changes in Eastern North ...
    Oct 9, 2024 · We find a general pattern of subsidence (causing accelerated relative sea-level rise) in the eastern United States region and uplift (causing ...Missing: bulge | Show results with:bulge
  77. [77]
    Contrasting Response of West and East Antarctic Ice Sheets to ...
    Jun 18, 2021 · This process, called glacial isostatic adjustment, is important to consider in ice sheet models because it can stabilize an ice sheet undergoing ...
  78. [78]
    An assessment of forward and inverse GIA solutions for Antarctica
    Over the ASE, the predicted mean uplift rates vary from 4.9 ± 2.0 and 5.5 ± 1.1 mm yr−1 in models such as G14 or RATES to 0.8 mm yr−1 as estimated by AGE‐1b ( ...
  79. [79]
    new glacial isostatic adjustment model for Antarctica: calibrated and ...
    The maximum range of uplift rates produced by the 17 earth models that lie within the 95 per cent confidence limit of the RSL data analysis is ∼2.7 mm yr−1 (Fig ...
  80. [80]
    The impact of 3-D Earth structure on far-field sea level following ...
    Dec 1, 2021 · 3-D models yield greater subsidence of the peripheral bulge. •. 1-D models overestimate sea level by up to ∼0.3 m in the far field.
  81. [81]
    Unprecedented mass gain over the Antarctic ice sheet between ...
    The AIS showed a record-breaking mass gain of 129.7 ± 69.6 Gt yr -1 between 2021 and 2022. During this period, the mass gain over the East AIS and Antarctic ...
  82. [82]
    Rapid early Holocene sea-level rise in Prydz Bay, East Antarctica
    Following deglaciation, as crustal rebound outpaced sea-level rise, they became isolated and transformed into freshwater, or brackish, lakes. The sediments in ...Missing: emergence | Show results with:emergence
  83. [83]
    A Multiple 1D Earth Approach (M1DEA) to account for lateral ...
    In this paper, we explore the impact of variations in Earth structure on local uplift rates for loads at various distances, with a focus on Antarctic ...<|separator|>
  84. [84]
    On Postglacial Sea Level | Geophysical Journal International
    Summary. An exact method is presented for calculating the changes in sea level that occur when ice and water masses are rearranged on the surface of elasti.
  85. [85]
    On post-glacial sea level: I. General theory - Oxford Academic
    Modern analyses of sea level changes due to glacial isostatic adjustment (GIA) are based on the classic sea level equation derived by Farrell & Clark (1976).
  86. [86]
    Postglacial variations in the level of the sea: Implications for climate ...
    Nov 1, 1998 · During each glacial stage, global sea level fell by ∼120 m on average, as extensive ice sheets formed and thickened on the surfaces of the continents.Missing: equation | Show results with:equation
  87. [87]
    SELEN4 (SELEN version 4.0): a Fortran program for solving ... - GMD
    Dec 4, 2019 · We present SELEN 4 (SealEveL EquatioN solver), an open-source program written in Fortran 90 that simulates the glacial isostatic adjustment (GIA) process.
  88. [88]
    The Antarctica component of postglacial rebound model ICE-6G_C ...
    The Southern Antarctic Peninsula is inferred to be rising at 2 mm yr–1, requiring there to be less Holocene ice loss there than in the prior model ICE-5G (VM2).
  89. [89]
    Glacial‐Isostatic Adjustment Models Using Geodynamically ...
    Oct 25, 2021 · We construct 18 3D Earth structures that are derived from seismic tomography models and are geodynamically constrained.
  90. [90]
    Limitations on the inversion for mantle viscosity from postglacial ...
    SUMMARY. Observations of postglacial rebound (PGR) can provide important constraints on mantle viscosity structure. In this study, we investigate how well.
  91. [91]
    Uncertainties of Glacial Isostatic Adjustment Model Predictions in ...
    Apr 17, 2020 · We quantify GIA prediction uncertainties of 250 1D and 3D glacial isostatic adjustment (GIA) models through comparisons with deglacial relative sea-level (RSL) ...
  92. [92]
    An analysis of anisotropic mantle viscosity, and its possible effects ...
    This paper presents a preliminary analysis of the possibility of using post-glacial rebound observations to place constraints on the anisotropy of upper-mantle ...Missing: estimates | Show results with:estimates
  93. [93]
    The influence of lateral Earth structure on inferences of global ice ...
    Aug 15, 2022 · We calculate gravitationally self-consistent sea-level change for each ice history and Earth model pairing (ICE-6G with VM2 and the VM2-based 3- ...
  94. [94]
    A global semi-empirical glacial isostatic adjustment (GIA) model ...
    Feb 14, 2020 · Nonetheless, GIA model uncertainties are still one of the main source of errors for, e.g., GRACE-based estimates of global mean ocean mass ...
  95. [95]
    Moment tensors, state of stress and their relation to post-glacial ...
    In the stable craton of northeastern Canada moderate seismicity as well as glacial isostatic adjustment (GIA) has been observed. We investigate five earthquakes ...Missing: glacio- | Show results with:glacio-
  96. [96]
    Postglacial rebound and fault instability in Fennoscandia
    There, the magnitude is predicted to be about 5 MPa at around 8 ka BP, about five times larger than that predicted for Laurentia. If these values of fault ...
  97. [97]
    Dependence of horizontal stress magnitude on load dimension in ...
    He concluded that for eastern Canada, the background stress is as important as the rebound stress in determining if postglacial rebound can trigger earthquakes.
  98. [98]
    Reactivation Potential of Intraplate Faults in the Western Quebec ...
    Jul 16, 2021 · This seismicity likely results mainly from the reactivation under the present-day tectonic stress field of inherited structures such as late ...
  99. [99]
    Intraplate Earthquakes and Postglacial Rebound in Eastern Canada ...
    Jul 3, 2019 · The causal relationship between postglacial rebound and intraplate earthquakes in Eastern Canada and Northern Europe is investigated with the ...
  100. [100]
    (PDF) Recent, tectonically induced, surficial stress-relief structures in ...
    Aug 10, 2025 · There has been a debate over whether recent faulting occurs as a result of tectonic stress or glacio-isostatic adjustment, or both (e.g., Adams, ...
  101. [101]
    Intraplate seismicity in northern Central Europe is induced by the ...
    Jul 1, 2015 · There is growing evidence that climate-induced melting of large ice sheets has been able to trigger fault reactivation and earthquakes ...
  102. [102]
    Time-variable strain and stress rates induced by Holocene glacial ...
    May 5, 2023 · Large-scale non-tectonic processes, (eg, ice loading) may induce time-variable deformation across wide areas, beyond the directly loaded region.1. Introduction · 2. Modelling Approach · 4. Regional Examples
  103. [103]
    Evidence for the release of long‐term tectonic strain stored in ...
    Jun 20, 2016 · We argue that the tectonically stable continental lithosphere can store elastic strain on long timescales, the release of which may be triggered ...
  104. [104]
    Effects of present‐day deglaciation in Iceland on mantle melt ...
    Jul 2, 2013 · Melting will occur everywhere the decompression lowers the solidus temperature below the local mantle temperature. Melts generated outside ...
  105. [105]
    An assessment of potential causal links between deglaciation and ...
    Apr 11, 2023 · It is critical to determine whether the decompression of crustal magma systems via deglaciation has resulted in enhanced eruption rates along volcanic arcs.
  106. [106]
    Tectonism and volcanism enhanced by deglaciation events in ...
    Nov 22, 2019 · Occurrence of volcanism and seismicity in Iceland is commonly related to rifting events. Subglacial volcanic events seem moreover to have been ...<|separator|>
  107. [107]
    https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2025GC012290
  108. [108]
    [PDF] A revised crustal stress orientation database for Canada
    These deviations from the overall stress pattern can be possibly explained by post-glacial ... large wavelengths (N1000 km) in the area of the Canadian Shield, ...
  109. [109]
    Crustal stresses in Canada - GeoScienceWorld
    Post-glacial near-surface stresses were dominated by the radial ... Within the Canadian Shield there are stress anomalies in the basement rocks ...
  110. [110]
    [PDF] Glacial isostatic adjustment strain rate - stress paradox in the ...
    Apr 19, 2023 · GIA explains uplift and deformation, but its stress perturbations don't drive observed seismicity, creating a strain rate-stress paradox due to ...Missing: deviatoric | Show results with:deviatoric
  111. [111]
    Global Mean Sea Level from TOPEX & Jason Altimetry
    A correction for Glacial Isostatic Adjustment (GIA) is also applied (increasing the rate of global mean sea level by 0.3 mm/year) to account for an increase in ...
  112. [112]
    Global sea level rise and glacial isostatic adjustment - ScienceDirect
    The ongoing global process of glacial isostatic adjustment (GIA) contributes significantly to present-day observed rates of secular sea level change.
  113. [113]
    What is glacial isostatic adjustment (GIA), and why do you correct for ...
    Removing known components of sea level change, such as GIA or the solid earth and ocean tides, reveals the remaining signals contained in the altimetry ...Missing: isolation | Show results with:isolation
  114. [114]
    Separating Geophysical Signals Using GRACE and High ...
    Nov 14, 2018 · We demonstrate the utility of the developed approach by separating present-day ice-mass changes and glacial isostatic adjustment (GIA) in ...
  115. [115]
    Glacial isostatic adjustment as a control on coastal processes
    Mar 9, 2017 · We conclude that regional variations in relative sea-level change driven by glacial isostatic adjustment should be considered when interpreting ...Missing: variance | Show results with:variance<|separator|>
  116. [116]
    Constraining models of glacial isostatic adjustment in eastern North ...
    Jun 15, 2024 · We constrain a set of optimal Earth-ice models for eastern Canada and northeastern US through 1D (radially varying viscosity structure) GIA modelling.
  117. [117]
    [PDF] The glacial isostatic adjustment signal at present day in northern ...
    Jun 15, 2018 · D1 and D3 have respective peak uplift rates of 9.8 and. 9.2 mm yr−1. ... an improved postglacial land uplift model over the Nordic-. Baltic ...<|separator|>
  118. [118]
    [PDF] Uncertainty Estimation in Regional Models of Long‐Term GIA Uplift ...
    In Scandinavia, all models predict peak uplift rates in the Gulf of Bothnia of up to 8 mm/yr, although the spatial location of peak uplift differs slightly.
  119. [119]
    Subsidence along the Atlantic Coast of North America: Insights from ...
    Mar 17, 2016 · The present-day subsidence rates in these areas are approximately double the long-term geologic rates, which has important implications for ...
  120. [120]
    Study: From NYC to D.C. and beyond, cities on the East Coast are ...
    Jan 2, 2024 · These groundbreaking new maps show that a large area of the East Coast is sinking at least 2 mm per year, with several areas along the mid- ...Missing: peripheral bulge
  121. [121]
    Sea-Level Fingerprints Due to Present-Day Water Mass ... - MDPI
    Nov 19, 2021 · Here, we attempt to detect the spatial trend patterns of the fingerprints associated with present-day land ice melt and terrestrial water mass changes.<|separator|>
  122. [122]
    Sea-level changes at the LGM from ice-dynamic reconstructions of ...
    Although our results provide unambiguous evidence that the post LGM rise of eustatic sea-level was very close to the widely supported estimate of 120 m, the ...<|separator|>
  123. [123]
    [PDF] Late Quaternary meltwater pulses and sea level change
    During the largest and most rapid of these, MWP1A, rates of GMSL rise reached approximately 4 m per century (Hanebuth et al., 2000; Peltier and. Page 7. 6.<|control11|><|separator|>
  124. [124]
    Why Is Sea Level Rising Faster in Some Places Along the US East ...
    Dec 19, 2018 · The research team found that post-glacial rebound accounted for most of the variation in sea level rise along the East Coast. But ...Missing: milestones | Show results with:milestones
  125. [125]
    A Global Perspective on Local Sea Level Changes - MDPI
    Glacial Isostatic Adjustment or GIA is the process where the earth mantle is recovering from the presence of an iceshelf that has since disappeared [124].
  126. [126]
    Time and Tide Gauges wait for no Voortman - RealClimate
    Sep 18, 2025 · Actual analyses of the tide gauge record show that acceleration is sea level rise is not only widespread, but it is increasingly clear: Analysis ...
  127. [127]
    Chapter 9: Ocean, Cryosphere and Sea Level Change
    GIA uncertainty is caused by uncertainty in the rheological structure of the solid Earth, which drives the longer-term viscous Earth deformation, as well as ...Missing: criticisms | Show results with:criticisms
  128. [128]
    Uncertainties of Glacial Isostatic Adjustment Model Predictions in ...
    Apr 17, 2020 · We quantify GIA prediction uncertainties of 250 1D and 3D glacial isostatic adjustment (GIA) models through comparisons with deglacial relative sea-level (RSL) ...
  129. [129]
    [PDF] Predicting Present-Day Rates of Glacial Isostatic Adjustment Using a ...
    NAD83 is the official civil reference frame in Canada and the US. ▫ Defined as a time-dependent 7 parameter transformation from ITRF. ▫ Different realizations ...Missing: corrections | Show results with:corrections
  130. [130]
    [PDF] NAD83v70VG: A new national crustal velocity model for Canada
    It incorporates GPS observations with the crustal uplift predictions of Glacial Isostatic Adjustment (GIA) and elastic rebound models, which are especially.
  131. [131]
    (PDF) An improved and extended GPS-derived 3D velocity field of ...
    Aug 7, 2025 · We present a new GPS-derived 3D velocity field for the Fennoscandia glacial isostatic adjustment (GIA) area. This new solution is based upon ∼ ...Missing: ETRF | Show results with:ETRF<|separator|>
  132. [132]
    [PDF] Geodetic Fixing of Tide Gauge Bench Marks
    leveling. The strategy for measuring the effects of glacial rebound is to: Establish networks of absolute gravity stations in the regions of maximum ...<|control11|><|separator|>
  133. [133]
    Deformation of the North American plate interior from a decade of ...
    Jun 10, 2006 · The principal strains show a band of approximately north-south shortening at a rate on the order of 10−9 yr−1 in the northeastern United States ...
  134. [134]
    The impact of late Cenozoic uplift and erosion on hydrocarbon ...
    Several enhancing effects of the uplift are documented, including the development of fracture permeability in reservoirs, the remigration of hydrocarbons to ...
  135. [135]
    Fault reactivation, leakage potential, and hydrocarbon column ...
    Aug 6, 2025 · Fault reactivation and hydrocarbon leakage in this area appears to be caused by three factors: (1) locally elevated pore pressure due to buoyant ...
  136. [136]
    [PDF] Implications of Glacial Isostatic Adjustment on Petroleum Reservoirs ...
    Existing models of Glacial Isostatic Adjustment are used in conjunction with present-day geophysical observations to simulate differential vertical crustal ...Missing: imaging | Show results with:imaging
  137. [137]
    Fault reactivation, leakage potential, and hydrocarbon column ...
    Fault reactivation and leakage are caused by elevated pore pressure, optimal fault orientation, and postglacial rebound. Leakage potential influences  ...
  138. [138]
    A hierarchical Constrained Bayesian (ConBay) approach to jointly ...
    A hierarchical Constrained Bayesian (ConBay) approach to jointly estimate water storage and Post-Glacial Rebound from GRACE(-FO) and GNSS data · Full Article ...
  139. [139]
    [PDF] Post-Glacial Isostatic Adjustment and Global Warming in Subarctic ...
    Nov 23, 2009 · ABSTRACT. When Rupert's Land and the North-Western Territory became a part of Canada as the Northwest Territories.<|separator|>
  140. [140]
    the science legacy of Canada's extended continental shelf mapping ...
    This paper highlights a number of segments of Canada's continental margins to showcase this scientific evidence and how it is applied in the UNCLOS context.
  141. [141]
    [PDF] Viability of UNCLOS amid Emerging Global Maritime Challenges
    UNCLOS adopted in 1982 a modern dispute settlement system, which was tailored to the needs of an efficient ocean governance as it was seen 40 years ago and ...