Assessment of troposphere mapping functions using three-dimensional ray-tracing

The troposphere delay is an important source of error for precise GNSS positioning due to its high correlation with the station height parameter. It has been demonstrated that errors in mapping functions can cause sub-annual biases as well as affect the repeatability of GNSS solutions, which is a particular concern for geophysical studies. Three-dimensional ray-tracing through numerical weather models (NWM) is an excellent approach for capturing the directional and daily variation of the tropospheric delay. Due to computational complexity, its use for positioning purposes is limited, but it is an excellent tool for evaluating current state-of-the-art mapping functions used for geodetic positioning. Many mapping functions have been recommended in the past such as the Niell Mapping Function (NMF), Vienna Mapping Function 1 (VMF1), and the Global Mapping Function (GMF), which have been adopted by most IGS analysis centers. A new Global Pressure Temperature model (GPT2) has also been developed, which has been shown to improve upon the original atmospheric model used for the GMF. Although the mapping functions mentioned above use the same functional formulation, they vary in terms of their atmospheric source and calibration approach. A homogeneous data set of three-dimensional ray-traced delays is used to evaluate all components of the mapping functions, including their underlying functional formulation, calibration, and compression method. Additionally, an alternative representation of the VMF1 is generated using the same atmospheric source as the truth data set to evaluate the differences in ray-tracing methods and their effect on the end mapping function. The results of this investigation continue to support the use of the VMF1 as the mapping function of choice when geodetic parameters are of interest. Further support for the GPT2 and GMF as reliable back-ups when the VMF1 is not available was found due to their high consistency with the NWM-derived mapping function. Additionally, a small latitude-dependent bias in station height was found in the current mapping functions. This bias was identified to be due to the assumption of a constant radius of the earth and was largest at the poles and at the equator. Finally, an alternative version of the VMF1 is introduced, namely the UNB-VMF1 which provides users with an independent NWM-derived mapping function to support geodetic positioning.     

Multifrequency algorithms for precise point positioning: MAP3

We present the new MAP3 algorithms to perform static precise point positioning (PPP) from multifrequency and multisystem GNSS observations. MAP3 represents a two-step strategy in which the least squares theory is applied twice to estimate smoothed pseudo-distances, initial phase ambiguities, and slant ionospheric delay first, and the absolute receiver position and its clock offset in a second adjustment. Unlike the classic PPP technique, in our new approach, the ionospheric-free linear combination is not used. The combination of signals from different satellite systems is accomplished by taking into account the receiver inter-system bias. MAP3 has been implemented in MATLAB and integrated within a complete PPP software developed on site and named PCube. We test the MAP3 performance numerically and contrast it with other external PPP programs. In general, MAP3 positioning accuracy with low-noise GPS dual-frequency observations is about 2.5 cm in 2-h observation periods, 1 cm in 10 h, and 7 mm after 1 day. This means an improvement in the accuracy in short observation periods of at least 7 mm with respect to the other PPP programs. The MAP3 convergence time is also analyzed and some results obtained from real triple-frequency GPS and GIOVE observations are presented.     

Use of the L2C signal for inversions of GPS radio occultation data in the neutral atmosphere

Results from processing FORMOSAT-3/COSMIC radio occultations (RO) with the new GPS L2C signal acquired both in phase locked loop (PLL) and open loop (OL) modes are presented. Analysis of L2P, L2C, and L1CA signals acquired in PLL mode shows that in the presence of strong ionospheric scintillation not only L2P tracking, but also L1CA tracking often fails, while L2C tracking is most stable. The use of L2C improves current RO processing in the neutral atmosphere mainly by increasing the number of processed occultations (due to significant reduction in the number of L2 tracking failures) and marginally by a reduction in noise in statistics. The latter is due to the combination of reduced L2C noise (compared to L2P) and increased L1CA noise in those occultations where L2P would have failed. This result suggests application of OL tracking for L1CA and L2C signals throughout an entire occultation to optimally acquire RO data. Two methods of concurrent processing of L1CA and L2C RO signals are considered. Based on testing of individual occultations, these methods allow: (1) reduction in uncertainty of bending angles retrieved by wave optics in the lower troposphere and (2) reduction in small-scale residual errors of the ionospheric correction in the stratosphere.     

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.     

Automated approaches for displaying spatially oriented time‐series data through image processing techniques

This study demonstrated the cartographic implications of automated image processing and computer graphics for the study of time‐series data. Automated statistical and image processing techniques were applied to a case study data set consisting of weekly Crop Moisture Index (CMI) values summarized at 174 state cooperative weather stations within Oklahoma for the time period between February and October, 1980. Computer generated isoline maps of the CMI values were interpolated and rescaled into a series of 32 grid matrices for input into a raster‐based ERDAS image processing software system. Principal Components Analysis (PCA) was used to develop graphic models that synthesized the multi‐temporal data into statistical dimensions that represented the most significant elements of CMI variability. Graphic models of the PCA statistical vectors were displayed individually, in conjunction with eigenvector loadings, and as composite images. Resultant images were analyzed statistically and graphically through the generated CMI grid matrices to ascertain the location, severity, and progression of drought represented in the CMI values. Traditional image processing techniques and devices were combined with the ERDAS software system to transform the multi‐temporal CMI data into multi‐dimensional images that represented the drought's spatial and temporal signature unobscured by redundant information.     

Mapping burns and natural reforestation using thematic Mapper data

Remote sensing techniques are specially suitable to detect and to map areas affected by forest fires. In this work, Landsat 5 Thematic Mapper (TM) data has been used to study a number of forest fires that occurred in the province of Valencia (Spain) and to monitor the vegetation regeneration over burnt areas.

A reference area (non‐burnt forest) was established to assess the change produced by fire. The radiance in the thermal band (10.4–12.5 μm) and the normalized difference in reflectance between near 1R (0.76–0.90 μm) and middle IR (2.08–2.35 μm) were the most suitable parameters to map burnt areas. This index can also be used for monitoring vegetation regeneration in burnt areas. About a month after the fire, the burns show temperatures of 5–6 °C higher than those found in the reference area, and the vegetation index shows negative values whereas the reference area values remain positive. The differences between the burns and the reference area for the vegetation index decrease with time as vegetation regenerates.     

Statewide groundwater‐vulnerability assessment in nebraska using the drastic/GIS model

The paper summarizes a technique for implementing the “DRASTIC” groundwater‐vulnerability model within the context of an automated raster‐based geographic information system. Discussion focuses on a methodological development and a statewide project completed recently in Nebraska. The final products, a comprehensive flow chart illustrating procedures and a map of calculated potential pollution hazard, are presented. The methodology can be executed successfully with minimal training and experience. Areas of Nebraska considered vulnerable to groundwater pollution are identified.     

Determination of turbidity patterns of the Zhujiang estuarine region,South China,using satellite images and a GIS approach

Landsat MSS images and SPOT HRV data were employed to map the changes in turbidity levels in the Zhujiang estuarine region, South China, during the dry season in the period 1973–1987 at low and high tides. Analysis of turbidity trends and changes in the spatial pattern of high turbidity class was carried out with a GIS software—IDRISI. It was concluded that with the use of OVERLAY and RECLASS functions in the GIS approach a large number of turbidity maps could be easily compared and the turbidity trend determined. The GIS approach further permitted evaluation of the importance of such factors as water depths, mean tidal differences, and water salinity to sedimentation in the study region.     

Comparative evaluation of Indian remote sensing satellite and landsat thematic mapper data for geological and geomorphological applications

In order to evaluate the potentials of IRS‐1A Linear Imaging Self‐scanning Sensor (LISS‐I) data for geological and geomorphological applications and also to compare the IRS‐1A LISS‐I data with Landsat Thematic Mapper (TM) data, a study has been attempted for parts of Uttar Pradesh and Madhya Pradesh in Northern India. The first four spectral bands of Landsat TM sensor data which are similar and close to IRS‐1A LISS‐I senor have been utilised for the comparative evaluation. Various techniques employed for both the data set to derive the required geology and geomorphology related information include (i) band combination (ii) spectral response analysis (iii) principal component analysis (iv) supervised classification techniques and (v) visual observation of various outputs generated by the above methods. The Optimum Index Factor (OIF) method adopted for selecting suitable band combinations showed similar OIF rankings for IRS‐1A LISS‐I data and Landsat TM data. It has been visually observed that the band combination 1, 3 & 4 offers relatively better feature display. The spectral responses derived for various major geologic rock units such as Deccan Trap, Vindhyan Formation, Bundelkhand Granite and for a few landcovers such as surface water bodies and black soil show striking similarity in pattern for both LISS‐I and TM. The Principal Component (PC) analysis of both data sets suggested that the total scene brightness tends to dominate in the first PC. The percentage information contributed by PCs 1&2 as also by PCs 1,2 & 3 in both the LISS‐I and TM are comparable. It was observed from the classified image generated by performing supervised classification with a maximum likelihood algorithm that major geomorphic landforms were clearly distinguishable. Thus the qualitative and quantitative evaluation of both IRS‐1A LISS‐I and Landsat TM data showed that significant similarities exist between them. The study also revealed that IRS‐1A LISS‐I data can be effectively used for deriving geology and geomorphology related details.     

Surface description for urban climate studies: A gis based methodology

In order to understand the nature of the urban climate, predict the effects of urbanization, or attempt to ameliorate some of the negative hydroclimatic effects of urbanization, it is necessary to have a good understanding of the role and significance of the urban surface. This paper presents a methodology which uses GIS to represent the characteristics and morphology of the urban surface, which can be used to describe a site objectively, model fluxes, or ensure spatial consistency between measured and modelled data, all of which can vary through time. The methodology is illustrated with respect to Chicago, Illinois. Surface data collected at three spatial scales were used to construct a georeferenced database which was linked to an objective, dynamic accessing system. Spatial variability of surface cover, derived hydroclimatic attributes, and modelled fluxes associated with changes in the urban environment are used to illustrate potential applications of the approach.