Geocenter coordinates estimated from GNSS data as viewed by perturbation theory

Time series of geocenter coordinates were determined with data of two global navigation satellite systems (GNSSs), namely the U.S. GPS (Global Positioning System) and the Russian GLONASS (Global’naya Nawigatsionnaya Sputnikowaya Sistema). The data was recorded in the years 2008–2011 by a global network of 92 permanently observing GPS/GLONASS receivers. Two types of daily solutions were generated independently for each GNSS, one including the estimation of geocenter coordinates and one without these parameters. A fair agreement for GPS and GLONASS was found in the geocenter x- and y-coordinate series. Our tests, however, clearly reveal artifacts in the z-component determined with the GLONASS data. Large periodic excursions in the GLONASS geocenter z-coordinates of about 40 cm peak-to-peak are related to the maximum elevation angles of the Sun above/below the orbital planes of the satellite system and thus have a period of about 4 months (third of a year). A detailed analysis revealed that the artifacts are almost uniquely governed by the differences of the estimates of direct solar radiation pressure (SRP) in the two solution series (with and without geocenter estimation). A simple formula is derived, describing the relation between the geocenter z-coordinate and the corresponding parameter of the SRP. The effect can be explained by first-order perturbation theory of celestial mechanics. The theory also predicts a heavy impact on the GNSS-derived geocenter if once-per-revolution SRP parameters are estimated in the direction of the satellite’s solar panel axis. Specific experiments using GPS observations revealed that this is indeed the case. Although the main focus of this article is on GNSS, the theory developed is applicable to all satellite observing techniques. We applied the theory to satellite laser ranging (SLR) solutions using LAGEOS. It turns out that the correlation between geocenter and SRP parameters is not a critical issue for the SLR solutions. The reasons are threefold: The direct SRP is about a factor of 30–40 smaller for typical geodetic SLR satellites than for GNSS satellites, allowing it in most cases to not solve for SRP parameters (ruling out the correlation between these parameters and the geocenter coordinates); the orbital arc length of 7 days (which is typically used in SLR analysis) contains more than 50 revolutions of the LAGEOS satellites as compared to about two revolutions of GNSS satellites for the daily arcs used in GNSS analysis; the orbit geometry is not as critical for LAGEOS as for GNSS satellites, because the elevation angle of the Sun w.r.t. the orbital plane is usually significantly changing over 7 days.

Time-variable gravity signal in Greenland revealed by high-low satellite-to-satellite tracking

[1] In the event of a termination of the Gravity Recovery and Climate Experiment (GRACE) mission before the launch of GRACE Follow-On (due for launch in 2017), high-low satellite-to-satellite tracking (hl-SST) will be the only dedicated observing system with global coverage available to measure the time-variable gravity field (TVG) on a monthly or even shorter time scale. Until recently, hl-SST TVG observations were of poor quality and hardly improved the performance of Satellite Laser Ranging observations. To date, they have been of only very limited usefulness to geophysical or environmental investigations. In this paper, we apply a thorough reprocessing strategy and a dedicated Kalman filter to Challenging Minisatellite Payload (CHAMP) data to demonstrate that it is possible to derive the very long-wavelength TVG features down to spatial scales of approximately 2000 km at the annual frequency and for multi-year trends. The results are validated against GRACE data and surface height changes from long-term GPS ground stations in Greenland. We find that the quality of the CHAMP solutions is sufficient to derive long-term trends and annual amplitudes of mass change over Greenland. We conclude that hl-SST is a viable source of information for TVG and can serve to some extent to bridge a possible gap between the end-of-life of GRACE and the availability of GRACE Follow-On.

SLR-derived terrestrial reference frame using the observations to LAGEOS-1/2, Starlette, Stella, and AJISAI

Currently, the contributions of Starlette, Stella, and AJISAI are not taken into account when defining the International Terrestrial Reference Frame (ITRF), despite the large amount of data collected in a long time-span. Consequently, the SLR-derived parameters and the SLR part of the ITRF are almost exclusively defined by LAGEOS-1 and LAGEOS-2. We investigate the potential of combining the observations to several SLR satellites with different orbital characteristics. Ten years of SLR data are homogeneously processed using the development version 5.3 of the Bernese GNSS Software. Special emphasis is put on orbit parameterization and the impact of LEO data on the estimation of the geocenter coordinates, Earth rotation parameters, Earth gravity field coefficients, and the station coordinates in one common adjustment procedure. We find that the parameters derived from the multi-satellite solutions are of better quality than those obtained in single satellite solutions or solutions based on the two LAGEOS satellites. A spectral analysis of the SLR network scale w.r.t. SLRF2008 shows that artifacts related to orbit perturbations in the LAGEOS-1/2 solutions, i.e., periods related to the draconitic years of the LAGEOS satellites, are greatly reduced in the combined solutions.

Multi-Decadal Record of Ice Dynamics on Daugaard Jensen Gletscher, East Greenland, from Satellite Imagery and Terrestrial Measurements

The history of ice velocity and calving front position of Daugaard Jensen Gletscher, a large outlet glacier in East Greenland, is reconstructed from field measurements, aerial photography and satellite imagery for the period 1950-2001. The calving terminus of the glacier has remained in approximately the same position over the past similar to 50 years. There is no evidence of a change in ice motion between 1968 and 2001, based on a comparison of velocities derived from terrestrial surveying and feature tracking using sequential satellite images. Estimates of flux near the entrance to the fjord vs snow accumulation in the interior catchment show that Daugaard Jensen Gletscher has a small negative mass balance. This result is consistent with other mass-balance estimates for the inland region of the glacier.

A New Velocity Map for Byrd Glacier, East Antarctica, from Sequential Aster Satellite Imagery

New ice-velocity measurements are obtained for the main trunk of Byrd Glacier, East Antarctica, using recently acquired Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. The velocities are derived from the application of a cross-correlation technique to sequential images acquired in 2000 and 2001. Images were co-registered and ortho-rectified with the aid of a digital elevation model (DEM) generated from ASTER stereo imagery. This paper outlines the process of DEM generation, image co-registration and correction, and the application of the cross-correlation technique to obtain ice velocities. Comparison of the new velocity map with earlier measurements of velocity from 1978 indicates that the glacier has undergone a substantial deceleration between observations. Portions of the glacier flowing at speeds of similar to 850 m a(-1) in 1978/79 were flowing at similar to 650 m a(-1) in 2000/01. The cause of this change in ice dynamics is not known, but the observation shows that East Antarctic outlet glaciers can undergo substantial changes on relatively short timescales.

Rapid Volume Loss from Two East Greenland Outlet Glaciers Quantified Using Repeat Stereo Satellite Imagery

The coastal portions of Kangerdlugssuaq and Helheim glaciers in southeast Greenland lost at least 51 +/- 8 km(-3) yr(-1) of ice between 2001-2006 due to thinning and retreat, according to an analysis of sequential digital elevation models (DEMs) derived from stereo ASTER satellite imagery. The dominant contribution to this ice loss was dynamic thinning caused by the acceleration in flow of both glaciers. Peak rates of change, including thinning rates of similar to 90 m yr(-1), coincided with the rapid increases in flow speed. Extrapolation of the measured data to the ice divides yields an estimated combined catchment volume loss of similar to 122 +/- 30 km(-3) yr(-1), which accounts for half the total mass loss from the ice sheet reported in recent studies. These catchment-wide volume losses contributed similar to 0.31 +/- 0.07 mm yr(-1) to global sea level rise over the 5-year observation period with the coastal regions alone contributing at least 0.1 +/- 0.02 mm yr(-1).

Satellite-Derived Variability in Chlorophyll, Wind Stress, Sea Surface Height, and Temperature in the Northern California Current System

Satellite-derived data provide the temporal means and seasonal and nonseasonal variability of four physical and biological parameters off Oregon and Washington ( 41 degrees - 48.5 degrees N). Eight years of data ( 1998 - 2005) are available for surface chlorophyll concentrations, sea surface temperature ( SST), and sea surface height, while six years of data ( 2000 - 2005) are available for surface wind stress. Strong cross-shelf and alongshore variability is apparent in the temporal mean and seasonal climatology of all four variables. Two latitudinal regions are identified and separated at 44 degrees - 46 degrees N, where the coastal ocean experiences a change in the direction of the mean alongshore wind stress, is influenced by topographic features, and has differing exposure to the Columbia River Plume. All these factors may play a part in defining the distinct regimes in the northern and southern regions. Nonseasonal signals account for similar to 60 - 75% of the dynamical variables. An empirical orthogonal function analysis shows stronger intra-annual variability for alongshore wind, coastal SST, and surface chlorophyll, with stronger interannual variability for surface height. Interannual variability can be caused by distant forcing from equatorial and basin-scale changes in circulation, or by more localized changes in regional winds, all of which can be found in the time series. Correlations are mostly as expected for upwelling systems on intra-annual timescales. Correlations of the interannual timescales are complicated by residual quasi-annual signals created by changes in the timing and strength of the seasonal cycles. Examination of the interannual time series, however, provides a convincing picture of the covariability of chlorophyll, surface temperature, and surface height, with some evidence of regional wind forcing.

Satellite-Detected Fluorescence Reveals Global Physiology of Ocean Phytoplankton

Phytoplankton photosynthesis links global ocean biology and climate-driven fluctuations in the physical environment. These interactions are largely expressed through changes in phytoplankton physiology, but physiological status has proven extremely challenging to characterize globally. Phytoplankton fluorescence does provide a rich source of physiological information long exploited in laboratory and field studies, and is now observed from space. Here we evaluate the physiological underpinnings of global variations in satellite-based phytoplankton chlorophyll fluorescence. The three dominant factors influencing fluorescence distributions are chlorophyll concentration, pigment packaging effects on light absorption, and light-dependent energy-quenching processes. After accounting for these three factors, resultant global distributions of quenching-corrected fluorescence quantum yields reveal a striking consistency with anticipated patterns of iron availability. High fluorescence quantum yields are typically found in low iron waters, while low quantum yields dominate regions where other environmental factors are most limiting to phytoplankton growth. Specific properties of photosynthetic membranes are discussed that provide a mechanistic view linking iron stress to satellite-detected fluorescence. Our results present satellite-based fluorescence as a valuable tool for evaluating nutrient stress predictions in ocean ecosystem models and give the first synoptic observational evidence that iron plays an important role in seasonal phytoplankton dynamics of the Indian Ocean. Satellite fluorescence may also provide a path for monitoring climate-phytoplankton physiology interactions and improving descriptions of phytoplankton light use efficiencies in ocean productivity models.