Nonlinear least-squares method via an isomorphic geometrical setup

The resolution of a nonlinear parametric adjustment model is addressed through an isomorphic geometrical setup with tensor structure and notation, represented by a u-dimensional “model surface” embedded in a flat n-dimensional “observational space”. Then observations correspond to the observational-space coordinates of the pointQ, theu initial parameters correspond to the model-surface coordinates of the “initial” pointP, and theu adjusted parameters correspond to the model-surface coordinates of the “least-squares” point . The least-squares criterion results in a minimum-distance property implying that the vector Q must be orthogonal to the model surface. The geometrical setup leads to the solution of modified normal equations, characterized by a positive-definite matrix. The latter contains second-order and, optionally, thirdorder partial derivatives of the observables with respect to the parameters. This approach significantly shortens the convergence process as compared to the standard (linearized) method.     

Seismic hazard monitoring in Kamchatka and its applications to the M = 6.6 earthquake of 6 October 1987

The earthquake of 6 October 1987 (M = 6.6), which occurred near the Shipunsky Cape, Kamchatka, was the largest crustal event in the vicinity of the main city of Kamchatka — Petropavlovsk-Kamchatsky — during the last three decades. It was followed by numerous aftershocks. This earthquake allowed us to test the effectiveness of the seismic hazard monitoring in Kamchatka, including the seismological, geodetic and hydrogeochemical surveys. The seismic survey provided the location and source nature of the main shock and aftershocks and the seismic environment of the main shock. The geodetic and hydrogeochemical surveys have yielded data on the response to earthquakes of the Earth's surface deformations, water level, and chemical elements concentration in the underground water. As a result, the following data were obtained:

u

The earthquake of 6 October had a seismic moment 4–10 E18 Nm, thrust type of faulting and the source volume of 20 × 20 × 10 km³. The maximum intensity was VI–VII (MSK-64 scale) and maximum acceleration 88 cm/s².

Before this event, a relative increase in the number of the upper mantle (depth more than 100 km) moderate magnitude earthquakes during 5 years and a one-year period of seismic quiescence for small shallow earthquakes, were recognized. Significant anomalies in HCO₃ and H₃BO₃ concentrations in the underground waters were observed in the wells a week before the main shock.

     

Low-low satellite-to-satellite tracking: a comparison between analytical linear orbit perturbation theory and numerical integration

Low-low satellite-to-satellite tracking (ll-SST) range-rate observations have been predicted by two methods: one based on a linear perturbation theory in combination with the Hill equations, and one based on solving the equations of motion of two low-flying satellites by numerical integration. The two methods produce almost equivalent Fourier spectra of the range-rate observations after properly taking into account a few resonant terms. For a typical GRACE-type configuration, where the two satellites trail each other at a distance of 300 km at an altitude of 460 km, and in the presence of the EGM96 gravity field model, complete to spherical harmonic degree and order 70, the agreement between the Fourier spectra is about 1 mm/s compared to a root-mean-square (RMS) value of more than 220 mm/s for the range-rate signal. The discrepancy of 1 mm/s can be reduced significantly when not taking into account perturbations caused by the J₂ term. Excluding the J₂ term, the agreement between the two methods improves to 0.4 mm/s compared to a RMS value of 6 mm/s for the range-rate signal. These values are 0.01 and 2.3 mm/s when ignoring the spectrum for frequencies below two cycles per orbital revolution, reducing the discrepancy even further to about 0.5% of the signal. The selected linear perturbation theory is thus capable of modeling gravity field induced range-rate observations with very high precision for a large part of the spectrum.     

Wet tropospheric effects on precise relative GPS height determination

Summary Considerable interest has been generated recently in the use of the Global Positioning System (GPS) for precise height determination. A major error source in these measurements is the propagation delay due to atmospheric water vapour. In order to achieve the high precisions required for such applications as absolute sea-level monitoring improvement of wet delay modelling is necessary. Results from a GPS campaign show a significant correlation (0.91) between the variability of the wet delay measured using a water vapour radiometer (WVR) at the Onsala site and the absolute value of the residual error in the height determination of a 134 km baseline from Onsala to Jönköping. This correlation indicates that the atmosphericvariability as inferred from the WVR data includes information on the quality of the GPS height estimate. During periods of high atmospheric activity, e.g., during the passage of a weather front, the use of a six-parameter gradient model reduces the spread for the vertical coordinate from 40 mm to 20 mm (with standard deviations of 17 mm and 9 mm respectively) over the 134 km baseline (less than 1 × 10^–7) using 8 hour data spans on 11 different days over a six month period.     

Land cover dynamics in Garhwal Himalayas — a case study of balkhila sub-watershed

The present paper focuses on land cover dynamics pattern in Balkhila sub-watershed situated in Garhwal Himalayas. This study contributes example from human shaped ecosystem in mountainous regions where fragmentation of natural resources is active. The remote sensing and GIS has been used to understand the land cover dynamics along the topography. The results report that the land cover dynamics is dependent on the aspect due to sun illumination. The altitude and slope are no more a barrier for resource extraction and the human activity zone is shifting towards higher altitudes and slopes. The changes are also defined along the road and settlements.     

An automated seamless mosaicing system of multi-charge coupled devices of panchromatic data

Panchromatic data of pixel resolution 5.8 m obtained from IRS-1C and IRS-1D satellites proved to be very useful for mapping purposes. One of the popular data product is the 70 km swath mosaic which is covered by a combination of 3 CCD line sensors, each with 4096 pixels. Each CCD-line sensor with different imaging times causes geometric problems of mosaicing three strips data together. In this paper, we propose the details of the design elements of system that caters to the need for accurate and automatic multi strip image registration without any second resampling of the data. The systematic geometric correction grid mapping is improved to facilitate accurate mosaicing by automatic image registration task that makes use of the overlap data within image strips and image registration is achieved up to sub-pixel level.     

Crop acreage estimation using ERS-1 SAR data

Synthetic Aperture Radar (SAR) data from the European Remote Sensing Satellite (ERS-1) acquired in July, October and November, 1992 covering the kharif season of the region were used separately and in combination to identify the major crops and for estimation of their acreage before harvest Separability indices were calculated for major cover types and it was found that single-date SAR data cannot be used for accurate identification of various crops. Multi-temporal colour composite facilitated better identification of crop types. Comparison of area estimates made with ERS-1 SAR and IRS-1B LISS II data showed that the commonly used digital data analysis techniques (per pixel classifiers) are not adequate for accurate estimation of crop acreage using SAR data.     

A review of algebraic constraints in terrestrial reference frame datum definition

 Station coordinates are combined with velocities estimated by space geodesy techniques to produce the International Terrestrial Reference System. The input is sets of coordinates and velocities calculated by International Earth Rotation Service analysis centers using space geodesy techniques. The working reference system of individual analysis centers is generally conventionally defined. However, the implications of such processing can have an effect on the resulting combined set. The problem of datum definition as a function of coordinate combinations is reviewed. In particular, the problem of minimum constraints is clearly emphasized and the reference system effect is defined. The goal is to build a process that could be used generally to remove uncertainties in the underlying coordinate system without disturbing the underlying information with additional unnecessary information. Received: 25 January 1999 / Accepted: 20 September 2000     

Solution and econometric estimation of spatial dynamic models in continuous space and continuous time

 This paper presents and demonstrates a general approach to solving spatial dynamic models in continuous space and continuous time that characterize the behaviour of intertemporally and interspatially optimizing agents and estimating from discrete data the parameters of such models. The approach involves the use of a projection method to solve the models and a quasi-Newton algorithm to update quasi-FIML parameter estimates. Received: 26 July 2000 / Accepted: 31 January 2001     

Corridor analysis in Rajaji-Corbett elephant reserve — A Remote sensing and GIS approach

This paper deals with corridor analysis in Rajaji-Corbett Elephant Reserve in the Shivaliks of Uttaranchal state. Efforts were made to detect changes in the state of forest covervis- à -vis the status of corridors during the three periods i.e. 1967. 1986 and 1998 using remote sensing and GIS. TheERDAS Imagine digital image processing andArcView GIS software packages were used for this purpose. Temporal satellite imagery and ground observations in the Rajaji-Corbett Elephant Reserve revealed forest loss, degradation and disturbances in the corridor areas, hindering elephant movement and restricting them to forests islands. Motichur-Chilla corridor, despite being a highly favoured habitat for the elephants was found to be highly threatened followed by Kotdwar and Ramnagar corridors. Construction of wide bridges across Kunao-Chilla Canal, recreation of corridors through reforestation, reduction of all kinds of pressures in the corridor areas and providing higher protection to corridors are recommended.