Calibrating the GOCE accelerations with star sensor data and a global gravity field model

A reliable and accurate gradiometer calibration is essential for the scientific return of the gravity field and steady-state ocean circulation explorer (GOCE) mission. This paper describes a new method for external calibration of the GOCE gradiometer accelerations. A global gravity field model in combination with star sensor quaternions is used to compute reference differential accelerations, which may be used to estimate various combinations of gradiometer scale factors, internal gradiometer misalignments and misalignments between star sensor and gradiometer. In many aspects, the new method is complementary to the GOCE in-flight calibration. In contrast to the in-flight calibration, which requires a satellite-shaking phase, the new method uses data from the nominal measurement phases. The results of a simulation study show that gradiometer scale factors can be estimated on a weekly basis with accuracies better than 2 × 10⁻³ for the ultrasensitive and 10⁻² for the less sensitive axes, which is compatible with the requirements of the gravity gradient error. Based on a 58-day data set, scale factors are found that can reduce the errors of the in-flight-calibrated measurements. The elements of the complete inverse calibration matrix, representing both the internal gradiometer misalignments and scale factors, can be estimated with accuracies in general better than 10⁻³.     

The ITG-Goce02 gravity field model from GOCE orbit and gradiometer data based on the short arc approach

In this contribution, we describe the global GOCE-only gravity field model ITG-Goce02 derived from 7.5 months of gradiometer and orbit data. This model represents an alternative to the official ESA products as it is computed completely independently, using a different processing strategy and a separate software package. Our model is derived using the short arc approach, which allows a very effective decorrelation of the highly correlated GOCE gradiometer and orbit data noise by introducing a full empirical covariance matrix for each arc, and gives the possibility to downweight ‘bad’ arcs. For the processing of the orbit data we rely on the integral equation approach instead of the energy integral method, which has been applied in several other GOCE models. An evaluation against high-resolution global gravity field models shows very similar differences of our model compared to the official GOCE results published by ESA (release 2), especially to the model derived by the time-wise approach. This conclusion is confirmed by comparison of the GOCE models to GPS/levelling and altimetry data.     

GOCE gradiometer: estimation of biases and scale factors of all six individual accelerometers by precise orbit determination

A method has been implemented and tested for estimating bias and scale factor parameters for all six individual accelerometers that will fly on-board of GOCE and together form the so-called gradiometer. The method is based on inclusion of the individual accelerometer observations in precise orbit determinations, opposed to the baseline method where so-called common-mode accelerometer observations are used. The method was tested using simulated data from a detailed GOCE system simulator. It was found that the observations taken by individual accelerometers need to be corrected for (1) local satellite gravity gradient (SGG), and (2) rotational terms caused by centrifugal and angular accelerations, due to the fact that they are not located in the satellite’s center of mass. For these corrections, use is made of a reference gravity field model. In addition, the rotational terms are derived from on-board star tracker observations. With a perfect a priori gravity field model and with the estimation of not only accelerometer biases but also accelerometer drifts, scale factors can be determined with an accuracy and stability better than 0.01 for two of the three axes of each accelerometer, the exception being the axis pointing along the long axis of the satellite (more or less coinciding with the flight direction) for which the scale factor estimates are unreliable. This axis coincides with the axis of drag-free control, which results in a small variance of the signal to be calibrated and thus an inaccurate determination of its scale factor in the presence of relatively large (colored) accelerometer observation errors. In the presence of gravity field model errors, it was found that still an accuracy and stability of about 0.015 can be obtained for the accelerometer scale factors by simultaneously estimating empirical accelerations.     

Optimum data weighting and error calibration for estimation of gravitational parameters

Summary A new approach has been developed for determining consistent satellite-tracking data weights in solutions for the satellite-only gravitational models. The method employs subset least-squares solutions of the satellite data contained within the complete solution and requires that the differences of the parameters of subset solutions and the complete solution to be in agreement with their error estimates by adjusting the data weights. GEM-T2 model was recently computed and adjusted through a direct application of this method. The estimated data weights are markedly smaller than the weights implied by the formal uncertainties of the measurements. Orbital arc tests as well as surface gravity comparisons show significant improvements for solutions when more realistic data weighting is achieved.     

Simulation of regional gravity field recovery from satellite gravity gradiometer data using collocation and FFT

When planning a satellite gravity gradiometer (SGG) mission, it is important to know the quality of the quantities to be recovered at ground level as a function of e.g. satellite altitude, data type and sampling rate, and signal variance and noise. This kind of knowledge may be provided either using the formal error estimates of wanted quantities using least-squares collocation (LSC) or by comparing simulated data at ground level with results computed by methods like LSC or Fast Fourier Transform (FFT). Results of a regional gravity field recovery in a 10^o×20^o area surrounding the Alps using LSC and FFT are reported. Data used as observations in satellite altitude (202 or161 km) and for comparison at ground level were generated using theOSU86F coefficient set, complete to degree 360. These observations are referred to points across simulated orbits. The simulated quantities were computed for a 45 days mission period and 4 s sampling. A covariance function which also included terms above degree 360 was used for prediction and error estimation. This had the effect that the formal error standard deviation for gravity anomalies were considerably larger than the standard deviations of predicted minus simulated quantities. This shows the importance of using data with frequency content above degree 360 in simulation studies. Using data at202 km altitude the standard deviation of the predicted minus simulated data was equal to8.3 mgal for gravity and0.33 m for geoid heights.     

Evaluation of the first GOCE static gravity field models using terrestrial gravity,vertical deflections and EGM2008 quasigeoid heights

Recently, four global geopotential models (GGMs) were computed and released based on the first 2 months of data collected by the Gravity field and steady-state Ocean Circulation Explorer (GOCE) dedicated satellite gravity field mission. Given that GOCE is a technologically complex mission and different processing strategies were applied to real space-collected GOCE data for the first time, evaluation of the new models is an important aspect. As a first assessment strategy, we use terrestrial gravity data over Switzerland and Australia and astrogeodetic vertical deflections over Europe and Australia as ground-truth data sets for GOCE model evaluation. We apply a spectral enhancement method (SEM) to the truncated GOCE GGMs to make their spectral content more comparable with the terrestrial data. The SEM utilises the high-degree bands of EGM2008 and residual terrain model data as a data source to widely bridge the spectral gap between the satellite and terrestrial data. Analysis of root mean square (RMS) errors is carried out as a function of (i) the GOCE GGM expansion degree and (ii) the four different GOCE GGMs. The RMS curves are also compared against those from EGM2008 and GRACE-based GGMs. As a second assessment strategy, we compare global grids of GOCE GGM and EGM2008 quasigeoid heights. In connection with EGM2008 error estimates, this allows location of regions where GOCE is likely to deliver improved knowledge on the Earth’s gravity field. Our ground truth data sets, together with the EGM2008 quasigeoid comparisons, signal clear improvements in the spectral band ~160–165 to ~180–185 in terms of spherical harmonic degrees for the GOCE-based GGMs, fairly independently of the individual GOCE model used. The results from both assessments together provide strong evidence that the first 2 months of GOCE observations improve the knowledge of the Earth’s static gravity field at spatial scales between ~125 and ~110 km, particularly over parts of Asia, Africa, South America and Antarctica, in comparison with the pre-GOCE-era.     

Ionospheric calibration of single frequency VLBI and GPS observations using dual GPS data

Summary The ionospheric effect is one of the main sources of error in Very Long Baseline Interferometry (VLBI) and Global Positioning System (GPS) high precision geodesy. Although the use of two frequencies allows the estimation of this effect, in some cases dual observations are not possible due to the available equipment or the type of observation. This paper presents the ionospheric calibration of single frequency VLBI and GPS observations based on the ionospheric electron content estimated from dual frequency GPS data. The ionospheric delays obtained with this procedure and the VLBI baseline length results have been compared with those obtained with dual frequency data. For the European geodetic VLBI baselines, both solutions agree at the 3–5 parts in 10^–9 level. The noise introduced by the GPS-based calibration is in the order of 3 cm for the VLBI observables and of 10 cm for the GPS observables.     

Estimating crop net primary production using national inventory data and MODIS-derived parameters

National estimates of spatially-resolved cropland net primary production (NPP) are needed for diagnostic and prognostic modeling of carbon sources, sinks, and net carbon flux between land and atmosphere. Cropland NPP estimates that correspond with existing cropland cover maps are needed to drive biogeochemical models at the local scale as well as national and continental scales. Existing satellite-based NPP products tend to underestimate NPP on croplands. An Agricultural Inventory-based Light Use Efficiency (AgI-LUE) framework was developed to estimate individual crop biophysical parameters for use in estimating crop-specific NPP over large multi-state regions. The method is documented here and evaluated for corn (Zea mays L.) and soybean (Glycine max L. Merr.) in Iowa and Illinois in 2006 and 2007. The method includes a crop-specific Enhanced Vegetation Index (EVI), shortwave radiation data estimated using the Mountain Climate Simulator (MTCLIM) algorithm, and crop-specific LUE per county. The combined aforementioned variables were used to generate spatially-resolved, crop-specific NPP that corresponds to the Cropland Data Layer (CDL) land cover product. Results from the modeling framework captured the spatial NPP gradient across croplands of Iowa and Illinois, and also represented the difference in NPP between years 2006 and 2007. Average corn and soybean NPP from AgI-LUE was 917 g C m⁻² yr⁻¹ and 409 g C m⁻² yr⁻¹, respectively. This was 2.4 and 1.1 times higher, respectively, for corn and soybean compared to the MOD17A3 NPP product. Site comparisons with flux tower data show AgI-LUE NPP in close agreement with tower-derived NPP, lower than inventory-based NPP, and higher than MOD17A3 NPP. The combination of new inputs and improved datasets enabled the development of spatially explicit and reliable NPP estimates for individual crops over large regional extents.     

基于天宫一号高光谱数据与Landsat TM数据的湿地利用变化监测

本文利用2013年4月18日获取的天宫一号高光谱成像仪数据和1985年5月1日获取的Landsat TM数据,对黄河入海口地区湿地的土地利用状况进行了监测研究.通过对比分类结果发现,该地区近30年来人类活动对湿地的影响强烈,滩涂开垦、水产养殖业占用自然水面现象严重.同时由于天宫一号高光谱成像仪数据具有更高的光谱和空间分辨率,其分类精度优于Landsat TM,更适合于湿地的监测.     

Harmonic analysis of total electron contents time series: methodology and results

In an attempt to model regular variations of the ionosphere, the least-squares harmonic estimation is applied to the time series of the total electron contents (TEC) provided by the JPL analysis center. Multivariate and modulated harmonic estimation spectra are introduced and estimated for the series to detect the regular and modulated dominant frequencies of the periodic patterns. Two significant periodic patterns are the diurnal and annual signals with periods of 24/n hours and 365.25/n days (n = 1, 2, …), which are the Fourier series decomposition of the regular daily and yearly periodic variations of the ionosphere. The spectrum shows a cluster of periods near 27 days, thereby indicating irregularities at this solar cycle period. A series of peaks, with periods close to the diurnal signal and its harmonics, are evident in the spectrum. In fact, the daily signal harmonics of ω _i  = 2πi are modulated with the annual signal harmonics of ω _j  = 2πj/365.25 as ω _ijM  = 2πi(1 ± j/365.25i). Among them, at low and midlatitudes, the largest variations belong to the diurnal signal modulated to the semiannual signal. Some preliminary results on the modulated part are presented. The maximum ranges of the modulated daily signal are ±15 TECU and ±6 TECU at high and low solar periods, respectively. A model consisting of purely harmonic functions plus modulated ones is capable of studying known regular anomalies of the ionosphere, which is currently in progress.     

Derivation of crop phenological parameters using multi-date SPOT-VGT-NDVI data: A case study for Punjab

Crop phenological parameters, such as the start and end time of the crop growth, the total length of the growing season, time of peak vegetation and rate of greening and senescence are important for planning crop management and crop diversification/intensification. Multi-temporal remote sensing data provides opportunity to characterize the crop phenology at regional level. This study was conducted during the kharif season of the year 2001–02 for Punjab. The ten-day Normalised Difference Vegetation Index (NDVI) composite products, with 1 km spatial resolution, available from the Vegetation sensor onboard SPOT4 were used for the study. Twenty-one temporal datasets from May 1, 2001 to November 21, 2001 were used. Logical modelling approach was followed to compute the minimum and maximum NDVI, the amplitude of NDVI, the threshold NDVI during sowing and harvest, the crop duration, integrated NDVI and skewness of profile. The analysis showed that before July beginning, in the whole of Punjab, sowing/planting was over. It was found that the crop emergence in the eastern part of the state started earlier than the western region. The maximum NDVI, which represented peak vegetative stage, was above 0.7 and occurred mostly during August. The duration of crops ranged between 90–140 days, with majority between 110–120 days. Total integrated NDVI in Punjab was generally above 60. Using principal component analysis and divergence analysis seven best metrics were selected for crop discrimination.     

Keeping it simple: Monitoring flood extent in large data-poor wetlands using MODIS SWIR data

Characterising inundation conditions for flood-pulsed wetlands is a critical first step towards assessment of flood risk as well as towards understanding hydrological dynamics that underlay their ecology and functioning. In this paper, we develop a series of inundation maps for the Okavango Delta, Botswana, based on the thresholding of the SWIR band (b7) MODIS MCD43A4 product. We show that in the Okavango Delta, SWIR is superior to other spectral bands or derived indices, and illustrate an innovative way of defining the spectral threshold used to separate inundated from dry land. The threshold is determined dynamically for each scene based on reflectances of training areas capturing end-members of the inundation spectrum. The method provides a very good accuracy and is suitable for automated processing.     

Absolute calibration of GPS antennas: laboratory results and comparison with field and robot techniques

A critical assessment of the accuracy of GPS antenna calibration is most effectively done by comparison between different calibration methods. We present new chamber calibrations of five different GPS receiver antenna types in an anechoic chamber and a comparison of an individual antenna calibrated by the absolute field calibration technique with robot mount of IfE/GEO++. The accuracy is described using standard error parameters which allow the characterization of the quality of different antennas. The results validate the absolute calibration methods at the 1-mm level and confirm the presence of significant variations in quality between antennas of different design. For the antenna pattern we directly use the measured phase variations and do not have to fit any functions for the chamber calibrations. We include the results of an earlier test made with a set of identical antennas calibrated at five different institutions: two using the absolute field technique with robot mount and three others applying the standard field calibration with reference antenna.     

融合时间序列环境卫星数据与物候特征的水稻种植区提取

获取高精度的区域水稻种植面积对于农业规划、配置与决策具有重要意义。区域尺度的水稻面积获取依赖于高时空分辨率影像,但受卫星回访周期和气候影响,难以获取足够时间序列的高时空分辨率影像,从而影响水稻种植面积遥感提取的精度。为此,提出适应于中国南方多雨云天气地区,基于国产环境卫星(HJ-1A/1B)与MODIS融合数据的水稻种植面积提取的新方法。以洞庭湖区为实验区,利用STARFM模型融合环境卫星NDVI数据与MODIS13Q1数据,获取时间序列的环境卫星NDVI数据,利用水稻关键期的NDVI数据结合物候特征参数对水稻种植区域进行提取。结果表明,该方法能有效提取区域水稻种植的面积,水稻种植面积提取的总体精度与Kappa系数分别达到91.71%与0.9024,分类结果明显优于仅采用多光谱影像或NDVI数据。该研究为中国南方多雨云天气地区水稻种植面积提取提供了有效的方法。     

Geoid and high resolution sea surface topography modelling in the mediterranean from gravimetry,altimetry and GOCE data: evaluation by simulation

The determination of local geoid models has traditionally been carried out on land and at sea using gravity anomaly and satellite altimetry data, while it will be aided by the data expected from satellite missions such as those from the Gravity field and steady-state ocean circulation explorer (GOCE). To assess the performance of heterogeneous data combination to local geoid determination, simulated data for the central Mediterranean Sea are analyzed. These data include marine and land gravity anomalies, altimetric sea surface heights, and GOCE observations processed with the space-wise approach. A spectral analysis of the aforementioned data shows their complementary character. GOCE data cover long wavelengths and account for the lack of such information from gravity anomalies. This is exploited for the estimation of local covariance function models, where it is seen that models computed with GOCE data and gravity anomaly empirical covariance functions perform better than models computed without GOCE data. The geoid is estimated by different data combinations and the results show that GOCE data improve the solutions for areas covered poorly with other data types, while also accounting for any long wavelength errors of the adopted reference model that exist even when the ground gravity data are dense. At sea, the altimetric data provide the dominant geoid information. However, the geoid accuracy is sensitive to orbit calibration errors and unmodeled sea surface topography (SST) effects. If such effects are present, the combination of GOCE and gravity anomaly data can improve the geoid accuracy. The present work also presents results from simulations for the recovery of the stationary SST, which show that the combination of geoid heights obtained from a spherical harmonic geopotential model derived from GOCE with satellite altimetry data can provide SST models with some centimeters of error. However, combining data from GOCE with gravity anomalies in a collocation approach can result in the estimation of a higher resolution geoid, more suitable for high resolution mean dynamic SST modeling. Such simulations can be performed toward the development and evaluation of SST recovery methods.     

GIS 优化 SRM 储流函数模型变量与参数方法研究

水文模型的结构确定后,模型参数的选择对水文模型整体性能和水文预报结果的好坏有着至关重要的影响。为使得模型的输出值能够尽可能地接近实际值,文中结合数字高程模型(DEM)和流域地质资料,使用GIS方法分别从优化数据源和加入地质因素影响分析优化模型参数两方面入手,对SRM储流函数模型进行优化,并以北海道沙流川流域为例对上述方法进行检验,结果说明了方法的可行性。     

Monitoring urbanization dynamics in India using DMSP/OLS night time lights and SPOT-VGT data

India is a rapidly urbanizing country and has experienced profound changes in the spatial structure of urban areas. This study endeavours to illuminate the process of urbanization in India using Defence Meteorological Satellites Program – Operational Linescan System (DMSP-OLS) night time lights (NTLs) and SPOT vegetation (VGT) dataset for the period 1998–2008. Satellite imagery of NTLs provides an efficient way to map urban areas at global and national scales. DMSP/OLS dataset however lacks continuity and comparability; hence the dataset was first intercalibrated using second order polynomial regression equation. The intercalibrated dataset along with SPOT-VGT dataset for the year 1998 and 2008 were subjected to a support vector machine (SVM) method to extract urban areas. SVM is semi-automated technique that overcomes the problems associated with the thresholding methods for NTLs data and hence enables for regional and national scale assessment of urbanization. The extracted urban areas were validated with Google Earth images and global urban extent maps. Spatial metrics were calculated and analyzed state-wise to understand the dynamism of urban areas in India. Significant changes in urban proportion were observed in Tamil Nadu, Punjab and Kerala while other states also showed a high degree of changes in area wise urban proportion.     

Monitoring aboveground forest biomass dynamics over three decades using Landsat time-series and single-date inventory data

Understanding forest biomass dynamics is crucial for carbon and environmental monitoring, especially in the context of climate change. In this study, we propose a robust approach for monitoring aboveground forest biomass (AGB) dynamics by combining Landsat time-series with single-date inventory data. We developed a Random Forest (RF) based kNN model to produce annual maps of AGB from 1988 to 2017 over 7.2 million ha of forests in Victoria, Australia. The model was internally evaluated using a bootstrapping technique. Predictions of AGB and its change were then independently evaluated using multi-temporal Lidar data (2008 and 2016). To understand how natural and anthropogenic processes impact forest AGB, we analysed trends in relation to the history of disturbance and recovery. Specifically, change metrics (e.g., AGB loss and gain, Years to Recovery - Y2R) were calculated at the pixel level to characterise the patterns of AGB change resulting from forest dynamics. The imputation model achieved a RMSE value of 132.9 Mg ha⁻¹ (RMSE% = 46.3%) and R² value of 0.56. Independent assessments of prediction maps in 2008 and 2016 using Lidar-based AGB data achieved relatively high accuracies, with a RMSE of 108.6 Mg ha⁻¹ and 135.9 Mg ha⁻¹ for 2008 and 2016, respectively. Annual validations of AGB maps using un-changed, homogenous Lidar plots suggest that our model is transferable through time (RMSE ranging from 109.65 Mg ha⁻¹ to 112.27 Mg ha⁻¹ and RMSE% ranging from 25.38% to 25.99%). In addition, changes in AGB values associated with forest disturbance and recovery (decrease and increase, respectively) were captured by predicted maps. AGB change metrics indicate that AGB loss and Y2R varied across bioregions and were highly dependent on levels of disturbance severity (i.e., a greater loss and longer recovery time were associated with a higher severity disturbance). On average, high severity fire burnt from 200 Mg ha⁻¹ to 550 Mg ha⁻¹ of AGB and required up to 15 years to recover while clear-fell logging caused a reduction in 250 Mg ha⁻¹ to 600 Mg ha⁻¹ of AGB and required nearly 20 years to recover. In addition, AGB within un-disturbed forests showed statistically significant but monotonic trends, suggesting a mild gradual drop over time across most bioregions. Our methods are designed to support forest managers and researchers in developing forest monitoring systems, especially in developing regions, where only a single date forestry inventory exists.