Consideration of time-correlated errors in a Kalman filter applicable to GNSS

An algorithm for considering time-correlated errors in a Kalman filter is presented. The algorithm differs from previous implementations in that it does not suffer from numerical problems; does not contain inherent time latency or require reinterpretation of Kalman filter parameters, and gives full consideration to additive white noise that is often still present but ignored in previous implementations. Simulation results indicate that the application of the new algorithm yields more realistic and therefore useful state and covariance information than the standard implementation. Results from a field test of the algorithm applied to the problem of kinematic differential GPS demonstrate that the algorithm provides slightly pessimistic covariance estimates whereas the standard Kalman filter provides optimistic covariance estimates.     

Comparing seven candidate mission configurations for temporal gravity field retrieval through full-scale numerical simulation

The goal of this contribution is to focus on improving the quality of gravity field models in the form of spherical harmonic representation via alternative configuration scenarios applied in future gravimetric satellite missions. We performed full-scale simulations of various mission scenarios within the frame work of the German joint research project “Concepts for future gravity field satellite missions” as part of the Geotechnologies Program, funded by the German Federal Ministry of Education and Research and the German Research Foundation. In contrast to most previous simulation studies including our own previous work, we extended the simulated time span from one to three consecutive months to improve the robustness of the assessed performance. New is that we performed simulations for seven dedicated satellite configurations in addition to the GRACE scenario, serving as a reference baseline. These scenarios include a “GRACE Follow-on” mission (with some modifications to the currently implemented GRACE-FO mission), and an in-line “Bender” mission, in addition to five mission scenarios that include additional cross-track and radial information. Our results clearly confirm the benefit of radial and cross-track measurement information compared to the GRACE along-track observable: the gravity fields recovered from the related alternative mission scenarios are superior in terms of error level and error isotropy. In fact, one of our main findings is that although the noise levels achievable with the particular configurations do vary between the simulated months, their order of performance remains the same. Our findings show also that the advanced pendulums provide the best performance of the investigated single formations, however an accuracy reduced by about 2–4 times in the important long-wavelength part of the spectrum (for spherical harmonic degrees ${<}50$ ), compared to the Bender mission, can be observed. Concerning state-of-the-art mission constraints, in particular the severe restriction of heterodyne lasers on maximum range-rates, only the moderate Pendulum and the Bender-mission are beneficial options, of course in addition to GRACE and GRACE-FO. Furthermore, a Bender-type constellation would result in the most accurate gravity field solution by a factor of about 12 at long wavelengths (up to degree/order 40) and by a factor of about 200 at short wavelengths (up to degree/order 120) compared to the present GRACE solution. Finally, we suggest the Pendulum and the Bender missions as candidate mission configurations depending on the available budget and technological progress.     

Reducing the draconitic errors in GNSS geodetic products

Systematic errors at harmonics of the GPS draconitic year have been found in diverse GPS-derived geodetic products like the geocenter $Z$ -component, station coordinates, $Y$ -pole rate and orbits (i.e. orbit overlaps). The GPS draconitic year is the repeat period of the GPS constellation w.r.t. the Sun which is about 351 days. Different error sources have been proposed which could generate these spurious signals at the draconitic harmonics. In this study, we focus on one of these error sources, namely the radiation pressure orbit modeling deficiencies. For this purpose, three GPS+GLONASS solutions of 8 years (2004–2011) were computed which differ only in the solar radiation pressure (SRP) and satellite attitude models. The models employed in the solutions are: (1) the CODE (5-parameter) radiation pressure model widely used within the International GNSS Service community, (2) the adjustable box-wing model for SRP impacting GPS (and GLONASS) satellites, and (3) the adjustable box-wing model upgraded to use non-nominal yaw attitude, specially for satellites in eclipse seasons. When comparing the first solution with the third one we achieved the following in the GNSS geodetic products. Orbits: the draconitic errors in the orbit overlaps are reduced for the GPS satellites in all the harmonics on average 46, 38 and 57 % for the radial, along-track and cross-track components, while for GLONASS satellites they are mainly reduced in the cross-track component by 39 %. Geocenter $Z$ -component: all the odd draconitic harmonics found when the CODE model is used show a very important reduction (almost disappearing with a 92 % average reduction) with the new radiation pressure models. Earth orientation parameters: the draconitic errors are reduced for the $X$ -pole rate and especially for the $Y$ -pole rate by 24 and 50 % respectively. Station coordinates: all the draconitic harmonics (except the 2nd harmonic in the North component) are reduced in the North, East and Height components, with average reductions of 41, 39 and 35 % respectively. This shows, that part of the draconitic errors currently found in GNSS geodetic products are definitely induced by the CODE radiation pressure orbit modeling deficiencies.     

Land Cover Mapping for Namdapha National Park (Arunachal Pradesh), India Using Harmonized Land Cover Legends

Since last few decades RS-GIS is playing vital role in studying and mapping spatiotemporal responses of land cover, however, as a matter of fact, the mapping outputs largely depend on the expert's/user's preferences because location specific and people specific land cover classification systems are adopted autonomously for image classification in GIS. This may actually lead to an ambiguous definition of a particular land cover type when such different maps are compared at global level. In 1993, FAO and UNEP started efforts for development of a software tool know as LCCS which is a comprehensive standardized tool capable of providing land cover characterization to all possible land cover types in the world regardless of spatial relevance, mapping scale, data collection method etc. Adding to the global efforts of land cover legend harmonization and mapping, this study presents development of harmonized land cover legends for Namdapha National Park located in north-eastern Indian Himalayan region using LCCS and subsequent mapping. The potential of Remote Sensing (RS) and Geographical Information Systems (GIS) in forest/land cover mapping is very well recognized. Therefore, adopting the developed harmonized legends for the study area, land cover mapping was done using RS-GIS approach.     

Multifrequency algorithms for precise point positioning: MAP3

We present the new MAP3 algorithms to perform static precise point positioning (PPP) from multifrequency and multisystem GNSS observations. MAP3 represents a two-step strategy in which the least squares theory is applied twice to estimate smoothed pseudo-distances, initial phase ambiguities, and slant ionospheric delay first, and the absolute receiver position and its clock offset in a second adjustment. Unlike the classic PPP technique, in our new approach, the ionospheric-free linear combination is not used. The combination of signals from different satellite systems is accomplished by taking into account the receiver inter-system bias. MAP3 has been implemented in MATLAB and integrated within a complete PPP software developed on site and named PCube. We test the MAP3 performance numerically and contrast it with other external PPP programs. In general, MAP3 positioning accuracy with low-noise GPS dual-frequency observations is about 2.5 cm in 2-h observation periods, 1 cm in 10 h, and 7 mm after 1 day. This means an improvement in the accuracy in short observation periods of at least 7 mm with respect to the other PPP programs. The MAP3 convergence time is also analyzed and some results obtained from real triple-frequency GPS and GIOVE observations are presented.     

Non-linear motions of Australian geodetic stations induced by non-tidal ocean loading and the passage of tropical cyclones

We investigate daily and sub-daily non-tidal oceanic and atmospheric loading (NTOAL) in the Australian region and put an upper bound on potential site motion examining the effects of tropical cyclone Yasi that crossed the Australian coast in January/February 2011. The dynamic nature of the ocean is important, particularly for northern Australia where the long-term scatter due to daily and sub-daily oceanic changes increases by 20–55 % compared to that estimated using the inverted barometer (IB) assumption. Correcting the daily Global Positioning System (GPS) time series for NTOAL employing either a dynamic ocean model or the IB assumption leads to a reduction of up to 52 % in the weighted scatter of daily coordinate estimates. Differences between the approaches are obscured by seasonal variations in the GPS precision along the northern coast. Two compensating signals during the cyclone require modelling at high spatial and temporal resolution: uplift induced by the atmospheric depression, and subsidence induced by storm surge. The latter dominates ( \(>\) 135 %) the combined net effect that reaches a maximum of 14 mm, and 10 mm near the closest GPS site TOW2. Here, 96 % of the displacement is reached within 15 h due to the rapid transit of cyclones and the quasi-linear nature of the coastline. Consequently, estimating sub-daily NTOAL is necessary to properly account for such a signal that can be 3.5 times larger than its daily-averaged value. We were unable to detect the deformation signal in 2-hourly GPS processing and show that seasonal noise in the Austral summer dominates and precludes GPS detection of the cyclone-related subsidence.     

Prediction of Land Cover Change Using Markov and Cellular Automata Models: Case of Al-Ain,UAE, 1992-2030

The UAE has witnessed rapid urban development and economic growth in recent years. With its ambitious vision to become one of the advanced nations by 2021, planners and policy-makers need to know the most likely direction of future urban development. In this study, remotely sensed imagery coupled with cellular automata models were used to predict land cover in Al Ain, the second largest city in the Emirate of Abu Dhabi. Markov and cellular automata models were used for 1992 and 2006 to predict land cover in 2012. Land Use and Land Cover maps for the study area were derived from 1992, 2006, and 2012 Landsat satellite images (TM, ETM⁺). The models achieved an overall accuracy of approximately 80 %. A Markov model was applied for 2006 and 2012 to predict land cover in 2030. The results conformed to the general trend of the Al Ain Master Plan 2030. This study demonstrates that remote sensing, with the availability of free Landsat data, is a viable technology that could be used to help in the prediction process especially in developing countries, where data availability is a problem.     

Generation of Vegetation Fraction and Surface Albedo Products Over India from Ocean Colour Monitor (OCM) Data Onboard Oceansat-2

The use of Local Area Coverage (LAC) data from Ocean Color Monitor (OCM) sensor of Oceansat-2 with its high radiometric resolution (12 bits/pixel) and 2-day repeat cycle for rapid monitoring of vegetation growth and estimating surface albedo for the Indian region is demonstrated in this study. For the vegetation monitoring, normalized difference vegetation index (NDVI) and vegetation fraction (VF) products were estimated by maximum value composite approach fortnightly and were resampled to 1 km. The surface albedo products were realized by converting narrow-band eight-band spectral reflectance OCM data to a) visible (300–700 nm) and b) broad band (300–3,000 nm) data. For validation, the derived products were compared with respective MODIS global products and found to be in good agreement.     

Monitoring of salt affected soils — A case study using aerial photographs,salyut‐7 space photographs,and landsat TM data

Soils of part of Ukai‐kakarapar Command area, Gujarat (India) have been mapped at 1:25, 000 scale using aerial photographs of December 1977. It was observed that about 36.3% of the area was affected by soil salinity/alkalinity. The test area has been remapped using Salyut‐7 space photographs taken during Indo‐Soviet joint flight in April, 1984. The area affected by soil salinity/alkalinity was found to be substantially higher (80.3%). The earlier mapping using aerial photographs was done when the soil surface was compartively moist (December 1977) as compared to date of Salyut‐7 photography (April 1984), when the soil surface was likely to be devoid of moisture and the salts moved to the surface. To have easy comparision with the map prepared by using aerial photographs, Landsat TM data of December, 1985 was used in which 45.7% of the total area was mapped as salt affected. The extent of area delineated using Landsat TM was higher than that of 1977 but much lesser than the area delineated using Salyut‐7 (MKF‐6M) photographs. This indicated that the increase in the extent of salt affected area in the map prepared using the MKF‐6M photographs might be partly due to actual increase in the salinity/alkalinity and partly due to the seasonal affects. Among the various bands of MKF‐6M, band ‐4 was found to be the best for delineating the salt affected soils. The boundaries were sharper in the FCC and band No.4 of MKF‐6M than in the aerial photographs.