Variations in Earth's 1D viscosity structure in different tectonic regimes

Earth’s long-wavelength geoid provides insights into the thermal, structural, and compositional evolution of the mantle. Historically, most estimates of mantle viscosity using the long-wavelength geoid have considered radial variations with depth in a symmetric Earth. Global estimates of this kind suggest an increase in viscosity from the upper mantle to lower mantle of roughly 2 -- 3 orders of magnitude. Using a spatio-spectral localization technique with the geoid, here we estimate a series of locally constrained viscosity-depth profiles covering two unique regions, the Pacific and Atlantic hemispheres, which show distinct rheological properties. The Pacific region exhibits the conventional Earth's 1D rheology with a factor of roughly 80-100 increase in viscosity occurring at transition zone depths (400 - 800 km). The Atlantic region in contrast does not show significant viscosity jumps with depth, and instead has a near uniform viscosity in the top 1000~km. The inferred viscosity variations between our two regions could be due to the prevalence of present-day subduction in the Pacific and the infrequence of slabs in the Atlantic, combined with a possible hydrated transition zone and mid-mantle of the Atlantic region by ancient subduction during recent supercontinent cycles. Rigid slab material within the top 800 km, with about 90\% Majoritic garnet in the form of subducted oceanic crust, coupled with unique regional mantle structures, may be generating a strong transition zone viscosity interface for the Pacific region. These effective lateral variations in mantle viscosity could play a role in the observed deformation differences between the Pacific and Atlantic hemispheres.

Regional mantle viscosity constraints for North America reveal upper mantle strength differences across the continent

We present regional constraints of mantle viscosity for North America using a local Bayesian joint inversion of mantle flow and glacial isostatic adjustment (GIA) models. Our localized mantle flow model uses new local geoid kernels created via spatio-spectral localization using Slepain basis functions, convolved with seismically derived mantle density to calculate and constrain against the regional free-air gravity field. The joint inversion with GIA uses two deglaciation of ice sheet models (GLAC1D-NA and ICE-6G-NA) and surface relative sea level data. We solve for the local 1D mantle viscosity structure for the entire North America (NA) region, the eastern region including Hudson Bay, and the western region of North America extending into the Pacific plate. Our results for the entire NA region show one order of magnitude viscosity jump at the 670 km boundary using a high seismic density scaling parameter (e.g., δlnp/δlnvs = 0.3). Seismic scaling parameter demonstrates significant influence on the resulting viscosity profile. However, when the NA region is further localized into eastern and western parts, the scaling factor becomes much less important for dictating the resulting upper mantle viscosity characteristics. Rather the respective local mantle density heterogeneities provide the dominate control on the upper mantle viscosity. We infer local 1D viscosity profiles that reflect the respective tectonic settings of each region's upper mantle, including a weak and shallow asthenosphere layer in the west, and deep sharp viscosity jumps in the eastern transition zone, below the suggested/proposed depth range of the eastern continental root.

Gravity constraints on structure of the East-West Antarctic lithospheric transition zone

<p><b>The differing structural evolution of cratonic East Antarctica and younger West Antarctica has resulted in contrasting lithospheric and asthenospheric mantle viscosities between the two regions. Combined with poor constraints on the upper mantle viscosity structure of the continent, estimates of surface uplift in Antarctica predicted from models of glacial isostatic adjustment (GIA) and observed by Global Satellite Navigation System (GNSS) contain large misfits. This thesis presents a gravity study ofthe lithospheric transition zone beneath the Taylor Valley, Antarctica, conducted to constrain the variation in lithological parameters such as viscosity and density of the upper mantle across this region.</b></p> <p>During this study 119 new gravity observations were collected in the ice-free regions of the Taylor Valley and amalgamated with 154 existing land-based gravity observations, analysed alongside aerogravity measurements of southern Victoria Land. Gravity data are used to construct 2D gravity models of the subsurface beneath this region. An eastward gradient in Bouguer anomalies of ~- 1.6 mGal/km is observed within the Taylor Valley. Models reveal thickening of the Moho from 23±5 km beneath the Ross Sea to 35±5 km in the Polar Plateau (dipping at 24.5±7.2°), and lithospheric mantle 100 km thicker in East Antarctica (~200±30 km) than West Antarctica (~90±30 km). </p> <p>Models of predicted surface uplift history are used to estimate an asthenospheric mantle viscosity of 2.1x1020 Pa.s at full surface recovery beneath the Ross Embayment, differing by ~14% from the viscosity at 50% recovery. The temperature contrast between lithospheric and asthenospheric mantle is estimated as ~400°C, equivalent to a viscosity that decreases by a factor of about 30 over the mantle boundary.</p> <p>Results demonstrate that the history of surface uplift in the study area may be complicated, resulting in observations of uplift, or subsidence, at GNSS stations. Future work should incorporate additional geophysical methods, such as seismicity and electrical resistivity, improving constraints on gravity models. A better understanding of the surface uplift (or subsidence) history in the Transantarctic Mountains is critical, with implications in reducing uncertainty in GIA models.</p>

Gravity constraints on structure of the East-West Antarctic lithospheric transition zone

<p><b>The differing structural evolution of cratonic East Antarctica and younger West Antarctica has resulted in contrasting lithospheric and asthenospheric mantle viscosities between the two regions. Combined with poor constraints on the upper mantle viscosity structure of the continent, estimates of surface uplift in Antarctica predicted from models of glacial isostatic adjustment (GIA) and observed by Global Satellite Navigation System (GNSS) contain large misfits. This thesis presents a gravity study ofthe lithospheric transition zone beneath the Taylor Valley, Antarctica, conducted to constrain the variation in lithological parameters such as viscosity and density of the upper mantle across this region.</b></p> <p>During this study 119 new gravity observations were collected in the ice-free regions of the Taylor Valley and amalgamated with 154 existing land-based gravity observations, analysed alongside aerogravity measurements of southern Victoria Land. Gravity data are used to construct 2D gravity models of the subsurface beneath this region. An eastward gradient in Bouguer anomalies of ~- 1.6 mGal/km is observed within the Taylor Valley. Models reveal thickening of the Moho from 23±5 km beneath the Ross Sea to 35±5 km in the Polar Plateau (dipping at 24.5±7.2°), and lithospheric mantle 100 km thicker in East Antarctica (~200±30 km) than West Antarctica (~90±30 km). </p> <p>Models of predicted surface uplift history are used to estimate an asthenospheric mantle viscosity of 2.1x1020 Pa.s at full surface recovery beneath the Ross Embayment, differing by ~14% from the viscosity at 50% recovery. The temperature contrast between lithospheric and asthenospheric mantle is estimated as ~400°C, equivalent to a viscosity that decreases by a factor of about 30 over the mantle boundary.</p> <p>Results demonstrate that the history of surface uplift in the study area may be complicated, resulting in observations of uplift, or subsidence, at GNSS stations. Future work should incorporate additional geophysical methods, such as seismicity and electrical resistivity, improving constraints on gravity models. A better understanding of the surface uplift (or subsidence) history in the Transantarctic Mountains is critical, with implications in reducing uncertainty in GIA models.</p>

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

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

Dislocation interactions in olivine control postseismic creep of the upper mantle

Changes in stress applied to mantle rocks, such as those imposed by earthquakes, commonly induce a period of transient creep, which is often modelled based on stress transfer among slip systems due to grain interactions. However, recent experiments have demonstrated that the accumulation of stresses among dislocations is the dominant cause of strain hardening in olivine at temperatures ≤600 °C, raising the question of whether the same process contributes to transient creep at higher temperatures. Here, we demonstrate that olivine samples deformed at 25 °C or 1150–1250 °C both preserve stress heterogeneities of ~1 GPa that are imparted by dislocations and have correlation lengths of ~1 μm. The similar stress distributions formed at these different temperatures indicate that accumulation of stresses among dislocations also provides a contribution to transient creep at high temperatures. The results motivate a new generation of models that capture these intragranular processes and may refine predictions of evolving mantle viscosity over the earthquake cycle.