Large-scale global surface mass variations inferred from GPS
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
We examined the sensitivity of GRACE Follow-On (GRACE-FO) laser ranging interferometer (LRI) measurements to sub-monthly time-variable gravity (TVG) signals caused by transient, high-frequency mass changes in the Earth system. GRACE-FO LRI provides complementary inter-satellite ranging measurements with higher precision over a wider range of frequencies than the baseline K-band microwave ranging system. The common approach for studying mass variation relies on the inverted TVG or mascon solutions over a period of, for example, one month or 10 days which are adversely affected by temporal aliasing and/or smoothing. In this article, we present the alternative along-orbit analysis methodology in terms of line-of-sight gravity difference (LGD) to fully exploit the higher precision LRI measurements for examination of sub-monthly mass changes. The discrepancy between “instantaneous” LGD LRI observations and monthly mean LGD (from Level-2 data) at satellite altitude indicates the sub-monthly gravitational variability not captured by monthly-mean solutions. In conjunction with the satellite ocean altimetry observations, high frequency non-tidal atmosphere and ocean models, and hydrology models, we show that the LGD LRI observations detect the high-frequency oceanic mass variability in the Argentine Basin and the Gulf of Carpentaria, and sub-monthly variations in surface (river) water in the Amazon Basin. We demonstrate the benefits gained from repeat ground track analysis of GRACE-FO LRI data in the case of the Amazon surface water flow. The along-orbit analysis methodology based on LGD LRI time series presented here is especially suitable for quantifying temporal and spatial evolution of extreme, rapidly changing mass variations.
Both surface and internal mass redistribution causes geocenter variations. The global surface mass redistribution can be measured or modeled and is considered as the primary contributor to the geocenter variations on seasonal time scale.
Since the launch of the Gravity Recovery And Climate Experiment (GRACE) in March 2002, the Earth's surface mass variations have been monitored with unprecedented accuracy and resolution. Compared to the classical spherical harmonic solutions, global high-resolution mascon solutions allows the retrieval of mass variations with higher spatial and temporal sampling (2 degrees and 10 days). We present here the validation of the GRACE global mascon solutions by comparing mass estimates to a set of about 100 ocean bottom pressure (OSP) records, and show that the forward modelling of continental hydrology prior to the inversion of the K-band range rate data allows better estimates of ocean mass variations. We also validate our GRACE results to OSP variations modelled by different state-of-the-art ocean general circulation models, including ECCO (Estimating the Circulation and Climate of the Ocean) and operational and reanalysis from the MERCATOR project.
Explore the source record for details and available documents.
Estimates of the overall mass balance and seasonal and inter-annual variations in the surface mass balance are obtainable from time-series of ice surface elevations measured by satellite altimetry. Beginning in 2001, NASA's ICESat laser altimeter and lidar mission will significantly improve the range accuracy, the orbit accuracy, and the spatial coverage for measurement of ice sheet elevations (to 86 S) , as compared to previous radar altimeters designed for ocean measurements The radar altimeters on Seasat and Geosat provided ice sheet measurements to 72 S, and on ERS-1 and ERS-2 to 81 S. Although radar altimetry has significant limitations in coverage (due to loss of tracking) and accuracy over sloping surfaces, information on ice-sheet surface-elevation changes has been derived for parts of Antarctica. Recently, the accuracy of the ice measurements by Seasat (3 months of 1978) and Geosat (1985 to 1989) have been improved by new calculations of the satellite orbit heights and other altimeter corrections. Residual orbit errors and inter-satellite biases are evaluated by crossover analysis and by global adjustments to an ocean surface derived from altimeter data. The standard deviation of the orbit error is less than 9 cm, and the long-term trend in the error appears to be less than 1 cm/yr. Orbit errors can be further reduced by adjustment to the ocean surface, but false signals of several cm/yr may be also introduced by the adjustments. These false signals are caused mainly by residual errors in the altimeter corrections over the ocean, and secondary by real changes in the ocean surface elevation. Maps of ice sheet elevation changes north of 72 S are derived from Seasat-Geosat crossovers and from 4.5 years of Geosat crossovers. A notable ice thinning rate of about 50 cm/yr is found at elevations below 2200 meters between 70 and 72 S to the East of the Amery ice shelf, in both the Seasat-Geosat and Geosat-Geosat time intervals Above 2200 meters, to the ridge around 95 E, the thinning is about 10 to 20 cm/yr. For the Seasat/Geosat interval, a thickening is noted on the area East of the ridge between about 95 E and 115 E. Preliminary results from analysis of ERS-1 and ERS-2 radar altimeter data are also described.
We present optical and IR surface photometry of M51 (NGC 5194) at B, V, R, I, J, K, and CO. These data are used to establish whether K-band images of spiral galaxies provide reliable maps of stellar surface mass density features such as massive spiral arms or bars. The main distorting agents in the mapping at shorter wavelengths are dust extinction and luminous young stars. From modeling the color changes across the main dust lanes in M51, we find the optical depths to be about 0.5 in the K band. For these optical depths the K-band flux is attenuated by only less than 10 percent even in the dust lanes. From monitoring the gravity-sensitive CO index across the spiral arms we find that young red supergiants do not significantly distort the K-band image. On this basis, we conclude that K-band images do trace the massive disk star population and allow to map azimuthal variation in the surface mass density of the stellar disk. In M51 we find the surface mass density contrast to range from 1.8 to 3, as expected from N-body simulations of the galaxy's tidal encounter with NGC 5195. The spiral arm amplitude in M51 clearly shows smooth, strong radial variations, with a maximum at about 130 arcsec and minima at 45 and 170 arcsec.
Accurate estimates of surface mass balance over the Greenland ice sheet (GrIS) would contribute to understanding the cause of recent changes and would help to better estimate the future contribution of the GrIS to sea-level rise. Given the limitations of in-situ measurement, modeling, and remote sensing, it is critical to explore the opportunity to merge the available data to better characterize the spatial and temporal variation of the GrIS surface mass balance (SMB). This work utilizes a particle batch smoother data assimilation technique that yields SMB estimates that benefit from the snow model Crocus and a 16-day albedo product derived from satellite remote sensing data. Comparison of the results against in-situ SMB measurements shows that the assimilation of the albedo product reduces the root mean square error (RMSE) of the posterior estimates of SMB by 51% and reduces bias by 95%.
Although geoid or surface gravity anomalies cannot be uniquely related to an interior distribution of mass, they can be related to a surface mass distribution. However, over horizontal distances greater than about 100 km, the condition of isostatic equilibrium above the asthenosphere is a good approximation and the total mass per unit column is zero. Thus the surface distribution of mass is also zero. For this case we show that the surface gravitational potential anomaly can be uniquely related to a surface dipole distribution of mass. Variations in the thickness of the crust and lithosphere can be expected to produce undulations in the geoid.
The NASA/German Gravity Recovery and Climate Experiment (GRACE) was launched in March 2002. Rather than looking downward, GRACE continuously monitors the locations of and precise distance between twin satellites which orbit in tandem about 200 km apart. Variations in mass near Earth's surface cause heterogeneities in its gravity field, which in turn affect the orbits of satellites. Thus scientists can use GRACE data to map Earth's gravity field with enough accuracy to discern month to month changes caused by ocean circulation and redistribution of water stored on and in the land. Other gravitational influences, such as atmospheric circulation, post-glacial rebound, and solid earth movements are either independently determined and removed or are negligible on a monthly to sub-decadal timescale. Despite its coarse spatial (>150,000 sq km at mid-latitudes) and temporal (approx monthly) resolutions, GRACE has enabled significant advancements in the oceanic, hydrologic, and cryospheric science, and has great potential for climate monitoring, because it is the only global observing system able to measure ocean bottom pressures, total terrestrial water storage, and ice mass changes. The best known GRACE results are estimates of Greenland and Antarctic ice sheet loss rates. Previously, scientists had estimated ice mass losses using ground and satellite based altimetry and surface mass balance estimates based on snowfall accumulation and glacier discharge. While such measurements are still very useful for their spatial detail, they are imperfectly correlated with large-scale ice mass changes, due to snow and ice compaction and incomplete spatial coverage. GRACE enables scientists to generate monthly time series of Greenland and Antarctic ice mass, which have confirmed the shrinking of the polar ice sheets, one of the most obvious and indisputable manifestations of climate change. Further, GRACE has located and quantified hot spots of ice loss in southeastern Greenland and western Antarctica. For 2002 to present, the rate of ice mass loss has been 200 to 300 GT/yr in Greenland and 70 to 210 GT/yr in Antarctica, and some scientists are suggesting that the rates are accelerating. Similarly, GRACE has been used to monitor mass changes in alpine glaciers. Tamisiea et al. first characterized glacier melt along the southern coast of Alaska, more recently estimated to be occurring at a rate of 84 GT/yr. Chen et al. estimated that Patagonian glaciers are melting at a rate of 28 GT/yr, and estimated that the high mountains of central Asia lose ice at a rate of 47 GT/yr. Tapley et al. and Wahr et al. presented the first GRACE based estimates of changes in column-integrated terrestrial water storage (TWS; the sum of ground-water, soil moisture, surface waters, snow, ice, and water stored in vegetation) at continental scales. Since then, dozens of studies have shown that GRACE based estimates of regional to continental scale TWS variations agree with independent information, and some innovative uses of GRACE data have been developed. Rodell et al. (2004) and Swenson and Wahr (2006) demonstrated that by combining GRACE derived terrestrial water storage changes with observations of precipitation and runoff in a river basin scale water budget, it was possible to produce new estimates of evapotranspiration and atmospheric moisture convergence, essential climate variables that are difficult to estimate accurately. Similarly, GRACE has been used to constrain estimates of global river discharge and the contribution of changes in TWS to sea level rise. Crowley et al. observed a negative correlation between interannual TWS anomalies in the Amazon and the Congo River basin. Yeh et al. and Rodell et al. estimated regionally averaged groundwater storage variations based on GRACE and auxiliary observations. Rodell et al. and Tiwari et al. applied that method to quantify massive groundwater depletion in northern India caused by over reliance on aquifers for irration, and Famiglietti et al. found a similar situation in California's Central Valley. Zaitchik et al. and Lo et al. described approaches to use GRACE to constrain hydrological models, enabling integration of GRACE data with other observations and achieving much higher spatial and temporal resolutions than GRACE alone. Such approaches are now supporting applications including drought and water resources monitoring. Oceanography has likewise benefitted from the independent nature of GRACE observations. One application is measurement of the mass component of sea level rise, which complements radar altimetry and in situ measurements. GRACE also measures ocean bottom pressures (OBP), which help to refine understanding and modeling of ocean circulation and the ocean's fresh water budget, among other things. For example, Hayakawa et al. showed that GRACE observes OBP patterns absent from the background models of oceanic variability. Morison et al. used GRACE to describe important decadal scale shifts in circulation and an ongoing trend of freshening of the western Arctic, important indicators of climate variability. The research of Song and Zlotnicki and Chambers and Willis on GRACE-derived ocean bottom pressures in the sub-polar gyre led to the discovery of an ENSO teleconnection and a long-term change in OBP in the North Pacific sub-polar gyre that was not predicted by an ocean model. Further, Chambers and Willis were able to identify an internal redistribution of mass between Atlantic and Pacific Oceans lasting at least six years, which was not predicted by ocean models and was the first direct evidence of sustained mass transport from one ocean basin to another on periods longer than a year. Boening et al. observed a record increase in OBP over part of the southeastern Pacific in late 2009 and early 2010, primarily caused by wind stress curl associated with a strong and persistent anticyclone and likely related to the concurrent Central Pacific El Nino. GRACE has far surpassed its 5-year design lifetime, but it will likely succumb to the aging of batteries and instrument systems sometime in the next few years. NASA has begun initial development of a follow-on to GRACE with very similar design, which could launch as soon as 2016 and would provide continuity in the data record while improving resolution slightly. Higher resolution time variable gravity missions are also on the drawing board .
Lunar surface formation, emphasizing density variations and mass concentrations
A model and simulation code have been developed for time-dependent axisymmetric disk accretion onto a compact object including for the first time the influence of an ordered magnetic field. The accretion rate and radiative luminosity of the disk are naturally coupled to the rate of outflow of energy and angular momentum in magnetically driven (+/- z) winds. The magnetic field of the wind is treated in a phenomenological way suggested by self-consistent wind solutions. The radial accretion speed u(r, t) of the disk matter is shown to be the sum of the usual viscous contribution and a magnetic contribution proportional to r(exp 3/2)(B(sub p exp 2))/sigma, where B(sub p)(r,t) is the poloidal field threading the disk and sigma(r,t) is the disk's surface mass density. An enhancement or variation in B(sub p) at a large radial distance leads to the formation of a soliton-like structure in the disk density, temperature, and B-field which propagates implosively inward. The implosion gives a burst in the power output in winds or jets and a simultaneous burst in the disk radiation. The model is pertinent to the formation of discrete fast-moving components in jets observed by very long baseline interferometry. These components appear to originate at times of optical outbursts of the active galactic nucleus.
Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.
The atmospheric angular momentum is closely related to variations in the Earth rotation. The atmospheric excitation function (AEF), or namely atmospheric effective angular momentum function, is introduced in studying the atmospheric excitation of the Earth's variable rotation. It may be separated into two portions, i.e, the "wind" terms due to the atmospheric motion relative to the mantle and the "pressure" terms due to the variations of atmospheric mass distribution evident through surface pressure changes. The AEF wind terms during the period of 1948-2004 are re-processed from the NCEP/NCAR (National Centers for Environmental Prediction-National Center for Atmospheric Research) reanalysis 6-hourly wind and pressure fields. Some previous calculations were approximate, in that the wind terms were integrated from an isobaric lower boundary of 1000 hPa. To consider the surface topography effect, however, the AEF is computed by integration using the winds from the Earth's surface to 10 hPa, the top atmospheric model level, instead of from 1000 hPa. For these two cases, only a minor difference, equivalent to approx. 0.004 milliseconds in length-of-day variation, exists with respect to the axial wind term. However, considerable differences, equivalent to 5-6 milliarcseconds in polar motion, are found regarding equatorial wind terms. We further compare the total equatorial AEF (with and without the topographic effect) with the polar motion excitation function (PMEF) during the period of 1980-2003. The equatorial AEF gets generally closer to the PMEF, and improved coherences are found between them when the topography effect is included. Keywords: Atmospheric angular momentum, Atmospheric excitation function, Earth rotation, Topography, Wind, Pressure.
The force of Earth's gravity field varies in proportion to the amount of mass near the surface. Spatial and temporal variations in the gravity field can be measured via their effects on the orbits of satellites. The Gravity Recovery and Climate Experiment (GRACE) is the first satellite mission dedicated to monitoring temporal variations in the gravity field. The monthly gravity anomaly maps that have been delivered by GRACE since 2002 are being used to infer changes in terrestrial water storage (the sum of groundwater, soil moisture, surface waters, and snow and ice), which are the primary source of gravity variability on monthly to decadal timescales after atmospheric and oceanic circulation effects have been removed. Other remote sensing techniques are unable to detect water below the first few centimeters of the land surface. Conventional ground based techniques can be used to monitor terrestrial water storage, but groundwater, soil moisture, and snow observation networks are sparse in most of the world, and the countries that do collect such data rarely are willing to share them. Thus GRACE is unique in its ability to provide global data on variations in the availability of fresh water, which is both vital to life on land and vulnerable to climate variability and mismanagement. This chapter describes the unique and challenging aspects of GRACE terrestrial water storage data, examples of how the data have been used for research and applications related to fresh water vulnerability and change, and prospects for continued contributions of satellite gravimetry to water resources science and policy.
Redistribution of mass near Earth's surface alters its rotation, gravity field, and geocenter location. Advanced techniques for measuring these geodetic variations now exist, but the ability to attribute the observed modes to individual Earth system processes has been hampered by a shortage of reliable global data on such processes, especially hydrospheric processes. To address one aspect of this deficiency, 17 yrs of monthly, global maps of vegetation biomass were produced by applying field-based relationships to satellite-derived vegetation type and leaf area index. The seasonal variability of biomass was estimated to be as large as 5 kg m(exp -2). Of this amount, approximately 4 kg m(exp -2) is due to vegetation water storage variations. The time series of maps was used to compute geodetic anomalies, which were then compared with existing geodetic observations as well as the estimated measurement sensitivity of the Gravity Recovery and Climate Experiment (GRACE). For gravity, the seasonal amplitude of biomass variations may be just within GRACE'S limits of detectability, but it is still an order of magnitude smaller than current observation uncertainty using the satellite-laser-ranging technique. The contribution of total biomass variations to seasonal polar motion amplitude is detectable in today's measurement, but it is obscured by contributions from various other sources, some of which are two orders of magnitude larger. The influence on the length of day is below current limits of detectability. Although the nonseasonal geodynamic signals show clear interannual variability, they are too small to be detected.