Bringing GPS into harsh environments for fully automated deformation monitoring

Engineering projects that require deformation monitoring frequently utilize geodetic sensors to measure displacements of target points located in the deformation zone. In situations where control stations and targets are separated by a kilometer or more, GPS can offer higher precision position updates at more frequent intervals than can normally be achieved using total station technology. For large-scale deformation projects requiring the highest precision, it is therefore advisable to use a combination of the two sensors. In response to the need for high precision, continuous GPS position updates in harsh deformation monitoring environments, a software has been developed that employs triple-differenced carrier-phase measurements in a delayed-state Kalman filter. Two data sets were analyzed to test the capabilities of the software. In the first test, a GPS antenna was displaced using a translation stage to mimic slow deformation. In the second test, data collected at a large open pit mine were processed. It was shown that the delayed-state Kalman filter developed could detect millimeter-level displacements of a GPS antenna. The actual precision attained depends upon the amount of process noise infused at each epoch to accommodate the antenna displacements. Higher process noise values result in quicker detection times, but at the same time increase the noise in the solutions. A slow, 25 mm displacement was detected within 30 min of the full displacement with sigma values in E, N and U of ±10 mm or better. The same displacement could also be detected in less than 5 h with sigma values in E, N and U of ±5 mm or better. The software works best for detecting long period deformations (e.g., 20 mm per day or less) for which sigma values of 1–2 mm are attained in all three solution components. It was also shown that the triple-differenced carrier-phase observation can be used to significantly reduce the effects of residual tropospheric delay that would normally plague double-differenced observations in harsh GPS environments.

A seismic hazard assessment and the results of detailed seismic zoning for urban territories of Sakhalin Island

This work briefly discusses the main features of probabilistic seismic hazard analysis (PSHA). Special attention is paid to the identification and quantification of uncertainties related to seismic source characteristics and seismic engineering models for prediction of strong ground motions. The principal seismic models and the results of PSHA application for detailed seismic zoning of urban territories in Sakhalin Island are presented.     

Estimating urban above ground biomass with multi-scale LiDAR

Background

Urban trees have long been valued for providing ecosystem services (mitigation of the “heat island” effect, suppression of air pollution, etc.); more recently the potential of urban forests to store significant above ground biomass (AGB) has also be recognised. However, urban areas pose particular challenges when assessing AGB due to plasticity of tree form, high species diversity as well as heterogeneous and complex land cover. Remote sensing, in particular light detection and ranging (LiDAR), provide a unique opportunity to assess urban AGB by directly measuring tree structure. In this study, terrestrial LiDAR measurements were used to derive new allometry for the London Borough of Camden, that incorporates the wide range of tree structures typical of an urban setting. Using a wall-to-wall airborne LiDAR dataset, individual trees were then identified across the Borough with a new individual tree detection (ITD) method. The new allometry was subsequently applied to the identified trees, generating a Borough-wide estimate of AGB.

Results

Camden has an estimated median AGB density of 51.6 Mg ha^–1 where maximum AGB density is found in pockets of woodland; terrestrial LiDAR-derived AGB estimates suggest these areas are comparable to temperate and tropical forest. Multiple linear regression of terrestrial LiDAR-derived maximum height and projected crown area explained 93% of variance in tree volume, highlighting the utility of these metrics to characterise diverse tree structure. Locally derived allometry provided accurate estimates of tree volume whereas a Borough-wide allometry tended to overestimate AGB in woodland areas. The new ITD method successfully identified individual trees; however, AGB was underestimated by ≤?25% when compared to terrestrial LiDAR, owing to the inability of ITD to resolve crown overlap. A Monte Carlo uncertainty analysis identified assigning wood density values as the largest source of uncertainty when estimating AGB.

Conclusion

Over the coming century global populations are predicted to become increasingly urbanised, leading to an unprecedented expansion of urban land cover. Urban areas will become more important as carbon sinks and effective tools to assess carbon densities in these areas are therefore required. Using multi-scale LiDAR presents an opportunity to achieve this, providing a spatially explicit map of urban forest structure and AGB.

     

Impact of different sampling rates on precise point positioning performance using online processing service

In this study, the effect of different sampling rates (i.e. observation recording interval) on the Precise Point Positioning (PPP) solutions in terms of accurac...     

An optimal determination of variable antenna height and power for in-orbit SAR operations

An optimal power and antenna height determination over a wide range of incidence angles is proposed for in-orbit operations of the spaceborne stripmap synthetic aperture radar (SAR) with variable antenna height and variable radiating element power. The antenna height and the average transmitted radio-frequency power of a radiating element for each incidence angle are determined so that the mapping area can be maximized with the constraints on the signal-to-noise ratio, ambiguities, and swath width, given the SAR geometry (an altitude and incidence angles) and radar wavelength. With the receive-chain power consumption of one transmit/receive module (P/sub r/=10 mW) and the SAR instrument power consumption without the phased array antenna (P/sub I/=5 W), the illustrative design example shows that the mapping area of the proposed SAR for the maximum incidence angle (/spl eta//sup max/=75/spl deg/) is about 269.1% of that resulted from utilizing only the variable antenna height. The proposed optimization would have a more beneficial effect on applications that need the high-incidence-angle operations.     

Theory of satellite ground-track crossovers

The fundamental geometry of satellite ground tracks and their crossover problem are investigated. For idealized nominal ground tracks, the geometry is governed by a few constant parameters whose variations lead to qualitative changes in the crossover solutions. On the basis that the theory to locate crossovers has not been studied in sufficient detail, such changes are described in regard to the number of crossover solutions in conjunction with their bifurcations. Employing the spinor algebra as a tool for establishing the ground-track crossing condition, numerical methodologies to locate crossovers appearing in general dual-satellite ground-track configurations are also presented. The methodologies are applied to precisely determined orbital ephemerides of the GEOSAT, ERS-1, and TOPEX/POSEIDON altimeter satellites. Received: 19 November 1996 / Accepted: 12 May 1997     

Estimation and improvement in the geolocation accuracy of rational polynomial coefficients with minimum GCPs using KOMPSAT-3A

ABSTRACT

In this paper, we propose a method to regenerate Rational Polynomial Coefficients (RPCs) using KOMPSAT-3A imagery and to reduce the geolocation error using minimum ground control points (GCPs). To estimate the new RPCs, the physical sensor model fitted to KOMPSAT-3A imagery was utilized and virtual GCPs over the study area were created. The size of the virtual grid used was 20x20x20. To remove the sensor-related errors in physical sensor model, three different image correction models (image coordinate translation model, shift and drift model, and affine transformation model) were additionally applied. We evaluated our proposed method in two areas within Korea, one in urban (Seoul) and one in rural (Goheung) areas. The results showed that there was a significant improvement after applying the suggested approach in the two areas. The image coordinate translation model is suggested in terms of GCP requirement and expected errors estimated from the error propagation analysis using Gauss–Markov Model (GMM).     

Relation of CME Speed and Magnetic Helicity in CME Source Regions on the Sun during the Early Phase of Solar Cycles 23 and 24

To investigate the relations between coronal mass ejection (CME) speed and magnetic field properties measured in the photospheric surface of CME source regions, we selected 22 disk CMEs in the rising and early maximum phases of the current Solar Cycle 24. For the CME speed, we used two-dimensional (2D) projected speed observed by the Large Angle and Spectroscopic Coronagraph onboard the Solar and Heliospheric Observatory (SOHO/LASCO), as well as a 3D speed calculated from the triangulation method using multi-point observations. Two magnetic parameters of CME source regions were considered: the average of magnetic helicity injection rate and the total unsigned magnetic flux. We then classified the selected CMEs into two groups, showing: i) a monotonically increasing pattern with one sign of helicity (group A: 16 CMEs) and ii) a pattern of significant helicity injection followed by its sign reversal (group B: 6 CMEs). We found that: 1) 3D speed generally shows better correlations with the magnetic parameters than the 2D speed for 22 CME events in Solar Cycle 24; 2) 2D speed and the magnetic parameters of 22 CME events in this solar cycle have lower values than those of 47 CME events in Solar Cycle 23; 3) all events of group B in Solar Cycle 24 occur only after the beginning of the maximum phase, a trend well consistent with that shown in Solar Cycle 23; 4) the 2D speed and the helicity parameter of group B events continue to increase in the declining phase of Solar Cycle 23, while those of group A events abruptly decrease in the same period. Our results indicate that the two CME groups have a different tendency in the solar cycle variations of CME speed and the helicity parameters. Active regions that show a complex helicity evolution pattern tend to appear in the maximum and declining phases, while active regions with a relatively simple helicity evolution pattern appear throughout the whole solar cycle.