Viscosity Structure Research Articles

Vertical Displacements and Sea‐Level Changes in Eastern North America Driven by Glacial Isostatic Adjustment: An Ensemble Modeling Approach

AbstractGlacial isostatic adjustment (GIA) describes the response of the solid Earth, oceans, and gravitational field to the spatio‐temporal evolution of ice sheets during a glacial cycle. Present‐day vertical displacements and sea‐level changes vary throughout eastern North America in response to the melting of the Laurentide Ice Sheet following the Last Glacial Maximum. We use the open‐source software SELEN4.0 (a SealEveL EquatioN solver) to investigate the influence of GIA on vertical land motions and sea‐level changes in eastern North America. Further, we evaluate the uncertainties associated with the lithospheric thickness and viscosity structure using an ensemble modeling approach (129,956 total simulations). We identify the best‐fitting rheological profiles by comparing modeled vertical displacements to vertical velocity rates derived from Global Positioning System (GPS). We find a general pattern of subsidence (causing accelerated relative sea‐level rise) in the eastern United States region and uplift (causing relative sea‐level fall) in the eastern Canada region consistent with previous studies for two tested ice sheet models (ICE‐6G(VM5a) and ICE‐7G(VM7)). Overall, we find lower rates of modeled vertical displacement using ICE‐6G(VM5a) compared with ICE‐7G(VM7), which produces lower residuals when compared with the GPS‐derived vertical velocity rates. Our ensemble analysis identifies adjustments to the nominal VM5a and VM7 viscosity models that improve fits to the GPS‐imaged vertical velocity rates throughout eastern North America and on the North American Atlantic Coast. The differences in our best‐fitting models for inland versus coastal regions highlight the importance of exploring lateral viscosity variations for GIA modeling throughout North America and elsewhere.

Detectability of low-viscosity zone along lithosphere–asthenosphere boundary beneath the Nankai Trough, Japan, based on high-fidelity viscoelastic simulation

After the 2011 Tohoku-oki earthquake, a seafloor observation system observed rapid landward deformation. These seafloor observation data are potentially explicable via modeling of a low-viscosity zone under the subducting plate, called the lithosphere–asthenosphere boundary (LAB). However, the effect of a low-viscosity LAB has not been empirically examined in past earthquakes because of the absence of seafloor observation data, which are expected to have high sensitivity in the detection of a low-viscosity LAB. The Nankai Trough, southwest Japan, is one location where seafloor geodetic stations are in operation that could detect a low-viscosity LAB after a future earthquake. Therefore, we quantitatively model how the low-viscosity of the LAB is expected to affect postseismic crustal deformation following an anticipated future large earthquake under the Nankai Trough. Calculating viscoelastic deformation using a model that takes into account the realistic crustal structure in the southwestern Japanese Archipelago, we examine the effects of LAB viscosity on surface deformation patterns. The results show for very low LAB viscosities (2.5×1017\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2.5 \ imes 10^{17}$$\\end{document} Pa s), rapid landward displacement as large as 38 cm/yr occurs in the offshore area, with uplift and subsidence patterns highly dependent on the viscosity structure. Especially in the first year after the earthquake, the difference in deformation between the cases with and without a low-viscosity LAB was much larger than the observational error level of the current seafloor observation system. For such low LAB viscosities, the probability of detection is therefore high. On the other hand, for moderately low LAB viscosities (2.5×1018\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2.5 \ imes 10^{18}$$\\end{document} Pa s), the deformation patterns in the cases with and without a low-viscosity LAB are similar. Indeed, the difference in displacement lies within the observational error level of the current observational network, so detection is expected to be difficult for a LAB of only moderately low-viscosity. However, such a LAB may be detected by improving observational accuracy in the Global Navigation Satellite System–acoustic observation technique by performing more frequent measurements. Additional detection potential lies in expected future accuracy improvements in vertical displacement estimates at the DONET stations.Graphical

Lateral and radial viscosity structure beneath Fennoscandia inferred from seismic and magnetotelluric observations

Fennoscandia is continuously uplifting in response to past deglaciation, termed glacial isostatic adjustment or GIA, and its mantle viscosity is well constrained from ice sheet and sea level data. Here, we compare those GIA-constrained viscosities for the Fennoscandian upper mantle with geophysically-constrained viscosities. We construct the upper mantle viscosity structure of Fennoscandia by inferring temperature and water content from seismic and magnetotelluric (MT) data. Using a 1-D MT model for Fennoscandian cratons together with a global seismic model, we infer an upper mantle viscosity (below 250 km) of ∼1021±2 Pa·s, which encompasses the GIA-constrained viscosities of 1020 − 1021 Pa·s. The GIA viscosities are better matched if the Fennoscandian upper mantle is a wet harzburgite or a dry pyrolite, where pyrolite is ∼10 times more viscous than harzburgite. Using the average temperatures and water contents for harzburgitic upper mantle, the GIA viscosities require 1–4 mm grain sizes indicating a diffusion creep regime. In northwestern Fennoscandia, where a high-resolution 2-D resistivity model is available, greater inferred mantle water content implies viscosities that are 10–100 times lower than those for the Fennoscandian Craton. Our work suggests that the combination of seismic and MT observations can improve upper mantle viscosity estimates, especially for regions with laterally-varying viscosity structures or where GIA constraints are not available. Although our method represents an important step forward, viscosity uncertainty can be further reduced by incorporating additional constraints on rock composition, grain size and mantle stress, as well as more accurate geophysical data, into the viscosity calculation.

Spatial migration of temporal earthquake clusters driven by the transfer of differential stress between neighbouring fault/shear-zone structures

Uncertainty concerning the processes responsible for slip-rate fluctuations associated with temporal clustering of surface faulting earthquakes is a fundamental, unresolved issue in tectonics, because strain-rates accommodated by fault/shear-zone structures are the key to understanding the viscosity structure of the crust and seismic hazard. We constrain the timing and amplitude of slip-rate fluctuations that occurred on three active normal faults in central Italy over a time period of 20–30 kyrs, using in situ36Cl cosmogenic dating of fault planes. We identify five periods of rapid slip on individual faults lasting a few millennia, separated time periods of up to 10 millennia with low or zero slip-rate. The rapid slip pulses migrated across the strike between the faults in two waves from SW to NE. We replicate this migration with a model where rapid slip induces changes in differential stress that drive changes in strain-rate on viscous shear zones that drive slip-rate variability on overlying brittle faults. Earthquakes increase the differential stress and strain-rate on underlying shear zones, which in turn accumulate strain, re-loading stress onto the overlying brittle fault. This positive feedback produces high strain-rate episodes containing several large magnitude surface faulting earthquakes (earthquake clusters), but also reduce the differential stress on the viscous portions of neighbouring fault/shear-zones slowing the occurrence of large-magnitude surface faulting earthquakes (earthquake anticlusters). Shear-zones on faults experiencing anticlusters continue to accumulate viscous strain at a lowered rate, and eventually this loads the overlying brittle fault to failure, initiating a period of rapid slip through the positive feedback process described above, and inducing lowered strain-rates onto neighbouring fault/shear-zones. We show that these patterns of differential stress change can replicate the measured earthquake clustering implied by the 36Cl data. The stress changes are related to the fault geometry in terms of distance and azimuth from the slipping structure, implying that (a) strain-rate and viscosity fluctuations for studies of continental rheology, and (b) slip-rates for seismic hazard purposes are to an extent predictable given knowledge of the fault system geometry.

Timescales of glacial isostatic adjustment in Greenland: is transient rheology required?

SUMMARY The possibility of a transient rheological response to ice age loading, first discussed in the literature of the 1980s, has received renewed attention. Transient behaviour across centennial to millennial timescales has been invoked to reconcile apparently contradictory inferences of steady-state (Maxwell) viscosity based on two distinct data sets from Greenland: Holocene sea-level curves and Global Navigation Satellite System (GNSS) derived modern crustal uplift data. To revisit this issue, we first compute depth-dependent Fréchet kernels using 1-D Maxwell viscoelastic Earth models and demonstrate that the mantle resolving power of the two Greenland data sets is highly distinct, reflecting the differing spatial scale of the associated surface loading: the sea-level records are sensitive to viscosity structure across the entire upper mantle while uplift rates associated with post-1000 CE fluctuations of the Greenland Ice Sheet have a dominant sensitivity to shallow asthenosphere viscosity. Guided by these results, we present forward models which demonstrate that a moderate low viscosity zone beneath the lithosphere in Maxwell Earth models provides a simple route to simultaneously reconciling both data sets by significantly increasing predictions of present-day uplift rates in Greenland whilst having negligible impact on predictions of Holocene relative sea-level curves from the region. Our analysis does not rule out the possibility of transient deformation, but it suggests that it is not required to simultaneously explain these two data sets. A definitive demonstration of transient behaviour requires that one account for the resolving power of the data sets in modelling the glacial isostatic adjustment process.

Constraining subducting slab viscosity with topography and gravity fields in free-surface mantle convection models

Subducted slabs provide the primary driving force for mantle convection and plate tectonics. The strength of the slab directly controls the transfer of forces between it and the lithospheric plates at the surface. However, slab viscosity constrained by surface topography and gravity data is significantly lower than that suggested by mineral physics laboratory experiments. It is unclear whether the free-slip top boundary condition used by many studies affects the inverted slab viscosity with gravity or geoid data. In this study, we constrain the subducting slab viscosity structure by comparing the computed topography and gravity in slab subduction models with a free surface with observations. The free-surface model results support relatively weak slabs (20–120 times more viscous than the upper mantle) at the bending region, consistent with previous studies with a free-slip top boundary. The viscosity of the slab below the bending region barely affects the surface topography and gravity field, and both strong and weak slabs fit the observed topography and gravity field, suggesting that extra independent observations are needed to constrain the deep slab viscosity. We investigate the comprehensive relations between subduction interface viscosity, surface topography and gravity anomaly, and trench motion. Models with trench advance have significantly low topography and gravity above the volcanic arc, contradicting subduction zone observations. Together with present trench motion observations and previous studies, we support the idea that the trench retreats under normal single-slab subduction conditions.

GIA imaging of 3-D mantle viscosity based on palaeo sea level observations - Part I: Sensitivity kernels for an Earth with laterally varying viscosity.

A key initial step in geophysical imaging is to devise an effective means of mapping the sensitivity of an observation to the model parameters, that is to compute its Fréchet derivatives or sensitivity kernel. In the absence of any simplifying assumptions and when faced with a large number of free parameters, the adjoint method can be an effective and efficient approach to calculating Fréchet derivatives and requires just two numerical simulations. In the Glacial Isostatic Adjustment problem, these consist of a forward simulation driven by changes in ice mass and an adjoint simulation driven by fictitious loads that are applied at the observation sites. The theoretical basis for this approach has seen considerable development over the last decade. Here, we present the final elements needed to image 3-D mantle viscosity using a dataset of palaeo sea-level observations. Developments include the calculation of viscosity Fréchet derivatives (i.e. sensitivity kernels) for relative sea-level observations, a modification to the numerical implementation of the forward and adjoint problem that permits application to 3-D viscosity structure, and a recalibration of initial sea level that ensures the forward simulation honours present-day topography. In the process of addressing these items, we build intuition concerning how absolute sea-level and relative sea-level observations sense Earth's viscosity structure and the physical processes involved. We discuss examples for potential observations located in the near field (Andenes, Norway), far field (Seychelles), and edge of the forebulge of the Laurentide ice sheet (Barbados). Examination of these kernels: (1) reveals why 1-D estimates of mantle viscosity from far-field relative sea-level observations can be biased; (2) hints at why an appropriate differential relative sea-level observation can provide a better constraint on local mantle viscosity and (3) demonstrates that sea-level observations have non-negligible 3-D sensitivity to deep mantle viscosity structure, which is counter to the intuition gained from 1-D radial viscosity Fréchet derivatives. Finally, we explore the influence of lateral variations in viscosity on relative sea-level observations in the Amundsen Sea Embayment and at Barbados. These predictions are based on a new global 3-D viscosity inference derived from the shear-wave speeds of GLAD-M25 and an inverse calibration scheme that ensures compatibility with certain fundamental geophysical observations. Use of the 3-D viscosity inference leads to: (1) generally greater complexity within the kernel; (2) an increase in sensitivity and presence of shorter length-scale features within lower viscosity regions; (3) a zeroing out of the sensitivity kernel within high-viscosity regions where elastic deformation dominates and (4) shifting of sensitivity at a given depth towards distal regions of weaker viscosity. The tools and intuition built here provide the necessary framework to explore inversions for 3-D mantle viscosity based on palaeo sea-level data.

The Mantle Viscosity Structure of Venus

AbstractThe long‐wavelength gravity and topography of Venus are dominated by mantle convective flows, and are hence sensitive to the planet's viscosity structure and mantle density anomalies. By modeling the dynamic gravity and topography signatures and by making use of a Bayesian inference approach, we investigate the viscosity structure of the Venusian mantle by constraining radial viscosity variations. We performed inversions under a wide range of model assumptions that consistently predicted the existence of a thin low‐viscosity zone in the uppermost mantle. The zone is about 235 km thick and has a viscosity reduction of 5–15 times with respect to the underlying mantle. Drawing a parallel with the Earth, the reduced viscosity could be a result of partial melting as suggested for the origin of the asthenosphere. These results support the interpretation that Venus is a geologically active world predominantly governed by ongoing magmatic processes.

Probabilistic Assessment of Antarctic Thermomechanical Structure: Impacts on Ice Sheet Stability

AbstractUncertainty in present‐day glacial isostatic adjustment (GIA) rates represents at least 44% of the total gravity‐based ice mass balance signal over Antarctica. Meanwhile, physical couplings between solid Earth, sea level and ice dynamics enhance the dependency of the spatiotemporally varying GIA signal on three‐dimensional variations in mantle rheology. Improved knowledge of thermomechanical mantle structure is therefore required to refine estimates of current and projected ice mass balance. Here, we present a Bayesian inverse method for self‐consistently mapping shear‐wave velocities from high‐resolution adjoint tomography into thermomechanical structure using calibrated parameterisations of anelasticity at seismic frequency. We constrain the model using regional geophysical data sets containing information on upper mantle temperature, attenuation and viscosity structure. Our treatment allows formal quantification of parameter covariances, and naturally permits propagation of material parameter uncertainties into thermomechanical structure estimates. We find that uncertainty in steady‐state viscosity structure at 150 km depth can be reduced by 4–5 orders of magnitude compared with a forward‐modeling approach neglecting covariance between viscoelastic parameters. By accounting for the dependence of apparent viscosity on loading timescale, we find good agreement between our estimates of mantle viscosity beneath West Antarctica, and those derived from satellite GPS. Direct access to temperature structure allows us to estimate lateral variations in lithosphere‐asthenosphere boundary (LAB) depth, geothermal heat flow (GHF), and associated uncertainties. We find evidence for shallow LAB depths (63 ± 13 km), and high GHF (76 ± 7 mW m−2) beneath West Antarctica that, combined with low asthenospheric viscosities, indicate a highly dynamic response to ice mass loss.

Influences of across-strike heterogeneous viscosity on the earthquake cycle in a three-dimensional strike-slip fault model

Geodetic observations have shown that there exist large differences in the viscosity of the deep lithosphere across many large strike-slip faults. Heterogeneity in lithospheric viscosity structure can influence the efficiency of stress transfer and thus may have a significant effect on the earthquake cycle. Until now, how the lateral viscosity variation across strike-slip faults affects the earthquake cycles is still not well understood. Here, we investigate the effects of across-strike viscosity variation on long-term earthquake behaviors with a three-dimensional strike-slip fault model. Our model is a quasi-static model which is controlled by the slip-weakening friction law and power-law rheology. By comparing with the reference case, we find that low viscosity on one side of the fault results in a smaller rupture area but with a higher Coulomb stress drop on the ruptured fault region. In addition, low viscosity also leads to a small Coulomb stress accumulation rate. These combined effects increase the earthquake recurrence interval by approximately 10% and the earthquake moments by about 30% when the low viscosity is related to a geothermal gradient of 30 ​K/km. In addition, across-strike viscosity variation causes asymmetric interseismic ground surface deformation rate. As the viscosity contrast increases, the difference in the interseismic ground surface deformation rate between the two sides of the fault gradually increases, although the asymmetric feature is not pronounced. This asymmetry of interseismic ground deformation rate across a strike-slip fault is supposed to result in asymmetric coseismic deformation if the long-term plate motion velocity is invariant. As a result, this kind of asymmetry of interseismic deformation may influence the evaluation of potential earthquake hazards along large strike-slip faults with lateral viscosity contrast.

Inference of the Timescale‐Dependent Apparent Viscosity Structure in the Upper Mantle Beneath Greenland

AbstractContemporary crustal uplift and relative sea level (RSL) change in Greenland is caused by the response of the solid Earth to ongoing and historical ice mass change. Glacial isostatic adjustment (GIA) models, which seek to match patterns of land surface displacement and RSL change, typically employ a linear Maxwell viscoelastic model for the Earth's mantle. In Greenland, however, upper mantle viscosities inferred from ice load changes and other geophysical phenomena occurring over a range of timescales vary by up to two orders of magnitude. Here, we use full‐spectrum rheological models to examine the influence of transient deformation within the Greenland upper mantle, which may account for these differing viscosity estimates. We use observations of shear wave velocity combined with constitutive rheological models to self‐consistently calculate mechanical properties including the apparent upper mantle viscosity and lithosphere thickness across a broad spectrum of frequencies. We find that the contribution of transient behavior is most significant over loading timescales of 102–103 years, which corresponds to the timeframe of ice mass loss over recent centuries. Predicted apparent lithosphere thicknesses are also in good agreement with inferences made across seismic, GIA, and flexural timescales. Our results indicate that full‐spectrum constitutive models that more fully capture broadband mantle relaxation provide a means of reconciling seemingly contradictory estimates of Greenland's upper mantle viscosity and lithosphere thickness made from observations spanning a range of timescales.

Weak upper-mantle base revealed by postseismic deformation of a deep earthquake.

Mantle viscosity plays a key role in the Earth's internal dynamics and thermal history. Geophysical inferences of the viscosity structure, however, have shown large variability depending on the types of observables used or the assumptions imposed1-3. Here, we study the mantle viscosity structure by using the postseismic deformation following a deep (approximately 560 km) earthquake located near the bottom of the upper mantle. We apply independent component analysis4 to geodetic time series to successfully detect and extract the postseismic deformation induced by the moment magnitude 8.2, 2018 Fiji earthquake. To search for the viscosity structure that can explain the detected signal, we perform forward viscoelastic relaxation modelling5,6 with a range of viscosity structures. We find that our observation requires a relatively thin (approximately 100 km), low-viscosity (1017 to 1018 Pa s) layer at the bottom of the mantle transition zone. Such a weak zone could explain the slab flattening7 and orphaning8 observed in numerous subduction zones, which are otherwise challenging to explain in the whole mantle convection regime. The low-viscosity layer may result from superplasticity9 induced by the postspinel transition, weak CaSiO3 perovskite10, high water content11 or dehydration melting12.

Exploring the Origin of Geoid Low and Topography High in West Antarctica: Insights from Density Anomalies and Mantle Convection Models

The deepest geoid low globally with respect to hydrostatic equilibrium is in the Ross Sea area. Nearby in West Antarctica is a residual topography high. Both are in a region with thin lithosphere, where a mantle plume has been suggested. Hence upper mantle viscosity could be regionally reduced, allowing for faster rebound than elsewhere upon melting of the West Antarctic Ice Sheet, one of the global climate system’s tipping elements. To study possible causes of the geoid low / topography high combination, we compute the effects of disk-shaped density anomalies. With -1% density anomaly and a global average radial viscosity structure, geoid low and topography high can be explained with disk radius about 10° and depth range ~150-650 km. Alternatively, there may be two separate disks somewhat laterally displaced, one just below the lithosphere and mainly causing a dynamic topography high and one below the transition zone causing the geoid low. If viscosity in the uppermost mantle is reduced by a factor 10 (from 50 to 350 km depth) to 100 (from 100 to 220 km), one shallow disk in the depth range 50-350 km would also be sufficient. In order to test the feasibility of such density models, we perform computations of a thermal plume that enters at the base of a cartesian box corresponding to a region in the upper mantle, as well as some whole-mantle thermal plume models, with ASPECT. These plume models have typically a narrow conduit and the plume tends to only become wider as it spreads beneath the lithosphere, typically shallower than ~300 km. These results are most consistent with the shallow disk model with reduced uppermost mantle viscosity, hence providing further support for such low viscosities beneath West Antarctica.

Investigating the Effect of Mantle Flow and Viscosity Structure on Surface Velocities in Alaska Using 3‐D Geodynamic Models

AbstractWe utilize 3‐D finite element geodynamic models, incorporating long‐term kinematic estimates of upper plate motion, to better understand the roles that viscosity structure and mantle tractions play in generating plate motions and continental interior deformation in Alaska. Surface deformation in the Pacific‐North American plate boundary zone in Alaska and northwest Canada is strongly influenced by the complex interactions between flat‐slab subduction, gravitational collapse and mantle tractions. Predictions of long‐term tectonic block motion derived from recent Global Positioning System datasets (GPS) reveal surface motions atypical of other continental convergent plate boundary zones. Specifically, in northern and northwestern Alaska, southeastward motion is observed directed back toward the plate boundary. Geodynamic models that incorporate a southeastward directed long‐wavelength mantle traction of ∼2.5–3.8 MPa best replicate surface velocities in Alaska. These mantle tractions, in conjunction with the collision of the Yakutat microplate, appear to drive the uplift and deformation in the Mackenzie Mountains. Furthermore, the extent of the northward motion in southern and central Alaska is controlled by the location of the leading edge of the Yakutat flat slab.

“Ghost Transient” Corrections to the Southern California GPS Velocity Field from San Andreas Fault Seismic Cycle Models

Abstract Seismic cycle models of the northern and southern San Andreas fault (SAF) were developed to estimate corrections for viscoelastic transients embedded in the present-day Global Positioning System (GPS) velocity field. These models incorporate a viscous shear zone embedded in either a layered or a 3D viscoelastic structure. In the southern SAF model, the 3D effective viscosity structure is based on the Southern California Earthquake Center (SCEC) community rheology model (CRM) and the SCEC community thermal model. Viscoelastic materials are represented with Burgers or Maxwell viscoelastic rheologies, and ranges of rheological parameters and earthquake chronologies are explored. For the northern and southern SAF models, surface velocity perturbations or “ghost transients” (GTs) of up to 3.5–6 mm/yr are obtained for a wide range of conditions. In all cases, the GT velocities are consistent with left-lateral rather than right-lateral motion across the SAF, and subtracting them from the present-day GPS velocity field increases the apparent SAF offset rate. For the southern SAF model, GT velocities are fairly insensitive to the timing of large SAF earthquakes prior to 1857, but they are sensitive to rheology (i.e., temperatures and strain rates assumed for computing CRM effective viscosities). For the northern SAF model, GT velocities are sensitive to both rheology (mantle, lower crust and shear zone effective viscosities, and Burgers Body parameters Δb and RB) and mean earthquake recurrence interval. For the reference northern and southern SAF models provided to U.S. Geological Survey National Seismic Hazard Mapping Project deformation modelers, the maximum GT velocities are 3.73 and 4.84 mm/yr, respectively.

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Abstract. Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ∼ 13 of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2 % for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

On the Relation Between Basal Erosion of the Lithosphere and Surface Heat Flux for Continental Plume Tracks

AbstractWhile hotspot tracks beneath thin oceanic lithosphere are visible as volcanic island chains, the plume‐lithosphere interaction for thick continental or cratonic lithosphere often remains hidden due to the lack of volcanism. To identify plume tracks with missing volcanism, we characterize the amplitude and timing of surface heat flux anomalies following a plume‐lithosphere interaction using mantle convection models. Our numerical results confirm an analytical relationship in which surface heat flux increases with the extent of lithosphere thinning, which is primarily controlled by the viscosity structure of the lower lithosphere and the asthenosphere. We find that lithosphere thinning is greatest when the plate is above the plume conduit, while the maximum heat flux anomaly occurs about 40–140 Myr later. Therefore, younger continental and cratonic plume tracks can be identified by observed lithosphere thinning, and older tracks by an increased surface heat flux, even if they lack extrusive magmatism.

Insight into the Viscosity–Structure Relationship of MnO–SiO2–MgO–Al2O3 Fused Submerged Arc Welding Flux

In this study, the influence of Al2O3 on the viscosity and structure of MnO–SiO2–MgO(–Al2O3) submerged arc welding fluxes was investigated. The results showed that the viscosity initially decreased and subsequently increased with increasing Al2O3 contents from 0 to 25 mass pct, which correlates well with the “V”-shaped variation trend of the activation energy. The enhanced (Q2 + Q3 + Al–O0)/(Q0 + Q1 + Al–O−) ratio indicated the continuous polymerization of the flux structure with higher Al2O3 contents, which was further confirmed by the distribution characteristics of various oxygen species. The variation trend may likely be attributed to the competing effects of the overall weakened bond energy and the enhanced degree of polymerization of the flux structure.

Preparation and properties of riboflavin-loaded sanxan microcapsules

Sanxan is a novel microbial polysaccharide approved by the Chinese government for use in the food industry. In this study, we evaluated the usability of sanxan as a drug carrier using riboflavin as a model drug. Sanxan microcapsules were prepared using an extrusion dripping method. The microcapsules were loaded with riboflavin and the release effect was evaluated by in vitro digestion. Microcapsules with a sphericity factor greater than 0.95 and regular appearance could be obtained when the concentration of sanxan was 0.4% and 0.06 M HCl was used as the fixing solution; the hardness of the microcapsules was 633.3 g and their elasticity was 0.44. In the simulated gastric fluid, 3.3% of riboflavin was released, while in the simulated intestinal fluid, this value was 47.4%; this happened when the loading ratio of riboflavin to sanxan in the microcapsules was 7 g g−1. Scanning electron microscopy showed that the viscosity structure of sanxan was an important factor for the slow release of riboflavin. These results revealed that sanxan can be used as a new polysaccharide material for oral delivery.

Influence of 3D Earth Structure on Glacial Isostatic Adjustment in the Russian Arctic

AbstractAnalyses of glacial isostatic adjustment (GIA) and deglacial relative sea‐level (RSL) change in the Russian Arctic deliver important insights into the Earth's viscosity structure and the deglaciation history of the Eurasian ice sheet complex. Here, we validate the 1D GIA models ICE‐6G_C (VM5a) and ICE‐7G_NA (VM7) and select new 3D GIA models in the Russian Arctic against a quality‐controlled deglacial RSL database of >500 sea‐level data points from 24 regions. Both 1D models correspond to the RSL data along the southern coast of the Barents Sea and Franz Josef Land from ∼11 ka BP to present but show notable misfits (>50 m at 10 ka BP) with the White Sea data. We find 3D model predictions of deglacial RSL resolve most of the misfits with the observed data for the White Sea while retaining comparable fits in other regions of the Russian Arctic. Our results further reveal: (a) RSL in the western Russian Arctic is sensitive to elastic lithosphere with lateral thickness variation and 3D viscosity structure in the upper mantle; and (b) RSL in the whole Russian Arctic is less sensitive to 3D viscosity structure in the lower mantle compared to the upper mantle. The 3D models reveal a compromise in the upper mantle between the background viscosity and scaling factor to best fit the RSL data, which needs to be considered in future 3D GIA studies.