Cited by 6Geological and Tectonic Studies in Latin AmericaCryospheric studies and observationsGeology and Paleoclimatology Research

Fast uplift in the southern Patagonian Andes due to long- and short-term deglaciation and the asthenospheric window underneath

Open access
Muller, Veleda A. P., Sternai, Pietro, Sue, Christian·2024·Solid Earth

Published version of a preprint

Abstract

Abstract. An asthenospheric window underneath much of the South American continent increases the heat flow in the southern Patagonian Andes where glacial–interglacial cycles drive the building and melting of the Patagonian Icefields since the latest Miocene. The Last Glacial Maximum (LGM) was reached ∼26 000 yr BP (years before present). Significant deglaciation onsets between 21 000 and 17 000 yr BP were subject to an acceleration since the Little Ice Age (LIA), which was ∼400 yr BP. Fast uplift rates of up to 41±3 mm yr−1 are measured by global navigation satellite system (GNSS) around the Southern Patagonian Icefield and are currently ascribed to post-LIA lithospheric rebound, but the possible longer-term post-LGM rebound is poorly constrained. These uplift rates, in addition, are 1 order of magnitude higher than those measured on other glaciated orogens (e.g. the European Alps), which raises questions about the role of the asthenospheric window in affecting the vertical surface displacement rates. Here, we perform geodynamic thermo-mechanical numerical modelling to estimate the surface uplift rates induced by post-LIA and post-LGM deglaciation, accounting for temperature-dependent rheologies and different thermal regimes in the asthenosphere. Our modelled maximum post-glacial rebound matches the observed uplift rate budget only when both post-LIA and post-LGM deglaciation are accounted for and only if a standard continental asthenospheric mantle potential temperature is increased by 150–200 °C. The asthenospheric window thus plays a key role in controlling the magnitude of presently observed uplift rates in the southern Patagonian Andes.

AI snapshot

Written from the paper’s full text.

TL;DR

Thermo-mechanical numerical models show that the exceptionally high GNSS uplift rates around the Southern Patagonian Icefield require both post-LIA and post-LGM deglaciation, plus a 150–200 °C asthenospheric thermal anomaly. The asthenospheric window lowers mantle viscosity and enhances post-glacial rebound.

Key findings

  • Modelled maximum post-glacial rebound matches the observed uplift-rate budget only when both post-LIA and post-LGM deglaciation are accounted for and only if a standard continental asthenospheric mantle potential temperature is increased by 150–200 °C.
  • Post-LIA deglaciation alone produces maximum uplift rates of ∼14–26 mm yr⁻¹ toward the end of LIA deglaciation with TA = 150–200 °C, but does not reach the observed maximum of 41 ± 3 mm yr⁻¹ even with a very low-viscosity asthenosphere.
  • Post-LGM experiments with a low-viscosity asthenosphere and deglaciation lasting 20 000 and 10 000 years reach up to ∼25 and ∼42 mm yr⁻¹ during the final stages, respectively; the slower post-LGM scenario may contribute up to 40 % to the present-day uplift-rate budget.
  • Asthenosphere viscosity ranges between 10²²–10¹⁹ Pa s for TA = 0–100 °C and between 10¹⁹–10¹⁶ Pa s for TA = 150–200 °C; TA > 200 °C is considered unlikely and would decrease viscosities to unrealistic values below 10¹⁶ Pa s.

Methods

  • Fully coupled thermo-mechanical, visco-elasto-plastic numerical geodynamic model with a 700 km × 120 km domain: 10 km sticky air, 30 km continental crust, 30 km lithospheric mantle, and 50 km asthenospheric mantle; finite differences with marker-in-cell, 51 × 61 nodes, 400 × 400 Lagrangian markers, and free-slip boundaries.
  • Deglaciation forcing uses a 2 km thick pseudo-ice cap: post-LGM Model set 1 imposes 75 % ice loss in 20 000 years (1500 m drop), Model set 2 imposes 95 % ice loss in 10 000 years (1900 m drop), and post-LIA imposes 10 % ice loss in 400 years (200 m drop); asthenospheric thermal anomalies TA up to 200 °C are added in 50 °C steps to a reference potential temperature of 1265 °C.
  • Modelled uplift rates are compared with GNSS-derived data from 31 GPS stations around the Southern Patagonian Icefield published in Lange et al. (2014); uplift rate is calculated as surface elevation change divided by the viscoelastic time step.

Limitations

  • The paper states that the modelling imposes linear and uniform ice loss instead of a realistic spatiotemporal melting pattern, and that information on melting velocities, ice thickness, and redistribution of surface masses is limited; GNSS, SRTM, and GRACE data also show discrepancies.
  • The paper states that it assumes a homogeneous lithosphere, neglects lateral viscosity variations in the asthenosphere, neglects longer-term tectonic deformation, keeps the pseudo-ice cap lateral extent fixed, and treats erosion as a negligible contribution.
  • Evident from the text: the study focuses on the magnitude rather than the spatial pattern of uplift, and the results depend on simplified deglaciation histories and the parameterized asthenospheric thermal anomaly.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.5194/se-15-387-2024

Cites

129 of 131 references verified0 citation errors · 2 not checkable

Statistics

Nothing to recompute.

Claims the paper makes

17 of 89 cited sources checked · 7 supported · 0 flagged · accuracy not yet measured

  • Premisecitation

    The last recorded glacial advance was the Little Ice Age with an apex around 1630 CE.

    “A glacial minimum must have been attained around 13 000 yr BP, but several glacial advances were recorded since that time, and the last recorded one was the Little Ice Age (LIA) with an apex around 1630 CE, well dated by terminal moraines around the present-day NPI and SPI (Ivins andJames, 1999, 2004;McCulloch et al., 2000;Glasser et al., 2004Glasser et al., , 2008Glasser et al., , 2011;;Davies”
    Quote in the paperNot a substantial claim1 of 3 sources checked
    • cites Ivins and James, 1999 source

      The Introduction identifies the Little Ice Age as the last of several Neoglacial advances. It does not date its apex to around 1630 CE: it describes the advance’s timing as uncertain and reports moraine-based maxima in the mid-18th and mid-19th centuries. Section 6.3.6 places the maximum in the late 19th century in the authors’ preferred model. (LLM, conf=0.88)

      Experimental — accuracy not yet measured
    • cites McCulloch et al., 2000 source
    • cites Glasser et al., 2004 source
  • Premisecitation

    A glacial minimum was attained around 13 000 yr BP.

    “A glacial minimum must have been attained around 13 000 yr BP, but several glacial advances were recorded since that time, and the last recorded one was the Little Ice Age (LIA) with an apex around 1630 CE, well dated by terminal moraines around the present-day NPI and SPI (Ivins andJames, 1999, 2004;McCulloch et al., 2000;Glasser et al., 2004Glasser et al., , 2008Glasser et al., , 2011;;Davies”
    Quote in the paperSubstantial claim1 of 3 sources checked
    • cites Ivins and James, 1999 source

      The Introduction (page unavailable) says Patagonia was substantially deglaciated by 13 ka, supporting markedly reduced glaciation around that time. It does not establish that a glacial minimum was attained specifically around 13,000 yr BP; it also says the icefields may have been substantially smaller throughout 13–5 ka. (LLM, conf=0.82)

      Experimental — accuracy not yet measured
    • cites McCulloch et al., 2000 source
    • cites Glasser et al., 2004 source
  • Premisecitation

    Tracing back the post-LGM or Holocene ice loss rate from current measurements is difficult considering that climate was at least 6 °C colder during the LGM.

    “Tracing back the post-LGM or Holocene ice loss rate from current measurements is difficult considering that climate was at least 6 °C colder during the LGM”
    Quote in the paperNot a substantial claim1 of 3 sources checked
    • cites Hulton et al., 2002 source
    • cites Sugden et al., 2002 source
    • cites Seltzer et al., 2021 source

      The Abstract (page unavailable) reports 5.8 ± 0.6 °C LGM cooling for low-altitude land between 45°S and 35°N. That supports roughly 6 °C cooling in the studied region, but not a claim that climate generally was at least 6 °C colder. The paper does not address whether current measurements make post-LGM or Holocene ice-loss rates difficult to reconstruct. (LLM, conf=0.88)

      Experimental — accuracy not yet measured
  • Premisecitation

    Increased asthenospheric temperatures beneath southern Patagonia are highly supported by geophysical data.

    “increased asthenospheric temperatures beneath southern Patagonia are highly supported by the geophysical data”
    Quote in the paperSubstantial claim1 of 3 sources checked
    • cites Russo et al., 2010 source
    • cites Russo et al., 2022

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Mark et al., 2022 source
  • Interpretationreasoning

    Lithospheric thinning due to the asthenospheric window underneath southern Patagonia affects the regional uplift rates.

    “Lithospheric thinning due to the asthenospheric window underneath southern Patagonia thus affects the regional uplift rates as previously suggested”
    Quote in the paperSubstantial claim2 of 3 sources checked
    • cites Ávila and Dávila, 2018, 2020

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Mark et al., 2022 source
    • cites Ben-Mansour et al., 2022 source

      The Introduction says the slab window produces lithospheric thickness anomalies, and Section 4.1 describes thermal and mechanical erosion of the mantle lithosphere. Our copy does not say that this thinning affects regional uplift rates; a reference title mentions implications for glacial isostatic adjustment, but that does not establish the claimed effect. (LLM, conf=0.84)

      Experimental — accuracy not yet measured
  • May combine two claims. 29
  • Classified as probably not a substantial claim. 10
  • The quote only partly states this claim. 2
  • May combine two claims.

    Premisecitation

    An asthenospheric window underneath much of the South American continent increases the heat flow in the southern Patagonian Andes, where glacial-interglacial cycles drive the building and melting of the Patagonian Icefields since the latest Miocene.

    “An asthenospheric window underneath much of the South American continent increases the heat flow in the southern Patagonian Andes where glacial-interglacial cycles drive the building and melting of the Patagonian Icefields since the latest Miocene.”
    Quote in the paperSubstantial claim1 of 1 sources checked
    • cites Ávila and Dávila, 2018 source

      Split into 2 atoms — 0 confirmed, 1 partially supported, 0 disputed, 0 misattributed, 1 unverifiable. Expand the claim to see per-atom verdicts.

      Experimental — accuracy not yet measured
  • Classified as probably not a substantial claim.

    Premisecitation

    Previous GIA studies infer a low asthenosphere viscosity (in the order of 10^18 Pa s) and a thin elastic lithosphere (∼35 km thick) to match the very high observed uplift rate budget.

    “To match the very high observed uplift rate budget, previous GIA studies infer a low asthenosphere viscosity (in the order of 10 18 Pa s) and a thin elastic lithosphere (∼ 35 km thick) (Ivins andJames, 1999, 2004;Klemann et al., 2007;Dietrich et al., 2010;Lange et al., 2014;Richter et al., 2016;Ávila and Dávila, 2020;Mark et al., 2022;Lenzano et al., 2023).”
    Quote in the paperSubstantial claim1 of 3 sources checked
    • cites Ivins and James, 1999 source

      Figure 11 (section 6.4.1; page not supplied) models Patagonian uplift with a 35 km lithosphere and mantle viscosities of 1–2 × 10^19 Pa s; Figure 13 considers a viscosity reduced by half. These are predicted responses, not parameters inferred by matching a very high observed uplift rate. The paper describes the motion as potentially geodetically detectable rather than reporting such an observational fit. (LLM, conf=0.88)

      Experimental — accuracy not yet measured
    • cites Klemann et al., 2007 source
    • cites Dietrich et al., 2010

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Classified as probably not a substantial claim.

    Premisecitation

    The contribution of post-LGM deglaciation to present-day rock uplift rate was marginally addressed.

    “In addition, the contribution of post-LGM deglaciation to present-day rock uplift rate was marginally addressed (Ivins andJames, 1999, 2004;Klemann et al., 2007).”
    Quote in the paperSubstantial claim1 of 2 sources checked
    • cites Ivins and James, 1999 source
    • cites Klemann et al., 2007 source
  • May combine two claims.

    Study design and conductown data

    The rheologic and thermal structure of the reference model gives a lithospheric elastic thickness, Te, of ∼30 km, comparable to previous estimates underneath the SPI.

    “The rheologic and thermal structure of the reference model gives a lithospheric elastic thickness, T e (sensu Burov and Diament, 1995), of ∼ 30 km, comparable to previous estimates underneath the SPI based on GIA models (Ivins and James, 1999;Dietrich et al., 2010;Lange et al., 2014), heat flow data (Ávila and Dávila, 2018), waveform inversion (Robertson Maurice et al., 2003), and low-temperature”
    Quote in the paperNot a substantial claim1 of 3 sources checked
    • cites Ivins and James, 1999 source
    • cites Dietrich et al., 2010

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Lange et al., 2014 source
  • May combine two claims.

    Premisecitation

    Low-temperature thermochronology studies suggest using erosion values between 0.1 and 1 mm yr-1 from 7 to 4 Ma, followed by a period of erosional quiescence (<0.1 mm yr-1) and a possible increase to 1 mm yr-1 in the last ∼2 Myr in the SPI region.

    “Fosdick et al. (2013), Herman and Brandon (2015), Fernandez et al. (2016), andMuller et al. (2023) suggest using values between 0.1 and 1 mm yr -1 from 7 to 4 Ma, followed by a period of erosional quiescence (< 0.1 mm yr -1 ) and a possible increase to 1 mm yr -1 in the last ∼ 2 Myr in the SPI region”
    Quote in the paperSubstantial claim1 of 3 sources checked
    • cites Fosdick et al., 2013

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Herman and Brandon, 2015 source
    • cites Muller et al., 2023 source

      The Conclusions (page unavailable) describe rapid exhumation at Fitz Roy and Torres del Paine between 7 and 4 Ma, subsequent erosional quiescence, and a partly recorded Quaternary acceleration. Section 5.2 (page unavailable) reports 7–4 Ma exhumation rates of 0.6–3 km/Myr (equivalent to mm/yr), not a uniform 0.1–1 mm/yr range; it gives no <0.1 mm/yr rate for the quieter period. It reports ~1 mm/yr Quaternary exhumation in parts of Torres del Paine, with onset varying from 3 to 0.5 Ma, rather than establishing that rate throughout the SPI region over the last ~2 Myr. (LLM, conf=0.88)

      Experimental — accuracy not yet measured
  • Classified as probably not a substantial claim.

    Study design and conductcitation

    Given the geologically short time window investigated here, deformation related to longer-term tectonic forces is neglected.

    “Given the geologically short time window investigated here, we neglect deformation related to longer-term tectonic forces (Breitsprecher and Thorkelson, 2009;Guillaume et al., 2013;Eagles and Scott, 2014;Muller et al., 2021).”
    Quote in the paperSubstantial claim1 of 3 sources checked
    • cites Breitsprecher and Thorkelson, 2009

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Guillaume et al., 2013 source
    • cites Eagles and Scott, 2014 source
  • May combine two claims.

    Findingcitation

    The asthenosphere viscosity from the models when TA > 100 °C is <10^19 Pa s, with the lowest viscosity value of 10^16 Pa s imposed where partial melting occurs.

    “the asthenosphere viscosity from our models when TA > 100 °C is < 10 19 Pa s, with the lowest viscosity value of 10 16 Pa s imposed where partial melting, supported by the regional Holocene volcanism (Stern and Kilian, 1996) and by geophysical data (e.g.shear wave velocity data by Mark et al., 2022), occurs.”
    Quote in the paperSubstantial claim1 of 2 sources checked
    • cites Stern and Kilian, 1996 source
    • cites Mark et al., 2022 source

      Section 3.1 reports anomalously low shear-wave velocities in the Patagonian slab window, so Mark et al. provide the geophysical data cited. Section 3.2 says thermal erosion likely requires melt or fluids and that melt may be supplied by upwelling mantle; it does not establish that partial melting occurs at the specific locations where the citing authors impose their lowest viscosity. The viscosity values are attributed to the citing authors’ own models, not to Mark et al. (LLM, conf=0.84)

      Experimental — accuracy not yet measured
  • May combine two claims.

    Findingown data

    In the post-LIA model set, the maximum uplift rate reaches ∼8.3 and ∼23 mm yr-1 during the first 100 years of deglaciation when TA is respectively 150 and 200 °C.

    “whereas it reaches ∼ 8.3 and ∼ 23 mm yr -1 during the same interval when TA is respectively 150 and 200 °C”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The asthenosphere viscosity ranges between 10^22-10^19 Pa s in simulations with TA equal to 0, 50, and 100 °C and between 10^19-10^16 Pa s in simulations with TA equal to 150 and 200 °C.

    “The asthenosphere viscosity ranges between 10 22 -10 19 Pa s in simulations with TA equal to 0 (reference model), 50, and 100 °C and between 10 19 -10 16 Pa s in simulations with TA equal to 150 and 200 °C (Fig. 3a-d).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingreasoning

    This results in predominant elastic deformation in the upper crust (below the ∼300 °C isotherm) and upper mantle lithosphere (below the ∼700 °C isotherm) and viscous deformation in the lower crust and lower lithospheric mantle and asthenosphere.

    “This results in predominant elastic deformation in the upper crust (below the ∼ 300 °C isotherm) and upper mantle lithosphere (below the ∼ 700 °C isotherm) and viscous deformation in the lower crust and lower lithospheric mantle and asthenosphere”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    After the end of the deglaciation, the maximum uplift rate takes longer than about 5000 years to re-equilibrate to 0 mm yr−1 when TA ≤ 100 °C, whereas it drops to 0 mm yr−1 almost immediately when TA is 150 or 200 °C.

    “After the end of the deglaciation, the maximum uplift rate takes longer than about 5000 years to re-equilibrate to 0 mm yr -1 when TA ≤ 100 °C, whereas it drops to 0 mm yr -1 almost immediately when TA is 150 or 200 °C (Fig. 7a).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In Model set 1 for post-LGM deglaciation, when TA equals 0 (reference model) the maximum uplift rate is <1 mm yr−1 during the first 5000 years of the deglaciation, increasing gradually up to 9.5 mm yr−1 in the later stages of the deglaciation (i.e. 20 000 years).

    “In Model set 1 for post-LGM deglaciation, when TA equals 0 (reference model) the maximum uplift rate is < 1 mm yr -1 during the first 5000 years of the deglaciation, increasing gradually up to 9.5 mm yr -1 in the later stages of the deglaciation (i.e.20 000 years, Fig. 4).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingreasoning

    Lithospheric warming due to increasing asthenospheric mantle potential temperature leads to a reduction in the lower lithosphere viscosity (from 10^22 to 10^20 Pa s), thereby decreasing the integrated lithospheric strength.

    “Lithospheric warming due to increasing asthenospheric mantle potential temperature also leads to a reduction in the lower lithosphere viscosity (from 10 22 to 10 20 Pa s), thereby decreasing the integrated lithospheric strength.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In the post-LIA model set, the maximum uplift rate is ∼1.4, ∼2.3, and ∼2.2 mm yr-1 during the first 100 years of deglaciation when TA is respectively 0, 50, and 100 °C.

    “In the post-LIA model set, the maximum uplift rate is ∼ 1.4, ∼ 2.3, and ∼ 2.2 mm yr -1 during the first 100 years of deglaciation when TA is respectively 0, 50, and 100 °C”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    When TA equals 50, 100, 150, and 200 °C, the maximum uplift rates can reach up to ∼2, ∼5, ∼12, and ∼15 mm yr−1, respectively, already in the first 1000 years of the deglaciation.

    “When TA equals 50, 100, 150, and 200 °C, the maximum uplift rates can reach up to ∼ 2, ∼ 5, ∼ 12, and ∼ 15 mm yr -1 , respectively, already in the first 1000 years of the deglaciation (Fig. 4a).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingreasoning

    When higher temperatures are imposed in the asthenospheric mantle, shallower 300 and 700 °C isotherms decrease Te and increase the isostatic surface uplift rates.

    “when we impose higher temperatures in the asthenospheric mantle, shallower 300 and 700 °C isotherms decreases T e and increases the isostatic surface uplift rates.”
    Quote in the paperSubstantial claim
  • The quote only partly states this claim.

    Findingown data

    The modelled maximum post-glacial rebound matches the observed uplift rate budget only if a standard continental asthenospheric mantle potential temperature is increased by 150-200 °C.

    “post-LIA and post-LGM deglaciation are accounted for and only if a standard continental asthenospheric mantle potential temperature is increased by 150-200 °C.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingcitation

    The deglaciation triggers uplift in the region covered by the melting pseudo-ice cap and subsidence in the neighbouring regions.

    “In agreement with the theory of lithospheric flexure (e.g.Turcotte and Schubert, 2002), the deglaciation triggers uplift in the region covered by the melting pseudo-ice cap and subsidence in the neighbouring regions (Figs.456).”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Turcotte and Schubert, 2002 source
  • May combine two claims.

    Findingown data

    The elastic thickness of the lithosphere (Te) varies between the simulations according to the imposed asthenospheric thermal anomaly, but it is generally lower than 30 km, resulting in a decoupled lithospheric rheology.

    “The elastic thickness of the lithosphere (T e ) varies between the simulations according to the imposed asthenospheric thermal anomaly, but it is generally lower than 30 km, resulting in a decoupled lithospheric rheology”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Burov and Diament, 1995 source
  • May combine two claims. The quote only partly states this claim.

    Findingown data

    Experiments that account for a low viscosity asthenosphere and long-term post-LGM deglaciation lasting for 20 000 and 10 000 years reach up to ∼25 and ∼42 mm yr-1 of uplift rate during the final stages of the deglaciation, respectively, comparable to present-day values.

    “Experiments that account for a low viscosity asthenosphere and long-term post-LGM deglaciation lasting for 20 000 and 10 000 years reach up to ∼ 25 and ∼ 42 mm yr -1 of an uplift rate during the final stages of the deglaciation”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In the faster post-LGM deglaciation scenario (Model set 2), the observed maximum uplift rate budget is attained in about 10 000 years of deglaciation, but only minor residual rebound could be observed today regardless of the amount of ice loss.

    “In the faster post-LGM deglaciation scenario (Model set 2) the observed maximum uplift rate budget is attained in about 10 000 years of deglaciation, but only minor residual rebound could be observed today regardless of the amount of ice loss”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisecitation

    The Chile triple junction at ∼46° S delimits the surface tip of the asthenospheric window, which opened during the last ∼16 Myr from south to north.

    “The Chile triple junction (CTJ) at ∼ 46° S delimits the surface tip of the asthenospheric window, which opened during the last ∼ 16 Myr from south to north (Ramos and Kay, 1992;Breitsprecher and Thorkelson, 2009).”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites Ramos and Kay, 1992 source
    • cites Breitsprecher and Thorkelson, 2009

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Premisecitation

    The Patagonian Ice Sheet covered the southern Patagonian Andes between ∼47 000 and ∼17 000 yr BP, extending from latitudes 38 to 55° S with an estimated area of ∼490 000 km2, volume of ∼550 000 km3, and average and maximum thickness of 1100 and 2500 m.

    “The Patagonian Ice Sheet covered the southern Patagonian Andes between ∼ 47 000 and ∼ 17 000 yr BP, extending from latitudes 38 to 55° S with an estimated area of ∼ 490 000 km 2 (Fig. 1a); volume of ∼ 550 000 km 3 ; and average and maximum thickness of 1100 and 2500 m, respectively, based on preserved glacial geomorphology, stratigraphy, palaeoecology, and geochronological data (Moreno et al., 199”
    Quote in the paperSubstantial claim0 of 3 sources checked
    • cites Moreno et al., 1999

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Hulton et al., 2002 source
    • cites Glasser et al., 2004 source
  • Classified as probably not a substantial claim.

    Premisecitation

    The high uplift rates in the northern part of the SPI are currently ascribed to lithospheric viscoelastic GIA following the LIA, which was responsible for an ice loss of 503 ± 101.1 km3 in the SPI.

    “Such outstandingly high uplift rates, especially in the northern part of the SPI, are currently ascribed to the lithospheric viscoelastic GIA following the LIA, which was responsible for an ice loss of 503 ± 101.1 km 3 in the SPI (Glasser et al., 2011).”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Glasser et al., 2011 source
  • May combine two claims.

    Premisecitation

    GNSS-measured data show ongoing vertical rock uplift rates between 18 ± 3 and 41 ± 3 mm yr−1 in the northern part (18-50.5°S) of the SPI, decreasing to values between 2 ± 6 and 17 ± 5 mm yr−1 in its southern part (50.5-51.5°S).

    “GNSS-measured data show ongoing vertical rock uplift rates between 18 ± 3 and 41 ± 3 mm yr -1 in the northern part (18-50.5°S) of the SPI (Fig. 1b), decreasing to val-ues between 2 ± 6 and 17 ± 5 mm yr -1 in its southern part (50.5-51.5°S) (Ivins and James, 2004;Dietrich et al., 2010;Lange et al., 2014;Richter et al., 2016;Lenzano et al., 2023).”
    Quote in the paperSubstantial claim0 of 3 sources checked
    • cites Ivins and James, 2004 source
    • cites Dietrich et al., 2010

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Lange et al., 2014 source
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisenone given

    The forcing of global cooling in increasing erosion rates during the Quaternary is debated, and not widely quantified in Patagonia nor worldwide.

    “The forcing of global cooling in increasing erosion rates during the Quaternary is, however, debated, and not widely quantified in Patagonia nor worldwide”
    Quote in the paperNot a substantial claim0 of 3 sources checked
    • cites Valla et al., 2012

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Champagnac et al., 2014

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Herman et al., 2013 source
  • Classified as probably not a substantial claim.

    Study design and conductnone given

    The simplification of a constant lateral extent of the pseudo-ice cap does not significantly affect the magnitude of post-glacial rebound.

    “Although this simplification may affect the inferred pattern of post-glacial rebound, it greatly facilitates the simulation of deglacial lithospheric unloading without significantly affecting the magnitude of post-glacial rebound, which is the main focus here.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Interpretationreasoning

    Supposing that these erosion rates still apply in the last ∼20 000 years, this would translate into 2-20 m of rocks eroded on average since the LGM, leading to local unloading of approximately 60-600 kPa if one assumes a crustal density of 3000 kg m-3.

    “Supposing that these erosion rates still apply in the last ∼ 20 000 years, this would translate into 2-20 m of rocks eroded on average since the LGM, leading to local unloading of approximately 60-600 kPa if one assumes a crustal density of 3000 kg m -3”
    Quote in the paperNot a substantial claim
    • The authors state this as an assumption: The erosion rates from low-temperature thermochronology still apply in the last ∼20 000 years.
    • The authors state this as an assumption: The crustal density is 3000 kg m-3.
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisereasoning

    Previous models have assumed simple deglaciation histories as well.

    “As a result, previous models have assumed simple deglaciation histories as well”
    Quote in the paperNot a substantial claim0 of 3 sources checked
    • cites Ivins and James, 1999, 2004

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Hulton et al., 2002 source
    • cites Klemann et al., 2007 source
  • May combine two claims.

    Study design and conductown data

    Overall, notwithstanding these limitations in the model, the fully coupled numerical thermomechanical geodynamic experiments provide realistic uplift rates that one can compare to current geodetic observations.

    “Overall, notwithstanding these limitations in the model, our fully coupled numerical thermomechanical geodynamic experiments provide realistic uplift rates (Figs.4567) that one can compare to current geodetic observations.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationreasoning

    Such stress change is approximately equivalent to the melting of about 6-60 m of ice, whereas the simulations simulate the melting of 200-1500 m of ice.

    “Such stress change is approximately equivalent to the melting of about 6-60 m of ice, whereas we simulate the melting of 200-1500 m of ice in our simulations.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Interpretationown data

    Overall, a warmer and less viscous asthenosphere generates a higher magnitude post-glacial rebound which takes much longer to re-equilibrate to 0 mm yr-1 after the end of the deglaciation than a cooler and more viscous asthenosphere.

    “Overall, a warmer and less viscous asthenosphere generates a highermagnitude post-glacial rebound which, however, takes much longer to re-equilibrate to 0 mm yr -1 after the end of the deglaciation than a cooler and more viscous asthenosphere.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    GNSS, SRTM, and GRACE data, constraining the net ice mass balance only during the last few decades, still show some discrepancies.

    “GNSS, SRTM, and GRACE data, constraining the net ice mass balance only during the last few decades, still showing some discrepancies”
    Quote in the paperNot a substantial claim0 of 3 sources checked
    • cites Aniya, 1996 source
    • cites Aniya et al., 1997 source
    • cites Rignot et al., 2003 source
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductreasoning

    Our modelling is simplistic in that we impose a linear and uniform ice loss instead of a more realistic ice sheet melting pattern in space and time.

    “Our modelling is simplistic in that we impose a linear and uniform ice loss instead of a more realistic ice sheet melting pattern in space and time”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Interpretationown data

    Asthenospheric thermal anomalies higher than 200 °C are unlikely and would decrease the asthenospheric viscosities to unrealistic values (less than 10^16 Pa s).

    “Asthenospheric thermal anomalies higher than 200 °C are unlikely and would decrease the asthenospheric viscosities to unrealistic values (less than 10 16 Pa s).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationreasoning

    Using 10^16 < v < 10^18 Pa s and λ = 200 km leads to ∼2000 < τr <∼200 000 years, a time range considerably longer than the post-LIA deglaciation and including full Pleistocene glacial-interglacial cycles.

    “Using 10 16 < v < 10 18 Pa s and λ = 200 km leads to ∼ 2000 < τ r <∼ 200 000 years, a time range considerably longer than the post-LIA deglaciation and including full Pleistocene glacial-interglacial cycles”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites Turcotte and Schubert, 2002 source
    • cites Ruddiman et al., 1986 source
  • May combine two claims.

    Interpretationreasoning

    It appears that post-LIA rebound alone cannot cover the entire budget of the observed uplift rates even with the highest tested TA, which points to a non-negligible contribution from post-LGM deglaciation.

    “it appears that post-LIA rebound alone cannot cover the entire budget of the observed uplift rates even with the highest tested TA, which points to a non-negligible contribution from post-LGM deglaciation.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationown data

    Results indicate that the outstanding observational budget of rock uplift in the SPI is matched only when accounting for higher-than-normal asthenospheric mantle temperatures, thereby highlighting the relevance of the regional asthenospheric window.

    “Results, therefore, indicate that the outstanding observational budget of rock uplift in the SPI is matched only when accounting for higher-thannormal asthenospheric mantle temperatures, thereby highlighting the relevance of the regional asthenospheric window.”
    Quote in the paperSubstantial claim
38 of 38
  • The authors state this as an assumption.

    The eroded material since the LGM is still in the transport zone and therefore does not significantly contribute to unloading the surface of the orogen.

    “it seems reasonable to assume that the eroded material is still in the transport zone and therefore does not significantly contribute to unloading the surface of the orogen.”
  • The authors state this as an assumption.

    The erosion rates from low-temperature thermochronology still apply in the last ∼20 000 years.

    “Supposing that these erosion rates still apply in the last ∼ 20 000 years”
  • The authors state this as an assumption.

    The crustal density is 3000 kg m-3.

    “if one assumes a crustal density of 3000 kg m -3”
  • The authors state this as an assumption.

    Long-term erosion rates contribute similarly to regional uplift rates as in the European Alps, generally a fraction of a mm yr-1, and this contribution is negligible in the context of the southern Patagonian Andes.

    “we assume a similar contribution to regional uplift rates (i.e.generally a fraction of a mm yr -1 ; Sternai et al., 2019); this is a negligible contribution in the context of the southern Patagonian Andes.”
  • The authors state this as an assumption.

    The lithosphere is homogeneous and lateral viscosity variations in the asthenosphere can be neglected.

    “We also assume a homogeneous lithosphere and neglect lateral viscosity variations in the asthenosphere”
  • The authors state this as an assumption.

    GNSS-measured rock uplift rates are mostly related to the deglaciation history and only marginally controlled by longer-term geodynamics.

    “we discuss our results assuming that GNSS-measured rock uplift rates are mostly related to the deglaciation history and only marginally controlled by the longer-term geodynamics”
6 of 6

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

  • Claimed 25 mm yr−1 uplift exceeds Model set 1 results

    Major

    The discussion says the 20,000-year post-LGM experiments reach about 25 mm yr−1 in their final stages. Table 2(a) reports a highest end-of-deglaciation rate of 16.26 mm yr−1, and none of its listed rates reaches 25 mm yr−1.

  • Stated range of simulated ice melt omits Model set 2

    Minor

    The discussion says the simulations melt 200-1500 m of ice, but the methods specify a 1900 m drop for Model set 2.

  • First-1,000-year uplift rates disagree with Table 2

    Minor

    For Model set 1, the text gives approximately 2, 5, and 15 mm yr−1 at thermal anomalies of 50, 100, and 200 °C within the first 1,000 years. Table 2(a) lists 0.56, 3.58, and 11.48 mm yr−1, respectively, at that time.

3 of 3

Cites

Flagged first, then by citations

Research works printed but not linked: 2

  1. Harvey, A. H.: Properties of Ice and Supercooled Water, in: CRC Handbook of Chemistry and Physics, 97th Edn., edited by: Haynes, W. M., Lide, D. R., and Bruno, T. J., CRC Press, Boca Raton, FL, ISBN 978-1-4987-5429-3, 2017.
  2. Ranalli, G.: Rheology of the Earth, Springer Science &amp;amp; Business Media, ISBN 0-412-54670-1, 1995.

Cited by

Sorted by citations

Authors
Muller, V. A. P., Sternai, P., & Sue, C.
Affiliations
  • University of Milano-Bicocca (IT)
  • Centre National de la Recherche Scientifique (FR)
Corresponding author
  • Sternai, Pietro
  • Sue, Christian
Venue
Solid Earth
Published
2024-04-03
Volume
15
Issue
4
Pages
387-404
ISSN
  • 1869-9510 (ISSN-L)
  • 1869-9529
OpenAlex
W4393901312
Cited by
6
References
124
Citations by year
2026: 1 · 2025: 3 · 2024: 2
FWCI
1.58
Type
journal article
Language
en
Open access
Locations
  • Solid Earth (publishedVersion, cc-by)
  • BOA (University of Milano-Bicocca) (submittedVersion, cc-by)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
Funders
Fondazione Cariplo (Grant n° 2022 – 1546_001)
Topics
  • Physical Sciences › Earth and Planetary Sciences › Geophysics › Geological and Tectonic Studies in Latin America
  • Physical Sciences › Earth and Planetary Sciences › Atmospheric Science › Cryospheric studies and observations
  • Physical Sciences › Earth and Planetary Sciences › Atmospheric Science › Geology and Paleoclimatology Research
Keywords
postglacial rebound, deglaciation, Southern Patagonian Icefield, uplift rate, thermomechanical modeling, mantle potential temperature, Little Ice Age, Last Glacial Maximum, glacial isostatic adjustment, Patagonian Icefields, temperature-dependent rheology, heat flow

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Claims read by deepseek-v4.1-flash under definition version 3 (extract-v3) on 2026-10-18.

Classified by jev-1.13-20260917 (gold-guideline-v0.2) on 2026-10-18.

Cited sources are checked by gpt-6-sol.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review