Reducing the draconitic errors in GNSS geodetic products

Systematic errors at harmonics of the GPS draconitic year have been found in diverse GPS-derived geodetic products like the geocenter $Z$ -component, station coordinates, $Y$ -pole rate and orbits (i.e. orbit overlaps). The GPS draconitic year is the repeat period of the GPS constellation w.r.t. the Sun which is about 351 days. Different error sources have been proposed which could generate these spurious signals at the draconitic harmonics. In this study, we focus on one of these error sources, namely the radiation pressure orbit modeling deficiencies. For this purpose, three GPS+GLONASS solutions of 8 years (2004–2011) were computed which differ only in the solar radiation pressure (SRP) and satellite attitude models. The models employed in the solutions are: (1) the CODE (5-parameter) radiation pressure model widely used within the International GNSS Service community, (2) the adjustable box-wing model for SRP impacting GPS (and GLONASS) satellites, and (3) the adjustable box-wing model upgraded to use non-nominal yaw attitude, specially for satellites in eclipse seasons. When comparing the first solution with the third one we achieved the following in the GNSS geodetic products. Orbits: the draconitic errors in the orbit overlaps are reduced for the GPS satellites in all the harmonics on average 46, 38 and 57 % for the radial, along-track and cross-track components, while for GLONASS satellites they are mainly reduced in the cross-track component by 39 %. Geocenter $Z$ -component: all the odd draconitic harmonics found when the CODE model is used show a very important reduction (almost disappearing with a 92 % average reduction) with the new radiation pressure models. Earth orientation parameters: the draconitic errors are reduced for the $X$ -pole rate and especially for the $Y$ -pole rate by 24 and 50 % respectively. Station coordinates: all the draconitic harmonics (except the 2nd harmonic in the North component) are reduced in the North, East and Height components, with average reductions of 41, 39 and 35 % respectively. This shows, that part of the draconitic errors currently found in GNSS geodetic products are definitely induced by the CODE radiation pressure orbit modeling deficiencies.     

The effects of frequency-dependent quasar variability on the celestial reference frame

We examine the relationship between source position stability and astrophysical properties of radio-loud quasars making up the International Celestial Reference Frame (ICRF2). Understanding this relationship is important for improving quasar selection and analysis strategies, and therefore reference frame stability. We construct flux density time series, known as light curves, for 95 of the most frequently observed ICRF2 quasars at both the 2.3 and 8.4 GHz geodetic very long baseline interferometry (VLBI) observing bands. Because the appearance of new quasar components corresponds to an increase in quasar flux density, these light curves alert us about potential changes in source structure before they appear in VLBI images. We test how source position stability depends on three astrophysical parameters: (1) flux density variability at X band; (2) time lag between flares in S and X bands; (3) spectral index root-mean-square (rms), defined as the variability in the ratio between S and X band flux densities. We find that the time lag between S and X band light curves provides a good indicator of position stability: sources with time lags $<$ 0.06 years are significantly more stable ( $>$ 20 % improvement in weighted rms) than sources with larger time lags. A similar improvement is obtained by observing sources with low $(<$ 0.12) spectral index variability. On the other hand, there is no strong dependence of source position stability on flux density variability in a single frequency band. These findings can be understood by interpreting the time lag between S and X band light curves as a measure of the size of the source structure. Monitoring of source flux density at multiple frequencies therefore appears to provide a useful probe of quasar structure on scales important to geodesy. The observed astrometric position of the brightest quasar component (the core) is known to depend on observing frequency. We show how multi-frequency flux density monitoring may allow the dependence on frequency of the relative core positions along the jet to be elucidated. Knowledge of the position–frequency relation has important implications for current and future geodetic VLBI programs, as well as the alignment between the radio and optical celestial reference frames.     

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.     

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.     

Mapping spiny aster infestations with QuickBird imagery

QuickBird satellite imagery acquired in June 2003 and September 2004 was evaluated for detecting the noxious weed spiny aster [Leucosyris spinosa (Benth.) Greene] on a south Texas, USA rangeland area. A subset of each of the satellite images representing a diversity of cover types was extracted and used as a study site. The satellite imagery had a spatial resolution of 2.8 m and contained 11-bit data. Unsupervised and supervised classification techniques were used to classify false colour composite (green, red, and near-infrared bands) images of the study site. Imagery acquired in June was superior to that obtained in September for distinguishing spiny aster infestations. This was attributed to differences in spiny aster phenology between the two dates. An unsupervised classification of the June image showed that spiny aster had producer's and user's accuracies of 90% and 93.1%, respectively, whereas a supervised classification of the June image had producer's and user's accuracies of 90% and 81.8%, respectively. These results indicate that high resolution satellite imagery coupled with image analysis techniques can be used successfully for detecting spiny aster infestations on rangelands.     

A global inventory of coral reef stressors based on satellite observed nighttime lights

In this article, we present a satellite-based approach to gather information about the threat to coral reefs worldwide. Three chosen reef stressors – development, gas flaring and heavily lit fishing boat activity – are analysed using nighttime lights data derived from the Defense Meteorological Satellite Program (DMSP) produced at the National Oceanic & Atmospheric Administration, National Geophysical Data Center (NOAA/NGDC). Nighttime lights represent a direct threat to coral reef ecosystems and are an excellent proxy measure for associated human-caused stressors. A lights proximity index (LPI) is calculated, measuring the distance of coral reef sites to each of the stressors and incorporating the stressor's intensity. Colourized maps visualize the results on a global scale. Area rankings clarify the effects of artificial night lighting on coral reefs on a regional scale. The results should be very useful for reef managers and for state administrations to implement coral reef conservation projects and for the scientific world to conduct further research.