Mass Redistribution Research Articles

Interannual variations of degree 2 from geodetic observations and surface processes

Summary Geodetic observations from space continuously record surface deformation and global mass redistribution with an increasing accuracy. In parallel, surficial processes (oceanic, atmospheric, and hydrological loading) are more and more precisely modeled.We propose a confrontation of the geodetic Global Positioning System (GPS) and gravity-field satellite laser ranging (SLR) observations at decadal and interannual time scales, in terms of resolution, correlation and comparison with surficial loading models. We focus on the largest global scale signals of degree 2. At interannual periods, surface deformations retrieved from GPS time-series do not exceed 0.8 mm. Our analysis does not reveal the presence of a dominant signal at a specific period, except perhaps for a signal of approximately 3 yr likely connected to the loading response to El Nio / Southern Oscillations. Contrary to the results of previous studies, we do not find in GPS time-series a clear 6-yr oscillation associated with a degree-2 order-2 pattern. Interannual variations in the degree-2 Stokes coefficients of the gravity field do not exceed 2 × 10−11. We do not detect a dominant gravity signal at one specific period but instead a broad spectrum of frequencies. The comparison between the degree 2 deformations built from GPS time-series with a prediction from SLR derived gravity variations reveals some correlations, though their differences remain important. This highlights the present day limitations of these techniques in their ability to characterize global scale interannual variations. Hydrological loading models show some correlations with both GPS and SLR signals, but we cannot firmly establish that continental hydrology is dominantly responsible for the observed variations. Given the current limits in the resolution of both gravity and surface deformation and in the modelling of surface processes, we conclude that it will be a challenge to retrieve a geodetic signal of sub-decadal period originating in the Earth’s core.

Tidal Geopotential Dependence on Earth Ellipticity and Seawater Density and Its Detection With the GRACE Follow‐On Laser Ranging Interferometer

AbstractOcean tides produce significant gravitational perturbations that affect near‐Earth orbiting spacecraft. The gravitational potential induced by tidal mass redistribution is routinely modeled for global gravity analysis and orbit determination, although generally by assuming a spherical Earth and a uniform seawater density. The inadequacy of these simplifications is here addressed. We have developed an accurate yet efficient algorithm to compute the ocean tidal geopotential, allowing for Earth's elliptical shape and variable seawater density. Using this new computation, we find that (1) the effect of ellipticity is several percent of the tide signal over mid to high‐latitude regions, which is comparable to elevation error in the state‐of‐the‐art ocean tide models; (2) the effect of seawater density variations on the potential is as large as 2–3 cm in water‐height equivalent, primarily in deep water where density increases 2%–3% from compressibility. Our analysis of new Gravity Recovery and Climate Experiment Follow‐On (GRACE‐FO) laser ranging interferometer measurements reveals evident errors when ellipticity and density variations are ignored. When accounted for, the GRACE‐FO residual tidal gravity perturbations are reduced by half, depending on the adopted tide model; only the remaining half likely represents actual model elevation error. The use of a spherical surface and a uniform seawater density is no longer tenable given the precision of gravity measurements from GRACE and GRACE‐FO satellites.

Self-consistent determination of the Earth’s GM, geocenter motion and figure axis orientation

The very low-degree Earth’s gravity coefficients, associated with the largest-scale mass redistribution in the Earth’s fluid envelope (atmosphere, oceans and continental hydrology), are the most poorly known. In particular, the first three degree geopotential terms are important, as they relate to intrinsic Earth’s mass references: gravitational coefficient (GM) of the Earth (degree 0), geocenter motion (degree 1), Earth’s figure axis orientation (degree 2). This paper presents a self-consistent determination of these three properties of the Earth. The main objective is to deal with the remaining sources of altimetry satellite orbit uncertainties affecting the fundamental record of sea surface height measurements. The analysis identifies the modeling errors, which should be mitigated when estimating the geocenter coordinates from Satellite Laser Ranging (SLR) observations. The long-term behavior of the degree-0 and -2 spherical harmonics is also observed over the 34-year period 1984–2017 from the long-time history of satellite laser tracking to geodetic spherical satellites. From the analysis of the evolution of these two coefficients, constraints regarding the Earth’s rheology and uncertainties in the value of GM could be inferred. Overall, the influence of the orbit characteristics, SLR station ranging/position biases and satellite signature effects, measurement modeling errors (tropospheric corrections, non-tidal deformations) are also discussed.

Nutation of the non-rigid Earth: Effect of the mass redistribution

In this research, we computed the nutation of the figure axis for a non-rigid Earth model due to the mass redistribution resulting from the lunisolar attraction on the deformable Earth, thus extending our previous work on the precessional motion. The basic Earth model is a two-layer structure composed of a fluid core and an anelastic mantle. We used the Hamiltonian approach, leading to closed-form analytical formulae that describe the nutations in longitude and obliquity of the figure axis as a sum of Poisson and Oppolzer terms. Those formulae were evaluated assuming different Earth rheologies by means of the Love number formalism. In particular, we first computed the effect using the standard model of the International Earth Rotation and Reference Systems Service Conventions (2010) solid tides, and then the Love numbers computed by Williams and Boggs, accounting for the complete oceanic tide contribution, which should provide more consistent and updated values for the nutations. The main amplitudes correspond to the 18.6 yr nutation component and reach 201 μas and −96 μas in the in-phase components in longitude and obliquity, respectively. The obtained values differ greatly from those considered in the current nutation model, IAU2000, of the International Astronomical Union (IAU) – and later similar studies – which includes this effect under the denomination of non-linear terms and derives its numerical contribution on the basis of the Sasao, Okubo, and Saito framework. The differences are significant and reach more than 30μas for some nutation amplitudes. They can be likely attributed to several factors: an incomplete modelling of the redistribution potential; a different treatment of the permanent tide; and the use of different oceanic tide models.

Temporal Variations of the Marine Geoid

AbstractThe effects of temporal changes in the marine geoid on estimates of the ocean dynamic topography are being investigated. Influences from mass redistribution due to changes of land hydrology, ice sheets, glacial isostatic adjustment (GIA), and ocean and atmospheric dynamics are considered, and the associated crustal deformation is included. The strongest signals are associated with the seasonal cycle caused by changes in terrestrial water storage and ice sheets as well as the redistribution of atmospheric mass. Second to this is the importance of an overall trend caused by GIA and decreasing ice sheets over Greenland and Antarctica. On long spatial scales, the amplitude of regional trends estimated for the geoid height has a sizable fraction of those from sea level anomaly (SLA) for the period 1993–2016, also after subtraction of steric height of the upper 1,000 m to analyze trends in deep ocean geostrophic currents. The estimated strong negative geoid height trend south of Greenland for the period 1993–2016 opposes changes in dynamic sea level for the same period thereby affecting past studies on changes of both the strength of the subpolar gyre based on SLA and the meridional overturning circulation on a section between Cape Farewell and Portugal applying ocean dynamic topography and hydrography. We conclude that temporal geoid height trends should be considered in studies of (multi)decadal trends in sea level and circulation on large spatial scales based on altimetry data referenced to a geoid field.

Geometrical effects of microchannel plates: Grazing incidence operation of time-of-flight mass spectrometry and comparison to standard carbon foil.

Direct detection of plasma composition in situ is integral to our understanding of ionospheric, magnetospheric, and heliospheric dynamics. Heavy, singly charged ions have been shown to contribute significantly to the mass and energy redistribution within the magnetosphere. However, most studies neglect the individual contributions of CNO group ions (instead attributing all effects to oxygen) due to the lack of mass separation of the spectrometer. In this paper, we investigate how the mass resolution of a traditional Time-of-Flight (ToF) mass spectrometer responds to different geometric applications of "grazing incidence" scattering. We apply this technique by utilizing standard Microchannel Plates (MCPs) to generate the "Start" signal, instead of the traditional carbon foil, reducing overall scattering via physically limiting plasma collisions to shallow angles of incidence. We considered only geometric differences, comparing LD 20 and LD 40 (length to diameter ratio) MCP effects on mass resolution observed by the ion composition and distribution function analyzer; traditional carbon foil spectra are contrasted against that produced by replacing one foil with these MCPs and subsequently illuminated with H+, H2 +, He+, N+, H2O+, N2 +, Ar+ ions over energies from ∼15 keV to 29 keV, similar to those seen in magnetospheric plasma. We found that among the MCP geometries, the increasing LD ratio corresponds to spectra that observe smaller FWHM and faster mean ToF, indicative of a significant drop in scattering. Additionally, the mass resolution is significantly better than when using the traditional carbon foil for masses above 4 amu. However, the geometry limits detection efficiencies due to the plasma's reduced cross section of collision with the surface that generates the "Start" signal.

A Novel Method for Changing the Dynamics of Slender Elements Using Sponge Particles Structures.

The paper concerns problems related to controlling the dynamic properties of beam-like elements. The parameters of the investigated system can be changed by external factors, resulting in partial changes in the system mass redistribution. It is assumed that it is possible to control the system dynamics by shaping the object frequency structure. The paper introduces the mathematical model of the investigated cantilever beam filled with a Sponge Particle Structure. The continuous model has been simplified to a discrete multi-degree of freedom system. The influence of the system parameters on its behavior is discussed in details. The possible applications of the presented concept are proposed. The spectral vibration analyses were carried out. Theoretical considerations enabled the use of the preliminary semi-active method for controlling the vibration frequencies through a mass redistribution. Experimental studies were carried out to verify the proposed mathematical model.

Preliminary Estimation and Validation of Polar Motion Excitation from Different Types of the GRACE and GRACE Follow-On Missions Data

The Gravity Recovery and Climate Experiment (GRACE) mission has provided global observations of temporal variations in the gravity field resulting from mass redistribution at the surface and within the Earth for the period 2002–2017. Although GRACE satellites are not able to realistically detect the second zonal parameter (ΔC20) of geopotential associated with the flattening of the Earth, they can accurately determine variations in degree-2 order-1 (ΔC21, ΔS21) coefficients that are proportional to variations in polar motion. Therefore, GRACE measurements are commonly exploited to interpret polar motion changes due to variations in the global mass redistribution, especially in the continental hydrosphere and cryosphere. Such impacts are usually examined by computing the so-called hydrological polar motion excitation (HAM) and cryospheric polar motion excitation (CAM), often analyzed together as HAM/CAM. The great success of the GRACE mission and the scientific robustness of its data contributed to the launch of its successor, GRACE Follow-On (GRACE-FO), which began in May 2018 and continues to the present. This study presents the first estimates of HAM/CAM computed from GRACE-FO data provided by three data centers: Center for Space Research (CSR), Jet Propulsion Laboratory (JPL), and GeoForschungsZentrum (GFZ). In this paper, the data series is computed using different types of GRACE/GRACE-FO data: ΔC21, ΔS21 coefficients of geopotential, gridded terrestrial water storage anomalies, and mascon solutions. We compare and evaluate different methods of HAM/CAM estimation and examine the compatibility between CSR, JPL, and GFZ data. We also validate different HAM/CAM estimations using precise geodetic measurements and geophysical models. Analysis of data from the first 19 months of GRACE-FO shows that the consistency between GRACE-FO-based HAM/CAM and observed hydrological/cryospheric signals in polar motion is similar to the consistency obtained for the initial period of the GRACE mission, worse than the consistency received for the best GRACE period, and higher than the consistency obtained for the terminal phase of the GRACE mission. In general, the current quality of HAM/CAM from GRACE Follow-On meets expectations. In the following months, after full calibration of the instruments, this accuracy is expected to increase.

Coastal sea level anomalies and associated trends from Jason satellite altimetry over 2002\u20132018

Climate-related sea level changes in the world coastal zones result from the superposition of the global mean rise due to ocean warming and land ice melt, regional changes caused by non-uniform ocean thermal expansion and salinity changes, and by the solid Earth response to current water mass redistribution and associated gravity change, plus small-scale coastal processes (e.g., shelf currents, wind & waves changes, fresh water input from rivers, etc.). So far, satellite altimetry has provided global gridded sea level time series up to 10–15 km to the coast only, preventing estimation of sea level changes very close to the coast. Here we present a 16-year-long (June 2002 to May 2018), high-resolution (20-Hz), along-track sea level dataset at monthly interval, together with associated sea level trends, at 429 coastal sites in six regions (Northeast Atlantic, Mediterranean Sea, Western Africa, North Indian Ocean, Southeast Asia and Australia). This new coastal sea level product is based on complete reprocessing of raw radar altimetry waveforms from the Jason-1, Jason-2 and Jason-3 missions.

GNSS characterization of hydrological loading in South and Southeast Asia

SUMMARY The elastic response of the lithosphere to surface mass redistributions produces geodetically measurable deformation of the Earth. This deformation is especially pronounced in South and Southeast Asia, where the annual monsoon produces large-amplitude hydrological loads. The Myanmar–India–Bangladesh–Bhutan (MIBB) network of about 20 continuously operating Global Navigation Satellite Systems (GNSS) stations, established in 2011, provides an opportunity to study the Earth's response to these loads. In this study, we use GRACE temporal gravity products as an estimate of long-wavelength surface water distribution and use this estimate in an elastic loading calculation. We compare the predicted vertical deformation from GRACE with that observed with GNSS. We find that elastic loading inferred from the GRACE gravity model is able to explain the phase and much of the peak-to-peak amplitude (typically 2–3 cm) of the vertical GNSS oscillations, especially in northeast India and central Myanmar. GRACE-based corrections reduce the RMS scatter of the GNSS data by 30–45% in these regions. However, this approach does not capture all of the seasonal deformation in central Bangladesh and southern Myanmar. We show by a synthetic test that local hydrological effects may explain discrepancies between the GNSS and GRACE signals in these places. Two independent hydrological loading models of water stored in soil, vegetation, snow, lakes and streams display phase lags compared to the GRACE and GNSS observations, perhaps indicating that groundwater contributes to the observed loading in addition to near-surface hydrology. The results of our calculations have implications for survey-mode GNSS measurements, which make up the majority of geodetic measurements in this region. By using the GNSS data together with estimates of hydrological loading from independent observations and models, we may be able to more accurately determine crustal motions caused by tectonic processes in South and Southeast Asia, while also improving our ability to monitor the annual monsoon and resulting water storage changes in the region.

Metamorphic Differentiation via Enhanced Dissolution along High Permeability Zones

AbstractMetamorphic differentiation, resulting in segregated mineral bands, is commonly recorded in metamorphic rocks. Despite the ubiquitous nature of compositionally layered metamorphic rocks, the processes that are responsible for metamorphic differentiation receive very little attention. Here, detailed petrography, quantitative mineral chemistry and bulk rock analyses are applied to investigate compositional variations and assemblage microstructure. Furthermore, thermodynamic modelling is applied to provide additional constraints on the P–T–XH2O conditions of assemblage formation and mass transfer. The studied outcrop, located within the Bergen arcs of southwestern Norway, preserves the hydration of anorthositic granulite at amphibolite-facies conditions. The amphibolite-facies hydration is expressed as both a statically hydrated amphibolite and a shear zone lithology, defined by the interlayering of amphibolite with leucocratic domains. Within the granulite, quartz-lined fractures surrounded by amphibolite-facies alteration haloes represent relics of initial fluid infiltration associated with brittle failure. The fracture assemblage (quartz + plagioclase + zoisite + kyanite ± muscovite ± biotite) is identical to that occurring within leucocratic domains of the shear zone. Consequently, the compositional layering of the shear zone lithology is linked to fluid infiltration along localized zones of high permeability that result from fracturing. Mass-balance calculations indicate that quartz-lined fractures and compositional differentiation of the shear zone resulted from mass redistribution internal to the shear zone rather than partial melting or precipitation of minerals from externally derived fluid. The process of internal fractionation within the shear zone is driven by enhanced dissolution along highly permeable fracture planes resulting in the loss of MgO, Fetot and K2O from the leucocratic domains. Elements dissolved in the fluid are then transported and ultimately either precipitated in comparatively impermeable amphibolite domains or removed from the system resulting in an overall mass loss. The mass transfer causing metamorphic differentiation of the shear zone is the result of coupled reaction and diffusion under differential stress. The mechanisms of mass redistribution observed within this shear zone provides further insight into the processes that facilitate mass transfer in the Earth’s crust.

Identification and target-pathway deconvolution of FFA4 agonists with anti-diabetic activity from Arnebia euchroma (Royle) Johnst

FFA4 is a novel therapeutic target for the treatment of metabolic diseases, such as type II diabetes. However, there are still few ligands with structural diversity, selectivity and high potency, and the signaling pathway downstream of FFA4 remains to be poorly characterized. In this study, a high performance liquid chromatography-corona charged aerosol detector (HPLC-CAD) combined with label-free dynamic mass redistribution (DMR) method was introduced to guide the discovery of FFA4 agonists from Arnebia euchroma (Royle) Johnst. Ten compounds were identified as FFA4 agonists and structure-activity relationship was obtained. Among them, shikonin displayed the most potent activity with pEC50 value of 6.02 ± 0.19. The activity of shikonin was confirmed by FLIPR (fluorometric imaging plate reader) assay. Signaling pathways of FFA4 were explored in HT-29 cells endogenously expressing FFA4 using shikonin and known FFA4 agonists α-linolenic acid (ALA) and TUG891. Multiple pathways included Gq/11-PLC-Ca2+-PKC, RohA, JNK, p38 MAPK, Gi/o and PI3K signaling but may not involve Gs signaling triggered by shikonin, ALA and TUG891. Besides, shikonin, TUG891 and ALA could induce ERK1/2 and AKT phosphorylation in HT-29 cells. Moreover, anti-diabetes effects of shikonin were evaluated on the glucose intolerance in diabetic db/db mice. Shikonin reduced plasma glucose level, suggesting that it had the potential in treatment of type II diabetes. The agonists identified in this study provided structure guidance for FFA4 drug design. This study was also useful for understanding FFA4 pharmacology and its biological function.

The deep magmatic cumulate roots of the Acadian orogen, eastern North America

Abstract The dearth of cumulate magmatic roots in accretionary orogens is a cornerstone of models that postulate redistribution of mass and energy within the crust for the genesis of intermediate to silicic magmatism. Likewise, the origin of the evolved Acadian (Devonian) plutonism in the New England Appalachians (northeastern USA) has long been explained by closed-system crustal melting due to the absence of associated coeval deep mafic counterparts. Here, we report the discovery of Acadian hydrous ultramafic cumulate rocks that formed by deep-seated (∼1.1 GPa) fractional crystallization processes from a mantle-derived parental melt (Connecticut, southern New England, USA). These rocks are the first of their kind identified in the Appalachian orogen, and one of only a handful of preserved deep subarc hydrous cumulates worldwide. We propose a genetic link between the studied rocks and the evolved coeval plutonism in central-southern New England, where the former represent the missing deep cumulate roots of the same magmatic arc. Our findings support the hypothesis that differentiation of mantle-derived hydrous magmas by fractional crystallization and assimilation processes in the deep crust is a fundamental process in the production of intermediate to silicic magmatism and the geochemical evolution of the continental crust.

Allosteric targeting of the FFA2 receptor (GPR43) restores responsiveness of desensitized human neutrophils.

The G protein‐coupled free fatty acid receptor 2 (FFA2R) is highly expressed on neutrophils and was previously described to regulate neutrophil activation. Allosteric targeting of G protein‐coupled receptors (GPCRs) is increasingly explored to create distinct pharmacology compared to endogenous, orthosteric ligands. The consequence of allosteric versus orthosteric FFA2R activation for neutrophil response, however, is currently largely elusive. Here, different FFA2R desensitization profiles in human neutrophils following allosteric or orthosteric activation are reported. Using a set of neutrophil functional assays to measure calcium flux, pERK1/2, chemotaxis, cellular degranulation, and oxidative burst together with holistic and pathway‐unbiased whole cell sensing based on dynamic mass redistribution, it is found that the synthetic positive allosteric modulator agonist 4‐CMTB potently activates neutrophils and simultaneously alters FFA2R responsiveness toward the endogenous, orthosteric agonist propionic acid (C3) after homologous and heterologous receptor desensitization. Stimulation with C3 or the hierarchically superior chemokine receptor activator IL‐8 led to strong FFA2R desensitization and rendered neutrophils unresponsive toward repeated stimulation with C3. In contrast, stimulation with allosteric 4‐CMTB engaged a distinct composition of signaling pathways as compared to orthosteric receptor activation and was able to activate neutrophils that underwent homologous and heterologous desensitization with C3 and IL‐8, respectively. Moreover, allosteric FFA2R activation could re‐sensitize FFA2 toward the endogenous agonist C3 after homologous and heterologous desensitization. Given the fact that receptor desensitization is critical in neutrophils to sense and adapt to their current environment, these findings are expected to be useful for the discovery of novel pharmacological mechanisms to modulate neutrophil responsiveness therapeutically.

Nanohole morphologies on photoresist surface produced by low-energy Ar+ ion bombardment under normal and near-normal incidence

Random nanohole morphologies on a photoresist surface were produced spontaneously by ion bombardment, from normal to near-normal incidence. At an optimized ion energy, the critical parameters (diameter and depth) and uniformity of these nanoholes can be tailored with multi-ion parameters, including ion fluence (i.e. bombardment time), ion flux, and incidence angle. The temporal evolution of the nanohole structures, at normal incidence, showed growth and stabilization regimes. The mass redistribution effect and the distribution of energy deposition led to the transition from nanoholes to nanoripples, which were deduced from the crater shape evolution at low incidence angles. The incidence angle-dependent evolution of the crater shape (or rim) on the resist surface is consistent with the crater function simulation, and confirms the ion bombardment-induced surface atomic currents and mass redistribution. Finally, X-ray photoelectron spectroscopy characterization confirmed ion bombardment-induced decomposition, in particular, graphitization in the surface layer of the resist.

Center of Pressure in the Paws of Clinically Sound Dogs in Comparison with Orthopedically Diseased Dogs

Simple SummaryThe analysis of the center of pressure is an important tool that allows conclusions to be drawn about the body balance of a patient. This balance can be altered by orthopedic or neurological diseases. To date, there are few data on the center of pressure in the paws of walking dogs. This study aimed to show the compensatory changes in the center of pressure parameters within the paw during stance-phase in the walking dog, with data being collected non-invasively by walking dogs with osteoarthritis and healthy dogs (control group) over a pressure platform. Differences in the center of pressure parameters were found between the affected and non-affected limbs of the diseased dogs, but also in comparison to the corresponding limbs of sound dogs. It was shown that dogs with osteoarthritis show different compensatory changes in the center of pressure parameters depending on whether a front or a hind limb are affected. This can reflect a compensatory redistribution of the body mass as well as compensatory changes of body balance. Based on these observations, a deeper investigation of the center of pressure within the paws of dogs would be of interest and a diagnostic application could become possible in the future.The center of pressure (COP) is recognized as a valuable tool for the assessment of orthopedic and neurologic disorders in humans. Relatively few studies are available in veterinary medicine, particularly concerning the COP in the individual paw. This study assessed the dynamic paw COP parameters during the stance phase of dogs with cox- or cubarthrosis (20 dogs each), as well as of 20 sound dogs. Data were obtained by walking over a pressure platform and analyzed within the diseased groups in comparison to the control group. Both diseased groups showed significant differences between the affected and non-affected limbs, but also in comparison to the reference limbs of sound animals. For coxathrosis, the primary increase was in the medio-lateral COP displacement and the COP area in both hind limbs. For cubarthrosis, the most prominent changes were an increase in the medio-lateral COP displacement in the ipsilateral hind limb and in the cranio-caudal COP displacement in the lame limb. Additionally, the COP area increased in both hind limbs. This can reflect a compensatory redistribution of the body mass as well as compensatory changes of body balance.

Changes of velocity of rotation of the Earth and its figure's deformation associated with them

It is shown, that continuous changes of angular velocity of rotation of the plastically-elastic Earth should cause the continuous coupled deformation of a crustal layer, redistribution of masses in a sub-crustal layer and the change of density associated with it, and also, as consequence of these phenomena, a polar pulsation of a figure when polar diameter of the Earth increases and decreases with time. The mechanism of occurrence of deformations of a body of a planet under action of a deforming (centrifugal) variable force is found, the tensors of deformations and pressure are written out, and on the basis of the rheological equations the equations of balance are deduced, and also calculation of the module of change of polar compression and radial displacements is made at real fluctuations of angular velocity of rotation of the Earth. The calculated values have given the quite real changes of compression and radial displacements of the Earth's crust and its other shells. The opposite process is also shown, namely: observed fluctuations of amplitude of the polar compression, leading to respective alterations of the moment of inertia of the Earth, quite correspond to real fluctuations of duration of a day.

Quantifying sediment mass redistribution from joint time-lapse gravimetry and photogrammetry surveys

Abstract. The accurate quantification of sediment mass redistribution is central to the study of surface processes, yet it remains a challenging task. Here we test a new combination of terrestrial gravity and drone photogrammetry methods to quantify sediment mass redistribution over a 1 km2 area. Gravity and photogrammetry are complementary methods. Indeed, gravity changes are sensitive to mass changes and to their location. Thus, by using photogrammetry data to constrain this location, the sediment mass can be properly estimated from the gravity data. We carried out three joint gravimetry–photogrammetry surveys, once a year in 2015, 2016 and 2017, over a 1 km2 area in southern Taiwan, featuring both a wide meander of the Laonong River and a slow landslide. We first removed the gravity changes from non-sediment effects, such as tides, groundwater, surface displacements and air pressure variations. Then, we inverted the density of the sediment with an attempt to distinguish the density of the landslide from the density of the river sediments. We eventually estimate an average loss of 3.7 ± 0.4 × 109 kg of sediment from 2015 to 2017 mostly due to the slow landslide. Although the gravity devices used in this study are expensive and need week-long surveys, new instrumentation currently being developed will enable dense and continuous measurements at lower cost, making the method that has been developed and tested in this study well-suited for the estimation of erosion, sediment transfer and deposition in landscapes.

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. The transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

Constrained Linear Deconvolution of GRACE Anomalies to Correct Spatial Leakage

Time-varying gravity observed by the Gravity Recovery and Climate Experiment (GRACE) satellites measures surface water and ice mass redistribution driven by weather and climate forcing and has emerged as one of the most important data types in measuring changes in Earth’s climate. However, spatial leakage of GRACE signals, especially in coastal areas, has been a recognized limitation in quantitatively assessing mass change. It is evident that larger terrestrial signals in coastal regions spread into the oceans and vice versa and various remedies have been developed to address this problem. An especially successful one has been Forward Modeling but it requires knowledge of geographical locations of mass change to be fully effective. In this study, we develop a new method to suppress leakage effects using a linear least squares operator applied to GRACE spherical harmonic data. The method is effectively a constrained deconvolution of smoothing inherent in GRACE data. It assumes that oceanic mass changes near the coast are negligible compared to terrestrial changes, with additional spatial regularization constraints. Some calibration of constraint weighting is required. We apply the method to estimate surface mass loads over Australia using both synthetic and real GRACE data. Leakage into the oceans is effectively suppressed and when compared with mascon solutions there is better performance over interior basins.