Scale factor calibration of a superconducting gravimeter at Esashi Station, Japan, using absolute gravity measurements

The scale factor of a superconducting gravimeter (SG) at the Esashi Earth Tides Station, Japan, was revised by repeating co-located absolute gravity measurements with an FG5 gravimeter. Although the calibration results from the absolute gravimeter (AG) show an apparent secular change in the scale factor of the SG (0.4% for the period 1993–2002), the relative scale factors, which are determined by tidal analysis with the response method, indicate that it has changed by no more than 0.01% during the above period. If the mean scale factor over the 10 years is adopted, a value of –56.082±0.029 Gal/V (1 Gal =10^–8 m s^–2) is obtained, which is about 0.4% smaller than that used in the global geodynamics project (GGP) database. Based on this newly determined scale factor, the tidal gravity factors at Esashi have been re-estimated. The observed tidal factors, corrected for the ocean tide effects with recent models, indicate that the theoretical gravity factors for an inelastic Earth model are more consistent with the observations than are those for an elastic model.     

The hierarchical tessellation model and its use in spatial analysis

The hierarchid tessellation model belongs to a class of spatial data models based on the recursive decomposition of space. The quadtree is one such tessellation and is characterized by square cells and a 1:4 decomposition ratio. To relax these constraints in the tessellation, a generalized hierarchical tessellation data model, called Adaptive Recursive Tessellations (ART), has been proposed. ART increases flexibility in the tessellation by the use of rectangular cells and variable decomposition ratios. In ART, users can specify cell sizes which are intuitively meaningful to their applications, or which can reflect the scales of data. ART is implemented in a data structure called Adaptive Recursive Run-Encoding (ARRE), which is a variant of two-dimensional run-encoding whose running path can vary with the different tessellation structures incorporated in an ART model. Given the recognition of the benefits of implementing statistical spatial analysis in GIS, the use of hierarchical tessellation models such as ART in spatial analysis is discussed. Three examples are introduced to show how ART can: (1) be applied to solve the quadrat size problem in quadrat analysis of point patterns; (2) act as the data model in the variable resolution block kriging technique for geostatistical data to reduce variation in kriging error; and (3) facilitate the evaluation of spatial autocorrelation for area data at multiple map resolutions via the construction of a connectivity matrix for calculating spatial autocorrelation indices based on ARRE.     

Delineating agricultural field boundaries from TM imagery using dyadic wavelet transforms

Agricultural field boundary information is important in the study of agricultural crops by remote sensing. In this paper, the dyadic wavelet transform is adopted to extract field boundaries directly from TM imagery. A field boundary is a set of locally connected pixels characterized by their abrupt spectral intensity variation (singularity) on the image. Preliminary extraction is made on the near-infrared Band 4 (0.76–0.90 μm) image. It is found that most of the boundaries are well characterized by the wavelet modulus maxima curves, particularly for fully vegetated fields, while better results are obtained from the NDVI image for cropped field boundaries. To extract additional boundaries while suppressing unwanted noisy edges, the wavelet transform is then applied to the Wetness image calculated by the Tasselled Cap Transform. A technique is implemented to integrate the edges extracted from each source. The combined results show that the majority of the field boundaries are delineated. The use of a single date image fails to delineate some of the fallow field boundaries and boundaries between two cropped fields with similar spectral properties. Although techniques have been applied to reduce noise, still great difficulties are encountered when trying to delineate edges that are field boundaries from all the detected singularities. Manual editing is still required to trace the remaining boundaries to obtain a complete boundary map. It is recommended that a multi-temporal data set be used to obtain a complete delineation with a fully automated algorithm.     

Response of GPS occultation signals to atmospheric gravity waves and retrieval of gravity wave parameters

We show that the amplitude of the Global Positioning System (GPS) signals in the radio occultation (RO) experiments is an indicator of the activity of the gravity waves (GW) in the atmosphere. The amplitude of the GPS RO signals is more sensitive to the atmospheric wave structures than is the phase. Early investigations used only the phase of the GPS occultation signals for statistical investigation of the GW activity in the height interval 10–40 km on a global scale. In this study, we use the polarization equations and Hilbert transform to find the 1-D GW radio image in the atmosphere by analyzing the amplitude of the RO signal. The radio image, also called the GW portrait, consists of the phase and amplitude of the GW as functions of height. We demonstrate the potential of this method using the amplitude data from GPS/Meteorology (GPS/MET) and satellite mission Challenge Mini-satellite Payload (CHAMP) RO events. The GW activity is nonuniformly distributed with the main contribution associated with the tropopause and the secondary maximums related to the GW breaking regions. Using our method we find the vertical profiles of the horizontal wind perturbations and its vertical gradient associated with the GW influence. The estimated values of the horizontal wind perturbations are in fairly good agreement with radiosonde data. The horizontal wind perturbations v(h) are ±1 to ±5 m s with vertical gradients dv/dh ±0.5 to ±15 m s km at height 10–40 km. The height dependence of the GW vertical wavelength was inferred through the differentiation of the GW phase. Analysis of this dependence using the dispersion relationship for the GW gives the estimation of the projection of the horizontal background wind velocity on the direction of the GW propagation. For the event considered, the magnitude of this projection changes between 1.5 and 10 m s at heights of 10–40 km. We conclude that the amplitude of the GPS occultation signals contain important information about the wave processes in the atmosphere on a global scale.     

Urban change detection related to earthquakes using an adaptive nonlinear mapping of high-resolution images

We propose a nonlinear mapping method for detecting geographical changes. In our system, two images photographed at different points of time are checked based on an iterative nonlinear mapping. The system automatically extracts feature changes of aerial imageries by computing distribution of the image matching score. We evaluate the change detection ability by using the receiver operating characteristic under misregistration cases, including rotational misalignments.     

Evaluation of Information Loss in Digital Elevation Models With Digital Photogrammetric Systems

Information loss is caused when a surface is sampled with a finite interval, such as in the production of a digital elevation model (DEM). This information loss can become the dominant part of the error in a DEM. The ability to quantify information loss enables guidance to be provided for an appropriate choice of grid interval and better accuracy assessment for the DEM. With the use of digital photogrammetric systems, evaluation of information loss has become much easier. This paper describes three methods of evaluating information loss. An example is given of the method which is most appropriate for use with a digital photogrammetric system, based on rock cliff surface data and the VirtuoZo system.     

Maximum-likelihood ambiguity resolution based on Bayesian principle

 Based on the Bayesian principle and the fact that GPS carrier-phase ambiguities are integers, the posterior distribution of the ambiguities and the position parameters is derived. This is then used to derive the maximum posterior likelihood solution of the ambiguities. The accuracy of the integer ambiguity solution and the position parameters is also studied according to the posterior distribution. It is found that the accuracy of the integer solution depends not only on the variance of the corresponding float ambiguity solution but also on its values. Received: 27 July 1999 / Accepted: 22 November 2000     

Geophysical signature of the Tunnunik impact structure,Northwest Territories,Canada

In 2011, the discovery of shatter cones confirmed the 28 km diameter Tunnunik complex impact structure, Northwest Territories, Canada. This study presents the first results of ground‐based electromagnetic, gravimetric, and magnetic surveys over this impact structure. Its central area is characterized by a ~10 km wide negative gravity anomaly of about 3 mGal amplitude, roughly corresponding to the area of shatter cones, and associated with a positive magnetic field anomaly of ~120 nT amplitude and 3 km wavelength. The latter correlates well with the location of the deepest uplifted strata, an impact‐tilted Proterozoic dolomite layer of the Shaler Supergroup exposed near the center of the structure and intruded by dolerite dykes. Locally, electromagnetic field data unveil a conductive superficial formation which corresponds to an 80–100 m thick sand layer covering the impact structure. Based on the measurements of magnetic properties of rock samples, we model the source of the magnetic anomaly as the magnetic sediments of the Shaler Supergroup combined with a core of uplifted crystalline basement with enhanced magnetization. More classically, the low gravity signature is attributed to a reduction in density measured on the brecciated target rocks and to the isolated sand formations. However, the present‐day fractured zone does not extend deeper than ~1 km in our model, indicating a possible 1.5 km of erosion since the time of impact, about 430 Ma ago.