Evaluation of gravitational curvatures of a tesseroid in spherical integral kernels

Proper understanding of how the Earth’s mass distributions and redistributions influence the Earth’s gravity field-related functionals is crucial for numerous applications in geodesy, geophysics and related geosciences. Calculations of the gravitational curvatures (GC) have been proposed in geodesy in recent years. In view of future satellite missions, the sixth-order developments of the gradients are becoming requisite. In this paper, a set of 3D integral GC formulas of a tesseroid mass body have been provided by spherical integral kernels in the spatial domain. Based on the Taylor series expansion approach, the numerical expressions of the 3D GC formulas are provided up to sixth order. Moreover, numerical experiments demonstrate the correctness of the 3D Taylor series approach for the GC formulas with order as high as sixth order. Analogous to other gravitational effects (e.g., gravitational potential, gravity vector, gravity gradient tensor), numerically it is found that there exist the very-near-area problem and polar singularity problem in the GC east–east–radial, north–north–radial and radial–radial–radial components in spatial domain, and compared to the other gravitational effects, the relative approximation errors of the GC components are larger due to not only the influence of the geocentric distance but also the influence of the latitude. This study shows that the magnitude of each term for the nonzero GC functionals by a grid resolution 15\(^{{\prime } }\,\times \) 15\(^{{\prime }}\) at GOCE satellite height can reach of about 10\(^{-16}\) m\(^{-1}\) s\(^{2}\) for zero order, 10\(^{-24 }\) or 10\(^{-23}\) m\(^{-1}\) s\(^{2}\) for second order, 10\(^{-29}\) m\(^{-1}\) s\(^{2}\) for fourth order and 10\(^{-35}\) or 10\(^{-34}\) m\(^{-1}\) s\(^{2}\) for sixth order, respectively.     

Determination of Phobos’ rotational parameters by an inertial frame bundle block adjustment

A functional model for a bundle block adjustment in the inertial reference frame was developed, implemented and tested. This approach enables the determination of rotation parameters of planetary bodies on the basis of photogrammetric observations. Tests with a self-consistent synthetic data set showed that the implementation converges reliably toward the expected values of the introduced unknown parameters of the adjustment, e.g., spin pole orientation, and that it can cope with typical observational errors in the data. We applied the model to a data set of Phobos using images from the Mars Express and the Viking mission. With Phobos being in a locked rotation, we computed a forced libration amplitude of \(1.14^\circ \pm 0.03^\circ \) together with a control point network of 685 points.     

What governs the presence of residual vegetation in boreal wildfires?

Wildfires are frequent boreal forest disturbances in Canada, and emulating their patterns with forest harvesting has emerged as a common forest management goal. Wildfires contain many patches of residual vegetation of various size, shape, and composition; understanding their characteristics provides insights for improved emulation criteria. We studied the occurrence of residual vegetation within eleven boreal wildfire events in a natural setting; fires ignited by lightning, no suppression efforts, and no prior anthropogenic disturbances. Relative importance of the measurable geo-environmental factors and their marginal effects on residual presence are studied using Random Forests. These factors included distance from natural firebreaks (wetland, bedrock and non-vegetated areas, and water), land cover, and topographic variables (elevation, slope, and ruggedness index). We present results at spatial resolutions ranging from four to 64 m while emphasizing four and 32 m since they mimic IKONOS- and Landsat-type images. Natural firebreak features, especially the proximity to wetlands, are among the most important variables that explain the likelihood residual occurrence. The majority of residual vegetation areas are concentrated within 100 m of wetlands. Topographic variables, typically important in rugged terrain, are less important in explaining the presence of residuals within our study fires.     

Generating GPS satellite fractional cycle bias for ambiguity-fixed precise point positioning

With the development of precise point positioning (PPP), the School of Geodesy and Geomatics (SGG) at Wuhan University is now routinely producing GPS satellite fractional cycle bias (FCB) products with open access for worldwide PPP users to conduct ambiguity-fixed PPP solution. We provide a brief theoretical background of PPP and present the strategies and models to compute the FCB products. The practical realization of the two-step (wide-lane and narrow-lane) FCB estimation scheme is described in detail. With GPS measurements taken in various situations, i.e., static, dynamic, and on low earth orbit (LEO) satellites, the quality of FCB estimation and the effectiveness of PPP ambiguity resolution (AR) are evaluated. The comparison with CNES FCBs indicated that our FCBs had a good consistency with the CNES ones. For wide-lane FCB, almost all the differences of the two products were within ±0.05 cycles. For narrow-lane FCB, 87.8 % of the differences were located between ±0.05 cycles, and 97.4 % of them were located between ±0.075 cycles. The experimental results showed that, compared with conventional ambiguity-float PPP, the averaged position RMS of static PPP can be improved from (3.6, 1.4, 3.6) to (2.0, 1.0, 2.7) centimeters for ambiguity-fixed PPP. The average accuracy improvement in the east, north, and up components reached 44.4, 28.6, and 25.0 %, respectively. A kinematic, ambiguity-fixed PPP test with observation of 80 min achieved a position accuracy of better than 5 cm at the one-sigma level in all three coordinate components. Compared with the results of ambiguity-float, kinematic PPP, the positioning biases of ambiguity-fixed PPP were improved by about 78.2, 20.8, and 65.1 % in east, north, and up. The RMS of LEO PPP test was improved by about 23.0, 37.0, and 43.0 % for GRACE-A and GRACE-B in radial, tangential, and normal directions when AR was applied to the same data set. These results demonstrated that the SGG FCB products can be produced with high quality for users anywhere around the world to carry out ambiguity-fixed PPP solutions.     

Assessing Performance of Tomo-SAR and Backscattering Coefficient for Hemi-Boreal Forest Aboveground Biomass Estimation

Tomo-SAR technique has been used for hemi-boreal forest height and further forest biomass estimation through allometric equation. Backscattering coefficient especially in longer wavelength (L- or P-band) is thought as a useful parameter for hemi-boreal forest biomass retrieval. The aim of this paper is to assess the performance of vertical backscattering power and backscattering coefficient for hemi-boreal forest aboveground biomass (AGB) estimation with airborne P-band data. The test site locates in southern Sweden called Remningstorp test site, and the in-situ forest AGB ranges from 14 t/ha to 245 t/ha at stand level. Multi-baseline P-band Pol-InSAR data in repeat-path mode collected during March and May in 2007 at Remningstorp test site was used. We found that the correlation coefficient (R) between backscattering coefficient of P-band HH polarization and the in-situ forest biomass reached 0.87. The R for P-band VV backscattering power at 5 m is 0.71 and 10 m is 0.72. Backscattering coefficient in HH polarization and vertical backscattering power at 5 m and 10 m were applied to construct a model for hemi-boreal forest AGB estimation by backward step-wise regression and cross-validation approach. The results showed that the estimated forest AGB ranges from 19 to 240 t/ha, and the constructed model obtained a higher R and smaller RMSE, the value of R is 0.91, RMSE is 30.43 t/ha at Remningstorp test site.     

Automatic Road Change Detection and GIS Updating from High Spatial Remotely-Sensed Imagery

This paper presents a tramework for road network change detectlon In order to upctate the Canadian National Topographic DataBase (NTDB). The methodology has been developed on the basis of road extraction from IRS-pan images (with a 5.8 m spatial resolution) by using a wavelet approach. The feature matching and conflation techniques are used to road change detection and updating. Elementary experiments have showed that the proposed framework could be used for developing an operational road database updating system.     

An intelligent agent based on virtual geographic environment system

On the basis of previous work, this paper designs an intelligent agent based on virtual geographic environment (VGE) system that is characterized by huge data, rapid computation, multi-user, multi-thread and intelligence and issues challenges to traditional GIS models and algorithms. The new advances in software and hardware technology lay a reliable basis for system design, development and application.     

Implementation of OGC web map service based on web service

OGC Web Map Service is one kind of OGC Portrayal Services belongs to OGC Web Service model and it provides multi-platform interoperability of spatial data set. This paper presents a method for implementing OGC Web Map Service based on Web Service technique and introduces the detailed process.     

Virtual reality and its application to slide annimation of landslide

A few researches adopt virtual reality in landslide. So it is significant to combine the physicsbased deformable method with the process simulation of a sliding landslide, this way the dynamic process of a going landslide can be simulated. The analysis of the dynamic law of landslide provides a foundation to forecast and make it realiable to foretell disasters. After recounting the condition of VR and landsliding body simulation, this paper puts forward simulating landslide by the physics-based deformable method. The method divides landsliding into deformable body and rigid one. The deformation of deformable body can be compute by DDA and FEM, and DDA fits the great deformation and FEM is good for the tiny deformation. This paper shows it is feasible to simulate the landslide animation by this method.