A note on observation decorrelation, variances of residuals, and redundancy numbers

A technique for processing correlated observations suitable for large, sparse, least-squares adjustments is reviewed. Correlated coordinate differences derived from the Global Positioning System are used as illustrative examples. However, the methods examined are suitable for all types of correlated observations. The computation of variances of residuals, redundancy numbers, and marginally detectable errors is considered for sparse systems.     

Combination of GPS and VLBI on the observation level during CONT11—common parameters,ties and inter-technique biases

Multi-technique space geodetic analysis software has been developed which allows to combine data on the observation level. In addition to local tie information, site-wise common parameters, i.e., troposphere and clocks, can be estimated with this software. Thus, it will be discussed how common parameters have to be estimated and where biases/offsets need to be taken into account. To test such a novel concept, Global Positioning System (GPS) and Very Long Baseline Interferometry (VLBI) data from the CONT11 campaign are being utilized. Since the VLBI baselines of this campaign extend over several thousands of kilometers, GPS data are processed in precise-point positioning mode and satellite orbits and clocks are kept fixed to the IGS final products. From the obtained results, it can be shown that the combination of space geodetic data on the observation level leads to a consistent improvement of station position repeatability as well as nuisance parameters like troposphere estimates. Furthermore, estimation of common parameters (troposphere or clocks) at co-located sites helps to improve the solution further and derive an utmost physically consistent model of the concerned parameters.     

Global positioning system carrier phase: Description and use

After removing the modulation from the Global Positioning System (GPS) signal (L ₁ orL ₂) a pure carrier signal remains. Suppose this carrier is continuously and precisely tracked by aGPS receiver. Furthermore, suppose the phase of the carrier is periodically measured and recorded (nearly simultaneously at two or more locations) with respect to receiver oscillators having the same nominal frequency as theGPS carrier. This paper first considers alternative modeling and processing approaches to these observational data for static operations. Then an approach to dynamic relative positioning using triple differences is presented. This approach should lend itself to performing centimeter accuracy relative surveys in seconds rather than hours. An approach to fixing cycle slips, automatically, is included.     

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.     

A new coarse-time GPS positioning algorithm using combined Doppler and code-phase measurements

A new coarse-time Global Positioning System (GPS) positioning algorithm based on the use of Doppler and code-phase measurements is proposed and described. The proposed method was demonstrated to be essential for reducing the time to first fix and the power consumption in a GPS receiver. Only 1 ms of data is required to obtain a positioning fix with accuracy comparable to that of the traditional GPS navigation algorithm. The algorithm is divided into two parts. In the first part, the Doppler measurement of the GPS signal is used to determine the coarse user position. With proper constraints, the required time accuracy for the Doppler measurements can be relaxed to be as long as 12 h. In the second part of the algorithm, the accurate user position is calculated by means of the 1 ms code-phase data. The traditional tracking process is no longer necessary in the proposed algorithm. Using the acquired 1-ms code-phase measurement, the positioning accuracy was determined to be approximately a few tens of meters in our experimental results. However, if the data length is extended to 10 ms, the positioning accuracy can be improved to within 10–20 m, which is similar to that of the traditional GPS positioning method. Various experiments were conducted to verify the usefulness of the proposed algorithm.     

Building Extraction from LIDAR Point Cloud Data Using Marked Point Process

This paper presents a new algorithm for building extraction from LIDAR (Light Detection and Ranging) point cloud data on the basis of a marked point process based building model. In this building model, the positions and geometries of buildings are modeled by a point process and its marks, respectively. The geometric marks for buildings include their length, width, direction, height. By Bayesian paradigm, a posterior distribution for the marked point process conditional on the LIDAR point cloud data is obtained. The Reversible Jump Markov Chain Monte Carlo (RJMCMC) based scheme is designed to simulate the posterior distribution. Finally, Maximum A Posteriori (MAP) strategy is used to obtain the optimal building detection. The proposed algorithm is tested by a set of LIDAR point cloud data. The results show its efficiency in complex residential environments.     

A complete VLBI delay model for deforming radio telescopes: the Effelsberg case

Deformations of radio telescopes used in geodetic and astrometric very long baseline interferometry (VLBI) observations belong to the class of systematic error sources which require correction in data analysis. In this paper we present a model for all path length variations in the geometrical optics of radio telescopes which are due to gravitational deformation. The Effelsberg 100 m radio telescope of the Max Planck Institute for Radio Astronomy, Bonn, Germany, has been surveyed by various terrestrial methods. Thus, all necessary information that is needed to model the path length variations is available. Additionally, a ray tracing program has been developed which uses as input the parameters of the measured deformations to produce an independent check of the theoretical model. In this program as well as in the theoretical model, the illumination function plays an important role because it serves as the weighting function for the individual path lengths depending on the distance from the optical axis. For the Effelsberg telescope, the biggest contribution to the total path length variations is the bending of the main beam located along the elevation axis which partly carries the weight of the paraboloid at its vertex. The difference in total path length is almost \(-\) 100 mm when comparing observations at 90 \(^\circ \) and at 0 \(^\circ \) elevation angle. The impact of the path length corrections is validated in a global VLBI analysis. The application of the correction model leads to a change in the vertical position of \(+120\)  mm. This is more than the maximum path length, but the effect can be explained by the shape of the correction function.     

A New Approach to 3D Dense LiDAR Data Classification in Urban Environment

Classification of Mobile Mapping LiDAR (Light Detection and Ranging) data is a challenge in the research community since the day when laser scanner system were integrated and mounted on vehicles for collection of 3D data in urban environment. The approach proposed here for classifying LiDAR data is analogous to the process followed for classifying data from satellite images. Pixel based and segmentation based methods have been employed in past for classifying images obtained from satellites. These methods were based on spectral properties of objects present in the images. But for Mobile mapping LiDAR data this approach has been applied and tested for the first time. The properties of this data are completely different from that of satellite images. So even if the basic approach remains the same, many changes have to be made in the entire classification process. The paper here aims to propose the basic procedure of using pixel-wise classification on dense 3D LiDAR data.     

Comparison of individual and combined zenith tropospheric delay estimations during CONT08 campaign

CONT campaigns are 2-week campaigns of continuous VLBI observations. The IERS working group on combination at the observation level uses these campaigns to study such combinations. In this work, combinations of DORIS, GPS, SLR, and VLBI technique measurements are studied during CONT08. We present different results concerning the use of common zenith tropospheric delay (ZTD) during the combination. We compare the ZTD obtained separately using each individual technique data processing, the combined ZTD, and the ZTD derived from a meteorological model. This resulted in a high level of consistency between each of these ZTD at a sub-centimeter level, a consistency which especially depends on the number of observations per estimated ZTD and the humidity level in the troposphere. We noted that GPS provides the main information about the combined ZTD, the other techniques providing complementary information when a lack of GPS observations occurs.     

Geoinformatic Modeling of Groundwater Resource Mapping of Shear Zone Regions- A Case Study on Attur Valley,Tamil Nadu,India

The general tendency of mapping groundwater resource using remote sensing and Geographic Information System (GIS) techniques involve assigning higher weightage to geomorphology. But this cannot be used as a thumb rule everywhere, especially an area where many ductile and brittle zones are prevalent. The influence of texture and structure of sheared rocks might play a control over retaining and permitting groundwater to flow. Attur valley is characterized by the presence of many shear zones and faults and hence the rocks are highly fissile within the shear zones. The present study tries to establish a new ranking and weightage scheme and hence a new spatial model for groundwater resource mapping in shear zone area like Attur Valley. This spatial model can be verified with field data such as water level data, pump test and resistivity data.     

Deformation analysis of the Lake Urmia causeway (LUC) embankments in northwest Iran: insights from multi-sensor interferometry synthetic aperture radar (InSAR) data and finite element modeling (FEM)

This paper presents deformation analysis of Lake Urmia causeway (LUC) embankments in northwest Iran using observations from interferometry synthetic aperture radar (InSAR) and finite element model (FEM) simulation. 58 SAR images including 10 ALOS, 30 Envisat and 18 TerraSAR-X are used to assess settlement of the embankments during 2003–2013. The interferometric dataset includes 140 differential interferograms which are processed using InSAR time series technique of small baseline subset approach. The results show a clear indication of large deformation on the embankments with peak amplitude of \(>\) 50 mm/year in 2003–2010, increasing to \(>\!\!80\)  mm/year in 2012–2013 in the line of sight (LOS) direction from ground to the satellite. 2D decomposition of InSAR observations from Envisat and ALOS satellites that overlap in the years 2007–2010 shows that the rate of the vertical settlement and horizontal motion is not uniform along the embankments; Both eastern and western embankments show significant vertical motion, while horizontal motion plays a more significant role in eastern embankment than western embankment. The InSAR results are then used to simulate deformation using FEM at two cross-sections at the distance of 4 and 9 km from the most western edge of the LUC for which detailed stratigraphy data are available. Results suggest that consolidation due to dissipation of excess pore pressure in embankments can satisfactory predict settlement of the LUC embankments. Our numerical modeling indicates that nearly half of the consolidation since the construction time of the causeway 30 years ago has been done.     

The airborne laser ranging system,its capabilities and applications

The Airborne Laser Ranging System is a proposed multibeam short pulse laser ranging system on board an aircraft. It simultaneously measures the distances between the aircraft and six laser retroreflectors (targets) deployed on the Earth's surface. Depending on the host aircraft and terrain characteristics, the system can interrogate hundreds of targets distributed over an area as large as 6×10⁴ sq. kilometers in a matter of hours. Potentially, a total of 1.3 million individual range measurements can be made in a six hour flight. The precision of these range measurements is approximately ±1 cm. These measurements are then used in a procedure which is basically an extension of trilateration techniques to derive the intersite vector between the laser ground targets. By repeating the estimation of the intersite vector, strain and strain rate errors can be estimated. These quantities are essential for crustal dynamic studies which include determination and monitoring of regional strain in the vicinity of active fault zones, land subsidence, and edifice building preceding volcanic eruptions.     

Further remarks on the altimetry-gravimetry problems

The two Altimetry-Gravimetry linearized boundary value problems, both characterized by gravity anomaly boundary data on the continents, while on the oceans either the gravity potential or its radial derivative is assumed to be known, are investigated in their modified version, in the spherical boundary approximation. A constant unknown bias is added to the data set on the oceans, while a vanishing zerodegree harmonic component is assumed as additional constraint. Existence, uniqueness and regularity results are obtained, it is proved that, in order to obtain regularization of the solution, not only regular data must be imposed, on the continents and on the oceans, but suitable constraints on the coast line are needed.     

Mapping of wastelands through remote sensing techniques and their reclamation—A case study

Wasteland map (1∶100,000) of Rewasa catchment (Sikar district) has been prepared using aerial photographs and Landsat TM imagery. Thematic Mapper data were helpful in identifying the types of wastelands and details could be derived from the aerial photographs. The types of wastelands identified are sands, gullied land, salt affected areas, and barren rocky area. Depending upon the nature of wasteland, suitable rehabilitation measures like plantations, afforestation have been proposed.     

Assessment of ECMWF-derived tropospheric delay models within the EUREF Permanent Network

The Global Positioning System (GPS) observations from the EUREF Permanent Network (EPN) are routinely analyzed by the EPN analysis centers using a tropospheric delay modeling based on standard pressure values, the Niell Mapping Functions (NMF), a cutoff angle of 3° and down-weighting of low elevation observations. We investigate the impact on EPN station heights and Zenith Total Delay (ZTD) estimates when changing to improved models recommended in the updated 2003 International Earth Rotation and Reference Systems Service (IERS) Conventions, which are the Vienna Mapping Functions 1 (VMF1) and zenith hydrostatic delays derived from numerical weather models, or the empirical Global Mapping Functions (GMF) and the empirical Global Pressure and Temperature (GPT) model. A 1-year Global Positioning System (GPS) data set of 50 regionally distributed EPN/IGS (International GNSS Service) stations is processed. The GPS analysis with cutoff elevation angles of 3, 5, and 10° revealed that changing to the new recommended models introduces biases in station heights in the northern part of Europe by 2–3 mm if the cutoff is lower than 5°. However, since large weather changes at synoptic time scales are not accounted for in the empirical models, repeatability of height and ZTD time series are improved with the use of a priori Zenith Hydrostatic Delays (ZHDs) derived from numerical weather models and VMF1. With a cutoff angle of 3°, the repeatability of station heights in the northern part of Europe is improved by 3–4 mm.     

Quantification of carbon stock to understand two different forest management regimes in Kayar Khola watershed,Chitwan, Nepal

The impact of forest management activities on the ability of forest ecosystems to sequester and store atmospheric carbon is of increasing scientific and social concern. This is because a quantitative understanding of how forest management enhances carbon storage is lacking in most forest management regimes. In this paper two forest regimes, government and community-managed, in Kayar Khola watershed, Chitwan, Nepal were evaluated based on field data, very high resolution (VHR) GeoEye-1 satellite image and airborne LiDAR data. Individual tree crowns were generated using multi-resolution segmentation, which was followed by two tree species classification (Shorea robusta and Other species). Species allometric equations were used in both forest regimes for above ground biomass (AGB) estimation, mapping and comparison. The image objects generated were classified per species and resulted in 70 and 82 % accuracy for community and government forests, respectively. Development of the relationship between crown projection area (CPA), height, and AGB resulted in accuracies of R² range from 0.62 to 0.81, and RMSE range from 10 to 25 % for Shorea robusta and other species respectively. The average carbon stock was found to be 244 and 140 tC/ha for community and government forests respectively. The synergistic use of optical and LiDAR data has been successful in this study in understanding the forest management systems.     

A demonstration of 1–2 parts in 107 accuracy using GPS

By interferometric analysis ofGPS phase observations made at Owens Valley, Mojave, and Mammoth Lakes, California, we determined the coordinate components of the71–245–313 km triangle of baselines connecting these sites. A separate determination was made on each of four days, April 1–4, 1985. The satellite ephemerides used in these determinations had been derived from observations on other baselines. Therms scatters of the four daily determinations of baseline vector components about their respective means ranged from a minimum of6 mm for the north component of the71-km baseline to a maximum of34 mm for the vertical component of the245-km baseline. To test accuracy, we compared the mean of ourGPS determinations of the245-km baseline between Owens Valley and Mojave with independent determinations by others using very-long-baseline interferometry(VLBI) and satellite laser ranging(SLR). TheGPS-VLBI difference was within 2 parts in10 ⁷ for every vector component. TheGPS-SLR difference was within6 parts in10 ⁸ in the horizontal coordinates, but83 mm in height.     

A note on palaeozoic bryozoa from the nayar phyllite formation,pauri garhwal,U.P.

This note records the discovery of bryozoan fossils collected from Nayar Phyllite Formation exposed around Satpuli, from a locality SSW of Ira (29° 53′32″: 78° 46′52″) in Pauri Garhwal district, U.P. On the basis of preliminary studies the fossil assemblage has been assigned to the families Phylloporinidae and Fenestellidae. The assemblage tentatively indicates Ordovician — early Silurian age for the formation.     

Corn Area Extraction by the Integration of MODIS-EVI Time Series Data and China’s Environment Satellite (HJ-1) Data

In this paper, we intent to use the remotely sensed MODerate resolution Imaging Spectroradiometer (MODIS) data and China’s Environment Satellite (HJ-1) data for extracting the corn cultivated area over a regional scale. The high resolution HJ-1 data was to extract corn distribution at a small scale class with Support Vector Machine (SVM). The mean Enhanced Vegetation Index (EVI) time series curve of corn from MODIS was derived for the reference area and validated in a larger area. The MODIS-EVI time series curve derived from the reference area instead of the MODIS-EVI time series curve derived from the study area after validation, which was taken as the standard MODIS-EVI time series curve in for generating a standard MODIS-EVI image of corn. The mean absolute distance (MAD) between the standard MODIS-EVI image of corn and the MODIS-EVI time series image was used to detect the maximum possible extent of corn distribution in the study area. The results showed that the overall accuracy of the method was 82.17 %, with commission and omission errors of 16.85 and 15.40 %, respectively; at the county level, the satellite-estimated corn area and statistical data were well correlated (R ²?=?0.85, N?=?50) for the whole Jilin Province. It indicated that the MODIS data integrated with higher spatial resolution of HJ-1 satellite data could be utilized to enhance the extraction accuracy of corn cultivated area at a larger scale.